http://www.digistats.net/z/index.php?%C5%FD%B7%D7%B3%D8 http://www.digistats.net/z/index.php?%C5%FD%B7%D7%B3%D8

正規分布の密度関数

  • x にに -5 から 5 まで 500 個のデータを入力。
     > x=seq(-5,5,length=500)
  • 平均が 0.5 , 1.0 , 1.5 という3つの正規分布を描く。そこで、i に3つの数値を入力。
     > i=c(0.5,1.0,1.5)
  • x の値を3つの平均に対応した正規分布の確率密度関数に入力。
     > y<-sapply(i,function(a){dnorm(x,mean=a,sd=1)})
  • y の値をプロット。
     > matplot(x,y,type="l",lty=1,col=1:3,ylab="density")
  • グラフに凡例を追加。凡例はグラフの x の目盛の -4.0 のところから書く。
     > legend(-4,0.4,paste("mean=",i,"sd=",1,sep=""),col=1:3,lty=1,bty="n")
  • 最後に目盛線を点線で追加。
     > abline(h=0,v=0,lty=3)


http://www.digistats.net/image/2008/03/sd1.gif


トップ   差分 バックアップ リロード   一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2010-02-19 (金) 23:59:41 (5172d)