[FLR-list] Problems with fwd()

Laurie lauriekell at googlemail.com
Wed Mar 21 18:06:54 CET 2012


This is wrong

target

Target
   year quantity min       val   max
1 2011        f  NA 0.4550558    NA
2 2012        f  NA 0.3200000    NA
3 2012    catch   0        NA 40000
4 2012        f  NA2.0000000     NA


     min         val         max
   1          NA     0.45506          NA
   2          NA     0.32000          NA
   3     0.00000          NA 40000.00000
   4          NA          NA2.00000


Try this! I have added an argument to fwd maxF that sets the upper bound on F

print("OK to here1")
       ## project MP to TAC year with the HCR and geom mean recr.
       c.=fwdControl(data.frame(val=c(rep(0.2,2),NA),min=c(NA,NA,0),max=c(NA,NA,40000),year=c(2011:2012,2012),quantity=c(rep("f",2),"catch")))
       MP<-fwd(MP,ctrl=c.,sr=list(model="mean",params=FLPar(1)),sr.residuals=SRrs,maxF=2)  ## HERE IS THE PROBLEM
   print("OK to here2")



Laureie
On 03/21/2012 05:15 PM, Luis Ridao wrote:
> here it comes ... sorry about it
>
> best,
> Luis
>
> On Wed, Mar 21, 2012 at 4:11 PM, Laurie <lauriekell at googlemail.com 
> <mailto:lauriekell at googlemail.com>> wrote:
>
>     There is no funcs.R file
>
>     Laurie
>
>
>     On 03/20/2012 01:09 PM, Luis Ridao wrote:
>>
>>     I have attached the input files and the script (the workspace is
>>     too large >60 MB)
>>
>>
>>     Best,
>>     Luis
>>
>>     On Tue, Mar 20, 2012 at 11:57 AM, Laurie
>>     <lauriekell at googlemail.com <mailto:lauriekell at googlemail.com>> wrote:
>>
>>         Send me the objects and code and I will try and run it.
>>
>>         Laurie
>>
>>
>>         On 03/20/2012 12:21 PM, Luis Ridao wrote:
>>>         FLR-help,
>>>
>>>         I'm having problems with a MSE.
>>>         I copy-paste the line which gives the error:
>>>
>>>         ...
>>>         ...
>>>         * print("OK to here1")*
>>>         *      ## project MP to TAC year with the HCR and geom mean
>>>         recr.*
>>>         *     
>>>         MP<-fwd(MP,ctrl=target,sr=list(model="mean",params=FLPar(1)),sr.residuals=SRrs)*
>>>         ...
>>>         ...
>>>
>>>         Year =2011
>>>         [1] "OK to here1"
>>>         *Error in .local(object, ctrl, ...) : 6*
>>>
>>>
>>>         My target is :
>>>         > *target*
>>>
>>>         Target
>>>           year min       val   max quantity season area unit spp
>>>         fleet metier rel.year
>>>         1 2011  NA 0.4550558    NA        f     NA   NA   NA  NA  
>>>          NA     NA       NA
>>>         2 2012  NA 0.3200000    NA        f     NA   NA   NA  NA  
>>>          NA     NA       NA
>>>         3 2012   0        NA 40000    catch     NA   NA   NA  NA  
>>>          NA     NA       NA
>>>         4 2012  NA 2.0000000 <tel:2.0000000>    NA        f     NA  
>>>         NA   NA  NA    NA     NA       NA
>>>           rel.season rel.area rel.unit
>>>         1         NA       NA       NA
>>>         2         NA       NA       NA
>>>         3         NA       NA       NA
>>>         4         NA       NA       NA
>>>             min        val        max
>>>           1         NA 4.5506e-01         NA
>>>           2         NA 3.2000e-01         NA
>>>           3 0.0000e+00         NA 4.0000e+04
>>>           4         NA         NA 2.0000e+00
>>>
>>>         the '*sr*' object :
>>>
>>>         > list(model="mean",params=FLPar(1))
>>>         $model
>>>         [1] "mean"
>>>
>>>         $params
>>>         An object of class "FLPar"
>>>         params
>>>         a
>>>         1
>>>         units:  NA
>>>
>>>         and the Sr residuals:
>>>
>>>         > *SRrs*
>>>         An object of class "FLQuant"
>>>         iters:  1000
>>>
>>>         , , unit = unique, season = all, area = unique
>>>
>>>              year
>>>         quant 2011     2012     2013
>>>           all 27355(0) 27355(0) 27355(0)
>>>
>>>         units:  NA
>>>
>>>         R version :               R.2.14.2
>>>         FLash Version:         2.5.0
>>>         FLAdvice Version:    1.0
>>>         FLXSA Version:       2.5
>>>         FLAssess Version:   2.5
>>>         FLCore Version:       2.5
>>>
>>>         Running on Ubuntu 10.04.04
>>>
>>>         Thanks in advance
>>>
>>>
>>>
>>>         _______________________________________________
>>>         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
>>
>>
>>         _______________________________________________
>>         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
>>
>>
>>
>>
>>     _______________________________________________
>>     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
>
>
>     _______________________________________________
>     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
>
>
>
>
> _______________________________________________
> flr-list mailing list
> flr-list at flr-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/flr-list/attachments/20120321/de6bd5e8/attachment.html>


More information about the flr-list mailing list