Description Usage Arguments Details Value Author(s) Examples
Reduce the size of Poisson-like counts by binomial thinning.
1 | thinCounts(x, prob=NULL, target.size=min(colSums(x)))
|
x |
numeric vector or array of non-negative integers. |
prob |
numeric scalar or vector of same length as |
target.size |
integer scale or vector of the same length as |
If prob
is not NULL
, then this function calls rbinom
with size=x
and prob=prob
to generate the new counts.
This is classic binomial thinning.
The new column sums are random, with expected values determined by prob
.
If prob
is NULL
, then this function does multinomial thinning of the counts to achieve specified column totals.
The default behavior is to thin the columns to have the same column sum, equal to the smallest column sum of x
.
If the elements of x
are Poisson, then binomial thinning produces new Poisson random variables with expected values reduced by factor prob
.
If the elements of each column of x
are multinomial, then multinomial thinning produces a new multinomial observation with a reduced sum.
A vector or array of the same dimensions as x
, with thinned counts.
Gordon Smyth
1 2 | x <- rpois(10,lambda=10)
thinCounts(x,prob=0.5)
|
Loading required package: limma
[1] 3 6 5 1 6 3 3 3 5 7
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.