Search
Question: R code for cn.mops -- if a chromosome does not have CNV ?
0
gravatar for Bogdan
11 months ago by
Bogdan480
Palo Alto, CA, USA
Bogdan480 wrote:

Dear Gunter,

please if I may add another 1-2 questions : after running cn.MOPS per chromosome on a PBS computer cluster, sometimes, in the output it says :

"No CNVs detected. Try changing "normalization", "priorImpact" or "thresholds".
Error in calcIntegerCopyNumbers(ref_analysis_norm0) :
  No CNV regions in result object. Rerun cn.mops with different parameters!
Calls: calcIntegerCopyNumbers -> calcIntegerCopyNumbers"

and of course, that is fine, as we do not aim to identify CNV on each chromosome ;) However, please could you advise :

-- is there a way to print a message when there are no detected CNV instead of receiving the message "Error in calcIntegerCopyNumbers(ref_analysis_norm0) : " , as it would help our pipeline

-- in the message it says "Try changing "normalization", "priorImpact" or "thresholds".". Shall we change any of these parameters though ?

thanks a lot !

 

 

 

 

ADD COMMENTlink modified 9 months ago by Günter Klambauer490 • written 11 months ago by Bogdan480
3
gravatar for Günter Klambauer
9 months ago by
Austria
Günter Klambauer490 wrote:

Hello again,

And apologies, again, I had not seen your message until now. I see that the error is problematic for your pipeline  - you can try to catch the error using "try()". However, I will change the behaviour from "error" to "warning" in the next version of cn.mops.

You don't need to change your parameters, if you have already adjusted them and they perform well.

Regards,

Günter

ADD COMMENTlink written 9 months ago by Günter Klambauer490
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 2.2.0
Traffic: 107 users visited in the last hour