Question: How to convert BED file of ATAC seq to Bigwig file
0
kuangliang128 • 0 wrote:
I want to convert such kind of BED file into Bigwig file for visualization on UCSC, and I have tried the "rtracklayer", the error message is as follows:
> library(rtracklayer) #bioconductor
>
> bed_loaded <- import(con="~/R/GSM3733316_E59_distal_femur_peaks_hg19.bed.gz", format="bed") #no need to unzip .gz
> # bed_loaded <- import.bed(con="~/R/GSM3733316_E59_distal_femur_peaks_hg19.bed.gz")
#if you unzip
>
> export.wig(object=bed_loaded, con="~/R/bed2wig.wig") Error in .local(object, con, format, ...) : The score must be numeric, without any NA's
I just start to use R, could someone help me? Thanks in advance!!
ADD COMMENT
• link
•
modified 4 weeks ago
by
GenoMax ♦ 96k
•
written
4 weeks ago by
kuangliang128 • 0
Converstion of BED with scores into bigWig