Nothing
jacobian <- function(func, nm, x, ...) {
arglist <- list(...)
arglist[[nm]] <- x
f <- do.call(func, arglist)
n <- length(x)
eps <- sqrt(.Machine$double.eps)
df <-matrix(NA, length(f), n)
newarglist <- arglist
for (i in 1:n) {
xx <- abs(x[i])
if(xx == 0)
delta <- eps
else
delta <- eps * xx
newx <- x
newx[i] <- newx[i] + delta
newarglist[[nm]] <- newx
df[,i] <- (do.call(func, newarglist) - f)/delta
}
df
}
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.