#当前文件路径
getwd()
#设置当前路径，注意转译
setwd("C://Users//Administrator//Desktop//R_test")
#导入数据
data<-read.csv("data.csv")
1、平稳性检验
单位根检验
library(tseries)
adf.test(data$gov_day)
adf.test(data$med_day)
adf.test(data$pub_day)
P值均大于0.05，所以两者都没有通过单位根检验（都是不平稳的时间序列），因此对两者进行一阶差分，然后再进行单位根检验。
adf.test(diff(data$gov_day))
adf.test(diff(data$med_day))
adf.test(diff(data$pub_day))
都通过了平稳性检验
gov_day1 = diff(data$gov_day)
med_day1 = diff(data$med_day)
pub_day1 = diff(data$pub_day)
#判断政府和媒体的格兰因果关系
sr.reg = lm(gov_day1~med_day1)
error = residuals(sr.reg)
plot(error,main = "gov_day1--med_day1")
adf.test(error)
summary(sr.reg)
残差通过单位根检验
协整检验通过，也就是确认两者具有协整关系
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格兰杰变化的原因
原文：https://www.cnblogs.com/chunfang/p/13571813.html