preprocessFunnorm in minfi
1
0
Entering edit mode
b.tawk • 0
@btawk-9390
Last seen 9.0 years ago

Hi everyone,

I am relatively new to R. I am working with 450K and trying to normalize my data using the minfi package.

Besides the QC, checking for bad probes, these are the following steps I use:

require(minfi)
require(minfiData)
base="XYZ"
targets <- read.450k.sheet(base)
RGset <- read.450k.exp(targets = targets)
MSET.fun<-preprocessFunnorm(RGset)

However, every time I try to run the "MSET.fun<-preprocessFunnorm(RGset)", I get a warning message in the end

"[preprocessFunnorm] Background and dye bias correction with noob 
[preprocessNoob] Using sample number 66 as reference level...
[preprocessFunnorm] Mapping to genome
[preprocessFunnorm] Quantile extraction
[preprocessFunnorm] Normalization
Warning message:
In .getSex(CN = CN, xIndex = xIndex, yIndex = yIndex, cutoff = cutoff) :
  An inconsistency was encountered while determining sex. One possibility is that only one sex is present. We recommend further checks, for example with the plotSex function."

However, using the plotSex function, predictedSex and using my clinical annotation, I know that I have males and females in this dataset.

I am not able to pinpoint where the error in the process is. It could be really obvious for everyone else. I would really appreciate any advice I get at this point.

Thanks a lot,

Bouchra

minfi 450k illuminahumanmethylation450k.db normalization • 2.8k views
ADD COMMENT
0
Entering edit mode

The issue isn't whether or not you have both genders, but whether or not the sex determination is consistent when predicted in two different ways. When you use plotSex() with the known sexes, does the plot make sense (e.g., are all the males blue and the females pink)?

ADD REPLY
0
Entering edit mode

Thank you for the reply! There is a small discrepancy in 2/100 samples (something that can be expected to happen in clinical databases) but otherwise, in all the other samples, the predictedSex matches the clinical annotation.

 

ADD REPLY
0
Entering edit mode
@james-w-macdonald-5106
Last seen 17 hours ago
United States

That clears it up then. The function is checking to see if everything looks OK, sees two samples that are questionable and gives you a warning to check them out. If it falls within your expectation, then you can just go ahead with the rest of the analysis.

ADD COMMENT
0
Entering edit mode

Super! Thank you very much!

Best,

B

ADD REPLY

Login before adding your answer.

Traffic: 867 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