全部版块 我的主页
› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛
请问大神，我这个函数编的有问题吗？
了空不了色
1448 9
2017-05-31
这是函数表达式，Y服从正态分布
这是我的程序，
eta<-0.01
sigma<-0.02
mu<-0.2
lamda<-300
n<-15
logl<-function(params){#eta,sigma,mu
f<--n/2*log(2*pi)-n*log(params[2])-n/2*log(params[1]+lamda)-
y^2/(2*params[2]^2*(params[1]+lamda))+lamda*params[3]*y/(params[2]^2*(params[1]+lamda))-
params[3]^2/(2*params[2]^2*(params[1]+lamda))*lamda^2
return(-f)
}
y<-seq(0,by=0,length=n)
for (i in 1:n) {
y<-rnorm(1,mu*lamda,sigma*sqrt(eta+lamda))
}
geyser.res<-nlminb(c(0.0001,0.0001,0.0001),logl,hessian=TURE,
lower=c(1e-20,1e-10,1e-10),
upper=c(Inf,Inf,Inf))
虽然这条也可以运行，但是eta总是估不出来，运行的话会触碰到它的下界，1e-20，然后就以1e-20
作为结果。但是Mu和sigma可以估出来 QQ图片20170531115900.jpg
二维码

扫码加我 拉你入群

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

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

全部回复
foozhencheng
2017-5-31 13:03:40
原问题是什么？
二维码

扫码加我 拉你入群

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

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

了空不了色
2017-5-31 15:33:37
foozhencheng 发表于 2017-5-31 13:03
原问题是什么？
原来就是要求图中那个函数的最大似然估计
二维码

扫码加我 拉你入群

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

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

foozhencheng
2017-6-2 10:20:28
之前手机一直没看清楚图，用电脑才看出来：模型有问题，不是代码的问题。期望和方差足以确定一个正态分布，第三个参数定不出来，属于identifiability problem～
二维码

扫码加我 拉你入群

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

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

了空不了色
2017-6-2 17:45:00
foozhencheng 发表于 2017-6-2 10:20
之前手机一直没看清楚图，用电脑才看出来：模型有问题，不是代码的问题。期望和方差足以确定一个正态分布， ...
但第三个参数就是方差的一部分啊？我 是eta这个参数定不出来
二维码

扫码加我 拉你入群

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

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

foozhencheng
2017-6-3 02:23:37
了空不了色 发表于 2017-6-2 17:45
但第三个参数就是方差的一部分啊？我 是eta这个参数定不出来
因为σ^2已经足以表达方差了，所以是初始建模的问题～
二维码

扫码加我 拉你入群

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

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

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

分享

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