R/lmFeatureDS.R

In isglobal-brge/dsOmics: DataSHIELD server site Omic functions

lmFeatureDS <- function(i, model, molecular.data, pheno.data){
  
  sel <- molecular.data$feature[i]
  
  mt <- as.formula(model)
  vars <- all.vars(mt)
  
  data.sel <- data.frame(as.numeric(t(molecular.data[which(molecular.data$feature==sel),-c(which(names(molecular.data)=='feature'))])), pheno.data[, vars])
  colnames(data.sel) <- c(sel, vars)
  rownames(data.sel) <- c()

  return(data.sel)
  
}
# ASSIGN FUNCTION