提取plm函数中Breusch-Godfrey测试的结果

我需要创建一个表，其中行中有面板数据回归的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问题吗？

以下是提取p值的方法：

---

是否有可能获得一个可复制的数据样本，其中包含dputheadyour_数据？我知道这会容易得多，但不幸的是，由于我工作的地方有保密协议，我无法提供。很抱歉

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