sciruby 0.1.3 → 0.2.1

data/data/r/man/Indometh.Rd

@@ -1,57 +0,0 @@
-% File src/library/datasets/man/Indometh.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2011 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{Indometh}
-\docType{data}
-\alias{Indometh}
-\title{Pharmacokinetics of Indomethacin}
-\description{
-  The \code{Indometh} data frame has 66 rows and 3 columns of data on
-  the pharmacokinetics of indometacin (or, older spelling,
-  \sQuote{indomethacin}).
-}
-\usage{Indometh}
-\format{
-  This object of class \code{c("nfnGroupedData", "nfGroupedData",
-    "groupedData", "data.frame")} containing the following columns:
- \describe{
-    \item{Subject}{
-      an ordered factor with containing the subject codes.  The
-      ordering is according to increasing maximum response.
-    }
-    \item{time}{
-      a numeric vector of times at which blood samples were drawn (hr).
-    }
-    \item{conc}{
-      a numeric vector of plasma concentrations of indometacin (mcg/ml).
-    }
-  }
-}
-\details{
-  Each of the six subjects were given an intravenous injection of
-  indometacin.
-  
-  This dataset was originally part of package \code{nlme}, and that has
-  methods (including for \code{[}, \code{as.data.frame}, \code{plot} and
-  \code{print}) for its grouped-data classes. 
-}
-\source{
-  Kwan, Breault, Umbenhauer, McMahon and Duggan (1976)
-  Kinetics of Indomethacin absorption, elimination, and
-    enterohepatic circulation in man.
-  \emph{Journal of Pharmacokinetics and Biopharmaceutics} \bold{4},
-  255--280.
-
-  Davidian, M. and Giltinan, D. M. (1995)
-  \emph{Nonlinear Models for Repeated Measurement Data},
-  Chapman & Hall (section 5.2.4, p. 129)
-
-  Pinheiro, J. C. and Bates, D. M. (2000) \emph{Mixed-effects Models in
-    S and S-PLUS}, Springer.
-}
-\seealso{
-  \code{\link{SSbiexp}} for models fitted to this dataset.
-}
-\keyword{datasets}

data/data/r/man/InsectSprays.Rd

@@ -1,45 +0,0 @@
-% File src/library/datasets/man/InsectSprays.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{InsectSprays}
-\docType{data}
-\alias{InsectSprays}
-\title{Effectiveness of Insect Sprays}
-\description{
-  The counts of insects in agricultural experimental units treated with
-  different insecticides.
-}
-\usage{InsectSprays}
-\format{
-  A data frame with 72 observations on 2 variables.
-  \tabular{rlll}{
-    [,1]  \tab count  \tab numeric  \tab Insect count\cr
-    [,2]  \tab spray  \tab factor   \tab The type of spray
-  }
-}
-\source{
-  Beall, G., (1942)
-  The Transformation of data from entomological field experiments,
-  \emph{Biometrika}, \bold{29}, 243--262.
-}
-\references{
-  McNeil, D. (1977) \emph{Interactive Data Analysis}.
-  New York: Wiley.
-}
-\examples{
-require(stats); require(graphics)
-boxplot(count ~ spray, data = InsectSprays,
-        xlab = "Type of spray", ylab = "Insect count",
-        main = "InsectSprays data", varwidth = TRUE, col = "lightgray")
-fm1 <- aov(count ~ spray, data = InsectSprays)
-summary(fm1)
-opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
-plot(fm1)
-fm2 <- aov(sqrt(count) ~ spray, data = InsectSprays)
-summary(fm2)
-plot(fm2)
-par(opar)
-}
-\keyword{datasets}

data/data/r/man/JohnsonJohnson.Rd

@@ -1,37 +0,0 @@
-% File src/library/datasets/man/JohnsonJohnson.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{JohnsonJohnson}
-\docType{data}
-\alias{JohnsonJohnson}
-\title{Quarterly Earnings per Johnson & Johnson Share}
-\description{
-  Quarterly earnings (dollars) per Johnson & Johnson share 1960--80.
-}
-\usage{JohnsonJohnson}
-\format{
-  A quarterly time series
-}
-\source{
-  Shumway, R. H. and Stoffer, D. S. (2000)
-  \emph{Time Series Analysis and its Applications}.
-  Second Edition.  Springer.  Example 1.1.
-}
-\examples{\donttest{
-require(stats); require(graphics)
-JJ <- log10(JohnsonJohnson)
-plot(JJ)
-## This example gives a possible-non-convergence warning on some
-## platforms, but does seem to converge on x86 Linux and Windows.
-(fit <- StructTS(JJ, type="BSM"))
-tsdiag(fit)
-sm <- tsSmooth(fit)
-plot(cbind(JJ, sm[, 1], sm[, 3]-0.5), plot.type = "single",
-     col = c("black", "green", "blue"))
-abline(h = -0.5, col = "grey60")
-
-monthplot(fit)
-}}
-\keyword{datasets}

data/data/r/man/LakeHuron.Rd

@@ -1,27 +0,0 @@
-% File src/library/datasets/man/LakeHuron.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{LakeHuron}
-\docType{data}
-\alias{LakeHuron}
-\title{Level of Lake Huron 1875--1972}
-\description{
-  Annual measurements of the level, in feet, of Lake Huron 1875--1972.
-}
-\usage{LakeHuron}
-\format{
-  A time series of length 98.
-}
-\source{
-  Brockwell, P. J. and Davis, R. A. (1991).
-  \emph{Time Series and Forecasting Methods}.
-  Second edition. Springer, New York. Series A, page 555.
-
-  Brockwell, P. J. and Davis, R. A. (1996).
-  \emph{Introduction to Time Series and Forecasting}.
-  Springer, New York.
-  Sections 5.1 and 7.6.
-}
-\keyword{datasets}

data/data/r/man/LifeCycleSavings.Rd

@@ -1,54 +0,0 @@
-% File src/library/datasets/man/LifeCycleSavings.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{LifeCycleSavings}
-\docType{data}
-\alias{LifeCycleSavings}
-\title{Intercountry Life-Cycle Savings Data}
-\description{
-  Data on the savings ratio 1960--1970.
-}
-\usage{LifeCycleSavings}
-\format{
-  A data frame with 50 observations on 5 variables.
-  \tabular{rlll}{
-    [,1]  \tab sr    \tab numeric  \tab aggregate personal savings \cr
-    [,2]  \tab pop15 \tab numeric  \tab \% of population under 15 \cr
-    [,3]  \tab pop75 \tab numeric  \tab \% of population over 75 \cr
-    [,4]  \tab dpi   \tab numeric  \tab real per-capita disposable
-    income \cr
-    [,5]  \tab ddpi  \tab numeric  \tab \% growth rate of dpi
-  }
-}
-\source{
-  The data were obtained from Belsley, Kuh and Welsch (1980).
-  They in turn obtained the data from Sterling (1977).
-}
-\details{
-  Under the life-cycle savings hypothesis as developed by Franco
-  Modigliani, the savings ratio (aggregate personal saving divided by
-  disposable income) is explained by per-capita disposable income, the
-  percentage rate of change in per-capita disposable income, and two
-  demographic variables: the percentage of population less than 15
-  years old and the percentage of the population over 75 years old.
-  The data are averaged over the decade 1960--1970 to remove the
-  business cycle or other short-term fluctuations. 
-}
-\references{
-  Sterling, Arnie (1977) Unpublished BS Thesis.
-  Massachusetts Institute of Technology.
-
-  Belsley, D. A., Kuh. E. and Welsch, R. E. (1980)
-  \emph{Regression Diagnostics}.
-  New York: Wiley.
-}
-\examples{
-require(stats); require(graphics)
-pairs(LifeCycleSavings, panel = panel.smooth,
-      main = "LifeCycleSavings data")
-fm1 <- lm(sr ~ pop15 + pop75 + dpi + ddpi, data = LifeCycleSavings)
-summary(fm1)
-}
-\keyword{datasets}

data/data/r/man/Loblolly.Rd

@@ -1,56 +0,0 @@
-% File src/library/datasets/man/Loblolly.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2011 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{Loblolly}
-\docType{data}
-\alias{Loblolly}
-\title{Growth of Loblolly pine trees}
-\description{
-  The \code{Loblolly} data frame has 84 rows and 3 columns of records of
-  the growth of Loblolly pine trees.
-}
-\usage{Loblolly}
-\format{
-  This object of class \code{c("nfnGroupedData", "nfGroupedData",
-    "groupedData", "data.frame")} containing the following columns:
-  \describe{
-    \item{height}{
-      a numeric vector of tree heights (ft).
-    }
-    \item{age}{
-      a numeric vector of tree ages (yr).
-    }
-    \item{Seed}{
-      an ordered factor indicating the seed source for the tree.
-      The ordering is according to increasing maximum height.
-    }
-  }
-}
-\details{   
-  This dataset was originally part of package \code{nlme}, and that has
-  methods (including for \code{[}, \code{as.data.frame}, \code{plot} and
-  \code{print}) for its grouped-data classes. 
-}
-\source{
-  Kung, F. H. (1986),
-  Fitting logistic growth curve with predetermined carrying capacity,
-  in \emph{Proceedings of the Statistical Computing Section,
-    American Statistical Association}, 340--343.
-
-  Pinheiro, J. C. and Bates, D. M. (2000)
-  \emph{Mixed-effects Models in S and S-PLUS}, Springer.
-}
-\examples{
-require(stats); require(graphics)
-plot(height ~ age, data = Loblolly, subset = Seed == 329,
-     xlab = "Tree age (yr)", las = 1,
-     ylab = "Tree height (ft)",
-     main = "Loblolly data and fitted curve (Seed 329 only)")
-fm1 <- nls(height ~ SSasymp(age, Asym, R0, lrc),
-           data = Loblolly, subset = Seed == 329)
-age <- seq(0, 30, length.out = 101)
-lines(age, predict(fm1, list(age = age)))
-}
-\keyword{datasets}

data/data/r/man/Nile.Rd

@@ -1,78 +0,0 @@
-% File src/library/datasets/man/Nile.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{Nile}
-\docType{data}
-\alias{Nile}
-\title{Flow of the River Nile}
-\usage{Nile}
-\description{
-  Measurements of the annual flow of the river Nile at Ashwan 1871--1970.
-}
-\format{
-  A time series of length 100.
-}
-\source{
-  Durbin, J. and Koopman, S. J. (2001) \emph{Time Series Analysis by
-    State Space Methods.}  Oxford University Press.
-  \url{http://www.ssfpack.com/DKbook.html}
-}
-\references{
-  Balke, N. S. (1993) Detecting level shifts in time series.
-  \emph{Journal of Business and Economic Statistics} \bold{11}, 81--92.
-
-  Cobb, G. W. (1978) The problem of the Nile: conditional solution to a
-  change-point problem.  \emph{Biometrika} \bold{65}, 243--51.
-}
-\examples{
-require(stats); require(graphics)
-par(mfrow = c(2,2))
-plot(Nile)
-acf(Nile)
-pacf(Nile)
-ar(Nile) # selects order 2
-cpgram(ar(Nile)$resid)
-par(mfrow = c(1,1))
-arima(Nile, c(2, 0, 0))
-
-## Now consider missing values, following Durbin & Koopman
-NileNA <- Nile
-NileNA[c(21:40, 61:80)] <- NA
-arima(NileNA, c(2, 0, 0))
-plot(NileNA)
-pred <-
-   predict(arima(window(NileNA, 1871, 1890), c(2,0,0)), n.ahead = 20)
-lines(pred$pred, lty = 3, col = "red")
-lines(pred$pred + 2*pred$se, lty=2, col="blue")
-lines(pred$pred - 2*pred$se, lty=2, col="blue")
-pred <-
-   predict(arima(window(NileNA, 1871, 1930), c(2,0,0)), n.ahead = 20)
-lines(pred$pred, lty = 3, col = "red")
-lines(pred$pred + 2*pred$se, lty=2, col="blue")
-lines(pred$pred - 2*pred$se, lty=2, col="blue")
-
-## Structural time series models
-par(mfrow = c(3, 1))
-plot(Nile)
-## local level model
-(fit <- StructTS(Nile, type = "level"))
-lines(fitted(fit), lty = 2)              # contemporaneous smoothing
-lines(tsSmooth(fit), lty = 2, col = 4)   # fixed-interval smoothing
-plot(residuals(fit)); abline(h = 0, lty = 3)
-## local trend model
-(fit2 <- StructTS(Nile, type = "trend")) ## constant trend fitted
-pred <- predict(fit, n.ahead = 30)
-## with 50\% confidence interval
-ts.plot(Nile, pred$pred,
-        pred$pred + 0.67*pred$se, pred$pred -0.67*pred$se)
-
-## Now consider missing values
-plot(NileNA)
-(fit3 <- StructTS(NileNA, type = "level"))
-lines(fitted(fit3), lty = 2)
-lines(tsSmooth(fit3), lty = 3)
-plot(residuals(fit3)); abline(h = 0, lty = 3)
-}
-\keyword{datasets}

data/data/r/man/Orange.Rd

@@ -1,57 +0,0 @@
-% File src/library/datasets/man/Orange.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2011 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{Orange}
-\docType{data}
-\alias{Orange}
-\title{Growth of Orange Trees}
-\description{
-  The \code{Orange} data frame has 35 rows and 3 columns of records of
-  the growth of orange trees.
-}
-\usage{Orange}
-\format{
-  This object of class \code{c("nfnGroupedData", "nfGroupedData",
-    "groupedData", "data.frame")} containing the following columns:
-  \describe{
-    \item{Tree}{
-      an ordered factor indicating the tree on which the measurement is
-      made.  The ordering is according to increasing maximum diameter.
-    }
-    \item{age}{
-      a numeric vector giving the age of the tree (days since 1968/12/31)
-    }
-    \item{circumference}{
-      a numeric vector of trunk circumferences (mm).  This is probably
-      \dQuote{circumference at breast height}, a standard measurement in
-      forestry.
-    }
-  }
-}
-\details{   
-  This dataset was originally part of package \code{nlme}, and that has
-  methods (including for \code{[}, \code{as.data.frame}, \code{plot} and
-  \code{print}) for its grouped-data classes. 
-}
-\source{
-  Draper, N. R. and Smith, H. (1998), \emph{Applied Regression Analysis
-    (3rd ed)}, Wiley (exercise 24.N).
-
-  Pinheiro, J. C. and Bates, D. M. (2000) \emph{Mixed-effects Models
-    in S and S-PLUS}, Springer.
-}
-\examples{
-require(stats); require(graphics)
-coplot(circumference ~ age | Tree, data = Orange, show.given = FALSE)
-fm1 <- nls(circumference ~ SSlogis(age, Asym, xmid, scal),
-           data = Orange, subset = Tree == 3)
-plot(circumference ~ age, data = Orange, subset = Tree == 3,
-     xlab = "Tree age (days since 1968/12/31)",
-     ylab = "Tree circumference (mm)", las = 1,
-     main = "Orange tree data and fitted model (Tree 3 only)")
-age <- seq(0, 1600, length.out = 101)
-lines(age, predict(fm1, list(age = age)))
-}
-\keyword{datasets}

data/data/r/man/OrchardSprays.Rd

@@ -1,62 +0,0 @@
-% File src/library/datasets/man/OrchardSprays.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{OrchardSprays}
-\docType{data}
-\alias{OrchardSprays}
-\title{Potency of Orchard Sprays}
-\description{
-  An experiment was conducted to assess the potency of various
-  constituents of orchard sprays in repelling honeybees, using a
-  Latin square design.
-}
-\usage{OrchardSprays}
-\format{
-  A data frame with 64 observations on 4 variables.
-  \tabular{rlll}{
-    [,1]  \tab rowpos    \tab numeric \tab Row of the design\cr
-    [,2]  \tab colpos    \tab numeric \tab Column of the design\cr
-    [,3]  \tab treatment \tab factor  \tab Treatment level\cr
-    [,4]  \tab decrease  \tab numeric \tab Response
-  }
-}
-\details{
-  Individual cells of dry comb were filled with measured amounts of lime
-  sulphur emulsion in sucrose solution.  Seven different concentrations
-  of lime sulphur ranging from a concentration of 1/100 to 1/1,562,500
-  in successive factors of 1/5 were used as well as a solution
-  containing no lime sulphur. 
-  
-  The responses for the different solutions were obtained by releasing
-  100 bees into the chamber for two hours, and then measuring the
-  decrease in volume of the solutions in the various cells.
-  
-  An \eqn{8 \times 8}{8 x 8} Latin square design was used and the
-  treatments were coded as follows:
-  \tabular{rl}{
-    A \tab highest level of lime sulphur\cr
-    B \tab next highest level of lime sulphur\cr
-    . \tab \cr
-    . \tab \cr
-    . \tab \cr
-    G \tab lowest level of lime sulphur\cr
-    H \tab no lime sulphur
-  }
-}
-\source{
-  Finney, D. J. (1947)
-  \emph{Probit Analysis}.
-  Cambridge.
-}
-\references{
-  McNeil, D. R. (1977)
-  \emph{Interactive Data Analysis}.
-  New York: Wiley.
-}
-\examples{
-require(graphics)
-pairs(OrchardSprays, main = "OrchardSprays data")
-}
-\keyword{datasets}

data/data/r/man/PlantGrowth.Rd

@@ -1,39 +0,0 @@
-% File src/library/datasets/man/PlantGrowth.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{PlantGrowth}
-\docType{data}
-\alias{PlantGrowth}
-\title{Results from an Experiment on Plant Growth}
-\description{
-  Results from an experiment to compare yields (as measured by dried
-  weight of plants) obtained under a control and two different treatment
-  conditions.
-}
-\usage{PlantGrowth}
-\format{
-  A data frame of 30 cases on 2 variables.
-
-  \tabular{rll}{
-    [, 1] \tab weight \tab numeric \cr
-    [, 2] \tab group  \tab factor
-  }
-
-  The levels of \code{group} are \sQuote{ctrl}, \sQuote{trt1}, and \sQuote{trt2}.
-}
-\examples{
-## One factor ANOVA example from Dobson's book, cf. Table 7.4:
-require(stats); require(graphics)
-boxplot(weight ~ group, data = PlantGrowth, main = "PlantGrowth data",
-        ylab = "Dried weight of plants", col = "lightgray",
-        notch = TRUE, varwidth = TRUE)
-anova(lm(weight ~ group, data = PlantGrowth))
-}
-\source{
-  Dobson, A. J. (1983)
-  \emph{An Introduction to Statistical Modelling}.
-  London: Chapman and Hall.
-}
-\keyword{datasets}

data/data/r/man/Puromycin.Rd

@@ -1,84 +0,0 @@
-% File src/library/datasets/man/Puromycin.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2010 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{Puromycin}
-\docType{data}
-\alias{Puromycin}
-\title{Reaction Velocity of an Enzymatic Reaction}
-\description{
-  The \code{Puromycin} data frame has 23 rows and 3 columns of the
-  reaction velocity versus substrate concentration in an enzymatic
-  reaction involving untreated cells or cells treated with Puromycin.
-}
-\usage{Puromycin}
-\format{
-  This data frame contains the following columns:
-  \describe{
-    \item{conc}{
-      a numeric vector of substrate concentrations (ppm)
-    }
-    \item{rate}{
-      a numeric vector of instantaneous reaction rates (counts/min/min)
-    }
-    \item{state}{
-      a factor with levels
-      \code{treated} 
-      \code{untreated} 
-    }
-  }
-}
-\details{
-  Data on the velocity of an enzymatic reaction were obtained
-  by Treloar (1974).  The number of counts per minute of radioactive
-  product from the reaction was measured as a function of substrate
-  concentration in parts per million (ppm) and from these counts the
-  initial rate (or velocity) of the reaction was calculated
-  (counts/min/min).  The experiment was conducted once with the enzyme
-  treated with Puromycin, and once with the enzyme untreated.
-}
-\source{
-  Bates, D.M. and Watts, D.G. (1988),
-  \emph{Nonlinear Regression Analysis and Its Applications},
-  Wiley, Appendix A1.3.
-
-  Treloar, M. A. (1974), \emph{Effects of Puromycin on
-    Galactosyltransferase in Golgi Membranes}, M.Sc. Thesis, U. of
-  Toronto. 
-}
-
-\seealso{
-  \code{\link{SSmicmen}} for other models fitted to this dataset.
-}
-\examples{
-require(stats); require(graphics)
-\testonly{options(show.nls.convergence=FALSE)}
-plot(rate ~ conc, data = Puromycin, las = 1,
-     xlab = "Substrate concentration (ppm)",
-     ylab = "Reaction velocity (counts/min/min)",
-     pch = as.integer(Puromycin$state),
-     col = as.integer(Puromycin$state),
-     main = "Puromycin data and fitted Michaelis-Menten curves")
-## simplest form of fitting the Michaelis-Menten model to these data
-fm1 <- nls(rate ~ Vm * conc/(K + conc), data = Puromycin,
-           subset = state == "treated",
-           start = c(Vm = 200, K = 0.05))
-fm2 <- nls(rate ~ Vm * conc/(K + conc), data = Puromycin,
-           subset = state == "untreated",
-           start = c(Vm = 160, K = 0.05))
-summary(fm1)
-summary(fm2)
-## add fitted lines to the plot
-conc <- seq(0, 1.2, length.out = 101)
-lines(conc, predict(fm1, list(conc = conc)), lty = 1, col = 1)
-lines(conc, predict(fm2, list(conc = conc)), lty = 2, col = 2)
-legend(0.8, 120, levels(Puromycin$state),
-       col = 1:2, lty = 1:2, pch = 1:2)
-
-## using partial linearity
-fm3 <- nls(rate ~ conc/(K + conc), data = Puromycin,
-           subset = state == "treated", start = c(K = 0.05),
-           algorithm = "plinear")
-}
-\keyword{datasets}

data/data/r/man/Theoph.Rd

@@ -1,84 +0,0 @@
-% File src/library/datasets/man/Theoph.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2011 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{Theoph}
-\docType{data}
-\alias{Theoph}
-\title{Pharmacokinetics of Theophylline}
-\description{
-  The \code{Theoph} data frame has 132 rows and 5 columns of data from
-  an experiment on the pharmacokinetics of theophylline.}
-\usage{Theoph}
-\format{
-  This object of class \code{c("nfnGroupedData", "nfGroupedData",
-    "groupedData", "data.frame")} containing the following columns:
-  \describe{
-    \item{Subject}{
-      an ordered factor with levels \code{1}, \dots, \code{12}
-      identifying the subject on whom the observation was made.  The
-      ordering is by increasing maximum concentration of theophylline
-      observed.
-    }
-    \item{Wt}{
-      weight of the subject (kg).
-    }
-    \item{Dose}{
-      dose of theophylline administered orally to the subject (mg/kg).
-    }
-    \item{Time}{
-      time since drug administration when the sample was drawn (hr).
-    }
-    \item{conc}{
-      theophylline concentration in the sample (mg/L).
-    }
-  }
-}
-\details{
-  Boeckmann, Sheiner and Beal (1994) report data from a study by
-  Dr. Robert Upton of the kinetics of the anti-asthmatic drug
-  theophylline.  Twelve subjects were given oral doses of theophylline
-  then serum concentrations were measured at 11 time points over the
-  next 25 hours.
-
-  These data are analyzed in Davidian and Giltinan (1995) and Pinheiro
-  and Bates (2000) using a two-compartment open pharmacokinetic model,
-  for which a self-starting model function, \code{SSfol}, is available.
-
-  This dataset was originally part of package \code{nlme}, and that has
-  methods (including for \code{[}, \code{as.data.frame}, \code{plot} and
-  \code{print}) for its grouped-data classes. 
-}
-\source{
-  Boeckmann, A. J., Sheiner, L. B. and Beal, S. L. (1994), \emph{NONMEM
-    Users Guide: Part V}, NONMEM Project Group, University of
-  California, San Francisco.
-    
-  Davidian, M. and Giltinan, D. M. (1995) \emph{Nonlinear Models for
-    Repeated Measurement Data}, Chapman & Hall (section 5.5, p. 145 and
-  section 6.6, p. 176)
-
-  Pinheiro, J. C. and Bates, D. M. (2000) \emph{Mixed-effects Models in
-    S and S-PLUS}, Springer (Appendix A.29)
-}
-\seealso{\code{\link{SSfol}}}
-\examples{
-require(stats); require(graphics)
-\testonly{options(show.nls.convergence=FALSE)}
-coplot(conc ~ Time | Subject, data = Theoph, show.given = FALSE)
-Theoph.4 <- subset(Theoph, Subject == 4)
-fm1 <- nls(conc ~ SSfol(Dose, Time, lKe, lKa, lCl),
-           data = Theoph.4)
-summary(fm1)
-plot(conc ~ Time, data = Theoph.4,
-     xlab = "Time since drug administration (hr)",
-     ylab = "Theophylline concentration (mg/L)",
-     main = "Observed concentrations and fitted model",
-     sub  = "Theophylline data - Subject 4 only",
-     las = 1, col = 4)
-xvals <- seq(0, par("usr")[2], length.out = 55)
-lines(xvals, predict(fm1, newdata = list(Time = xvals)),
-      col = 4)
-}
-\keyword{datasets}