I have been looking at dmpFinder function, and I was wondering if the function could work with beta values (as a matrix) and phenotypes (as a vector, for example race) without having the methylset or methylation object. My only worry is that how the function would be able to match the correct phenotypes to correct sample id beta values. Since the sample IDs are columns in the beta matrix, I am not sure what to do. Any suggestions would be appreciated. Thank you!
Hi James, Thank you a lot for responding. Hmm... so do you think if have matrix average_beta below:
and the phenotype dataset below pheno
where both ids are matching in column and row, would dmpfiner give me a correct answer? Thank you! Neyousha
Yes, but you should probably use M-values instead of betas (convert using
logit2
).However,
dmpFinder
is just a way to allow people to uselmFit
fromlimma
to fit a model without having to figure out how to uselmFit
. If you care to fit more than just one coefficient you will have to uselmFit
directly.