Question

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

Answer 1

这是提取p-value的方法：

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

从plus函数中提取Breusch Godfrey测试的结果

1 个答案: