plotGC {alpine} | R Documentation |
Plots smooth curves of the log fragment rate over fragment GC content.
plotGC(fitpar, model, col, lty, ylim, gc.range = NULL, return.type = 0)
fitpar |
a list of the output of fitBiasModels over samples |
model |
the name of one of the models |
col |
a vector of colors |
lty |
a vector of line types |
ylim |
the y limits for the plot |
gc.range |
a numeric of length two, the range of the fragment GC content. By default, [.2,.8] for plotting and [0,1] for returning a matrix |
return.type |
a numeric, either 0: make a plot, 1: skip the plot and return a matrix of log fragment rate, 2: skip the plot and return a matrix of probabilities |
Either plot, or if return.type
is 1 or 2, a matrix
# fitpar was fit using identical code # as found in the vignette, except with # 25 genes, and with fragment size in 80-350 bp data(preprocessedData) perf <- rep(1:2, each=2) plotGC(fitpar, "all", col=perf)