ks.test(X, "pweibull")这个行不行
用的是Kolmogorov-Smirnov检验,在dgof包里面
比如——
> set.seed(1)
> S<-rweibull(1000,1)
> ks.test(S, "pweibull",1)
One-sample Kolmogorov-Smirnov test
data: S
D = 0.0244, p-value = 0.5928
alternative hypothesis: two-sided
结论: D值很小, p-value>0.05,不能拒绝原假设,所以数据集S符合shape=1的weibull分布!