R: Mean of Functional Data

mean.fd {fda} | R Documentation

## Mean of Functional Data

### Description

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

### Usage

    ## S3 method for class '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

    ##
    ## 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))

[Package fda version 2.4.8 Index]