I am using an R package to perform an EWAS. The package uses
limma to fit the model for each CpG site. I check the code for
lmfit, specifically the function to fit a robust linear model (
mrlm). I saw that essentially this function uses
MASS in a for loop over the features (e.g., genes). Is there a way in
limma, without having to modify the package, to perform this step in paralle? Thank you.