GetMeanCurve {fdapace} | R Documentation |
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 NA s 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))
[Package
fdapace version 0.5.9
Index]