在R中的模型中获取具有分类变量的方差分解表

Question

我正在使用代码运行多元回归：

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

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

当我在模型中使用分类变量时，我能够获得方差分解表的输出，但我无法使用colldiag函数生成方差分解表（VD表）的输出当我在模型中有分类变量时，打包perturb。

R抛出以下错误：

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

是否有替代R代码用于为分类变量生成方差分解表？

Answer 1

这篇文章现在已经有一年了。但我也遇到了这个问题，我发布了答案以防万一其他人可以使用它。

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

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