「統計学関連なんでもあり」の過去ログ--- 048

No.22378　グラフ　　【課長】　2017/05/31(Wed) 22:58

ポアソン分布や二項分布が正規分布で近似できるところをグラフで視覚的に見てみたくて，Rでグラフを書こうとしたのですが，うまくゆきません。
正規分布のグラフは
curve(dnorm(x,mean=1,sd=1))
などでうまくゆくのですが，ポアソン分布などはうまくゆきません
curve(dpois(x,lambda=2))などとすると謎の図形が出てきてエラーも50くらい出てきます。
これまで統計をうわべだけ学んでいて基礎から学んだらわからないことだらけであることに気づきました。程度の低い質問ですが，どうかご教示いただければと思います。

No.22379　Re: グラフ　　【青木繁伸】　2017/05/31(Wed) 23:57

参考までに

# 二項分布
n = 20 # サンプルサイズ
p = 0.3 # 母比率
y = dbinom(0:n, n, prob=p)
x = barplot(y, space=0, ylim=c(0, max(y)*1.1))
text(x, -0.01, 0:n, xpd=TRUE)
Mean = n*p
SD = sqrt(n*p*(1-p))
y = dnorm(0:n, mean=Mean, sd=SD)
lines(0:n+0.5, y, col="red")
points(0:n+0.5, y, pch=19, col="red")

# ポアソン分布
n = 20 # サンプルサイズ
lambda = 10 # ポアソン定数
y = dpois(0:n, lambda)
x = barplot(y, space=0, ylim=c(0, max(y)*1.1))
text(x, -0.01, 0:n, xpd=TRUE)
Mean = lambda
SD = sqrt(lambda)
y = dnorm(0:n, mean=Mean, sd=SD)
lines(0:n+0.5, y, col="red")
points(0:n+0.5, y, pch=19, col="red")

No.22380　Re: グラフ　　【課長】　2017/06/01(Thu) 03:29

青木先生，いつもこんな未熟者のためにご丁寧な回答ありがとうございます。いただいたプログラムをまさに解剖するかのように分析して学習しようと思います。統計の基礎をやりなおすうえで，Rは強力な武器になると思ってますので，これからももしかしたら初歩的な質問をさせてもらうこともあるかもしれませんが，なるべく十分調べてからの質問にしようと思います。ありがとうございます。