https://github.com/cran/fda
Tip revision: 229f7206c9196b18bb88d6ee7c3c775e8b28e5d3 authored by J. O. Ramsay on 01 June 2009, 00:00:00 UTC
version 2.1.3
version 2.1.3
Tip revision: 229f720
refinery.Rd
\name{refinery}
\alias{refinery}
\title{
Reflux and tray level in a refinery
}
\description{
194 observations on reflux and "tray 47 level" in a distallation
column in an oil refinery.
}
\format{
A data.frame with the following components:
\describe{
\item{Time}{
observation time 0:193
}
\item{Reflux}{
reflux flow centered on the mean of the first 60 observations
}
\item{Tray47}{
tray 47 level centered on the mean of the first 60
observations
}
}
}
\source{
Ramsay, James O., and Silverman, Bernard W. (2006), \emph{Functional
Data Analysis, 2nd ed.}, Springer, New York, p. 4, Figure 1.4, and
chapter 17.
}
%\seealso{\code{\link{eval.penalty}}}
\examples{
attach(refinery)
# allow space for an axis on the right
op <- par(mar=c(5, 4, 4, 5)+0.1)
# plot uval
plot(Time, Reflux, type="l", bty="n")
# add yval
y.u <- diff(range(Tray47))/diff(range(Reflux))
u0 <- min(Reflux)
y0 <- min(Tray47)
lines(Time, u0+(Tray47-y0)/y.u, lty=3, lwd=1.5, col="red")
y.tick <- pretty(range(Tray47))
axis(4, at=u0+(y.tick)/y.u, labels=y.tick, col="red", lty=3,
lwd=1.5)
# restore previous plot margins
par(op)
detach(refinery)
}
% docclass is function
\keyword{datasets}