Displaying 1 result from an estimated 1 matches for "rationalcoef".
2012 May 11
0
NLS sensitivity to start= values or poles in data range
...yroot(c(1,B1,B2,B3))
cat((poles), '\n')
# Run the nls using above values as starting points
fullmodel <- {CTE ~ (A0 * K^0 + A1 * K^1 + A2 * K^2 + A3 * K^3)/
(1 + B1 * K^1 + B2 * K^2 + B3 * K^3)}
rationalfit <- nls(fullmodel,start=nlsstart)
summary(rationalfit)
rationalcoef <- as.list(coef(rationalfit))
attach(rationalcoef)
curve((A0*x^0+A1*x^1+A2*x^2+A3*x^3)/(1+B1*x^1+B2*x^2+B3*x^3),
from = min(K),
to = max(K),
add=TRUE, lwd=2)
legend(600,5,c("Starting Values","Initial Value Fit","Full Model Fit","Data"),
lt...