R语言_格兰因果检验

R语⾔_格兰因果检验

#当前⽂件路径

getwd()

#设置当前路径，注意转译

setwd("C://Users//Administrator//Desktop//R_test")

#导⼊数据

data<-read.csv("data.csv")

1、平稳性检验

library(tseries)

刀豆蛋白

计量st(data$med_day)分频器设计

P值均⼤于0.05，所以两者都没有通过单位根检验（都是不平稳的时间序列），因此对两者进⾏⼀阶差分，然后再进⾏单位根检验。st(diff(data$gov_day))

双性同体

都通过了平稳性检验

gov_day1 = diff(data$gov_day)

med_day1 = diff(data$med_day)

pub_day1 = diff(data$pub_day)

#判断政府和媒体的格兰因果关系

< = lm(gov_day1~med_day1)

error = )

自然基金plot(error,main = "gov_day1--med_day1")

用友华表cell插件)

残差通过单位根检验

协整检验通过，也就是确认两者具有协整关系

grangertest(gov_day1~med_day1,order=2)

#P值⼩于0.05，通过检验，也就是拒绝了gov_day1不是引起med_day1格兰杰变化的原因即gov_day1是引起med_day1格兰杰变化的原因grangertest(med_day1~gov_day1,order=2)

#P值⼤于0.05，没有通过检验，也就是接受了med_day1不是引起gov_day1格兰杰变化的原因即med_day1不是引起gov_day1格兰杰变化的原因

本文发布于:2024-09-20 23:36:20，感谢您对本站的认可！

标签：检验单位根媒体

留言与评论（共有 0 条评论）