Sample from the Dirichlet distribution

Sample from the Dirichlet distribution

Usage

rdirichlet(n, alpha, k = NULL)

Arguments

n

Number of samples (rows) to draw

alpha

Concentration parameters; either length 1 or length > 1 If length 1, assumes symmetric Dirichlet; k must not be null

k

Number of components (columns); ignored if length(alpha) > 1

Value

A matrix of n rows and length(alpha) or k columns

See also

Other generators: draw_corpus(), journal_specific(), peak_alpha()

Examples

  rdirichlet(10, .1, 5)
#>               [,1]         [,2]         [,3]         [,4]         [,5]
#>  [1,] 6.270234e-09 9.928948e-01 1.658366e-13 7.104752e-03 4.315166e-07
#>  [2,] 4.104904e-01 3.834932e-07 7.583033e-05 5.890177e-01 4.157263e-04
#>  [3,] 7.097726e-09 1.671456e-03 4.583121e-05 1.687170e-01 8.295657e-01
#>  [4,] 8.918237e-03 1.179325e-09 8.420814e-01 1.490003e-01 7.582908e-12
#>  [5,] 9.940553e-01 1.666257e-06 1.796612e-04 3.893606e-09 5.763321e-03
#>  [6,] 6.849839e-09 5.208332e-02 9.471826e-01 4.054764e-04 3.285600e-04
#>  [7,] 8.518659e-04 5.407577e-01 1.236959e-01 8.770145e-04 3.338175e-01
#>  [8,] 1.727927e-05 3.197111e-04 9.053879e-01 7.236306e-02 2.191208e-02
#>  [9,] 6.955948e-04 9.793753e-01 1.020914e-02 9.719905e-03 6.984616e-08
#> [10,] 9.673816e-02 9.031504e-01 5.762476e-16 1.241913e-06 1.101958e-04
  rdirichlet(10, c(.8, .1, .1))
#>            [,1]         [,2]         [,3]
#>  [1,] 0.9686239 2.418037e-18 3.137609e-02
#>  [2,] 0.2608741 7.391259e-01 6.199701e-12
#>  [3,] 0.7793189 4.552116e-06 2.206766e-01
#>  [4,] 0.9851283 1.484478e-02 2.696271e-05
#>  [5,] 0.9997213 2.787471e-04 1.426758e-09
#>  [6,] 0.7472511 5.481167e-03 2.472677e-01
#>  [7,] 0.4447568 5.552328e-01 1.043757e-05
#>  [8,] 0.9347345 1.948347e-03 6.331712e-02
#>  [9,] 0.6735244 2.935070e-07 3.264753e-01
#> [10,] 0.7809019 9.235622e-04 2.181746e-01