Refit SSRHE-Style Hessian-Energy Regression

Reuses the operator from fit.ssrhe.hessian.regression to fit new responses or new fixed penalty weights without rebuilding local PCA neighborhoods or Hessian-energy matrices.

Usage

# S3 method for class 'ssrhe.hessian.fit'
refit(
  object,
  y = NULL,
  lambda1 = object$lambda$lambda1,
  lambda2 = object$lambda$lambda2,
  weights = NULL,
  ridge = object$lambda$ridge,
  verbose = FALSE,
  ...
)

Arguments

object

A "ssrhe.hessian.fit" object.

y

New numeric response vector or matrix with one row per vertex. If NULL, the original response is reused.

lambda1

Nonnegative Hessian-energy penalty multiplier for \(f^\top B f = \|Af\|_2^2\).

lambda2

Nonnegative supplemental-stabilizer multiplier for \(f^\top B_S f\). If positive, stabilizer must be TRUE. Currently supported only for derivative.order = 2L.

weights

Optional nonnegative observation weights. May be NULL, a vector of length nrow(X), or a matrix with the same dimensions as y.

ridge

Nonnegative diagonal ridge added to the linear system for numerical stabilization.

verbose

Logical. If TRUE, print a short native construction message.

...

Additional arguments are not currently supported.

Value

A list of class "ssrhe.hessian.refit".

See also

refit for supported fitted objects and migration.

Examples

X <- matrix(seq(0, 1, length.out = 12), ncol = 1)
fit <- fit.ssrhe.hessian.regression(
  X, X[, 1]^2, k = 6L, tangent.dim = 1L, lambda1 = 0.1
)
refit(fit, y = X[, 1]^3)
#> SSRHE Hessian regression refit
#>   responses: 1 
#>   lambda1: 0.1 
#>   lambda2: 0