发帖
雷达卡
入门级别计量课的上机作业...
自行编写语句实现Arch（3）
关于其似然函数部分：
logLikFun=function(param)
{
t=length(y)
c=param[1]
a1=param[2]
a2=param[3]
a3=param[4]
x1<-y[1:(length(y)-3)]
xx1<-x1^2
x2<-y[2:(length(y)-2)]
xx2<-x2^2
x3<-y[3:(length(y)-1)]
xx3<-x3^2
sigma2= c+a1*xx1+a2*xx2+a3*xx3
loglik=-(t-3)*log(2*pi)/2-sum((((y[-1:-3])^2))/sigma2)-sum(log(sigma2))/2
return(loglik)
}
相当于默认y1,y2,y3给定，求y4...yn的条件分布 直接认为这是log-likelihood
但是跑出来结果跟软件结果不一样
求问：
是否y1,y2,y3也要算在分布内？
采用什么分布？
(即其中sigma1，sigma2，sigma3怎么取？？)
（上课讲ar的时候前3项是直接取期望得到的分布的参数，但这个地方直接期望就是默认为0了吧）
默认假设为最基本的
应该是正态吧
京公网安备 11010802022788号
