prior_unif_dirichlet {saturnin}

Uniform prior counts for the hyper-Dirichlet distribution.

Description

The function generates a (r,r,p,p)-array filled with uniform counts for the hyper-Dirichlet distribution used as prior in lweights_multinomial when there are p variables with r levels. Neq is the equivalent prior sample size.

Usage

prior_unif_dirichlet(p, r, Neq = 0.5 * r^2)

Arguments

p

Number of variables.

r

Number of levels.

Neq

Equivalent sample size.

Value

prior

A (r,r,p,p)-array containings counts.

Author(s)

Loïc Schwaller

References

This package implements the method described in the paper "Bayesian Inference of Graphical Model Structures Using Trees" by L. Schwaller, S. Robin, M. Stumpf, 2015 (submitted and availavable on arXiv).

Examples

library('saturnin')
p <- 100
r <- 10
prior <- prior_unif_dirichlet(p,r)