The Misra1d
data frame has 14 rows and 2 columns.
This is the same data as Misra1a
but a different model is fit.
This data frame contains the following columns:
A numeric vector of volume values.
A numeric vector of pressure values.
These data are the result of a NIST study regarding dental research in monomolecular adsorption. The response variable is volume, and the predictor variable is pressure.
Misra, D., NIST (1978). Dental Research Monomolecular Adsorption Study.
1 2 3 4 5 6 7 8 9 10 11 12 13 14  Try < function(expr) if (!inherits(val < try(expr), "tryerror")) val
plot(y ~ x, data = Misra1d)
Try(fm1 < nls(y ~ b1*b2*x*((1+b2*x)**(1)), data = Misra1d, trace = TRUE,
start = c(b1 = 500, b2 = 0.0001) ))
Try(fm1a < nls(y ~ b1*b2*x*((1+b2*x)**(1)), data = Misra1d, trace = TRUE,
alg = "port", start = c(b1 = 500, b2 = 0.0001) ))
Try(fm2 < nls(y ~ b1*b2*x*((1+b2*x)**(1)), data = Misra1d, trace = TRUE,
start = c(b1 = 450, b2 = 0.0003) ))
Try(fm2a < nls(y ~ b1*b2*x*((1+b2*x)**(1)), data = Misra1d, trace = TRUE,
alg = "port", start = c(b1 = 450, b2 = 0.0003) ))
Try(fm3 < nls(y ~ b2*x*((1+b2*x)**(1)), data = Misra1d, trace = TRUE,
start = c(b2 = 0.0001), algorithm = "plinear" ))
Try(fm4 < nls(y ~ b2*x*((1+b2*x)**(1)), data = Misra1d, trace = TRUE,
start = c(b2 = 0.0005), algorithm = "plinear" ))

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.