获得方差分解表与模型分类变量R中

model = lm(Price ~ factor(Make) + Cylinder + Liter + Doors + Cruise + Sound +leather, data=kuiper)

我没有任何错误的summary(model)输出。

我能当我在模型中使用分类变量得到方差分解表的输出，但我不能够生产使用colldiag功能包perturb当方差分解表（VD表）输出我在模型中有分类变量。

R的投掷以下错误：

print(colldiag(model)) 
Error in v^2 : non-numeric argument to binary operator

是否有用于产生分类变量的方差分解表中的任何备选方案R代码？

当您使用从非默认包的功能，请说明您正在使用的包。一些函数名称出现在多个包中。你在使用'perturb'软件包吗？ – 2014-09-04 10:30:44

是的，先生，我正在使用perturb包。当模型中使用分类变量时，我非常需要方差分解表的输出。请向我推荐任何替代代码。 – 2014-09-05 05:14:50

这篇文章现在大约一岁了。但我也有这个问题，我发布的答案，以防其他人可以使用它。

传递一个model.matrix对象似乎在我的情况下工作：

cd <- colldiag(model.matrix(model)) 
print(cd)

2015-08-24 12:15:25 koohyar

回答