Problems using intansv...
Entering edit mode
Phil Ross ▴ 10
Last seen 7.7 years ago
Hello, Within my lab we’ve recently started trying to incorporate different CNV detection programs into our pipeline and found the intansv ( package which is great for integrating the results of different CNV detection programs, combining them into one output dataset, and then using a GFF file to annotate the CNVs that were called by each. However, I’m consistently getting an error that states the following… Error in `$<`(`*tmp*`, "parent", value = c("PF3D7_0100100", : replacement has 6685 rows, data has 7245 The tool isn’t very complicated and I’m not sure why I consistently get this error no matter how I import the GFF file or which cnvcalls output file I use. Here is the code I’m currently using… ## LOAD LIBRARIES library(intansv) library(GenomicRanges) library(rtracklayer) ## IMPORT GFF FILE ### WARNING: PlasmoDB gff file comes improperly formatted...on first commented line change tab to space between 'gff-version' and '3' gffAnnotation <- import.gff("/Users/philross/Desktop/PlasmoDB-9.3_Pfalciparum3D7.gff", asRangedData=FALSE) ## IMPORT CNVNATOR CALLS ### WARNING: Make sure the file doesn't contain a header. If it does insert '#' as the first character on the first line structuralVariation <- readCnvnator(dataDir="/Users/philross/Desktop/cnvnator", method="CNVnator") ## SV ANNOTATION sa <- llply(structuralVariation, svAnnotation, genomeAnnotation=gffAnnotation) The error occurs while running the svAnnotation function seen on the last line. While trying to sort through the debugger I notice that the error occurs between this two calls within the svAnnotation function itself: structuralVariationGeneAnnoNna[structuralVariationGeneAnnoNna$anno tation != "gene", ]$parent <- (genomeAnnotation[as.numeric(structuralVar iationGeneAnnoNna[structuralVariationGeneAnnoNna$annotation != "gene", ]$subject)]$Parent)@unlistData structuralVariationGeneAnnoNa$parent <- NA So it’s the second line there that’s causing the error which I don’t understand…Any ideas? Am I not importing the gene annotation file correctly? Not using arguments somewhere where I should be? Perhaps it’s something wrong with the way the GFF file is formatted? Maintainer: "Wen Yao <>” Any help would be greatly appreciated! Best, Phil Ross [[alternative HTML version deleted]]
annotate annotate • 1.0k views
Entering edit mode

Hi Phill,

I am using intansv and am running into a similar problem.  I was wondering if you figured out how to get around it?



Entering edit mode
akv22 • 0
Last seen 5.9 years ago

I have the same problem. Any solutions yet?


Login before adding your answer.

Traffic: 237 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6