I have a .wig file with three columns (Columns: Pos, Positive Strand, Negative Strand). I am trying to use import.wig into R, however I get the error "data length [1545] is not a sub-multiple or multiple of the number of rows [2]". When opening the wig in R, the columns are shifted and incorrectly filled. How can I import this strand information and maintain the integrity of the columns?
•
written
by
k.smolyar •
1 answer
Perhaps you could convert with BEDOPS wig2bed and import the resulting BED file into R via usual means:
$ wig2bed < in.wig > out.bed
If it is a correctly-formatted Wiggle file, then it should convert correctly. If it is not a correctly-formatted Wiggle file, then some awk swiss-knifing may be needed, as well as a bug report to IGV.
Log in to answer this question.
-
Large number of down-regulated genes compares to upregulated genes
written by singlecell_bio •
-
Importing/formatting .wig file with strand information
written by k.smolyar •
-
How to annotate peaks without strand information
written by ashafik •
-
How to use scanbamparam and CIGAR from Rsamtools to count specific reads?
written by MAPK
-
Compare two columns of same file as paired rows
written by waqasnayab
-
wig file from fasta
written by boczniak767
-
Mapping WGBS probes to CGs (TCGA bed files)
written by sands3 •
-
Determining variant strand from dbSNP entrez document summary
written by sbrown8 •
-
How To Obtain Position, Strand Information And Whether A Sequence Is Transcriptionally Active Or Not
written by roll
-
Reading Bam Files With Rsamtools: Understanding Results For Negative Strand
written by cristianregep •
How was the wig file generated?
The .wig file was generated using igvtools, exporting "counts" coverage information from a BAM file into the three columns.