如何为混合模型生成残差与预测变量的图？ [关闭]

Question

我的混合模型如下：

model <- lme(Cost~1+Units, random=~1+Units|Factory, method="ML", data=A)

有人告诉我应用下面的代码来绘制残差与拟合值的关系，并且有效：

plot(fitted(model), resid(model))

但是，当我尝试绘制残差与预测变量Units的关系时，它显示了一条错误消息。

> plot(Units(model), resid(model))
Error: could not find function "Units"
Error in plot(Units(model), resid(model)) : 
error in evaluating the argument 'x' in selecting a method for function 'plot'

我应该怎么做才能解决这个问题？ 谢谢。

Answer 1

我将扩展@Didzis的注释（这是正确的），因此您将真正了解正在发生的事情。 Units是数据中的变量，而不是诸如residuals或fitted values类的某种变量的特殊名称（尽管通常，单位确实具有这种含义...），因此没有任何现成的功能可以为您提供这些功能。 在构建模型时，您已经在代码的第一行中使用了Units：在这里，您使用了带有值A参数data ，该参数告诉lme函数您在建模中使用的变量位于称为A的数据框内。

您可以使用带有要提取的变量名称的美元符号$来访问数据框中的变量（可以使用函数names(A)检查名称）。 因此，您可以通过使用命令A$Units来获得Units变量，因此可以使用命令来绘制Units变量与残差的关系

plot(A$Units,resid(model))