求大神帮助，在用R做极大似然时出现了Error in objective(.par, ...) :

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

LZ是 R语言新手...以前没接触过编程，所以都是边查边写，好多地方不知道对不对..求各位大神指点！
要求4个参数的极大似然值，这是似然方程：

我编的代码是这样的，其实并没有信心感觉他是对的

> LL<-function(params,ut){
+ m=param[1]
+ n=param[2]
+ l=param[3]
+ k=param[4]
+ f<-exp(-ut/(gamma(1)*exp(2*(param[1]+(param[2]+param[3])*log(ut))/(1+param[2]+param[3]-param[4]))))
+ ll=sum(log(f))
+ return(-log(ll))
+ }

然后在用nlminb函数的时候出现了以下的情况：
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)))

求各位指点问题出现在那里！谢谢谢谢！万分感谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Objective object 求大神帮助 Error 极大似然 objective Error

你的param是list还是vector? 如果是list试试：m=param[[1]]

bbslover 发表于 2016-5-29 18:25
你的param是list还是vector? 如果是list试试：m=param[[1]]

是vector....不过问题解决了，直接写的方程而不是用param这个方式来写的...

jsjyyjc 发表于 2016-5-31 16:33
是vector....不过问题解决了，直接写的方程而不是用param这个方式来写的...

请问你是怎么解决的？我也遇到了同样的问题