线性回归与In R中的群

Question

线性回归与In R中的群我想在R中使用lm()功能。我的数据是一个年度时间序列，其中一个字段用于年份(22年)，另一个字段用于州(50个州)。我想拟合每个状态的回归，这样在最后我有一个lm响应向量。我可以想象，对每个状态执行for循环，然后在循环中进行回归，并将每个回归的结果添加到向量中。然而，这似乎不太像R。在SAS中，我将执行‘by’语句，而在SQL中，我将执行‘GROUP BY’。做这件事的方式是什么？

Answer 1

下面是一种使用普雷一揽子：

d <- data.frame(
  state = rep(c('NY', 'CA'), 10),
  year = rep(1:10, 2),
  response= rnorm(20))library(plyr)# Break up d by state, then fit the specified model to each piece and# return a listmodels 
  <- dlply(d, "state", function(df) 
  lm(response ~ year, data = df))# Apply coef to each model and return a data frameldply(models, coef)# Print the summary of each 
  modell_ply(models, summary, .print = TRUE)