如何在R中繪制混合模型的因子水平的子集

Question

我用lme做了一個混合模型，其中有兩個隨機因子F1（6個級別）和F2（4個級別）。

MiModel<-lme(iv~d1+d1_id,list(Fact1=~1+d1, 
                              Fact2=~-1+d1),
             data=MiData,method="REML")

我想在一些不同的圖中繪制殘差與預測值的關系，每個F2水平對應一個。

我嘗試了這個：

plot(Model,resid(.)~fitted(.)|factor(F1)+factor(F2))

但是，我得到了一個6 x 4的大圖。

是否可以做一些類似的事情，但將圖形分成4個不同的圖，每個F2級一個？

提前致謝。

Answer 1

謝謝您的回答，但是您提出的選項不是我所需要的，因為我得到了四個圖，每個因子F2一張，所有F1在一起。

最終，我決定解決這個問題。 我不知道這是否是最好的，但它可以工作...

for (in in levels(factor(F2)){
NewData<-subset(MiData,F2==i)
NewData$fitt<-predict(MiModel,newdata=NewData)
xyplot(x~fitt|factor(F1), data=NewData)
}

現在我免費使用秤，但是在x和y上卻差不多。

再次感謝