什么是模型矩阵/设计矩阵

Question

我对stats::model.matrixR中的函数之以鼻。在描述中它说它将创建一个设计矩阵。它给了我奇怪的行数，既不对应于我的数据中的观察数，也不对应于我的模型中的参数数。

什么是设计矩阵/模型矩阵？

这是我的用法：

M03b <- glmer(APMs ~ PrePost + Gf + eyeFRF + (1|content) + (eyeFRF|ID),    data=mlmData, family=binomial("logit"))
X <- model.matrix(M03b)

它给了我2895x4矩阵 ID有105个级别，内容为28，所以行数对我来说意义不大。也许缺少价值是一个问题？

Answer 1

在回归模型中，以矩阵向量形式写为

Y = X * B + e,

矩阵X是设计矩阵，而Y是因变量的观测向量，B是响应系数的向量（每个解释变量一个），e是一个向量，其中包含模型误差项对于各种变量的值观察。在设计矩阵中，每一列都是一个解释变量的观察向量。

因此，X的大小必须使得它是一个X m矩阵，而B是amxp矩阵，这使得仅当X的列数等于B的行数时才定义乘积XB。米

给定X的2895x4暗，您应该能够确定B i M03b具有4个响应系数。