fda (version 2.4.7)

mean.fd: Mean of Functional Data

Description

Evaluate the mean of a set of functions in a functional data object.

Usage

# S3 method for fd
mean(x, ...)

Arguments

x

a functional data object.

...

Other arguments to match the generic function for 'mean'

Value

a functional data object with a single replication that contains the mean of the functions in the object fd.

See Also

stddev.fd, var.fd, sum.fd, center.fd mean

Examples

Run this code
# NOT RUN {
##
## 1.  univeriate:  lip motion
##
liptime  <- seq(0,1,.02)
liprange <- c(0,1)

#  -------------  create the fd object -----------------
#       use 31 order 6 splines so we can look at acceleration

nbasis <- 51
norder <- 6
lipbasis <- create.bspline.basis(liprange, nbasis, norder)

#  ------------  apply some light smoothing to this object  -------

lipLfdobj   <- int2Lfd(4)
lipLambda   <- 1e-12
lipfdPar <- fdPar(lipbasis, lipLfdobj, lipLambda)

lipfd <- smooth.basis(liptime, lip, lipfdPar)$fd
names(lipfd$fdnames) = c("Normalized time", "Replications", "mm")

lipmeanfd <- mean.fd(lipfd)
plot(lipmeanfd)

##
## 2.  Trivariate:  CanadianWeather
##
dayrng <- c(0, 365)

nbasis <- 51
norder <- 6

weatherBasis <- create.fourier.basis(dayrng, nbasis)

weather.fd <- smooth.basis(day.5, CanadianWeather$dailyAv,
                           weatherBasis)$fd

str(weather.fd.mean <- mean.fd(weather.fd))

# }

Run the code above in your browser using DataCamp Workspace