Извлечь AI C из FitARMA в R

Question

Я применяю FitARMA функцию из пакета FitARMA к определенному времени ser ie и получаю следующий результат:

> model <- FitARMA(ts, c(1,0,1))
> model
ARIMA(1,0,1)
length of series = 1593 ,  number of parameters = 3
loglikelihood = 5113 ,  aic = -10220 ,  bic =  -10203.9

Я хочу извлечь ai c к переменной. Однако в деталях модели нет ни ai c (экран с деталями модели ), ни информации о нем в документации к пакету.

Есть ли возможность сделать что-то вроде model_aic <- model$aic так как я хочу сделать для l oop для различных p, q порядков ARMA, поэтому я хотел бы извлечь ai c в переменную вместо того, чтобы набирать ее из консоли вручную?

Answer 1

Одним из способов является создание функции, которая вычисляет AI C модели FitARMA

library(FitARMA)
model <- FitARMA(AirPassengers, c(1,0,1))
model
ARIMA(1,0,1)
length of series = 144 ,  number of parameters = 3
loglikelihood = -496.55 ,  aic = 999.1 ,  bic =  1008

AICFitARMA <- function(model){
  k <- nrow(coef.FitARMA(model))
  AIC <- 2 * k - 2 * model$loglikelihood
  return(AIC)
}

AICFitARMA(model)
[1] 999.0944