Conditional posterior distribution of the latents Y

Share:

Description

This function simulates form the conditional posterior distribution of the latents Y.

Usage

1
fcondYXA(x, distr, Tau, J, sigma)

Details

For internal use.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
## The function is currently defined as
function (x, distr = 1, Tau, J, sigma) 
{
    K <- matrix(NA, nrow = length(Tau), ncol = length(x))
    for (i in seq(Tau)) {
        K[i, ] <- dk(x, distr = distr, mu = Tau[i], sigma = sigma) * 
            J[i]
    }
    pK <- prop.table(K, margin = 2)
    y <- apply(pK, 2, function(x) sample(Tau, size = 1, prob = x))
    return(y)
  }

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.