急求：如何用R进行多元的极大似然估计？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

log(1-theta+theta/sum(prod(a+b/i)))+n*log(1-theta+theta/sum(prod(a+b/i))-1+theta)+sum((n-i+1)*(a+b/i))是对数似然函数，请问怎么模拟出抽取10各，100个，1000个数据时候对三个参数的估计？

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：极大似然估计 似然估计 极大似然 如何用 Theta

Sunzhuowei789 发表于 2018-4-1 17:25
log(1-theta+theta/sum(prod(a+b/i)))+n*log(1-theta+theta/sum(prod(a+b/i))-1+theta)+sum((n-i+1)*(a+b ...

用optim函数

中国梦丶 发表于 2018-4-2 22:07
用optim函数

n<-10;
LL<-function(params,data)
{
    theta<-params[1]
    a<-params[2]
    b<-params[3]
ll<-log(1-theta+theta/sum(prod(a+b/i)))+n*log(1-theta+theta/sum(prod(a+b/i))-1+theta)+sum((n-i+1)*(a+b/i))
return(-ll)
}
geyser.res<-nlminb(c(0.5,0.5,0.2),LL,data=rnbinom(n,2,0.5),lower=c(-2,-2,0),upper=c(2,2,1))
geyser.res$par

可以帮忙修改一下吗？

Sunzhuowei789 发表于 2018-4-5 15:27
n

i是什么