请选择 进入手机版 | 继续访问电脑版
楼主: moonalone
780 2

[问答] 在编写极大似然法时出现的一个问题 [推广有奖]

  • 1关注
  • 0粉丝

本科生

20%

还不是VIP/贵宾

-

威望
0
论坛币
208 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
1058 点
帖子
19
精华
0
在线时间
140 小时
注册时间
2014-7-9
最后登录
2019-12-14

moonalone 学生认证  发表于 2017-6-22 16:50:26 |显示全部楼层 |坛友微信交流群

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

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

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
本人想先用标准正态分布函数,模拟出一些值。
然后把这些值放入正态分布的函数中 用极大似然法去估计这些值的期望和方差,最后看看是不是标准正态的。
可是代码写出来之后,发现估计出来的数值不对,望大神帮帮忙 看看问题出现在哪里
a<- rnorm(100,0,1)
LL<- function(params,data){
f<- 100*log(params[2])-sum(((data-params[1])^2)/(2*(params[2])^2))
ll<- f
return(-ll)
}
geyser.res<-nlminb(c(0.1,0.1),LL,data=a)
geyser.res$par
结果:
> a<- rnorm(100,0,1)
> LL<- function(params,data){
+ f<- 100*log(params[2])-sum(((data-params[1])^2)/(2*(params[2])^2))
+ ll<- f
+ return(-ll)
+ }
> geyser.res<-nlminb(c(0.1,0.1),LL,data=a)
> geyser.res$par
[1] -11228394205 2435854546936


二维码

扫码加我 拉你入群

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

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

关键词:极大似然法 极大似然 标准正态分布 正态分布函数 正态分布 极大似然法 编程

qoiqpwqr 发表于 2017-6-24 20:29:23 |显示全部楼层 |坛友微信交流群
f <- 那行的100前面少了个负号

使用道具

moonalone 学生认证  发表于 2017-6-25 12:26:58 |显示全部楼层 |坛友微信交流群
qoiqpwqr 发表于 2017-6-24 20:29
f
感谢感谢 我查了好久

使用道具

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

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

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

GMT+8, 2024-3-29 07:05