Fits 2D P-splines with anisotropic penalties in both directions.
pspline.2d(
fdataobj,
nbasis.s = 10,
nbasis.t = 10,
lambda.s = 1,
lambda.t = 1,
order = 2,
lambda.select = FALSE,
criterion = c("GCV", "AIC", "BIC")
)A list of class "pspline.2d" similar to pspline().
A 2D fdata object.
Number of B-spline basis functions in s direction.
Number of B-spline basis functions in t direction.
Smoothing parameter in s direction.
Smoothing parameter in t direction.
Order of the difference penalty (default 2).
Logical. If TRUE, select lambdas automatically.
Criterion for selection: "GCV", "AIC", or "BIC".
The 2D penalty uses Kronecker product structure: $$P = \lambda_s (I_t \otimes P_s) + \lambda_t (P_t \otimes I_s)$$