mean.fd | R Documentation |
Evaluate the mean of a set of functions in a functional data object.
## S3 method for class 'fd'
mean(x, ...)
x |
a functional data object. |
... |
Other arguments to match the generic function for 'mean' |
a functional data object with a single replication
that contains the mean of the functions in the object fd
.
stddev.fd
,
var.fd
,
sum.fd
,
center.fd
mean
oldpar <- par(no.readonly=TRUE)
##
## 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))
par(oldpar)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.