Formal class VRanges versus Large VRanges
0
0
Entering edit mode
@guillaumedachy-11994
Last seen 22 months ago
Brussels

Hi everyone, 

I am using the SomaticSignatures package with my "usual" script but I get today an error thath I don't manage to fix. 

By importing VCF files with the readVcfAsVRanges function, I get two different objects of different classes incompatible with each other (I can't build a vector to go futher in the analysis). 

VCF are output of Mutect2 pipeline. 

Here is an example of the code I use: 

vr_HNC7 = readVcfAsVRanges("C:/path/to/HNC7.vcf", "GenomeA")
mcols(vr_HNC7)$sample = "HNC7"

vr_HNC9 = readVcfAsVRanges("C:/path/to/HNC9.vcf", "GenomeA")
mcols(vr_HNC9)$sample = "HNC9"

This give me 2 different objects : a "Formal class VRanges" for HNC7 and a "Large VRanges" for HNC9. The only difference is the number of variant (and so the dimensions of the matrix), respectively 50x114 and 167x114. 

When I try to make a vector :

vr_HNCtot = c(vr_HNC2, vr_HNC4, vr_HNC7, vr_HNC9, vr_HNC10, vr_HNC12)
Error in rbind(...) : 
  number of columns of matrices must match (see arg 5)

Is there a way to get the same kind of object to be able to merge them in a vector and to make the SomaticSignatures analysis ?? 

Thanks a lot !!

Cheers

Guillaume

somaticsignatures vranges readvcfasvranges • 347 views
ADD COMMENT

Login before adding your answer.

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