plotPriors {baySeq} | R Documentation |
This function plots the density of the log values estimated for the mean rate in the data used to estimate a prior distribution for data under the assumption of a Negative Binomial distribution. This function is useful for looking for bimodality of the distributions, and thus determining whether we should try and identify data with no true expression.
plotPriors(cD, group, par = 1)
cD |
|
group |
Which group should we plot the priors for? In general,
should be the group that defines non-differentially expressed
data. Can be defined either as the number of the element in
|
par |
The parameter of the prior that will be plotted. |
If the plot of the data appears bimodal, then it may be sensible to
try and look for data with no true expression by using the option
nullPosts = TRUE
in getLikelihoods.NB
.
Plotting function.
Thomas J. Hardcastle
getPriors.NB
, getLikelihoods.NB
# We load in a `countData' object containing the estimated priors (see `getPriors'). data(CDPriors) plotPriors(CDPriors, group = "NDE", par = 1)