14/02/25 17:06:50.33 .net
質問です。
optim関数を実行した際に、「初期パラメータで関数を表示出来ません」とエラーが出てしまいます。
関数の式は間違ってないと思うのですが、原因がわからないため、ご教授お願いいたします。
コードはこちらです。
sai <- function(x)
{
T<-length(x); q<-3; eps<-diff(log(x))
ts.plot(eps)
Mean <- mean(x); Var <- var(x); S = 1e-6
parm <- c(0.1,0.1*Var )
archsai<-function(par){
for( t in 2:T){
alpha<-c(par[1])
omega<-par[2]
if(t<=1){
htt<-omega
}else {
htt <- (omega+alpha[1]*eps[t-1]^2)
ht <- sqrt(htt)
A <- -(log(htt))^2*(1/2)
B <- -(((eps[t])^2)/htt)*(1/2)
ln <- A+B
}
}
return(ln)
}
optim(parm,archsai,control = list(fnscale=-1) )
}