Hi, After normalized my data using methylumi, I used limma to find different methylated genes. If I used the pvals as the weights, is it reasonable? I do like this. w<-log(pvals(mldat),0.01) fit1 <- lmFit(exprs(mldat.norm), dm,weights=w) Thanks.