用尺度参数拟合学生分布

时间:2013-05-31 00:55:08

标签: r

我有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)

0 个答案:

没有答案