全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
4804 3
2021-03-23
BIC信息准则的R语言代码怎么搞,我想用BIC信息准则给ARMA模型定阶,请完整讲解,非常感谢
二维码

扫码加我 拉你入群

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

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

全部回复
2021-3-23 13:39:30
求解答啊> 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)
※像这样后面什么都没有是可以画图,但是图都是分开的,怎么合到一起去呢,急求,写论文需要,谢谢了
二维码

扫码加我 拉你入群

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

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

2021-3-23 23:54:22
看你的描述,问题很可能是出在你读取数据上面,你可以用RStudio的工具去读取,同时检查是不是有NA值。
二维码

扫码加我 拉你入群

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

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

2021-3-26 19:24:01
楼主你好,不知道你用的是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

附件列表
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群