r - 複数の予測子を持つ混合モデルでランダムな切片と勾配をプロットする方法は?

Question

複数の予測子がある場合、混合モデルのランダム切片または勾配をプロットすることは可能ですか?

1 つの予測子を使用すると、次のようになります。

#generate one response, two predictors and one factor (random effect)
resp<-runif(100,1, 100)
pred1<-c(resp[1:50]+rnorm(50, -10, 10),resp[1:50]+rnorm(50, 20, 5))
pred2<-resp+rnorm(100, -10, 10)
RF1<-gl(2, 50)

#gamm
library(mgcv)
mod<-gamm(resp ~ pred1, random=list(RF1=~1))
plot(pred1, resp, type="n")
for (i in ranef(mod$lme)[[1]]) {
abline(fixef(mod$lme)[1]+i, fixef(mod$lme)[2])
}

#lmer
library(lme4)
mod<-lmer(resp ~ pred1 + (1|RF1))
plot(pred1, resp, type="n")
for (i in ranef(mod)[[1]][,1]) {
abline(fixef(mod)[1]+i, fixef(mod)[2])
}

しかし、代わりに次のようなモデルがある場合はどうなるでしょうか?:

mod<-gamm(resp ~ pred1 + pred2, random=list(RF1=~1))

またはlmerで

mod<-lmer(resp ~ pred1 + pred2 + (1|RF1))

プロットしている変数のすべての係数または係数のみを考慮する必要がありますか?

ありがとう