请问R软件里如何得到ARIMA模型的BIC值

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位大侠，我刚刚学习使用R软件来处理时间序列数据，有诸多不懂，还请多多指导。谢谢！
请问R软件里如何得到ARIMA模型的BIC值。以下例为例，结果只有AIC而无BIC值？该如何得到BIC值呢？

> p=0
> q=1
> x=as.ts(gz19)
> d=0
> arima(x,order=c(p,d,q))
Call:
arima(x = x, order = c(p, d, q))
Coefficients:
          ma1  intercept
      -1.0000   106.9203
s.e.   0.2231     0.4263
sigma^2 estimated as 120.8:  log likelihood = -74.01,  aic = 154.01
>

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：ARIMA模型 ARIMA MA模型 Rim ima 模型 软件 ARIMA BIC

问题2：AICC值又该如何得到呢？有无已有函数？还是得自己编程？

问题3，eviews同R计算aic的公式有什么不同呢？
问题4，R软件里是否有时间序列专用的包（库），在哪里可以找到？哪里又有详细介绍包里有哪些东东如何使用呢？

#######arima
library(forecast)
Arima(lh, order = c(1,0,1))
Series: lh
ARIMA(1,0,1) with non-zero mean
Coefficients:
         ar1     ma1  intercept
      0.4522  0.1982     2.4101
s.e.  0.1769  0.1705     0.1358
sigma^2 estimated as 0.1923:  log likelihood=-28.76
AIC=65.52   AICc=66.45   BIC=73.01
#####
ps:
aic=65.52
npar=4
nstar=48
bic = aic + npar*(log(nstar) - 2)
bic  #   73.0048
#####lm model
     可以用公式算,也可以用既有 AIC, BIC function
AIC = n + n * log(2 * pi) + n * log(rss/n) + 2 * (P+1)
BIC = n + n * log(2 * pi) + n * log(rss/n) + log(n) * (p+1)

lm <- lm(Fertility ~ . , data = swiss)
n <- nrow(swiss)
rss <- sum(resid(lm) ^2)

#calculate AIC
n + n*log(2*pi) + n * log(rss / n ) + 2 * (6+1)
#[326.0716]
#call AIC function
AIC(lm)
#[326.0716]

#calculate BIC
n + n * log(2*pi) + n*log(rss/n) +  log(n) * (6+1)
#[339.0226]
#call AIC function to calculate BIC
AIC(lm, k=log(n))
#[339.0226]
#call BIC function from "nlme" package
library(nlme)
BIC(lm)
#[339.0226]

问题3  eviews同R计算aic的公式,只是小不同而已.

^_^，谢谢LS.

Pi是什么意思？

Pi什么意思????

pi = π = 3.141593

我想知道 p q d 是如何确定的？谢谢解答！！！！

看看，留个脚印先！