将字符串替换为{plm}模型

Question

本教程展示了如何将字符串替换为lm模型调用。我想做同样的事情，但是使用plm进行面板回归。目标是建立一个参数组合的网格，然后使用一个apply调用运行几个回归，每次在适当的变量中替换。

library(plm)    
#make grid of parameters to loop over; select just the first row to reproduce the error
parameters <- expand.grid(dv="mpg", x1=c("cyl", "gear"), x2=c("hp", "wt"), stringsAsFactors = FALSE)
row <- parameters[1,]

#lm works
lm(substitute(y ~ x*z, list(y=as.name(row[[1]]), x=as.name(row[[2]]), z=as.name(row[[3]]))), data=mtcars)

#plm does not
plm(substitute(y ~ x*z, list(y=as.name(row[[1]]), x=as.name(row[[2]]), z=as.name(row[[3]]))), data=mtcars)

这产生了：

错误:inherits(对象，“公式”)不是真

解决方案可能包括将变量名粘贴到plm中，而不是使用substitute。或者我可以避免为此目的使用plm吗？

Answer 1

plm需要面板数据。您可以使用来自Produc的?plm示例数据集。

library(plm)
parameters <- expand.grid(dv="pcap", 
                          x1=c("hwy", "water"), 
                          x2=c("util", "pc"), 
                          stringsAsFactors = FALSE)
row <- parameters[1,]

# I assign the formula to another object
# to avoid clutter in calling plm()
formula.plm <- substitute(y ~ x  * z, 
                   list(y = as.name(row[[1]]), 
                        x = as.name(row[[2]]), 
                        z = as.name(row[[3]]))) 

plm(eval(formula.plm), data = Produc)