Displaying 1 result from an estimated 1 matches for "fittb".
Did you mean:
fitb
2012 Aug 23
1
NLS bi exponential Fit
...p*exp(a*t) + q*exp(b*t)+rnorm(t,sd=0.3*(p*
exp(a*t) + q*exp(b*t)))
fittA <- nls(y~cbind(exp(a*t), exp(b*t)),
algorithm="plinear",start=list(a=-.1, b=-0.2), data=list(y=y, t=t),
trace=FALSE);fittA
# a b .lin1 .lin2
# -0.003074 -2.777 4512 -2399
fittB <- nls(y~cbind(exp(a*t), exp(b*t)),
algorithm="plinear",start=list(a=-.1, b=-0.3), data=list(y=y, t=t),
trace=FALSE);fittB
# a b .lin1 .lin2
# -0.02248 -0.04684 2414.86017 2052.96601
but
1 - the initial condition is very sensitive, is there...