[R] nls confidence intervals
Paul, David A
paulda at BATTELLE.ORG
Thu Aug 14 17:48:19 CEST 2003
You can use the well-known Taylor series approximation to the
variance of an arbitrary function:
Var( f(X) ) ~= Sum( s[i]^2*D2[i] ) + 2*Sum( Sum( s[i,j]*D[i]*D[j] ) )
where D2[i] is the second partial derivative of f(x) with respect
to the ith parameter and D[j] is the first partial derivative of f(x)
with respect to the jth parameter. The indices on the summations
for 2*Sum( Sum( ... ) ) are i=1 to (p-1) and j>i, respectively,
where p denotes the total number of parameters in the model. Also,
s[i]^2 denotes the ith diagonal element of the variance-covariance
matrix for the model, and s[i,j] denotes an off-diagonal element
of the same matrix. You should be able to use vcov( ) to extract
the variance-covariance matrix of your fitted model.
This approximation will estimate the functional form of the variance
of f(X). To get the approximate variance of f(X) for a specific
value of X simply plug in X=x.
After this, you will need to add the estimated variance of the
residuals and take the square root to obtain the standard error
used in calculating prediction intervals. I have used this approach
for some highly nonlinear functions in the past, and the approximation
is only "good" when it is reasonable to assume that E(f(x)) ~= f(E(x)).
[This assumption is present in the derivation of the Taylor
approximation.] When this is not a reasonable assumption, the
approximation can be horrible. In other words, the more locally linear
your nonlinear function is, the better the approximation will work.
Hope this helps,
david paul
-----Original Message-----
From: Enrique Portilla [mailto:portillae at marlab.ac.uk]
Sent: Thursday, August 14, 2003 9:28 AM
To: R-help at stat.math.ethz.ch
Subject: [R] nls confidence intervals
Hi,
Does anyone know how to compute the confidence prediction intervals for a
nonlinear least squares models (nls)?
I was trying to use the function 'predict' as I usually do for other models
fitting (glm, lm, gams...), but it seems that se.fit, and interval
computation is not implemented for the nls...
Cheers
Enrique
~~~~~~~~~~~~~~~~~~~~~~~~~~~
Fisheries Research Services,
Marine Laboratory,
Victoria Road,
Torry,
Aberdeen, UK.
Tel. 44 (0) 1224 295314
~~~~~~~~~~~~~~~~~~~~~~~~~~~
______________________________________________
R-help at stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
More information about the R-help
mailing list