全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2661 3
2016-05-29
LZ是 R语言新手...以前没接触过编程,所以都是边查边写,好多地方不知道对不对..求各位大神指点!
要求4个参数的极大似然值,这是似然方程:

捕获.PNG


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

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

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


附件列表
捕获.PNG

原图尺寸 77.47 KB

捕获.PNG

捕获.PNG

原图尺寸 26.17 KB

捕获.PNG

捕获.PNG

原图尺寸 8.34 KB

捕获.PNG

二维码

扫码加我 拉你入群

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

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

全部回复
2016-5-29 18:25:17
你的param是list还是vector? 如果是list试试:m=param[[1]]
二维码

扫码加我 拉你入群

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

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

2016-5-31 16:33:27
bbslover 发表于 2016-5-29 18:25
你的param是list还是vector? 如果是list试试:m=param[[1]]
是vector....不过问题解决了,直接写的方程而不是用param这个方式来写的...
二维码

扫码加我 拉你入群

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

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

2017-9-7 17:34:02
jsjyyjc 发表于 2016-5-31 16:33
是vector....不过问题解决了,直接写的方程而不是用param这个方式来写的...
请问你是怎么解决的?我也遇到了同样的问题
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群