ZikV 发表于 2023-4-1 18:55 
谢谢楼主,有没有相关的代码可以分享下呢?邮箱:
install.packages("mgc")
library(mgc)
# 定义数据
n <- 100
x <- rnorm(n)
y <- 0.5*x^2 + rnorm(n)
# 计算格兰杰因果指数
gc_res <- mgc::mgc(x, y, nperm = 1000, family = "x2")
gc_res$tstat # TVAL值
gc_res$observed.Tn # Tn值
其中,nperm 表示置换检验的次数,family 表示检验所采用的分布,这里选择了 x2 分布