Nothing
sheets <-
function (B, DD, yy, pp, lambda, smooth, nb, center, types)
{
nterms = length(nb)
m = length(yy)
np = length(pp)
lala <- matrix(c(rep(lambda, nterms), rep(lambda, nterms)),
nrow = nterms, ncol = 2, dimnames = list(1:nterms, c("curve",
"sheet")))
med = which(pp == 0.5)
if (smooth == "gcv") {
acv.min = nlminb(start = lala, objective = acv.sheets,
yy = yy, B = B, pp = pp, DD = DD, nb = nb, center = center,
lower = 0, upper = 10000)
min.lambda = matrix(abs(acv.min$par), ncol = 2)
lala[, 1] <- min.lambda[, 1]
lala[, 2] <- min.lambda[, 2]
ynp <- rep(yy, np)
ps <- rep(pp, each = m)
w <- runif(m * np)
p2f.new <- pspfit2d.new(B, DD, ps, ynp, w, lala[, 1],
lala[, 2], center)
}
else {
ynp <- rep(yy, np)
ps <- rep(pp, each = m)
w <- runif(m * np)
p2f.new <- pspfit2d.new(B, DD, ps, ynp, w, lala[, 1],
lala[, 2], center)
}
curves = p2f.new$curves
coefficients <- p2f.new$coef
diag.hat = matrix(p2f.new$hatma, ncol = np)
result = list(coefficients, lala, diag.hat)
result
}
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.