r - 在 R w/PLYR 中提取组回归系数

Question

我正在尝试对数据集中的每个邮政编码运行回归并将系数保存到数据框中，但我遇到了问题。

每当我运行下面的代码时，我都会得到一个名为“系数”的数据框，其中包含每个邮政编码，但每个邮政编码的截距和系数等于简单回归 lm(Sealed$hhincome ~密封 $square_footage)。

当我按照下面链接中 Ranmath 的示例中所示运行代码时，一切都按预期工作。在使用 STATA 多年后，我是 R 的新手，所以非常感谢任何帮助:)

R extract regression coefficients from multiply regression via lapply command

library(plyr)
Sealed <- read.csv("~/Desktop/SEALED.csv")

x <- function(df) {
      lm(Sealed$hhincome ~ Sealed$square_footage)
}

regressions <- dlply(Sealed, .(Sealed$zipcode), x)
coefficients <- ldply(regressions, coef)

Answer 1

因为 dlply 接受一个 ... 参数，允许将额外的参数传递给函数，您可以使事情变得更简单:

dlply(Sealed,.(zipcode),lm,formula=hhincome~square_footage)

lm 的前两个参数是formula 和data。由于此处指定了 formula，lm 将选取下一个给定的参数(Sealed 的相关邮政编码特定 block )作为 数据参数...