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,] 1.532803e-12 2.169466e-06 7.406266e-01 2.593712e-01 2.624708e-08
#>  [2,] 9.977698e-01 2.041237e-03 7.547927e-05 8.240927e-07 1.126776e-04
#>  [3,] 9.046657e-05 9.149266e-12 1.369167e-17 6.126945e-06 9.999034e-01
#>  [4,] 1.644043e-02 6.617479e-04 1.413637e-03 3.338899e-04 9.811503e-01
#>  [5,] 7.659434e-03 9.048060e-03 3.313666e-01 6.519259e-01 8.077258e-10
#>  [6,] 1.145890e-04 7.392255e-02 3.694376e-07 9.259611e-01 1.395315e-06
#>  [7,] 7.163768e-09 2.523227e-07 9.276730e-01 7.220412e-02 1.226443e-04
#>  [8,] 9.785355e-01 1.145658e-03 4.768537e-03 1.554207e-02 8.257981e-06
#>  [9,] 2.071258e-03 9.659243e-01 1.175474e-02 2.024974e-02 1.213172e-11
#> [10,] 9.647006e-02 2.556884e-03 1.855148e-02 5.378448e-01 3.445768e-01
  rdirichlet(10, c(.8, .1, .1))
#>             [,1]         [,2]         [,3]
#>  [1,] 0.89489435 1.051047e-01 9.937702e-07
#>  [2,] 0.99999997 2.546939e-08 5.208616e-16
#>  [3,] 0.99995797 4.202902e-05 2.243478e-16
#>  [4,] 0.86415746 1.231873e-01 1.265520e-02
#>  [5,] 0.05493031 6.767110e-01 2.683587e-01
#>  [6,] 0.14377540 8.483355e-01 7.889068e-03
#>  [7,] 0.97931985 1.905010e-02 1.630049e-03
#>  [8,] 0.27276990 4.473032e-02 6.824998e-01
#>  [9,] 0.85135935 1.486406e-01 5.826725e-09
#> [10,] 0.96300713 8.924995e-03 2.806787e-02