Simples script em R que exibe histograma de posições sequenciais de uma determinada massa de dados e uma plotagem usando curva gaussiana e curva de densidade

loop_hist.r

#!/usr/bin/Rscript

emg <- read.csv("emg_s1.csv", header = FALSE, sep = "|", dec = ",")

View(emg)
attach(emg)
par(mfrow = c(4, 5))

for(data in 1:20) {
  position <- paste(c("V", data), collapse = "")
  header <- paste(c("Histograma de ", position), collapse = "")
  hist(emg[[position]], main = header, xlab = header)
}

for(data in 1:5) {
  position <- paste(c("V", data), collapse = "")
  header <- paste(c("Histograma de ", position), collapse = "")
  hist(emg[[position]], main = header, xlab = header, prob = TRUE)

  media <- mean(emg[[position]])
  desvio <- sd(emg[[position]])

  curve(dnorm(x, mean = media, sd = desvio), add=TRUE, col = "red", xlab="Miles Per Gallon")
  lines(density(emg[[position]]), col = "blue", xlab="Miles Per Gallon")
}