paircomp {hyperdirichlet}R Documentation

Hyperdirichlet distribution for repeated Bernoulli trials

Description

Hyperdirichlet distribution with nonzero coefficients for terms consisting of a pair of ps.

Usage

paircomp(x)

Arguments

x A square matrix. If players i and j play, with i<j, then player j wins x[i,j] trials and player i wins x[j,i] trials

Details

The PDF is

propto prod_{i<j} frac{ p_i^{M[j,i]} p_j^{M[i,j]} }{ (p_i+p_j)^{M[i,j]+M[j,i]} }

Value

Returns a hyperdirichlet distribution (without normalizing factor) of dimension nrow(x).

Note

Elements on the diagonal of x are silently ignored.

Author(s)

Robin K. S. Hankin

See Also

justpairs

Examples

paircomp(matrix(1:9,6,6))

[Package hyperdirichlet version 1.1-7 Index]