我正在尝试使用glmnet
包中的glmnet
来运行LASSO回归。
我使用以下命令:
library(glmnet)
glmnet(a,b,family="binomial",alpha=1)
我收到了错误:
> Error in if (!all(o)) { : missing value where TRUE/FALSE needed
a
是一个矩阵,带有数值。
b
是一个以因子为值的向量。
但是,b
有一些缺失值。我怀疑这可能是造成错误的原因。但是,我没有看到在glmnet文档中排除NA
的选项。
答案 0 :(得分:4)
由于glmnet
不接受带公式的完整数据框(因此没有na.omit),但使用单独的响应和预测矩阵,您必须在b
中找到哪些值缺少,然后将预测矩阵子集以排除这些行。
library(glmnet)
set.seed(123)
a <- matrix(rnorm(100*20),100,20)
b <- as.factor(sample(0:1,100,replace = TRUE))
b[10] <- NA
na_index <- is.na(b)
res <- glmnet(a[!na_index, ], b[!na_index], family = "binomial", alpha = 1)