我已使用lme包中的nlme函数创建了一个双因素生物实验数据的线性混合效应模型，其中两个因素都是分类的：
model<-lme(responsevariable~factor1*factor2, data=dataset)
我试图在该模型上执行ANOVA之前绘制该模型的分位数 - 分位数图以评估正态性。
我尝试使用qqnorm函数：
qqnorm(model)
但是，在执行上述命令时，我收到以下错误消息：
＆＃34;＆＃34; xyplot＆＃34;没有适用的方法适用于类的对象 ＆＃34; c（＆＃34; double＆＃34;，＆＃34;数字＆＃34;）＆＃34;＆＃34;
我正在寻找关于此错误意味着什么的建议，如果有办法解决它以使qqnorm能够在我的模型上工作，或者如果没有，可以在R中使用哪些其他函数创建我的模型的分位数 - 分位数图。