synRNASeqNet (version 1.0)

entropyShrink: James-Stein Shrinkage Entropy Estimate

Description

Computing the James-Stein Shrinkage Entropy Estimate of cellCounts.

Usage

entropyShrink(cellCounts, unit = unit, shrinkageTarget = shrinkageTarget)

Arguments

cellCounts

an integer vector (or matrix) representing the number of times each particular count is obtained.

unit

the unit in which entropy is measured. One of "bit" (log2, default), "ban" (log10) or "nat" (natural units).

shrinkageTarget

shrinkage target frequencies. If not specified (default) it is estimated in a James-Stein-type fashion (uniform distribution).

Value

The entropyShrink function returns the value of the entropy of that gene H(X) (or pair of genes H(X,Y)).

References

James W., Stein C. (1961). Estimation with Quadratic Loss. Proceedings of the Fourth Berkeley Symposium on Mathematical Statistics and Probability, vol. 1 pp. 361-379.

See Also

entropyML, entropyMM, entropyBayes, entropyCS

Examples

simData <- simulatedData(p = 50, n = 100, mu = 100, sigma = 0.25,
                        ppower = 0.73, noise = FALSE)
cellCounts <- table(simData$counts[1, ])
eShrink <- entropyShrink(cellCounts, unit = "nat", shrinkageTarget = NULL)