GetMeanCurve: Mean Curve

In fdapace: Functional Data Analysis and Empirical Dynamics

Mean Curve

Description

Mean curve calculation for dense or sparse functional data.

Usage

GetMeanCurve(Ly, Lt, optns = list())

Arguments

Ly	A list of n vectors containing the observed values for each individual. Missing values specified by NAs are supported for dense case (dataType='Dense').
Lt	A list of n vectors containing the observation time points for each individual corresponding to y. Each vector should be sorted in ascending order.
optns	A list of options control parameters specified by list(name=value). See ‘Details’. Available control options are userBwMu The bandwidth value for the smoothed mean function (using 'CV' or 'GCV'); positive numeric - default: determine automatically based on 'methodBwMu' methodBwMu The bandwidth choice method for the mean function; 'GMeanAndGCV' (the geometric mean of the GCV bandwidth and the minimum bandwidth),'CV','GCV' - default: 5% of the support dataType The type of design we have (usually distinguishing between sparse or dense functional data); 'Sparse', 'Dense', 'DenseWithMV', 'p>>n' - default: determine automatically based on 'IsRegular' plot Plot mean curve; logical - default: FALSE kernel Smoothing kernel choice, common for mu and covariance; "rect", "gauss", "epan", "gausvar", "quar" - default: "gauss"; dense data are assumed noise-less so no smoothing is performed. kFoldMuCov The number of folds to be used for mean and covariance smoothing. Default: 10 methodMuCovEst The method to estimate the mean and covariance in the case of dense functional data; 'cross-sectional', 'smooth' - default: 'cross-sectional' numBins The number of bins to bin the data into; positive integer > 10, default: NULL useBinnedData Should the data be binned? 'FORCE' (Enforce the # of bins), 'AUTO' (Select the # of bins automatically), 'OFF' (Do not bin) - default: 'AUTO' userMu The user-defined smoothed mean function; list of two numerical vector 't' and 'mu' of equal size, 't' must cover the support defined 'Ly' - default: NULL useBW1SE Pick the largest bandwidth such that CV-error is within one Standard Error from the minimum CV-error, relevant only if methodBwMu ='CV' and/or methodBwCov ='CV'; logical - default: FALSE

Value

A list containing the following fields:

mu	A vector of length nWorkGrid containing the mean function estimate.
workGrid	A vector of length nWorkGrid. The internal regular grid on which the mean estimation is carried out.
bwMu	The selected (or user specified) bandwidth for smoothing the mean function.
optns	A list of actually-used options relevant to the mean function calculation.

Examples

set.seed(1)
n <- 20
pts <- seq(0, 1, by=0.025)
sampWiener <- Wiener(n, pts)
mu = sin(2*pi*pts)
sampWiener <- Sparsify(t(t(sampWiener) + mu), pts, 10)
res = GetMeanCurve(Ly = sampWiener$Ly, Lt = sampWiener$Lt, optns = list(plot = TRUE))

fdapace documentation built on July 3, 2024, 5:08 p.m.