[FLR-list] fwd(), catch differs from targeted catch

[Simon Fischer (Cefas)]

Hi Laurie,

Thank you very much for your solution, works perfectly!

Regards,
Simon

From: flr-list-bounces at lists.r-forge.r-project.org [mailto:flr-list-bounces at lists.r-forge.r-project.org] On Behalf Of laurie
Sent: 26 February 2016 15:32
To: flr-list at lists.r-forge.r-project.org
Subject: Re: [FLR-list] fwd(), catch differs from targeted catch


There is a problem with indexing ages, if the ages go from 1, then it is OK

load("/home/laurie/Desktop/temp/data.RData")
ctrl <- fwdControl(data.frame(year = 0, quantity = "catch", val = 111111111))
stock=stock[-1]
res <- fwd(stock, ctrl = ctrl, sr = srr)
catch(res)[,"0"]

i.e. any 1st age other than 1 causes a problem.

Laurie

On 26/02/16 13:25, Simon Fischer (Cefas) wrote:
Dear all,

I am working on mse simulations with hypothetical stocks. The stocks are created from life history parameters which are fished at 2Fmsy until an equilibrium is reached. This state is then used as starting point for the simulations.
I use fwd() to project the stock (FLStock) and use catch as a target. Somehow, for some stocks the catch predicted by fwd() differs considerably from the targeted catch (the actual catch is higher than the targeted catch).

Did anyone else encounter this behaviour or knows how to avoid it?

I use the current versions of R and FLR. Attached to this email is a minimum reproducible example that shows this behaviour, the stock object and my sessionInfo().

Thanks in advance,
Simon

Minimal reproducible example:
### load FLR
library(FLash)
library(FLAssess)

### load data
load("data.RData")
### contains FLStock "stock" and FLSR "srr"

### extend stock for 1 year
stock <- stf(object = stock, nyears = 1)

### create fwdControl object
ctrl <- fwdControl(data.frame(year = 1, quantity = "catch", val = 111697431))

### fwd prediction
res <- fwd(stock, ctrl = ctrl, sr = srr)

### check predicted catch
catch(res)[,ac(1)]
### 115308236 and NOT 111697431 as specified in the control object?


#> sessionInfo()
#R version 3.2.3 (2015-12-10)
#Platform: i386-w64-mingw32/i386 (32-bit)
#Running under: Windows 7 x64 (build 7601) Service Pack 1
#
#locale:
#[1] LC_COLLATE=English_United Kingdom.1252
#[2] LC_CTYPE=English_United Kingdom.1252
#[3] LC_MONETARY=English_United Kingdom.1252
#[4] LC_NUMERIC=C
#[5] LC_TIME=English_United Kingdom.1252
#
#attached base packages:
#[1] stats     graphics  grDevices utils     datasets  methods   base
#
#other attached packages:
#[1] FLAssess_2.5.20130716 FLash_2.5.2           FLCore_2.5.20160107
#[4] Matrix_1.2-3          MASS_7.3-45           lattice_0.20-33
#
#loaded via a namespace (and not attached):
#[1] grid_3.2.3   stats4_3.2.3





This email and any attachments are intended for the named recipient only. Its unauthorised use, distribution, disclosure, storage or copying is not permitted. If you have received it in error, please destroy all copies and notify the sender. In messages of a non-business nature, the views and opinions expressed are the author's own and do not necessarily reflect those of Cefas. Communications on Cefas' computer systems may be monitored and/or recorded to secure the effective operation of the system and for other lawful purposes.






_______________________________________________

flr-list mailing list

flr-list at flr-project.org<mailto:flr-list at flr-project.org>

https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list
This email and any attachments are intended for the named recipient only. Its unauthorised use, distribution, disclosure, storage or copying is not permitted.
If you have received it in error, please destroy all copies and notify the sender. In messages of a non-business nature, the views and opinions expressed are the author's own
and do not necessarily reflect those of Cefas. 
Communications on Cefas’ computer systems may be monitored and/or recorded to secure the effective operation of the system and for other lawful purposes.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/flr-list/attachments/20160226/e116c5ea/attachment.html>