MLE估计遇见Missing Value的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

最近做Ender课本第三章第4节中，PPI的GARCH估计
用MLE估计Engle的衰减权重的ARCH（4）模型
自己写的程序如下：
series pi=log(ppi)-log(ppi(-1))
'set sample
sample s1 1960:1 1960:1
sample s2 1960:2 2002:1
smpl s2

' declare coef vectors to use in ARCH likelihood
coef(1) mu = .1
coef(1) alpha = .1
coef(2) beta = .1
coef(1) omega = .1
coef(1) theta=.1


' get starting values from ARMA(1,(1,4))
equation eq_arma.ls pi c pi(-1) ma(1) ma(4)
mu(1) = eq_arma.c(1)
alpha(1)=eq_arma.c(2)
beta(1)=eq_arma.c(3)
beta(2)=eq_arma.c(4)
omega(1)=eq_arma.@se^2

' set presample values of expressions in logl
smpl s1
series sig2=omega(1)
series res=0

' set up ARCH likelihood
logl ll1
ll1.append @logl logl
ll1.append res=pi-mu(1)-alpha(1)*pi(-1)-beta(1)*res(-1)-beta(2)*res(-4)
ll1.append sig2 = omega(1)+theta(1)*(0.4*res(-1)^2+0.3*res(-2)^2+0.2*res(-3)^2+0.1*res(-4)^2)
ll1.append z = res/@sqrt(sig2)
ll1.append logl =  log(@dnorm(z)) - log(sig2)/2  

' estimate and display results
smpl s2
ll1.ml(showopts, m=1000, c=1e-5)
show ll1.output

在运行后提示：Missing values in @logl series at current coefficients at observation 1960Q2 in "Do_LL1.ml"
百思不得其解
望各位能施以援手，不胜感激，谢谢

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：missing value SSIN Miss alue value MLE missing 遇见

参考过别人的问题后，已解决问题
主要是样本区间取值的问题
有MA(4)而忘记将初始区间取足够值，导致做ML估计的时候无法取初值
改正：
sample s1 1960Q1 1960Q4
sample s2 1961Q1 2002Q1

似懂非懂。。。。

你好，我也有个极大似然估计的问题想请教一下，我的qq是32962219，希望和你取得联系

MLE是什么估计？

你好  我能加你QQ吗  也有一些极大似然估计的问题想请教  QQ:393745164

为什么复制你的代码   做不出来呢  可以发一份完整版的吗   QQ 2972581670