如何在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图一个？<​​/ p>

提前致谢。

Answer 1

感谢您的回答，但您提出的选项并非我所需要的，因为我得到了四个图，每个因子F2一个，但所有F1一起。

最后我决定用bucle来解决它。我不知道它是否是最好的但是有效......

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方面类似。

再次感谢