export.wig: Export values in WIG format

Description Usage Arguments Value Author(s) References Examples

Description

Export coverage/intensity values in WIG format, compatible with UCSC genome browser, IGB and others.

Usage

1
export.wig(data, name, chrom = "", filepath = name)

Arguments

data

Coverage/intensity values (numeric vector)

name

Name of the track

chrom

Information about chromosome if not inferrable from data (only for numeric vectors)

filepath

Filepath where to save the object. Chromosome name and "wig" extension will be automatically added

Value

(none)

Author(s)

Oscar Flores oflores@mmb.pcb.ub.es

References

WIG format specification: http://genome.ucsc.edu/FAQ/FAQformat#format6

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
# Load data
data(nucleosome_htseq)
cover <- coverage.rpm(nucleosome_htseq)

# Create wig file
export.wig(cover, name="example_track")

# This would create the file "example_track.chr1.wig" with:

# track type=wiggle_0 name="example_track"
# fixedStep chrom=chr1 start=1 step=1
# 55.55247
# 55.55247
# 55.55247
# 277.7623
# 388.8673
# ...

gthar/nucleR documentation built on May 28, 2019, 4:37 p.m.