Составить несколько объектов FitDist в одном и том же сюжете с разными цветами?

Question

У меня есть список fitdist объекты, которые я сохранил, используя этот кусок кода.

norm_dist_res <- list()
for(i in 1:10)
{
  x <- 1+(8000*(i-1))
  y <- 8000*i
  print (x)
  print(y)
  norm_dist_res[[i]] = norm_dist_res[[i]] <- fitdist(data=as.vector(g_all_p$data[x:y,]), distr="norm")

}

Есть ли способ построить все нормальные распределения, извлеченные из fittest с другим цветом, чтобы показать, как распределяются данные?

Или вообще как визуализировать несколько нормальных распределений?

Answer 1

Вы оцениваете параметры нормального распределения, поэтому просто нанесите на график плотности.

## Don't no what g_all_p is, so simplifying the data
library(fitdistrplus)
norm_dist_res <- list()
for(i in 1:10)
{
  norm_dist_res[[i]] = norm_dist_res[[i]] <- fitdist(data=rnorm(10), distr="norm")

}

Тогда просто заговор, используя lines и изменение цвета

x = seq(-5, 5, length.out=100)
plot(x, type="n", ylim=c(0, 1), xlim=range(x))
for(i in 1:10) {
  est = norm_dist_res[[i]]$estimate
  lines(x, dnorm(x, est[1], est[2]), col="grey90")
}

Получить