[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]ガンマ関数の値の計算(「入門 統計解析 -医学・自然科学編」(東京図書)pp.120-121)(2025.01.24)
- [R]グループに対するダミー変数(性別)(「44の例題で学ぶ計量経済学」(オーム社)pp.231-232)(2023.11.28)
- [R]重回帰モデルにおける仮説検定(「44の例題で学ぶ計量経済学」(オーム社)pp.205-206)(2023.11.27)
- [R]重回帰モデルにおけるt検定(「44の例題で学ぶ計量経済学」(オーム社)pp.185-188)(2023.11.26)
- [R]平均勤続年数と所定内賃金(「線形回帰分析」(朝倉書店)pp.116-118)(2023.11.06)
« [R]n年後の日付を得る | トップページ | [R]λ=1のときのポアソン分布(「統計解析のはなし」(東京図書)pp.93-94) »

コメント