model.matrix生成的行数少于原始data.frame

Question

为什么模型矩阵必须与数据帧具有相同的行数？

mergem = model.matrix(as.formula(paste(response, '~ .')), data=mergef)
dim(mergef)
# [1] 115562     71
dim(mergem)
# [1] 66786   973

我尝试在文档中查找提示但找不到任何内容。提前谢谢。

Answer 1

好吧，如果一行中包含NA，则该行（默认情况下）被删除：

d <- data.frame(x=c(1,1,2), y=c(2,2,4), z=c(4,NA,8))
m <- model.matrix(x ~ ., data=d)

nrow(d) # 3
nrow(m) # 2

此行为由选项“na.action”控制：

options(na.action="na.fail")
m <- model.matrix(x ~ ., data=d) # Error: missing values in object