Question: callLOH output file NA values
0
gravatar for twtoal
4 months ago by
twtoal0
twtoal0 wrote:

I updated to a newer version of PureCN and ran into a problem in my post-processing of PureCN output from callLOH. The problem was caused by NA values in the "M" column. My previous PureCN version never had NA values in the column. The PureCN doc file doesn't mention NA in its description of the "M" column. Is there a problem, or is NA expected and just means PureCN is unable to estimate M?

I'd like to suggest whenever possible, the doc files should be updated to indicate changes such as this, and furthermore, since they may have an effect on the user's downstream processing, they should be mentioned in the release notes. Thanks.

purecn na values • 99 views
ADD COMMENTlink modified 4 months ago by markus.riester110 • written 4 months ago by twtoal0
Answer: callLOH output file NA values
0
gravatar for markus.riester
4 months ago by
markus.riester110 wrote:

Yes, the callLOH function changed and I actually mention this in the NEWS file ( http://bioconductor.org/packages/release/bioc/news/PureCN/NEWS ).

LOH calling requires heterozygous SNPs, so NA means no SNPs. The change of defaults was necessary because some users were not aware that segments without SNPs were missing from this output.

ADD COMMENTlink modified 4 months ago • written 4 months ago by markus.riester110

Ok, thanks. It didn't occur to me that that change might introduce NAs in result. Since I don't call callLOH() directly, but use PureCN.R to do it for me, I'm not real tuned in to those function names. What might be helpful is to mention output filenames in the release notes, e.g.:

callLOH (PureCN.R output file *_loh.csv) now reports all segments, even if there are no informative SNPs (indicated by NA in some output columns) since some users were not aware that segments are missing from this output. Use keep.no.snp.segments = FALSE to restore old behaviour.
ADD REPLYlink written 4 months ago by twtoal0

It looks like the keep.no.snp.segments = FALSE argument to callLOH() has no equivalent option argument in PureCN.R ?

ADD REPLYlink written 4 months ago by twtoal0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 145 users visited in the last hour