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 npentane (psia).
partial pressure of isopentane (psia).
reaction rate for isomerization of npentane to isopentane (1/hr).
These are data on the reaction rate of catalytic isomerization of npentane 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. 391396.
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 = "npentane partial pressure (psia)",
ylab = expression(plain("Reaction rate (hr")^{1}*plain(")")))
xyplot(rate ~ iso.pen, data = Isom, type = c("g", "p"),
xlab = "isopentane partial pressure (psia)",
ylab = expression(plain("Reaction rate (hr")^{1}*plain(")")))
## Note  the model is misstated 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.