Description Usage Arguments Details Value Author(s) References See Also Examples
These functions estimate, via maximum likelihood methods, the parameters of the underlying distributions specified in the 'densityFunction' slot of the countData object. A special case is maintained for historical reasons; getPriors.NB estimates parameters for a negative binomial distribution using quasimaximumlikelihood methods.
1 2 3 4 5 
cD 
A 
samplesize 
How large a sample should be taken in estimating the priors? 
samplingSubset 
If given, the priors will be sampled only from the subset specified. 
consensus 
If TRUE, creates a consensus distribution rather than a separate distribution for each member of the groups structure in the ‘cD’ object. See Details. 
cl 
A SNOW cluster object. 
verbose 
Should status messages be displayed? Defaults to TRUE. 
equalDispersions 
Should we assume equal dispersions of data
across all groups in the 
estimation 
Defaults to "QL", indicating quasilikelihood estimation of priors. Currently, the only other possibilities are "ML", a maximumlikelihood method, and "edgeR", the moderated dispersion estimates produced by the 'edgeR' package. See Details. 
zeroML 
Should parameters from zero data (rows that within a group are all zeros) be estimated using maximum likelihood methods (which will result in zeros in the parameters? See Details. 
... 
Additional parameters to be passed to the

These functions empirically estimate prior parameters for the distributions used in estimating posterior likelihoods of each count belonging to a particular group.
For priors estimated for the negative binomial methods, three options
are available. Differences in the options focus on the way in which
the dispersion is estimated for the data. In simulation studies,
quasilikelihood methods ('estimation = "QL"'
) performed best
and so these are used by default. Alternatives are maximumlikelihood
methods ('estimation = "ML"'
), and the 'edgeR' packages
moderated dispersion estimates ('estimation = "edgeR"'
).
The priors estimated for the negative binomial methods
('getPriors.NB'
) may assume that the dispersion of data for a
given row is identical for all group structures defined in
'[email protected]'
('equalDispersions = TRUE'
). Alternatively,
the dispersions may be estimated individually for each group structure
('equalDispersions = FALSE'
). Unless there is a strong reason
for believing that the data are differently dispersed between groups,
'equalDispersions = TRUE'
is recommended. If 'estimation
= "edgeR"'
then this parameter is ignored and dispersion is assumed
identical for all group structures.
If all counts in a given row for a given group are zero, then maximum
and quasilikelihood estimation methods will result in a zero
parameter for the mean. In analyses where segment length is a factor, this makes
it hard to differentiate between (for example) a region which contains
no reads but is only ten bases long and one which likewise contains no
reads but is ten megabases long. If 'zeroML'
is FALSE
,
therefore, the dispersion is set to 1 and the mean estimated as the
value that leaves the likelihood of zero data at fifty percent.
If ‘consensus = TRUE’, then a consensus distribution is created and used for each group in the 'cD' object. This allows faster computation of the priors and likelihoods, but with some degradation of accuracy.
A 'cluster' object is recommended in order to estimate the priors for the negative binomial distribution. Passing NULL to this variable will cause the function to run in nonparallel mode.
A countData
object.
Thomas J. Hardcastle
Hardcastle T.J., and Kelly, K. baySeq: Empirical Bayesian Methods For Identifying Differential Expression In Sequence Count Data. BMC Bioinformatics (2010)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29  # See vignette for more examples.
# If we do not wish to parallelise the functions we set the cluster
# object to NULL.
cl < NULL
# Alternatively, if we have the 'snow' package installed we
# can parallelise the functions. This will usually (not always) offer
# significant performance gain.
## Not run: try(library(snow))
## Not run: try(cl < makeCluster(4, "SOCK"))
# load test data
data(simData)
# Create a {countData} object from test data.
replicates < c("simA", "simA", "simA", "simA", "simA", "simB", "simB", "simB", "simB", "simB")
groups < list(NDE = c(1,1,1,1,1,1,1,1,1,1), DE = c(1,1,1,1,1,2,2,2,2,2))
CD < new("countData", data = simData, replicates = replicates, groups = groups)
#estimate library sizes for countData object
libsizes(CD) < getLibsizes(CD)
# Get priors for negative binomial method
CDPriors < getPriors.NB(CD, samplesize = 10^5, estimation = "QL", cl = cl)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.