[降雪量のデータが】:玩具RコードIは、降雪観測の数を有する
x <- c(98.044, 107.696, 146.050, 102.870, 131.318, 170.434, 84.836, 154.686,
162.814, 101.854, 103.378, 16.256)
と私はそれに知られている標準偏差の正規分布に従うと言われました25.4しかし未知の平均mu
。ベイジアンフォーミュラを使ってmu
の推論をしなければなりません。
これは、次のmu
mean of snow | 50.8 | 76.2 | 101.6 | 127.0 | 152.4 | 177.8
---------------------------------------------------------------
probability | 0.1 | 0.15 | 0.25 |0.25 | 0.15 | 0.1
---------------------------------------------------------------
の情報では、私がこれまで試してみましたものですが、post
についての最終ラインが正常に動作しません。結果のプロットは水平線を与えます。
library(LearnBayes)
midpts <- c(seq(50.8, 177.8, 30))
prob <- c(0.1, 0.15, 0.25, 0.25, 0.15, 0.1)
p <- seq(50, 180, length = 40000)
histp <- histprior(p, midpts, prob)
plot(p, histp, type = "l")
# posterior density
post <- round(histp * dnorm(x, 115, 42)/sum(histp * dnorm(x, 115, 42)), 3)
plot(p, post, type = "l")