# NOT RUN {
## 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)
}
# }
Run the code above in your browser using DataLab