Nothing
skip_if(debug_mode)
order <- 4
deg <- order - 1
a <- sqrt(2)
x <- 0:100 / 100
# Create desired knots:
knots <- make_knots(9, deg = deg)
# Create the b-spline basis using splines
old <- splines::splineDesign(knots, x, deg + 1, derivs = 0L, outer.ok = TRUE)
# Create the b-spline basis using splines2
new <- as.matrix(make_basis_matrix(x, knots, deg))
dimnames(new) <- NULL
expect_equal(old, new)
kstep <- 0.1
order <- 3
deg <- order - 1
a <- sqrt(2)
pstep <- 0.01
# Create desired knots:
knots <- make_knots(9, deg = deg)
# Create the b-spline basis using splines
old <- splines::splineDesign(knots, x, deg + 1, derivs = 0L, outer.ok = TRUE)
# Create the b-spline basis using splines2
new <- as.matrix(make_basis_matrix(x, knots, deg))
dimnames(new) <- NULL
expect_equal(old, new)
# Create knot sequence
order <- 2
deg <- order - 1
n_inner <- 3
mu <- 0.3
intercept <- TRUE
knots <- make_knots(n_inner, mu = mu, deg = deg)
K <- length(knots)
B <- splines2_basis(x, knots, deg, intercept = intercept)
# There should be K-order splines
expect_true(dim(B)[2] == K - order)
# Periodic Case
B <- splines2_basis(x, knots, deg, periodic = TRUE, intercept = intercept)
expect_true(dim(B)[2] == K - 2 * deg - 1)
order <- 4
deg <- order - 1
knots <- make_knots(n_inner, mu = mu, deg = deg)
K <- length(knots)
B <- splines2_basis(x, knots, deg, intercept = intercept)
expect_true(!all((B[1, ] - B[nrow(B), ]) == 0))
# There should be K-order splines
expect_true(dim(B)[2] == K - order)
# Periodic Case
B <- splines2_basis(x, knots, deg, periodic = TRUE, intercept = intercept)
expect_true(dim(B)[2] == K - 2 * deg - 1)
expect_true(all((B[1, ] - B[nrow(B), ]) == 0))
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.