Nothing
ins <-
function (x, df = NULL, knots = NULL, intercept = FALSE, Boundary.knots = range(x),
from = 0, weight.fun = NULL, integrand.fun = NULL, ...) {
if (!is.null(weight.fun) && !is.function(weight.fun))
stop("'weight.fun' must be a function.\n")
ns.x <- if (is.null(knots)) {
splines::ns(x, df = df, intercept = intercept, Boundary.knots = Boundary.knots)
} else {
splines::ns(x, knots = knots, intercept = intercept, Boundary.knots = Boundary.knots)
}
kn <- attr(ns.x, "knots")
Bkn <- attr(ns.x, "Boundary.knots")
GK <- gaussKronrod(15)
wk <- GK$wk
sk <- GK$sk
P1 <- c(x + from) / 2
P2 <- c(x - from) / 2
st <- outer(P2, sk) + P1
out <- vector("list", 15L)
for (i in seq_len(15)) {
basis <- splines::ns(st[, i], knots = kn, Boundary.knots = Bkn,
intercept = intercept)
if (!is.null(integrand.fun))
basis <- integrand.fun(basis)
out[[i]] <- wk[i] * basis
if (!is.null(weight.fun)) {
ww <- c(weight.fun(st[, i], x, ...))
out[[i]] <- out[[i]] * ifelse(is.finite(ww), ww, 0)
}
}
out <- P2 * Reduce("+", out)
attr(out, "from") <- from
attr(out, "weight.fun") <- weight.fun
attr(out, "class") <- c("ins", "basis", "matrix")
out
}
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.