Nothing
"cdfgov" <-
function(x, para) {
if(! are.pargov.valid(para)) return()
U <- para$para[1]
A <- para$para[2]
B <- para$para[3]
lo <- U
hi <- U+A
"afunc" <- function(x, f) return(x - quagov(f, para))
f <- sapply(1:length(x), function(i) {
tmp <- NULL
try(tmp <- uniroot(afunc, lower=0, upper=1, x=x[i]), silent=TRUE)
ifelse(is.null(tmp), return(NA), return(tmp$root)) } )
f[x < lo] <- 0
f[x > hi] <- 1
names(f) <- NULL
f[! is.finite(f)] <- NA
return(f)
}
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.