This function calculates the density of observation(s) under a Mallows-Binomial distribution.

dmb(rankings, ratings, p, theta, M, pi0 = NULL, log = FALSE)

Arguments

rankings

A matrix of rankings, potentially with attribute "assignments" to signify separate reviewer assignments. One ranking per row.

ratings

A matrix of ratings, one row per judge and one column per object.

p

A vector specifying the underlying object qualities. All values between be between 0 and 1, inclusive.

theta

A numeric entry specifying the Mallows scale parameter.

M

Numeric specifying maximum (=worst quality) integer rating.

pi0

A vector specifying the consensus (modal probability) ranking; should be used only for tie-breaking equal values in p.

log

A boolean indicating if the log likelihood should be returned.

Value

A numeric value indicating the (log) likelihood of rankings and ratings under a Mallows distribution.

Examples

data(ToyData1)
dmb(rankings=ToyData1$rankings,ratings=ToyData1$ratings,p=c(.2,.5,.7),theta=1,M=ToyData1$M)
#> [1] 1.481917e-32
dmb(rankings=ToyData1$rankings,ratings=ToyData1$ratings,p=c(.25,.25,.7),theta=1,M=ToyData1$M,
pi0=c(1,2,3),log=TRUE)
#> [1] -57.95753