Question

我尝试了多种方法，包括qqmath，lattice densityplot()以及许多小组功能，例如panel.mathdensity和panel.densityplot。但是，我无法让他们做我想让他们做的事情。

关于此主题的互联网搜索会生成结果，这些结果关注于R中的基础图或不在一个面板中绘制两个分布。我可以使用R基础图形，但是，我还想为每一对绘制几个分布对和一个面板。

书籍“RGraphics”和“Lattice：多变量数据可视化与R”在这方面也无法启发我。

数据通常如下所示：

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

虽然我有更多使用lattice的经验，但我对ggplot或lattice套餐持开放态度。

如果您需要更多信息来帮助我，请告诉我。

Answer 1

您只提供了一个数据对象，但要求构造两个数据对象，因此我将尝试使用不同参数的两个相同分布的项目。当使用lattice或ggplot2时，需要构造一个正确排列的数据帧，这通常是数据的“长”格式。 “group”参数指定使用不同颜色在同一面板中绘图。

require(lattice)
?lattice  # essential reading
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")

如何使用R在一个面板中绘制贝叶斯先验和后验分布？

1 个答案: