请问大神，我这个函数编的有问题吗？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

这是函数表达式，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可以估出来

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：正态分布 函数表 表达式 ETA

原问题是什么？

foozhencheng 发表于 2017-5-31 13:03
原问题是什么？

原来就是要求图中那个函数的最大似然估计

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

foozhencheng 发表于 2017-6-2 10:20
之前手机一直没看清楚图，用电脑才看出来：模型有问题，不是代码的问题。期望和方差足以确定一个正态分布， ...

但第三个参数就是方差的一部分啊？我 是eta这个参数定不出来

了空不了色 发表于 2017-6-2 17:45
但第三个参数就是方差的一部分啊？我 是eta这个参数定不出来

因为σ^2已经足以表达方差了，所以是初始建模的问题～

foozhencheng 发表于 2017-6-3 02:23
因为σ^2已经足以表达方差了，所以是初始建模的问题～

不是很明白，为什么σ^2足以表达方差，可以解释一下吗

foozhencheng 发表于 2017-6-3 02:23
因为σ^2已经足以表达方差了，所以是初始建模的问题～

是因为整个方差可以用某个参数的平方来表示，所以σ^2足够了可以包含后面eta+300这一部分吗

了空不了色 发表于 2017-6-3 08:57
是因为整个方差可以用某个参数的平方来表示，所以σ^2足够了可以包含后面eta+300这一部分吗

整个分布就两个参数，均值和方差，因此你未知参数中只可能确定两个。

foozhencheng 发表于 2017-6-3 09:25
整个分布就两个参数，均值和方差，因此你未知参数中只可能确定两个。

谢谢，我懂你的意思了