[FLR-list] Initialising an FLSR

Mark Payne mpa at aqua.dtu.dk
Fri Oct 18 14:51:06 CEST 2013


Hi,

Cheers for the reply. I have tried that, but it doesn't actually speed it up! eg

#Without initialisation
system.time(fmle(nsher))  
 user  system elapsed 
  8.444   0.084   8.928 

#With correct initial guesses
system.time(fmle(nsher, start=params(nsher)))
  user  system elapsed 
  8.121   0.060   8.675 

But there may be something wierd with the dispatching - I get this warning to start with

> system.time(fmle(nsher, start=params(nsher)))
Note: method with signature ‘FLModel#FLPar’ chosen for function ‘fmle’,
 target signature ‘FLSR#FLPar’.
 "FLSR#ANY" would also be valid

Mark


________________________________________
Fra: flr-list-bounces at r-forge.wu-wien.ac.at [flr-list-bounces at r-forge.wu-wien.ac.at] på vegne af Iago MOSQUEIRA (JRC) [iago.mosqueira-sanchez at jrc.ec.europa.eu]
Sendt: 18. oktober 2013 14:43
Til: flr-list at flr-project.org
Emne: Re: [FLR-list] Initialising an FLSR

On 10/18/2013 02:27 PM, Mark Payne wrote:
> Hi,
>
> I have an FLSR object that I have fitted. I now wish to use these parameters as the initial conditions for fitting a second FLSR, based on similar but slightly different data. What's the best way to do this?
>
> The @initial slot seems to need a function, which I wasn't expecting... Can I feed the FLPar object to fmle for example?
>
> Cheers,
>
> Mark

Hi,

This should do it, use the start argument

library(FLCore)

data(nsher)

params(nsher)

out <- fmle(nsher, start=params(nsher))



Iago


--
Dr Iago Mosqueira

European Commission – Joint Research Center
Institute for the Protection and Security of the Citizen (IPSC)
Maritime Affairs Unit FISHREG – Scientific Support to Fisheries
TP 051, Via Enrico Fermi 2749
I-21027 Ispra (VA), Italy

Office : +39 0332 785413
Fax: +39 0332 789658

iago.mosqueira at jrc.ec.europa.eu

http://fishreg.jrc.ec.europa.eu/home

_______________________________________________
flr-list mailing list
flr-list at flr-project.org
https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list


More information about the flr-list mailing list