rdirichlet: Generate random samples from dirichlet distribution.

Description Usage Arguments Value Author(s) Examples

Description

Generate random samples from dirichlet distribution.

Usage

1
rdirichlet(n, parms)

Arguments

n

number of vectors to be drawn

parms

vector of parameters of dirichlet distribution

Value

n draws from the dirichlet density with parameters in parms

Author(s)

Kate Cowles

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
## The function is currently defined as
function (n, parms) 
{
# generate n random vectors from dirichlet
# rejection envelope

    l <- length(parms)
    x <- matrix(rgamma(l * n, parms), ncol = l, byrow = TRUE)
    sm <- x %*% rep(1, l)
    return(x/as.vector(sm))
  }

CARrampsOcl documentation built on May 2, 2019, 3:27 a.m.