獲得方差分解表與模型分類變量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

回答