我使用fitdist
包中的函数fitdistrplus
来获得适合我的数据的最佳分布并绘制ppcomp
图形,我使用示例代码?fitdistrplus
替换我的数据和代码。
data(groundbeef)
serving <- groundbeef$serving
fitW <- fitdist(serving, "weibull")
fitg <- fitdist(serving, "gamma")
fitln <- fitdist(serving, "lnorm")
ppcomp(list(fitW, fitg, fitln), legendtext=c("Weibull", "gamma", "lognormal"))
到目前为止,一切都很好!但是看图片的左下角,有一些空间,或者轴不是从零开始的!
所以我用谷歌搜索并找到了两种方法:一种方法在基础图中,xaxs
并被yaxs
使用。
set.seed(1234)
x <- runif(100, min=0, max=100)
y <- runif(100, min=0, max=100)
plot(x,y,xaxs="i",yaxs="i",xlim=c(0,100),ylim=c(0,100))
ggplot2
中的另一种方法expand=c(0,0)
被使用。
df <- data.frame(x=x,y=y)
ggplot(df,aes(x,y)) + geom_point() + scale_x_continuous(expand = c(0,0)) + scale_y_continuous(expand = c(0,0))
所以我尝试用这两种方法来画图ppcomp
,
ppcomp(list(fitW, fitg, fitln), legendtext=c("Weibull", "gamma", "lognormal"),xaxs="i",yaxs="i")
但发生错误:
Error in legend(x = xlegend, y = ylegend, bty = "n", legend = legendtext, :
unused arguments (xaxs = "i", yaxs = "i")
那么,我应该怎么做才能在没有错误的情况下完成目标,或者什么是正确的代码?