# NOT RUN {
## The function is currently defined as
function(xleft, xright, censor_code_filters, distr, Tauy, Tauz,
J) {
pJ <- J / sum(J)
K <- matrix(NA, nrow = length(Tauy), ncol = length(xleft))
for (i in seq(Tauy)) {
K[i, ] <- dkcens2(
xleft = xleft, xright = xright, c_code_filters = censor_code_filters,
distr = distr, mu = Tauy[i], sigma = Tauz[i]
)
}
fcondXA2cens <- apply(K, 2, function(x) sum(x * pJ))
return(fcondXA2cens)
}
# }
Run the code above in your browser using DataLab