lasso报错：#Error in nonzeroCoef(beta) : 没有“nonzeroCoef“这个函数

最新推荐文章于 2024-10-04 13:44:40 发布

hx2024

最新推荐文章于 2024-10-04 13:44:40 发布

阅读量112

点赞数

文章标签： r语言

本文链接：https://blog.csdn.net/hx2024/article/details/134223301

版权

文章介绍了R包glmnets中的nonzeroCoef函数，用于处理线性模型中系数的非零检测，支持bystep参数控制按步骤返回哪些变量在每个步骤中是非零的。该函数适用于稀疏矩阵并返回变量索引。

摘要由CSDN通过智能技术生成

#Error in nonzeroCoef(beta) : 没有"nonzeroCoef"这个函数

##手动添加

nonzeroCoef = function (beta, bystep = FALSE) 
{
  ### bystep = FALSE means which variables were ever nonzero
  ### bystep = TRUE means which variables are nonzero for each step
  nr=nrow(beta)
  if (nr == 1) {#degenerate case
    if (bystep) 
      apply(beta, 2, function(x) if (abs(x) > 0) 
        1
        else NULL)
    else {
      if (any(abs(beta) > 0)) 
        1
      else NULL
    }
  }
  else {
    beta=abs(beta)>0 # this is sparse
    which=seq(nr)
    ones=rep(1,ncol(beta))
    nz=as.vector((beta%*%ones)>0)
    which=which[nz]
    if (bystep) {
      if(length(which)>0){
        beta=as.matrix(beta[which,,drop=FALSE])
        nzel = function(x, which) if (any(x)) 
          which[x]
        else NULL
        which=apply(beta, 2, nzel, which)
        if(!is.list(which))which=data.frame(which)# apply can return a matrix!!
        which
      }
      else{
        dn=dimnames(beta)[[2]]
        which=vector("list",length(dn))
        names(which)=dn
        which
      }
      
    }
    else which
  }
}

glmnet source: R/nonzeroCoef.R (rdrr.io)