Doubt about WL in cn.mops
1
0
Entering edit mode
Cased • 0
@casedugr
Last seen 10 months ago
University of Granada

Hello:

I am Daniel Castillo Secilla, from University of Granada (Spain). I am using cn.mops for detecting CNVs but I have a question about how the paremeters WL (Windows Lentgh) works. When I used samples from the reference human genome hg19 I used a WL very low (5000) because I was trying different sizes, and the resultant plot was right, but then, I tried the same with samples from the current reference genome hg38 and, with this WL the plot has a lot of artifacts as you can see in grch38 with Artifacts image in the link below. So, I setted the WL with the value by default again (25000) and the artifacts of the plot were removed as you can see at figure grch38 without Artifacts image in the link below.


So I don't understand why these artifacts appear when the WL is too small and I hope that you can bring me some light in this situation. 

grch38 with Artifacts --> https://ibb.co/fJpf4o 

grch38 without Artifacts --> https://ibb.co/bA1bRd

I am looking forward to hearing from you.

Best regards.

dnacopy dnaseq genomic cn.mops cnv • 596 views
ADD COMMENT
1
Entering edit mode

see below
 

ADD REPLY
1
Entering edit mode
@gunter-klambauer-5426
Last seen 8 months ago
Austria

Hi,

This is an intrinsic problem of the read count data. If you have a small WL, the counts are low and you will easily get large fold-changes: Imagine that most samples of a readcount of 1 and one sample has a read count of 2. Then this would be a strong signal in terms of fold-changes (2-fold change).

Therefore, I suggest both in the paper and the manual to adjust the WL such that on average 100 reads are in a window. This depends of course on your coverage. The function "getReadCountsFromBAM" suggests and appropriate window length if you leave out that parameter.

Regards,
Günter

ADD COMMENT
0
Entering edit mode

Hi  Mr. Günter:

This information is very useful for me. Thank you very much for your quickly response! I understand your explanation and now I can continue my research without this doubt.

Best regards.

Dani.

ADD REPLY

Login before adding your answer.

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