面板数据平稳性与单位根检验报错求助：punitroot参数问题

错误原因分析 #

你遇到的unused arguments (test = "llc", lags = "AIC")错误，核心是对plm包中punitroot函数的参数使用有误：

• test参数的选项中没有"llc"，LLC检验对应的参数值是"levinlin"
• lags参数不支持直接传入"AIC"，需要以列表形式指定滞后选择方法与最大滞后阶数
• 另外，punitroot要求输入的是面板数据格式（pdata.frame），而非普通数据框的列。

修正后的plm包实现 #

首先确保数据转换为面板格式，再调用正确参数：

# 加载plm包
library(plm)

# 将普通数据框转换为面板数据（替换index中的变量为你的个体/时间标识）
data_3_p <- pdata.frame(data_3, index = c("个体ID", "年份"))

# 运行LLC单位根检验（对应test="levinlin"）
llc_test_std <- punitroot(
  data_3_p$y, 
  test = "levinlin", 
  lags = list(method = "AIC", max = 5) # 用AIC选择滞后，最大阶数设为5，可按需调整
)

# 查看检验结果
summary(llc_test_std)

如果要做IPS检验，只需将test参数改为"im"即可：

ips_test_std <- punitroot(data_3_p$y, test = "im", lags = list(method = "AIC", max = 5))
summary(ips_test_std)

替代方案：使用urca包 #

urca包的ur.llc/ur.ips函数参数更贴近你最初的预期，支持直接传入lags="AIC"，但需要将数据转换为个体为列、时间为行的矩阵：

# 加载urca包
library(urca)

# 将面板数据转换为矩阵格式
y_matrix <- t(as.matrix(data_3_p$y))

# LLC检验（type可选择"none"/"drift"/"trend"，按需调整）
llc_test <- ur.llc(y_matrix, type = "drift", lags = "AIC")
summary(llc_test)

# IPS检验
ips_test <- ur.ips(y_matrix, type = "drift", lags = "AIC")
summary(ips_test)

其他注意事项 #

1. 数据平衡性：LLC、IPS检验要求数据为平衡面板，如果你的数据是不平衡面板，建议先整理为平衡面板，或改用punitroot中的"madwu"（Mad-Wu检验）或"hadri"（Hadri异方差稳健检验）。
2. R版本：你的R 4.2.3版本完全适配当前主流面板分析包（plm、urca），无需升级版本。
3. 标准化无用：平稳性检验不需要对数据做标准化，这一步不会解决参数错误，也不影响检验结果。

内容的提问来源于stack exchange，提问作者Mafalda Pinto