View source: R/plotHistModels.R
plotHistModels | R Documentation |
Plot the model of an isoform pair
plotHistModels(deltaVal, dmix, plot.title = NULL, cex.main.size = NULL,
fit.line = FALSE, group.color = c("darkgreen", "darkorange",
"darkmagenta", "darkorchid"), lwd = 1)
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 |
None
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.