Nothing
qF <-
function(q, iter, Xpred, icov, P, Beta, Sigma) {
toInvert <- function(x, q, iter, Xpred, icov, P, Beta, Sigma) {
return(sum(P[iter,]*pnorm(x, mean = Xpred[icov,]%*%t(Beta[iter,,]), sd = Sigma[iter,])) - q)
}
res <- uniroot(toInvert, interval = c(-10^15,10^15), q, iter, Xpred, icov, P, Beta, Sigma)$root
return(res)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.