在 plm 函数中从 Breusch Godfrey 测试中提取结果答案

【问题标题】：Extract results from Breusch Godfrey Test in a plm function在 plm 函数中从 Breusch Godfrey 测试中提取结果
【发布时间】：2024-01-06 23:00:02
【问题描述】：

我需要创建一个表，其中行中包含面板数据回归的 id，列中包含模型的统计信息。我正在尝试从 Breusch Godfrey 测试中提取结果（统计、df、p 值）以插入表中。一个评论是我正在使用相同的命令运行许多回归。我试过做

GodfreyallModelsResultsF <- lapply(allModelsList, function(x) 
  bgtest(plm(x, data=data, model="within")))

GodfreyallModelsResultsF2 <- as.data.frame(
  matrix(unlist(GodfreyallModelsResultsF,use.names=TRUE), 
         nrow=length(unlist(GodfreyallModelsResultsF[1]))))

GodfreyallModelsResultsF <- t(GodfreyallModelsResultsF)
colnames(GodfreyallModelsResultsF) <- c("BPTest","df","name","Pvalue","code")

但它没有用。最后一点是，我能够使用此代码成功地从 Breusch Pagan 测试中提取。有人可以帮我解决 B-G 问题吗？

【问题讨论】：

请问dput(head(your_data))能否获得可重现的数据样本？
我知道这会容易得多，但很遗憾，由于我工作的地方有保密协议，我无法提供。对不起

标签： r regression panel-data

【解决方案1】：

这里是p-value的提取方式：

model <- plm(x, data = data, model = "within")
as.numeric(pbgtest(model)[4])

【讨论】：