首頁(yè) 猿問(wèn) 線性回歸與In R中的群

線性回歸與In R中的群

R語(yǔ)言

慕碼人8056858 2019-06-21 15:24:35

線性回歸與In R中的群我想在R中使用lm()功能。我的數(shù)據(jù)是一個(gè)年度時(shí)間序列，其中一個(gè)字段用于年份(22年)，另一個(gè)字段用于州(50個(gè)州)。我想擬合每個(gè)狀態(tài)的回歸，這樣在最后我有一個(gè)lm響應(yīng)向量。我可以想象，對(duì)每個(gè)狀態(tài)執(zhí)行for循環(huán)，然后在循環(huán)中進(jìn)行回歸，并將每個(gè)回歸的結(jié)果添加到向量中。然而，這似乎不太像R。在SAS中，我將執(zhí)行‘by’語(yǔ)句，而在SQL中，我將執(zhí)行‘GROUP BY’。做這件事的方式是什么？

查看完整描述

3 回答

慕容森

TA貢獻(xiàn)1853條經(jīng)驗(yàn) 獲得超18個(gè)贊

下面是一種使用普雷一攬子：

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)

反對(duì) 回復(fù) 2019-06-21