minfi sample limit
1
0
Entering edit mode
igor ▴ 50
@igor
Last seen 18 months ago
United States

I recently started using minfi for analyzing a large set of 450k methylation arrays. For a small trial run, everything worked as expected, but once you get over around 200 samples, there are some problems. Some functions like getBeta() work as expected. However, qcReport() and densityPlot() skip some samples. I don't get any errors. With density plots, it's hard to notice that any are missing, but with the bean plots, certain samples display the sample name, but the graph is blank. Changing the group labels changes which samples show up, so it's not a problem with the samples or the input object. Has anyone else run into this problem? Now I am concerned that other functions may be affected, but I am just not noticing it.

minfi illumina illimina 450k methylation methylation • 1.8k views
ADD COMMENT
0
Entering edit mode
@kasper-daniel-hansen-2979
Last seen 18 months ago
United States
This sounds really weird. I suspect your sample sheet is messed up. Kasper On Thu, Mar 17, 2016 at 10:59 AM, igor [bioc] <noreply@bioconductor.org> wrote: > Activity on a post you are following on support.bioconductor.org > > User igor <https: support.bioconductor.org="" u="" 7184=""/> wrote Question: > minfi sample limit <https: support.bioconductor.org="" p="" 79762=""/>: > > I recently started using minfi for analyzing a large set of 450k > methylation arrays. For a small trial run, everything worked as expected, > but once you get over around 200 samples, there are some problems. Some > functions like getBeta() work as expected. However, qcReport() and > densityPlot() skip some samples. I don't get any errors. With density > plots, it's hard to notice that any are missing, but with the bean plots, > certain samples display the sample name, but the graph is blank. Changing > the group labels changes which samples show up, so it's not a problem with > the samples or the input object. Has anyone else run into this problem? Now > I am concerned that other functions may be affected, but I am just not > noticing it. > > ------------------------------ > > Post tags: minfi, illumina, illimina 450k methylation, methylation > > You may reply via email or visit minfi sample limit >
ADD COMMENT
0
Entering edit mode

If the sample sheet is messed up, then subsetting it to problematic samples should produce problems again. However, if I subset the targets data frame to 20 samples that got skipped in a set of 200 and just use that, everything is fine.

ADD REPLY
0
Entering edit mode
It is not clear that what you say is true. Do you have any duplicates in the columns identifying the individual samples, for example the Basename column? On Thu, Mar 17, 2016 at 11:33 AM, igor [bioc] <noreply@bioconductor.org> wrote: > Activity on a post you are following on support.bioconductor.org > > User igor <https: support.bioconductor.org="" u="" 7184=""/> wrote Comment: minfi > sample limit <https: support.bioconductor.org="" p="" 79762="" #79764="">: > > If the sample sheet is messed up, then subsetting it to problematic > samples should produce problems again. However, if I subset the targets > data frame to 20 samples that got skipped in a set of 200 and just use > that, everything is fine. > > ------------------------------ > > Post tags: minfi, illumina, illimina 450k methylation, methylation > > You may reply via email or visit > C: minfi sample limit >
ADD REPLY
0
Entering edit mode

I checked both Basename and sample name. Both of them only have unique values.

ADD REPLY

Login before adding your answer.

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