Question

有人可以详细解释这段代码如何工作吗？

require(lattice)
?lattice  # essential reading

 data <- dgamma(seq(from=0.00001,to=0.01,by=0.00001),shape = .1, scale = .01)

dfrm <- data.frame(dgam = data, param="s.1.01")
dfrm <- rbind(dfrm, data.frame(dgam =
                                 dgamma( seq(from=0.00001,to=0.01,by=0.00001), 
                                 shape = .2, scale = .01), 
                               param="s.2.01") )
dfrm <- cbind( dfrm, X.val=seq(from=0.00001,to=0.01,by=0.00001) )
str(dfrm)
#'data.frame':  2000 obs. of  3 variables:
# $ dgam : num  5263 2817 1954 1507 1231 ...
# $ param: Factor w/ 2 levels "s.1.01","s.2.01": 1 1 1 1 1 1 1 1 1 1 ...
# $ X.val: num  1e-05 2e-05 3e-05 4e-05 5e-05 6e-05 7e-05 8e-05 9e-05 1e-04 ...
xyplot( dgam ~ X.val , 
        group=param, 
        data=dfrm, type="l")

例如，在哪个部分指定先验和后验？

我已经在How to plot Bayesian prior and posterior distributions in one panel using R?

中找到了答案中的代码。

Answer 1

该答案仅说明了如何在同一图中绘制两条曲线，但完全不涉及在给定先验和数据的情况下计算真实的后验。即，它仅绘制了两个shape=0.1和shape=0.2的Gamma分布，并且在数据框中通过包含因子列param来区分它们，其中"s.1.01"和{{1} }表示各自的分布。

使用R在一个面板中绘制贝叶斯先验和后验分布的解释

1 个答案: