運行中的R

Question

我在下面的函數使用mapply Logistic模型：

logModels <- function(data_idx_list,data) 

{ 
    x <- lapply(data_idx_list, function(m) 
     sparse.model.matrix(~.,data = data[m,])) 
    y <- lapply(data_idx_list, function(m) data$earlyR[m]) 

    logM <- mapply(function(x,y) { 
      cv.glmnet(x=x,y=y,family="binomial",alpha=0) 
}, x,y) 

return(logM) 

} 

其中x包含數據的5個樣品，並y是我的因變量。當我嘗試使用mapply我收到以下錯誤：

Error in glmnet(x, y, weights = weights, offset = offset, lambda = lambda, : number of observations in y (1) not equal to the number of rows of x (100000)

但是，當我運行模式如下原理：

lm1 = cv.glmnet(x=x[[1]],y=y,family="binomial",alpha=0) 

所以，我認爲我有我在接近x的方式中遇到了一些問題，例如我的mapply函數。

您的幫助將不勝感激。

Answer 1

如果我正確理解你的問題，你只需要包裝y成list：

logM <- mapply(function(x, y) cv.glmnet(x, y, family = "binomial", alpha = 0), 
       x, list(y)) 

的y否則單個元素傳遞給您的GLM功能的連續調用。

另外，寫這個的更簡潔的方式：

logM <- mapply(cv.glmnet, x = x, y = list(y), family = "binomial", alpha = 0))