r - distr包-如何在一个窗口中绘制两个图？

Question

我使用 distr 来形成以下发行版:

library(distr)

G1 <- Gammad(shape=1.64, scale=0.766)
G2<- Gammad(shape=0.243, scale=4.414)

现在为了比较这两个分布，我需要将它们绘制在一个窗口中，但我不知道如何做。我尝试了 ggplot，但显然它不适用于 gamma 函数。

Answer 1

使用ggplot

您可以使用stat_function

例如

# data that defines the range of x values you are interested in
DF <-data.frame(x = c(1,8))
ggplot(DF, aes(x=x)) + 
  stat_function(fun = d(G1), aes(colour = 'G1')) + 
  stat_function(fun = d(G2), aes(colour = 'G2')) + 
  scale_colour_manual('Distribution', 
           values = setNames(c('red', 'blue'), c('G1', 'G2')))

使用基础

distr::plot 的帮助文件显示了如何组合绘图。

您需要自己设置mfrow(或mfcol)，然后在绘图调用中设置mfColRow =FALSE。

例如:

par(mfrow=c(2,3))
plot(G1, mfColRow=F)
plot(G2, mfColRow=F)