plotHistModels: Plot the model of an isoform pair

In nghiavtr/ISOP: Isoform-Level Expression Patterns in Single-cell RNA-sequencing Data

Plot the model of an isoform pair

Description

Plot the model of an isoform pair

Usage

plotHistModels(deltaVal, dmix, plot.title = NULL, cex.main.size = NULL,
  fit.line = FALSE, group.color = c("darkgreen", "darkorange",
  "darkmagenta", "darkorchid"), lwd = 1)

Arguments

deltaVal	A vector of difference between two isoforms of the isoform pair
dmix	A mixture model from doGaussianMixtureModel function
plot.title	Title of the plot
cex.main.size	Font size of the plot
fit.line	An option for plotting the fitting line of distribution
group.color	A vector of colors for the cell groups
lwd	The thickness of lines

Value

None

Examples

data(isoformDataSample)
#preprocessing
isoformDataSample=ifelse(isoformDataSample <= 3,0,isoformDataSample)
isoformDataSample=isoformDataSample[which(rowSums(isoformDataSample)>0),]
#tranform read count dataset to log scale
isoformDataSample=ifelse(isoformDataSample==0,0,log2(isoformDataSample))
#now data is ready
tbreak=round(sqrt(ncol(isoformDataSample)))
x0=isoformDataSample[1,]
x1=isoformDataSample[2,]
model.res=doMixtureModelIsoformPair(x0,x1,tbreak=tbreak)
dmix=model.res$dmix.list[[model.res$AIC.min.Id]];
deltaVal=x0-x1
plotHistModels(deltaVal=deltaVal,dmix=dmix)

nghiavtr/ISOP documentation built on April 21, 2023, 3:57 p.m.