如何同时从grangertest（）函数从该对象提取F和Pvalue？

Question

Model 1: DLPFc ~ Lags(DLPFc, 1:1) + Lags(IPs, 1:1)
Model 2: DLPFc ~ Lags(DLPFc, 1:1)
  Res.Df Df      F    Pr(>F)    
1   4597                        
2   4598 -1 88.522 < 2.2e-16 ***

Answer 1

library(lmtest)
data(ChickEgg)
(out <- grangertest(egg ~ chicken, order = 3, data = ChickEgg))
#Granger causality test

#Model 1: egg ~ Lags(egg, 1:3) + Lags(chicken, 1:3)
#Model 2: egg ~ Lags(egg, 1:3)
#  Res.Df Df      F Pr(>F)
#1     44                 
#2     47 -3 0.5916 0.6238

提取F和p值

out[, c("F", "Pr(>F)")]
#       F Pr(>F)
#1              
#2 0.5916 0.6238

如果要排除第一行空白，请执行out[-1, c("F", "Pr(>F)")]

---

如果我们不需要新对象并立即获取值，我们可以

unlist(grangertest(egg ~ chicken, order = 3, data = ChickEgg)[-1, c("F", "Pr(>F)")])
#        F    Pr(>F) 
#0.5916153 0.6237862

哪个返回命名向量。