fda (version 1.1.4)

plot.fd: Plot a Functional Data Object

Description

Functional data observations, or a derivative of them, are plotted. These may be either plotted simultaneously, as matplot does for multivariate data, or one by one with a mouse click to move from one plot to another. The function also accepts the other plot specification arguments that the regular plot does.

Usage

plot.fd(x, Lfdobj=0, href=TRUE, nx=201,
        titles=NULL, xlab=xlabel, ylab=ylabel,
        xlim=rangex, ylim=rangey, ask=FALSE, ...)

Arguments

x

a functional data object to be plotted.

Lfdobj

either a nonnegative integer or a linear differential operator object. If present, the derivative or the value of applying the operator is plotted rather than the functions themselves.

href

a logical variable: I if TRUE, the observations are all plotted simultaneously. If not, they are plotted one by one, with a mouse click to move to the next plot.

nx

the number of points to use to define the plot. The default is usually enough, but for a highly variable function more may be required.

titles

a vector of strings for identifying curves

xlab

a label for the horizontal axis.

ylab

a label for the vertical axis.

xlim

a vector of length 2 containing axis limits for the horizontal axis.

ylim

a vector of length 2 containing axis limits for the vertical axis.

ask

a logical value: if TRUE, each curve is shown separately, and the plot advances with a mouse click

...

additional plotting arguments that can be used with function plot

Value

• none

Side Effects

a plot of the functional observations

Details

Note that for multivariate data, a suitable array must first be defined using the par function.

See Also

lines.fd, plotfit.fd

Examples

daytime   <- (1:365)-0.5
dayrange  <- c(0,365)
dayperiod <- 365
nbasis     <- 65
dayrange  <- c(0,365)
daybasis65 <- create.fourier.basis(dayrange, nbasis, dayperiod)
harmaccelLfd <- vec2Lfd(c(0,(2*pi/365)^2,0), dayrange)
harmfdPar     <- fdPar(daybasis65, harmaccelLfd, 1e5)
daytempfd <- data2fd(daily$tempav, daytime, daybasis65,
                     argnames=list("Day", "Station", "Deg C"))

#  plot all the temperature functions for the monthly weather data
#plot(daytempfd, main="Temperature Functions")
#  plot each temperature function in turn, advancing by a click
#plot(daytempfd, matplt=FALSE)