R/unstandardize.R

In lichen-lab/SICS: Phylogeny-Constrained Regularization for Sparse Generalized Linear Models

unstandardizeX <- function(b, center, scale, y,family="gaussian") {
  beta <- matrix(0, nrow=nrow(b), ncol=ncol(b))
  if(family=="gaussian"){
  	   	beta[-1,] <- b[-1,] / scale *sd(y)
  	   	beta[1,] <- (mean(y) - crossprod(center, beta[-1,,drop=FALSE])) #*sd(y)
  }else if(family=="binomial") {
  		beta[-1,] <- b[-1,] / scale
  		beta[1,] <- b[1,]- crossprod(center, beta[-1,,drop=FALSE])
  }
  beta
}