1

我正在尝试使用最大似然和模拟数据来估计 GARCH (1,1) 模型。这就是我得到的：

library(fGarch)
set.seed(1)
garch11<-garchSpec(model = list())
x<-garchSim(garch11, n = 1000)
y <- t(x)
r <- y[1, ]

### Calculate Residuals
CalcResiduals <- function(theta, r)
{
  n <- length(r)
  omega<-theta[1]
  alpha11<-theta[2]
  beta11<-theta[3]
  sigma.sqs <- vector(length = n)
  sigma.sqs[1] <- 0.02
  for (i in 1:(n-1)){
    sigma.sqs[i+1] <- omega + alpha11*(r[i]^2) + beta11*sigma.sqs[i]
  }
  return(list(et=r, ht=sigma.sqs))
}

###Calculate the log-likelihood

GarchLogl <- function(theta, r){
  res <- CalcResiduals(theta,r)
  sigma.sqs <- res$ht
  r <- res$et
  return(-sum(dnorm(r[-1], mean = 0, sd = sqrt(sigma.sqs[-1]), log = TRUE)))
}

fit2 <- nlm(GarchLogl, # function call 
            p = rep(1,3), # initial values = 1 for all parameters
            hessian = FALSE, # also return the hessian matrix
            r = r , # data to be used
            iterlim = 500) # maximum iteration

不幸的是，我收到以下错误消息并且没有结果

有 50 个或更多警告（使用 warnings() 查看前 50 个） 1：在 sqrt(sigma.sqs[-1]) 中：产生了 NaN

2: 在 nlm(GarchLogl, p = rep(1, 3), hessian = FALSE, data <- r, ... : NA/Inf durch größte positive Zahl ersetzt

你知道我的代码有什么问题吗？非常感谢！

4

0 回答 0