[FLR-list] Simulate for a4a not working

[Piera Carpi]

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



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