« [R]n年後の日付を得る | トップページ | [R]λ=1のときのポアソン分布(「統計解析のはなし」(東京図書)pp.93-94) »

2023年7月31日 (月)

[R]1年間に馬に蹴り殺された兵士の数とポアソン分布(「統計解析のはなし」(東京図書)pp.94-95)

> # データの読み込み
> dtf <- read.csv("data/table2_2_7.csv", header = TRUE)
> # 死亡した兵士数 x
> x <- dtf$x
> # 1年間に兵士がx人死亡した軍団の数 n
> n <- dtf$n
> # 度数の合計 200
> nn <- sum(n)
> # 相対度数 x/200
> rf <- n / nn
> # 死亡した兵士数の平均 x~
> xm <- sum(x * rf)
> # ポアソン分のパラメーターλは平均を採用
> lam <- xm
> # 階級ごとの確率 P(X=x)
> pp <- lam ^ x / factorial(x) * exp(-lam)
> # 計算結果の表示
> print(data.frame(x, n, 相対度数 = rf, ポアソン分布 = pp))
x n 相対度数 ポアソン分布
1 0 109 0.545 0.5433508691
2 1 65 0.325 0.3314440301
3 2 22 0.110 0.1010904292
4 3 3 0.015 0.0205550539
5 4 1 0.005 0.0031346457
6 5 0 0.000 0.0003824268

« [R]n年後の日付を得る | トップページ | [R]λ=1のときのポアソン分布(「統計解析のはなし」(東京図書)pp.93-94) »

R(本の計算を再現)」カテゴリの記事

コメント

コメントを書く

コメントは記事投稿者が公開するまで表示されません。

(ウェブ上には掲載しません)

« [R]n年後の日付を得る | トップページ | [R]λ=1のときのポアソン分布(「統計解析のはなし」(東京図書)pp.93-94) »

無料ブログはココログ

■■

■■■