plotLFRDistribution-methods: Plot the LFR distribution

Description Usage Arguments Value See Also Examples

Description

Plot the distribution of the LFR across the genome and across all the wells. Suitable to be call on genomic region of less than 2MB length. Segments represent fragments, the y-axis is the wells and the x-axis the genomic location. The thickness of the segment is proportional to the averagecoverage. The blue segment is an illustration of a fagment of 100kbp (if Region <5Mb) or 1Mb if (region >5Mb) with and average coverage of 1.

Usage

1
plotLFRDistribution(object, chrom, region = NULL, minLength = 10000)

Arguments

chrom

the chromosome to be considered

region

The genomic region to plot. For best results consider a region of less than 2MB.

minLength

Conside ronly the LFR having a length greater or equal to minLength.

x

an object of class LFRset

Value

See Also

plot, codesummary.LFRset, getLFRset, getWells

Examples

1
2
3
4
5
6
7
cat('\n Loading an LFRset data frame ')
LFR_11152Tissue_10k_file<-system.file('extdata',
       'D1000_10k_allfragments.tsv.gz',package = 'DigitalPicoTools')
D1000_10k_11152Tissue_LFR=read.table(file=LFR_11152Tissue_10k_file)
D1000_10k_11152Tissue_LFRset=LFRset(Table=D1000_10k_11152Tissue_LFR)
cat('\n Plotting the Fragment Distribution ')
plotLFRDistribution(D1000_10k_11152Tissue_LFRset,'chr2',region='chr22:22000000-23000000')

chedonat/DigitalPicoTools documentation built on May 13, 2019, 3:39 p.m.