BIC信息准则的R语言代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

BIC信息准则的R语言代码怎么搞，我想用BIC信息准则给ARMA模型定阶，请完整讲解，非常感谢

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：BIC R语言 ARMA模型定阶 arma模型 非常感谢 R语言

求解答啊> a<-read.table("D:/数据/江苏省教育经费.csv",sep=",",header=T)
> x<-ts(a$kilimeter,start=2013)
Error in ts(a$kilimeter, start = 2013) :
  'ts'对象至少必需有一个或多个观察量

※※这一步我总是做不对呀，$这个符号后面加东西就根本不出来东西
> x<-ts(a,start=2013)
> plot(x)
※像这样后面什么都没有是可以画图，但是图都是分开的，怎么合到一起去呢，急求，写论文需要，谢谢了

看你的描述，问题很可能是出在你读取数据上面，你可以用RStudio的工具去读取，同时检查是不是有NA值。

楼主你好，不知道你用的是R里面的哪个程序包，如果参数估计之后，你尝试去提取计算对象，但是没有看到bic，那可能就是这个程序包没有提供方法给你计算这个信息准则。比如TSA这个包对于ARMA参数估计后，就是只提供aic的结果，因为它只计算aic。那就得自己尝试写bic的计算代码，不过这个并不难，我把两个信息准则AIC和BIC的源码附在下面，like是你的loglikelihood极大似然估计值，Total是你参数估计ARMA模型时的观测值总数目，n是你算ARMA时需要进行参数估计的变量个数，代进去下面的代码进行计算，得到的就是BIC的结果。数学公式参照自R语言程序包rugarch程序说明论文的第28页，原作者是Alexios Ghalanos，网上有些信息准则的回答会有出入，以里面的公式为准即可。说明论文也附在下面。
# AIC (Akaike Information Criterion)
Index_AIC <- (-2*llike)/Total + (2*n)/Total

# BIC (Bayesian Information Criterion)
Index_BIC <- (-2*llike)/Total + (n*log(Total))/Total