[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