Description Usage Details Examples
This function simulates form the conditional posterior distribution of the latents (Y,Z).
1 | fcondYZXA(x, distr, Tauy, Tauz, J)
|
For internal use.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ## The function is currently defined as
function (x, distr = 1, Tauy, Tauz, J)
{
K <- matrix(NA, nrow = length(Tauy), ncol = length(x))
for (i in seq(Tauy)) {
K[i, ] <- dk(x, distr = distr, mu = Tauy[i], sigma = Tauz[i]) *
J[i]
}
if (any(is.na(K)))
print(K, Tauy, Tauz, J)
pK <- prop.table(K, margin = 2)
j <- apply(pK, 2, function(x) sample(length(Tauy), size = 1,
prob = x))
return(matrix(c(y = Tauy[j], z = Tauz[j]), nrow = length(x),
ncol = 2))
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.