Search
Question: PureCN Error in xg[1, ] : incorrect number of dimensions
0
gravatar for twtoal
12 weeks ago by
twtoal0
twtoal0 wrote:

 

I'm running PureCN (development version V1.11.8) on a new project for the first time, and some errors occur, one of which is:

WARN [2018-06-20 13:37:46] Could not impute mapping bias for all variants. Did you use calculateMappingBiasVcf?

INFO [2018-06-20 13:37:46] Sample sex: F

INFO [2018-06-20 13:37:46] Segmenting data...

Error in xg[1, ] : incorrect number of dimensions

Calls: runAbsoluteCN ... <Anonymous> -> .calculateMappingBias -> .adjustEmpBayes

In addition: Warning message:

In .bcfHeaderAsSimpleList(header) :

  duplicate keys in header will be forced to unique rownames

Execution halted

How can I diagnose this problem?

 

ADD COMMENTlink modified 12 weeks ago • written 12 weeks ago by twtoal0

I'll fix the crash, but the problem is very likely your input data. The warning means that for some of the variants in your VCF, PureCN could not find any neighboring SNPs in the mapping database. There is currently no size restriction, so neighborhood is the whole chromosome. This means either your input VCF does not contain many variants after filtering (can you post the log-file?), or your mapping bias database is essentially empty.

 

ADD REPLYlink written 12 weeks ago by markus.riester110
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: 224 users visited in the last hour