R/basis.R

Defines functions get_gram_matrix make_design get_knots

# Get inner and outer knots.
get_knots <- function(knots, type){
  deg <- ifelse(type == "smooth", 1, 0)
  knotl = 1 / ( length(knots) - 1 )  # mean length of a knot
  c(rep(-knotl, deg), knots, rep(1+knotl, deg))
}


# Design matrix of mean basis.
make_design <- function(t, knots, type) {
  deg <- ifelse(type == "smooth", 1, 0)
  knots <- get_knots(knots, type)
  design_mat <- splines::splineDesign(knots = knots, x = t, ord = deg + 1)
  design_mat
}


# Get b-spline Gram matrix.
get_gram_matrix <- function(knots, type){
  knots = get_knots(knots, type)
  if( type == "polygon"){
    diag(diff(knots))
  } else {
    osb_smooth = orthogonalsplinebasis::SplineBasis(knots, order = 2)  # degree = order - 1
    orthogonalsplinebasis::GramMatrix(osb_smooth)
  }
}

Try the elastes package in your browser

Any scripts or data that you put into this service are public.

elastes documentation built on July 4, 2022, 1:06 a.m.