关于极大似然估计参数的问题： 无法用初始参数来评估函数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

da=read.csv("C:\\Users\\Administrator\\Desktop\\论文\\数据\\INTC英特尔-table.csv")
KK=function(theta)
{
x=theta[1]
y=theta[2]
z=theta[3]
option1=0
for(i in 1:9130)
{
res=(1+x[i]*(da[i,4]-z[i])/y[i])^(-1/x[i])
option1=option1+res
}
option2=1
q=1
while(q<9130)
{
if(da[q,4]>=0.1)
{
con=((1+x[q]*(da[q,4]-z[q])/y[q])^((-1-x[q])/x[q]))/y[q]
option2=option2*con
}
q=q+1
}
return(-option1*option2)
}
optim(par=0,KK)

Error in optim(par = 0, KK) : 无法用初始参数来评估函数
此外: Warning message:
In optim(par = 0, KK) : 用Nelder-Mead方法来算一维优化不很可靠：
用"Brent"或直接用optimize()

问题：需要估计的参数为x[1:9130],y[1:9130],z[1:9130],da[,4]是对数收益率序列。这样的话optim里的初始值应该取多少呢？
报错里的无法用初始参数来评估函数应该怎么修改。。。
另外，function里面用
x=theta[1]
y=theta[2]
z=theta[3]
定义三个待估参数对吗？
求大神拯救

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：极大似然估计 似然估计 极大似然 function OPTIMIZE 评估

KK函数里的参数是三个变量，在optim里给的初始值也要是三个，比如
optim(c(0, 1, 1), KK)