cvfit 错误 {<data.frame> 的尚未实现的方法}

Question

我正在使用 glmnet 和 cv.glmnet 拟合多项式模型：

fit <- glmnet(x, y, family = "multinomial")
cvfit <- cv.glmnet(x, y, family = "multinomial")

我的变量是

x:
 $ Q1: int [1:1000] 3 3 5 3 3 4 3 3 4 4 ...
 $ Q2: int [1:1000] 4 3 5 5 5 5 5 5 4 4 ...
 $ Q3: int [1:1000] 3 3 4 3 5 5 4 5 5 5 ...
 $ Q4: int [1:1000] 3 5 5 4 5 4 4 4 1 4 ...
 $ Q5: int [1:1000] 2 4 5 3 5 3 4 4 3 5 ...
 $ Q6: int [1:1000] 2 1 4 3 4 3 3 4 3 3 ...
 $ Gender : Factor w/ 2 levels "Female","Male": 1 1 1 1 2  

  

y: Factor w/ 7 levels "A", "B", "C", "D", "E"

但是，我只在使用 cv.glmnet 时收到此错误，但在使用 glmnet 时不会：

Error in cbind2(1, newx) %*% (nbeta[[i]]) : 
  not-yet-implemented method for <data.frame> %*% <dgCMatrix>

希望有人能解释一下原因，谢谢！

Answer 1

您需要使用以下示例确保您的独立或外生x 是矩阵而不是数据框：

library(glmnet)

x = data.frame(matrix(sample(1:5,5000,replace=TRUE),ncol=5))
x$Gender = factor(sample(c("Female","Male"),1000,replace=TRUE))
y = factor(sample(LETTERS[1:7],1000,replace=TRUE))

我们使用model.matrix()，~.作为公式，表示数据框中的所有变量都是自变量，我删除了intercept列，因为它将包含在cv.glmnet中，你可以看到你的因子@987654325 @ 是二进制编码的：

xmat = model.matrix(~ .,data=x)[,-1]
head(xmat)
      X1 X2 X3 X4 X5 GenderMale
1  3  1  1  1  1          0
2  5  2  4  3  2          1
3  1  3  2  5  4          1
4  5  5  1  4  2          0
5  3  3  3  1  2          0
6  1  1  3  1  4          1

运行这个：

cv.glmnet(xmat,y,family="multinomial")

Call:  cv.glmnet(x = xmat, y = y, family = "multinomial") 

Measure: Multinomial Deviance 

     Lambda Index Measure       SE Nonzero
min 0.02215     1     3.9 0.004419       0
1se 0.02215     1     3.9 0.004419       0