Displaying 1 result from an estimated 1 matches for "model001epr".
2012 Apr 18
1
Add covariate in nlme?
...use the following code to group my data:
Parameterg <- groupedData( result ~ time | Batch,
data = Batchdata,
labels = list( x = "Time", y = "analysis")
)
and then uses the nlme function to fit a nonlinear mixed model that includes
Process as a fixed covariate:
nlme.model001epr <- nlme(result ~ A0 * exp(- ( exp(A1) * exp(-Ea /
(0.0083144*TEMP.K)) * exp(eps)) * time),
data = Parameterg,
fixed=list(A0+Ea~1,A1~Process),
random=eps~1,
start=c(93, 92, 34.5,37),
control=list(msVerbose=TRUE, maxIter = 200),
verbose=TRUE, method="REML",
na.action=n...