Question: PureCN Error in xg[1, ] : incorrect number of dimensions
0
gravatar for twtoal
15 months 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?

 

purecn • 223 views
ADD COMMENTlink modified 15 months ago • written 15 months 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 15 months 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 16.09
Traffic: 320 users visited in the last hour