如何将预测变量制作成矩阵?
if:formula =喜欢〜水分+甜味
预测变量是水分和甜味。现在我想用矩阵计算
来计算系数b=solve(t(X)%*%X)%*%t(X)%*%Y
如何将X定义为预测变量的矩阵?我需要在R中执行此操作。
答案 0 :(得分:1)
(你应该包括一个最小的可重复的例子,所以我们不必提供一个)
请参阅?model.matrix
这里我为R
附带的数据集计算X'X x <- model.matrix(dist~speed,cars)
crossprod(x)
(Intercept) speed
(Intercept) 50 770
speed 770 13228