[FLR-list] Simulate for a4a not working

[Ernesto Jardim]

Hi Piera,

Sorry for the late reply. I was away last week. Usually that error is 
related with FLIndices and FLStocks that are not fully compatible. Can 
you send me those objects ?

Best

EJ

On 13/02/16 16:59, Piera Carpi wrote:
> Dear all,
>
> I’m experimenting an issue with the function “simulate” for a4a.
>
> Whenever I try to use it on windows, it gives me the following error:
>
> *Error in X %*% b : non-conformable arguments*
>
> Here is the code:
>
> #########
> > fit <- sca(stock = stk, srmodel=rmodel, fmodel = fmod, qmodel=qmod, 
> indices = ids, fit = "assessment")
> > stks <- stk + simulate(fit, 1000)
> Error in X %*% b : non-conformable arguments
> > summary(fit)
>
> "An object of class "a4aFitSA"
>
> Name: MySTK
> Description: Imported from a VPA file. ( 
> ~/Assessment_2015/Data/AllData/stk//stk7514.ndx ).  Thu Nov 26 
> 15:20:03 2015 + FLSAM: stk_7514
> Range:  min max pgroup minyear maxyear minfbar maxfbar
>  0 4 4 1975 2014 1 2
> Quant: age
>
> stock.n       : [ 5 40 1 1 1 1 ], units =  NA
> harvest       : [ 5 40 1 1 1 1 ], units =  f
> catch.n       : [ 5 40 1 1 1 1 ], units =  NA"
>
>
> #########
> Session Info:
> > 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  LC_CTYPE=English_United 
> Kingdom.1252
> [3] LC_MONETARY=English_United Kingdom.1252 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] FLSAM_1.0             FLBRP_2.5.20140822  reshape2_1.4.1
>  [4] plyr_1.8.3            ggplotFL_2.5.20160119 ggplot2_2.0.0
>  [7] FLAssess_2.5.20130716 FLash_2.5.2 FLa4a_0.700.1
> [10] coda_0.18-1           copula_0.999-14 triangle_0.10
> [13] mgcv_1.8-9            nlme_3.1-122  FLCore_2.5.20160107
> [16] Matrix_1.2-3          MASS_7.3-45 lattice_0.20-33
>
> loaded via a namespace (and not attached):
>  [1] Rcpp_0.12.3      magrittr_1.5 munsell_0.4.2    ADGofTest_0.3
>  [5] colorspace_1.2-6 pspline_1.0-17 stringr_1.0.0    tools_3.2.3
>  [9] grid_3.2.3       gtable_0.1.2 gridExtra_2.0.0  stringi_1.0-1
> [13] scales_0.3.0     gsl_1.9-10.1 stats4_3.2.3     stabledist_0.7-0
> [17] mvtnorm_1.0-5
>
> #########
>
> The same code is working on linux.
>
> Any idea what’s going on?
>
> Thanks,
>
> Piera
>
> ----------------------------------
> Piera Carpi, Ph.D.
> Istituto di Scienze del Mare (CNR-ISMAR)
> Sezione Pesca Marittima-Ancona
> Largo Fiera della Pesca
> 60125, Ancona - ITALY
> Tel: 071-2078858
> Fax: 071-55313
> e-mail: piera.carpi at an.ismar.cnr.it <mailto:piera.carpi at an.ismar.cnr.it>
>
>
>
>
>
> _______________________________________________
> 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/20160214/148dc94f/attachment.html>