正确显示4个面板分布图的图例，并在R中显示fitdistrplus

Question

我在R中使用fitdistrplus软件包，并试图绘制一个漂亮的4面板分布图，如作者发表的论文(J. of Statistical，vol.64，第4期)：

​

​

但是，当我使用他们在论文中显示的相同代码时，与不同颜色编码的线条/曲线相关联的文本值只是任意放置，您无法分辨哪种颜色对应于哪个系列：

par(mfrow = c(2, 2))
plot.legend <- c("Weibull", "lognormal", "gamma")
denscomp(list(fw, fln, fg), legendtext = plot.legend)
qqcomp(list(fw, fln, fg), legendtext = plot.legend)
cdfcomp(list(fw, fln, fg), legendtext = plot.legend)
ppcomp(list(fw, fln, fg), legendtext = plot.legend)

​

​

我也尝试过使用, xlegend = "bottomright"选项。

Answer 1

如果您希望图例在每个情节中可见，您可以手动操作该位置。下面提供了一些示例代码。此外，使用cex缩放图例将避免图例在很大程度上重叠于绘图。准确的图例位置将基于您的数据。

library(fitdistrplus)
x <- exp(rnorm(100,0,1))

fw <- fitdist(x,distr="weibull")
fln <- fitdist(x,distr="lnorm")
fg <- fitdist(x,distr="gamma")


par(mfrow = c(2, 2))
plot.legend <- c("Weibull", "lognormal", "gamma")
denscomp(list(fw, fln, fg)
         , legendtext = plot.legend
         ,xlegend=10
         ,ylegend=0.5
         ,cex=0.7)
qqcomp(list(fw, fln, fg)
       , legendtext = plot.legend
       ,xlegend=8
       ,ylegend=6
       ,cex=0.7)
cdfcomp(list(fw, fln, fg)
        , legendtext = plot.legend
        ,xlegend=8
        ,ylegend=0.3
        ,cex=0.7)
ppcomp(list(fw, fln, fg)
       , legendtext = plot.legend
       ,xlegend=0.7
       ,ylegend=0.3
       ,cex=0.7)

​

​