cics_unif_explicit_smooth {ICSsmoothing}

Smooth given data set by k-component uniform clamped interpolating spline (UcICS).

Description

cics_unif_explicit_smooth constructs the uniform clamped interpolating spline with k components that smoothes given data set {(xx[i],yy[i]), i=1..length(xx)}.

Usage

cics_unif_explicit_smooth(
  xx,
  yy,
  k,
  clrs = c("blue", "red"),
  d,
  xlab = NULL,
  ylab = NULL,
  title = NULL,
  plotTF = TRUE
)

Arguments

xx

a vector of data set's x-coordinates (that are in increasing order).

yy

a vector of datanvidi set's y-coordinates.

k

a chosen number of components of smoothing UcICS (integer \ge 2).

clrs

a vector of colours that are used alternately to plot the graph of spline's components.

d

a vector (optional parameter) that contains two values of derivative, in the first and the last computed node. If missing, values of derivative are estimated by given linear regression model. If present, their contribution is removed from linear model and only function values are estimated.

xlab

a title (optional parameter) for the x axis.

ylab

a title (optional parameter) for the y axis.

title

a title (optional parameter) for the plot.

plotTF

a boolean value (optional parameter), if TRUE then plot.

Value

a list with components

nodes

vector of equidistant nodes, based on which we construct the smoothing spline.

est_spline_coeffs

4-element array of (k)x(k+3) matrices, whereas element in i-th row and j-th of l-th matrix contains coefficient by x^{l-1} of cubic polynomial, which is in i-th row and j-th column of matrix B from smoothing spline's explicit form

S=B.\gamma.

est_spline_polynomials

list of string representations of smoothing UcICS.

est_gamma

vector of estimated smoothing spline's coefficients (function values and exterior-node derivatives).

aux_BF

A basis function of the spline

aux_tridiag_inverse

An inverse of the tridiagonal matrix used for spline derivatives construction

aux_M

An estimation matrix used to compute est_gamma

Examples

cp <- cics_unif_explicit_smooth(
xx = CERN$x,
  yy = CERN$y,
  k = 19, #23,
  d = c(1, 0),
  xlab = "X axis",
  ylab = "Y axis"
)