Hi,
I still seem to be getting errors from trying to run my altered R script, any
advice?
Thanks
Jess
Model1A = function(meth_matrix,exposure, X1, X2, X3, batch) {
+
+ mod = lm(methcol ~ exposure+X1+X2+X3+batch, data = meth_matrix)
+
+
+ res=coef(summary(mod))[2,]
+
+
+ }>
>
> ##Run model1A####
>
> system.time(indiv.results <-mclapply(setNames(seq_len(ncol(cord_betas)),
dimnames(cord_betas)[[2]]), Model1A, meth_matrix=cord_betas,
exposure=filtered_pheno$Number3Mth, X1=covariates$k032, X2=covariates$k021,
X3=covariates$kz029, batch=pdata.B1221.cord$BCDPlate))
user system elapsed
0.044 0.526 0.272
Warning message:
In mclapply(setNames(seq_len(ncol(cord_betas)), dimnames(cord_betas)[[2]]), :
all scheduled cores encountered errors in user code
>
[[alternative HTML version deleted]]