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