This function performs linear regression with functional responses and scalar predictors by (1) fitting a separate linear model at each point along the function, and then (2) smoothing the resulting coefficients to obtain coefficient functions.
fosr2s(
Y,
X,
argvals = seq(0, 1, , ncol(Y)),
nbasis = 15,
norder = 4,
pen.order = norder - 2,
basistype = "bspline"
)
An object of class fosr
, which is a list with the following
elements:
object of class "{fd}"
representing the
estimated coefficient functions. Its main components are a basis and a
matrix of coefficients with respect to that basis.
\(d\times p\) matrix of coefficient estimates from
regressing on X
separately at each point along the function.
\(d\times p\) matrix of standard errors of the raw coefficient estimates.
\(n\times d\) matrix of fitted values.
\(d\times p\) matrix of coefficient function
estimates, obtained by smoothing the columns of raw.coef
.
\(d\times p\) matrix of coefficient function standard errors.
points at which the coefficient functions are evaluated.
smoothing parameters (chosen by REML) used to smooth the \(p\) coefficient functions with respect to the supplied basis.
the functional responses, given as an \(n\times d\) matrix.
\(n\times p\) model matrix, whose columns represent scalar predictors. Should ordinarily include a column of 1s.
the \(d\) argument values at which the functional responses are evaluated, and at which the coefficient functions will be evaluated.
number of basis functions used to represent the coefficient functions.
norder of the spline basis, when basistype="bspline"
(the default, 4, gives cubic splines).
order of derivative penalty.
type of basis used. The basis is created by an appropriate
constructor function from the fda package; see basisfd. Only "bspline"
and "fourier"
are
supported.
Philip Reiss phil.reiss@nyumc.org and Lan Huo
Unlike {fosr}
and {pffr}
, which obtain smooth
coefficient functions by minimizing a penalized criterion, this function
introduces smoothing only as a second step. The idea was proposed by Fan
and Zhang (2000), who employed local polynomials rather than roughness
penalization for the smoothing step.
Fan, J., and Zhang, J.-T. (2000). Two-step estimation of functional linear models with applications to longitudinal data. Journal of the Royal Statistical Society, Series B, 62(2), 303--322.
{fosr}
, {pffr}