全部版块 我的主页
› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛
（10论坛币求助）极大似然估计，我的Rcodes哪一行有问题呀
xiazaihao1990
1341 9
2020-04-02

似然方程如图所示： （参数只有一个是k,  xi是样本观察值（xi有五个观察值3,1,6,8,7））
11111.jpg

以下是我的R Codes, 但运行不出来 , 10个论坛币求答案，哭着感谢

x<-c(3,1,6,8,7)

likelifn<- function(x, k){

  return(prod((k^x)*exp(-k)/prod(x)))

}  

    loglikelifn<-function(x,par){

    return(-sum(log(likelifn(x, k = par[1]))) )  

    }


optimize(par, loglikelifn, x)


二维码

扫码加我 拉你入群

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

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

全部回复
mj谢
2020-4-2 20:19:33
f=function(k){
x=c(3,1,6,8,7)
return(prod(k^x)*exp(-5*k)/prod(x))
}
optimize(f,c(-100,100),maximum=TRUE,tol=0.0001)
你发的代码第三行漏了*5。第二个函数好像没有必要，或者直接第一个函数改成log形式的。
二维码

扫码加我 拉你入群

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

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

rosenbloog
2020-4-2 23:48:25
复制代码


二维码

扫码加我 拉你入群

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

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

xiazaihao1990
2020-4-2 23:57:05
rosenbloog 发表于 2020-4-2 23:48
请问老哥， 最后一行的c(0,120)起的啥作用啊
二维码

扫码加我 拉你入群

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

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

xiazaihao1990
2020-4-3 01:55:32
mj谢 发表于 2020-4-2 20:19
f=function(k){
x=c(3,1,6,8,7)
return(prod(k^x)*exp(-5*k)/prod(x))
老哥请问

c(-100,100),maximum=TRUE 这俩在这里啥作用啊
二维码

扫码加我 拉你入群

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

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

rosenbloog
2020-4-3 07:54:42
xiazaihao1990 发表于 2020-4-2 23:57
请问老哥， 最后一行的c(0,120)起的啥作用啊
自己看帮助？这是参数k的取值范围。
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
栏目导航
热门文章
推荐文章

说点什么

分享

微信QQ空间QQ微博
扫码加好友，拉您进群
各岗位、行业、专业交流群