[FLR-list] Hockey stick
Luis Ridao
luridao at gmail.com
Tue Mar 20 17:45:30 CET 2012
FLR-help,
I'm fitting a Hockey-stick function with following:
fsai.segreg <- fsaiSR
model(fsai.segreg) <- segreg()
fsai.segreg <- fmle(fsai.segreg)
and the parameter estimates are as follows:
> params(fsai.segreg)
An object of class "FLPar"
params
a b
0.63387 45.31551
units: NA
However I want to use the parameter estimates from ADMB
which are something like :
a = 0.6060456
b = 56.51466
How can i set the ADMB-estimates in the 'fsai.segreg' object and recalculate
the residuals and log-likelihood based on these estimates?
I do something like:
> model(fsai.segreg3) <- segreg(FLpar(a = exp(-0.5008), b = exp(4.0345)))
Error in segreg(FLpar(a = exp(-0.5008), b = exp(4.0345))) :
and :
> model(fsai.segreg3) <- segreg(params=list(FLpar(a = exp(-0.5008), b =
exp(4.0345))))
Error in segreg(params = list(FLpar(a = exp(-0.5008), b = exp(4.0345)))) :
unused argument(s) (params = list(FLpar(a = exp(-0.5008), b =
exp(4.0345))))
Thanks in advance
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/flr-list/attachments/20120320/319539e3/attachment-0001.html>
More information about the flr-list
mailing list