Description Format Details Source Examples
The Isom
data frame has 24 rows and 4 columns from an
isomerization experiment.
This data frame contains the following columns:
partial pressure of hydrogen (psia).
partial pressure of n-pentane (psia).
partial pressure of isopentane (psia).
reaction rate for isomerization of n-pentane to isopentane (1/hr).
These are data on the reaction rate of catalytic isomerization of n-pentane to isopentane.
Bates and Watts (1998), Nonlinear Regression Analysis and Its Applications, Wiley (Appendix A1.5).
Originally from Carr (1960), Industrial and Engineering Chemistry, 52, pp. 391-396.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | str(Isom)
splom(Isom, main = "Isom data")
xyplot(rate ~ hyd, data = Isom, type = c("g", "p"),
xlab = "Hydrogen partial pressure (psia)",
ylab = expression(plain("Reaction rate (hr")^{-1}*plain(")")))
xyplot(rate ~ n.pent, data = Isom, type = c("g", "p"),
xlab = "n-pentane partial pressure (psia)",
ylab = expression(plain("Reaction rate (hr")^{-1}*plain(")")))
xyplot(rate ~ iso.pen, data = Isom, type = c("g", "p"),
xlab = "iso-pentane partial pressure (psia)",
ylab = expression(plain("Reaction rate (hr")^{-1}*plain(")")))
## Note - the model is mis-stated on page 272, x2 and x3 are reversed
fm1 <-
nls(rate ~ b3*(n.pent - iso.pen/1.632)/(1+b2*hyd+b3*n.pent+b4*iso.pen),
data = Isom, start = c(b2 = 0.1, b3 = 0.1, b4 = 0.1),
algorithm = "plinear", trace = TRUE)
summary(fm1) # compare to Table 2.2, page 56
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.