export.wig: Export values in WIG format

In gthar/nucleR: Nucleosome positioning package for R

Description

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

Usage

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 [e-mail]

References

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

Examples

# 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.