0

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データがどのように分布しているかを示すために、抽出されたすべての正規分布を別の色でプロットする方法はありますか?

または、一般的に、複数の正規分布を視覚化する方法は?

4

1 に答える 1

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")
}

取得するため

ここに画像の説明を入力

于 2016-04-19T11:52:40.583 に答える