r - 使用 Lattice(或其他东西)在 R 中绘制 lme4 的回归结果

Question

感谢 previous answer，我已经使用 lme4 拟合了回归。 .现在我对每个州都有一个回归拟合，我想使用 lattice 为每个州绘制 QQ 图。我还想以格子格式绘制每个状态的误差图。如何使用 lme4 回归的结果制作格子图？

下面是一个使用两种状态的简单示例(是的，我喜欢一个很好的头韵)。我想制作一个由对象制成的两个面板格子适合。

library(lme4)
d <- data.frame(state=rep(c('NY', 'CA'), c(10, 10)), year=rep(1:10, 2), response=c(rnorm(10), rnorm(10)))
fits <- lmList(response ~ year | state, data=d)

Answer 1

我推荐使用更通用的 plyr 包，而不是使用 lmList。

library(plyr)

d <- data.frame(
 state = rep(c('NY', 'CA'), c(10, 10)), 
 year = rep(1:10, 2), 
 response = c(rnorm(10), rnorm(10))
)

# Create a list of models
# dlply = data frame -> list
models <- dlply(d, ~ state, function(df) { 
  lm(response ~ year, data = df)
})

# Extract the coefficients in a useful form
# ldply = list -> data frame
ldply(models, coef)

# We can get the predictions in a similar way, but we need
# to cast to a data frame so the numbers come out as rows,
# not columns.
predictions <- ldply(models, as.data.frame(predict))

predictions 是一个常规的 R 数据框，因此很容易绘制。