Draw binary indicators given true class memberships

Usage

draw_indicators(X, rho)

Arguments

X: Integer vector of length n. True latent class (1-indexed).
rho: T x K matrix. rho[t, k] = P(Y_k = 1 | X = t).

Value

An n x K integer matrix of 0/1 values.

Examples

rho <- make_rho(0.9)
X   <- draw_classes(50, c(1/3, 1/3, 1/3))
draw_indicators(X, rho)
#>       [,1] [,2] [,3] [,4] [,5] [,6]
#>  [1,]    0    0    0    0    0    0
#>  [2,]    0    1    1    1    0    0
#>  [3,]    1    1    1    0    1    1
#>  [4,]    0    0    0    0    0    0
#>  [5,]    0    1    1    0    0    0
#>  [6,]    1    1    1    0    1    1
#>  [7,]    1    1    0    0    0    0
#>  [8,]    1    0    1    1    1    0
#>  [9,]    0    1    1    0    0    1
#> [10,]    1    1    1    1    1    1
#> [11,]    0    0    0    0    0    0
#> [12,]    1    1    1    0    0    0
#> [13,]    1    1    1    1    1    0
#> [14,]    1    1    1    1    1    0
#> [15,]    0    0    0    1    0    0
#> [16,]    0    1    0    0    0    0
#> [17,]    0    0    0    0    0    0
#> [18,]    1    1    1    1    0    1
#> [19,]    1    1    1    1    1    1
#> [20,]    0    0    0    0    0    0
#> [21,]    1    0    1    1    1    1
#> [22,]    0    0    0    0    0    0
#> [23,]    1    0    0    1    0    0
#> [24,]    0    0    0    0    0    0
#> [25,]    1    1    0    1    1    1
#> [26,]    0    0    0    0    0    0
#> [27,]    0    0    0    0    0    0
#> [28,]    1    1    1    1    1    1
#> [29,]    1    1    1    1    1    1
#> [30,]    0    0    0    0    0    0
#> [31,]    1    1    0    1    1    1
#> [32,]    0    0    0    0    0    0
#> [33,]    1    0    0    1    0    0
#> [34,]    1    0    1    0    0    0
#> [35,]    1    1    1    1    0    0
#> [36,]    1    1    1    1    1    1
#> [37,]    1    1    1    0    0    0
#> [38,]    1    1    1    0    0    0
#> [39,]    1    0    1    1    1    0
#> [40,]    0    0    0    0    1    0
#> [41,]    1    0    0    0    0    1
#> [42,]    1    1    1    0    0    0
#> [43,]    1    1    1    1    1    0
#> [44,]    0    1    1    0    0    0
#> [45,]    0    1    1    1    1    1
#> [46,]    0    1    1    1    1    1
#> [47,]    1    0    0    0    0    0
#> [48,]    1    1    1    1    1    1
#> [49,]    1    1    1    1    0    1
#> [50,]    1    1    1    1    1    1