入门级别计量课的上机作业...
自行编写语句实现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了吧）


默认假设为最基本的
应该是正态吧