求大神帮助，怎么求参数的初值！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

小白刚接触R不久，要做一个用极大似然估计求参数的作业，在用nlminb时候发现参数的初值不会求。。。求大神指点！顺便请指点下错误之处！谢谢啦！作业是这个样子的：

[url=]entropy-hw1.pdf[/url]
entropy-hw1.pdf (153.98 KB)

2016-4-8 17:08:59 上传

数据：[url=] 1.txt (377 Bytes)

2016-4-8 17:30:31 上传

[/url]

代码：
> ut=read.table("1.txt")

> gg=log(ut)
> View(gg)
> acf(gg)
> pacf(gg)
> Box.test(gg,lag=15,type="Ljung")
> library("fGarch")

> m1=garchFit(~garch(2,0),data=gg,trace=F)
> summary(m1)

> m2=garchFit(~garch(1,0),data=gg,trace=F)
> summary(m2)

> plot(m1)
> plot(m2)

> bt=var(gg)



> LL<-function(params,ut){
+ f<-exp(-(ut)*(exp(params[1]+params[2]*abs(log((ut)/(bt)))+params[3]*log((ut)/(bt))+params[4]*log(bt)))^(-2))
+ ll<-sum(log(f1))
+ return(-log(ll))
+ }

mp<-optim(par=c(),LL,data=ut)
求问四个的参数的初值怎么求呢？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：求大神帮助 garchFit Summary Entropy Library library

初始值不是自己设的吗？

jinkelazzz 发表于 2016-4-8 19:42
初始值不是自己设的吗？

我也纳闷啊...自己设置的初始值带进去就提示出错

Error in objective(.par, ...) :
  unused argument (data = list(V1 = c(1.01157, 0.95575, 0.88368, 0.89658, 1.04106, 1.08622, 0.83943, 0.96935, 1.01142, 1.04429, 1.02087, 0.90694, 1.09974, 1.04155, 0.93418, 0.96656, 1.03199, 1.04394, 0.96934, 1.02265, 1.06791, 1.09918, 0.99796, 1.04255, 1.13939, 0.97173, 0.96149, 0.99774, 0.77768, 1.04982, 0.84358, 0.81743, 1.02211, 0.98401, 0.99872)))