Silverman / Ramsay log-spline density estimator (legacy R backend)

Original fda::density.fd algorithm, ported verbatim. Use the unified density_mpl() (default backend: "rust") unless you need bit-for-bit compatibility with prior results.

Usage

density_mpl_legacy(
  x,
  WfdParobj,
  conv = 1e-04,
  iterlim = 20,
  active = 1:nbasis,
  dbglev = 0,
  ...
)

Arguments

x

Observation vector or two-column matrix.

WfdParobj

An fdPar object specifying the basis and penalty.

conv

Convergence tolerance on the objective. Default 1e-4.

iterlim

Maximum Fisher-scoring iterations. Default 20.

active

Indices of coefficients that may be updated. Default all.

dbglev

Verbosity for the line-search trace. Default 0.

...

Currently unused; accepted for forward-compatibility.

Value

A list with Wfdobj, C, Flist, iternum, iterhist.