view: Viewing Data with Allele-specific Copy Number

In falcon: Finding Allele-Specific Copy Number in Next-Generation Sequencing Data

Description

This function generates three plots: The first plots the A-allele frequencies of the case (black) sample overlayed onto those of the control (gray) sample; the second plots the relative depth of the case over control adjusted by the ratio of total mapped reads, i.e. P*(read count in tumor)/(read count in normal), where P=(total reads mapped in normal)/(total reads mapped in tumor); the third plots the estimated parent-specific DNA copy numbers.

Usage

 view(output, pos=NULL, rdep=NULL, plot="all", ...) 

Arguments

output	The whole output from calling function "getASCN".
pos	A vector of the base positions for the SNPs. If this information is not provided, the x-axis of the plots will simply be the SNP ordering. If this information is provided, the x-axis of the plots will be the position information.
rdep	The relative depth of the case sample over the control sample. If it is not specified (NULL), then the value median(AT+BT)/median(AN+BN) will be used.
plot	This argument determines what to plot. By default, this function gives all three plots described above ("all"). You can also plot each one individually if you set this argument to either of "Afreq", "RelativeCoverage" or "ASCN".
...	Arguments from plot can be passed along.

See Also

getASCN

Examples

data(Example) 
cn = getASCN(readMatrix, tauhat=tauhat)
view(cn)

# to view the plot for only showing A-allele frequency of the case (black) sample overlayed 
# onto those of the control (gray) sample
par(mfrow=c(1,1))
view(cn, plot="Afreq")

# to view the relative depth of the case over control adjusted by the ratio of total mapped 
# reads in fixed size bins
par(mfrow=c(1,1))
view(cn, plot="RelativeCoverage")

# to view the estimated allele-specific DNA copy numbers
par(mfrow=c(1,1))
view(cn, plot="ASCN")

falcon documentation built on May 2, 2019, 3:30 p.m.