gene level multiple linear regression, correcting for total number of mutations

mlr.v2(mut, e, nmut, covariates = 1)