fda (version 6.1.4)

phaseplanePlot: Phase-plane plot

Description

Plot acceleration (or Ldfobj2) vs. velocity (or Lfdobj1) of a function data object.

Usage

phaseplanePlot(evalarg, fdobj, Lfdobj1=1, Lfdobj2=2,
        lty=c("longdash", "solid"),
      labels=list(evalarg=seq(evalarg[1], max(evalarg), length=13),
             labels=fda::monthLetters),
      abline=list(h=0, v=0, lty=2), xlab="Velocity",
      ylab="Acceleration", returnMatrix=FALSE, ...)

Value

Invisibly returns a matrix with two columns containg the points plotted.

Arguments

evalarg

a vector of argument values at which the functional data object is to be evaluated.

Defaults to a sequence of 181 points in the range specified by fdobj[["basis"]][["rangeval"]].

If(length(evalarg) == 1)it is replaced by seq(evalarg[1], evalarg[1]+1, length=181).

If(length(evalarg) == 2)it is replaced by seq(evalarg[1], evalarg[2], length=181).

fdobj

a functional data object to be evaluated.

Lfdobj1

either a nonnegative integer or a linear differential operator object. The points plotted on the horizontal axis are eval.fd(evalarg, fdobj, Lfdobj1). By default, this is the velocity.

Lfdobj2

either a nonnegative integer or a linear differential operator object. The points plotted on the vertical axis are eval.fd(evalarg, fdobj, Lfdobj2). By default, this is the acceleration.

lty

line types for the first and second halves of the plot.

labels

a list of length two:

evalarg = a numeric vector of 'evalarg' values to be labeled.

labels = a character vector of labels, replicated to the same length as labels[["evalarg"]] in case it's not of the same length.

abline

arguments to a call to abline.

xlab

x axis label

ylab

y axis label

returnMatrix

logical: If TRUE, a two-dimensional is returned using a special class from the Matrix package.

...

optional arguments passed to plot.

See Also

plot, eval.fd plot.fd nondurables

Examples

goodsbasis <- create.bspline.basis(rangeval=c(1919,2000),
                                   nbasis=161, norder=8)
LfdobjNonDur <- int2Lfd(4)
argvals = seq(1919,2000,len=length(nondurables))
logNondurSm <- smooth.basisPar(argvals,
                y=log10(nondurables), fdobj=goodsbasis,
                Lfdobj=LfdobjNonDur, lambda=1e-11)
oldpar <- par(no.readonly=TRUE)
phaseplanePlot(1964, logNondurSm$fd)
par(oldpar)