r - 一連の線形モデルからの予測値

Question

こんにちは: tidyr と dplyr を使用して構築されたデータ フレームに一連の線形モデルがあります。以下のようです。newdata の固定セットを使用して各モデルから予測値を生成するにはどうすればよいですか? 実際には 10 個の従属変数がありますが、独立変数は 2 つしかありません

#random data
x1<-rnorm(100, mean=10, sd=5)
x2<-rnorm(100, mean=5, sd=2 )
y1<-rnorm(100, mean=5, sd=1)
y2<-rnorm(100, mean=3, sd=1)
#create test data farame
df<-data.frame(y1, y2, x1, x2)
#create models
df%>%
  gather(dv, value, y1, y2, -x1,-x2) %>%
  group_by(dv)%>%
  do(mod=lm(value~x1+x2, data=.))