我有x和y值以下,如你所见,x大部分都是负数,基本上我只有我观察数据的PDF左侧。
我必须将其与学生分布相匹配,并找出自由度和比例参数。
问题是,估计的分布将具有非常小的方差(即小规模参数)。所以当我使用下面的方法来拟合分布时,无论我设置什么初始值,nls都不会收敛。
我在下面的代码中使用了额外的参数c,因为我使用它来扩展分布:dt(x/a,df)
。因此,为了保存概率,我不可避免地要对输出进行定时但是不变。我相信这个额外的参数导致收敛性差,但我不知道如何更好地适应分布。
我已经找了配件包,但是这些包需要完整的分发,而我只有它的左侧。
x y
1 -0.0050 0.000000
2 -0.0045 26.723019
3 -0.0040 28.557704
4 -0.0035 41.085068
5 -0.0030 66.258445
6 -0.0025 81.129807
7 -0.0020 83.751611
8 -0.0015 130.378353
9 -0.0010 157.806018
10 -0.0005 201.505657
11 0.0000 949.650354
12 0.0005 193.721270
dat<-data.frame(x=x,y=y)
res<-nls( y~(dt(x/a,df)*c), dat,
start=list(a=0.000201, df=0.9, c=2104), trace = TRUE)