我接受我要遵循的值的分布是一个左截尾weibull分布。我知道使用ptrunc命令的这个发行版的参数a、形状和规模:因此,我希望ks.testks.test(myvalues,"pweibull",scale=b,shape=c) #for normal weibullks.test(myvalues,ptrunc(x,"weibull
我试图通过求解两个联立方程来估计威布尔分布的参数: F(Q1)=0.25和F(Q3)=0.75。代数上,我知道答案应该是:形状参数a=0.8038,缩放参数b=1889,但是无论我把什么作为我的起始值p,我得到的答案是非常不同的。为什么下面的工作不工作呢?Q1<-401# function: |F(Q1)+F(Q3)-1| so perfect fit should = 0 abs(pweibull(Q1,params[1],params[2])+pweib