distr包 - 如何在一個窗口中繪製兩個圖？

Question

我使用distr，形成如下分佈：

library(distr) 

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

現在

，以這兩個分佈，我需要他們繪製在一個窗口比較，但我不知道怎麼辦。我試過ggplot，但顯然它不適用於伽瑪函數。

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)