【R言語】統計解析フリーソフトR 第5章【GNU R】
■ このスレッドは過去ログ倉庫に格納されています
0234132人目の素数さん
2014/02/22(土) 14:58:46.83X<-c(1,2,3,4,5,6,7,8,9,12,13,24)
T<-length(X) #データ数
q<-3 #次数
eps<-X #diff(log(X))#収益率
ty<-function(alpha,omega){
for( t in 1:T){
alpha[1]
alpha[2]
alpha[3]
omega
for( j in 1:q ){
AA <- alpha[j]*eps[t+3-j]*eps[t+3-j] + AA
A<-(log(omega+AA)^2)*(-1/2)
print(eps[abs(t+3-j)])
}
}
}
optimize(c(2,1,2,3),ty, maximum=TRUE)
■ このスレッドは過去ログ倉庫に格納されています