0
私は20万以上の値を持つ大きなデータを持っています。プライバシーとコードの再現性のために、私はmydataを使用して置き換えます。凡例パッケージのfitdistrplusでの凡例R
set.seed(1234)
mydata <- rlnorm(28000000,3.14,1.3)
は、私はパッケージfitdistrplus
で機能fitdist
を選びだしているので、知られているディストリビューションは、最高のmydata
に合わせた見つけたいです。
library(fitdistrplus)
fit.lnorm <- fitdist(mydata,"lnorm")
fit.weibull <- fitdist(mydata, "weibull")
fit.gamma <- fitdist(mydata, "gamma", lower = c(0, 0))
fit.exp <- fitdist(mydata,"exp")
はその後、私は私が最高のフィット分布を選択するためP-Pのプロットを描画するppcomp
関数を使用します。
library(RColorBrewer)
tiff("./pplot.tiff",res = 300,compression = "lzw",height = 6,width = 10,units = "in",pointsize = 12)
ppcomp(list(fit.lnorm,fit.weibull, fit.gamma,fit.exp), fitcol = brewer.pal(9,"Set1")[1:4],legendtext = c("lnorm","weibull", "gamma","exp"))
dev.off()
は絶対に、対数正規最高mydata
収まるが、プロットのlegend
を見て、異なる色のライン注釈が欠落している、唯一のテキスト注釈のショー、私は何をすべきでしょうか?
私はいくつかの値を持つデータセットを試してみました。大きなデータが問題につながります。私は伝説を完璧にするために何をすべきですか?
'ppcomp'の' legend'コマンドに問題があり、おそらく 'lty'を指定する必要があります。 –
@RichardTelfordええ、あなたは正しいのですか? –