rnaseqcomp (version 1.2.2)

plot2TX: Estimate And Plot Transcript Proportion Difference

Description

For any compared two replicates in each cell line, the proportion of one transcript for genes that only include two annotated transcripts can be different even flipped. This function estimates and plots the proportion difference stratefied by detrended logsignal. Means of absolute difference will be reported for three levels of detrened logsignals. Average is used when multiple two-replicate comparisons included.

Usage

plot2TX(dat, genes, step = 0.5, thresholds = c(1, 6), plotcell = 1, ...)

Arguments

dat

A rnaseqcomp S4 class object.

genes

A vector of gene names corresponding to quantified transcripts. Note that length(genes) should equal to nrow(dat@quantData[[1]]).

step

A number specifying the resolution on detrended logsignal for calculation and plotting the proportion difference. (default: 0.5)

thresholds

A vector of two numbers define cutoffs for three levels of detreded log signals, where one number summary will be generated. (default: c(1, 6))

plotcell

1 or 2 indicating which cell line will be plotted. If values other than 1 and 2, both cell lines will be plotted. This value won't affect estimation for both cell lines. (default: 1)

...

Parameters for base function plot.

Value

• plot
  2TX plots of quantification pipelines for selected cell line by plotcell.
• 2TX
  A matrix of mean proportion difference. Valuesa are based on averageing two cell lines.

Examples

data(simdata)
condInfo <- factor(simdata$samp$condition)
repInfo <- factor(simdata$samp$replicate)
evaluationFeature <- rep(TRUE, nrow(simdata$meta))
calibrationFeature <- simdata$meta$house & simdata$meta$chr == 'chr1'
unitReference <- 1
dat <- signalCalibrate(simdata$quant, condInfo, repInfo, evaluationFeature,
calibrationFeature, unitReference, calibrationFeature2 = calibrationFeature)
plot2TX(dat,genes=simdata$meta$gene)