楼主: xiazaihao1990
1095 9

[问答] (10论坛币求助)极大似然估计,我的Rcodes哪一行有问题呀 [推广有奖]

  • 0关注
  • 0粉丝

本科生

5%

还不是VIP/贵宾

-

威望
0
论坛币
819 个
通用积分
1.0030
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
607 点
帖子
36
精华
0
在线时间
64 小时
注册时间
2017-12-30
最后登录
2022-5-28

相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币

似然方程如图所示: (参数只有一个是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 |只看作者 |坛友微信交流群
  1. x_vec <- c(3,1,6,8,7)

  2. ll <- function(k,x){
  3.   - sum( log( (k^x)*exp(-k) / sapply(x, factorial) ) )
  4. }

  5. optimize(ll,c(0,120),maximum=FALSE,tol=0.0001,x=x_vec)
复制代码


使用道具

板凳
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的取值范围。

使用道具

7
rosenbloog 发表于 2020-4-3 08:25:21 |只看作者 |坛友微信交流群
另外,建议参考这篇:https://bbs.pinggu.org/thread-5758285-1-1.html

  1. ll <- function(k, dat=x_vec){
  2.   loglik = sapply(k, function(x) -sum(log((x^x_vec)*exp(-x) / sapply(x_vec, factorial))))
  3.   return(loglik)
  4. }
  5. curve(ll, 1, 15)
复制代码


已有 1 人评分学术水平 热心指数 收起 理由
xiazaihao1990 + 1 + 1 热心帮助其他会员

总评分: 学术水平 + 1  热心指数 + 1   查看全部评分

使用道具

8
xiazaihao1990 学生认证  发表于 2020-4-3 09:05:08 |只看作者 |坛友微信交流群
rosenbloog 发表于 2020-4-3 08:25
另外,建议参考这篇:https://bbs.pinggu.org/thread-5758285-1-1.html
感谢大哥

使用道具

9
mj谢 发表于 2020-4-3 10:50:44 |只看作者 |坛友微信交流群
xiazaihao1990 发表于 2020-4-3 01:55
老哥请问

c(-100,100),maximum=TRUE 这俩在这里啥作用啊
应该你解决了∠( ᐛ 」∠)_
第一个是参数范围,若无要求,尽量取范围大些。第二个表示取目标函数的极大值,如果原函数中取了负号,那么改成FALSE取极小值。
已有 1 人评分学术水平 收起 理由
xiazaihao1990 + 1 精彩帖子

总评分: 学术水平 + 1   查看全部评分

使用道具

10
xiazaihao1990 学生认证  发表于 2020-4-3 19:50:13 |只看作者 |坛友微信交流群
mj谢 发表于 2020-4-3 10:50
应该你解决了∠( ᐛ 」∠)_
第一个是参数范围,若无要求,尽量取范围大些。第二个表示取目标函数的 ...
嗯嗯,这个解释看懂了, 跪谢

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-25 07:05