Displaying 1 result from an estimated 1 matches for "meth_matrix".
Did you mean:
metamatrix
2014 Jul 12
0
lm and 450k data
Hi,
I'm quite new to R and currently trying to use lm to fit linear models but I am currently stuck my code is as follows:
Model1 = function(meth_matrix,exposure, X1, X2, X3, batch) {
mod = lm(meth_matrix[, methcol]~exposure+X1+X2+X3+batch)
res = summary(mod)$coef[2,]
system.time(indiv.results <- mclapply(setNames(ncol=473864(cord_betas)), dimnames(cord_betas)[[2]]), Model1, meth_matrix=cord_betas, exposure=filtered_pheno, X1=covariate...