gene level multiple linear regression

mlr(mut, e, covariates = rep(1, length(e)))