refinery {fda}R Documentation

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:

Time
observation time 0:193
Reflux
reflux flow centered on the mean of the first 60 observations
Tray47
tray 47 level centered on the mean of the first 60 observations

Source

Ramsay, James O., and Silverman, Bernard W. (2006), Functional Data Analysis, 2nd ed., Springer, New York, p. 4, Figure 1.4, and chapter 17.

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)

[Package fda version 2.1.1 Index]