Question: PureCN Error in xg[1, ] : incorrect number of dimensions
gravatar for twtoal
27 days ago by
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 27 days ago • written 27 days 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 27 days ago by markus.riester90
Please log in to add an answer.


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