Question: Problems using intansv...
5.0 years ago by
Phil Ross • 10
Phil Ross • 10 wrote:
Hello, Within my lab weâve recently started trying to incorporate different CNV detection programs into our pipeline and found the intansv (http://www.bioconductor.org/packages/2.13/bioc/html/intansv.html) 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 `$<-.data.frame`(`*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 <firstname.lastname@example.org>â Any help would be greatly appreciated! Best, Phil Ross [[alternative HTML version deleted]]
ADD COMMENT • link •