Probe filtering beta-values from rgSet
0
0
Entering edit mode
5.8 years ago
moldach ▴ 130

I've see many NA's for beta-values in my rgSet object (either because some probe/CpGs were specifically set NA or because of transformation from intensities to beta-values).

beta <- minfi::getBeta(rgSet)
sumis.na(beta)) # 16,928

I'm trying use the DNAmArray R packages probeFiltering() function to obtain filtered beta-values but I end up getting more missing beta-values which seems odd.

rgSetF <- DNAmArray::probeFiltering(rgSet)
betaF <- minfi::getBeta(rgSetF)
sumis.na(betaF)) # 477,470? Why is there more now?

The rgSet data is ~9Gb so I cannot share it via filedropper but I could probably share it with you via dropbox. I would have to check with my supervisor first but it should be okay if I de-identify samples first. I don't actually know how I would do this to an rgSet object so could someone provided a code example of how to do this please.

R DNAmArray minfi rgSet CpGs • 2.0k views
ADD COMMENT

Login before adding your answer.

Traffic: 2929 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6