SepVPA-methods           package:FLAssess           R Documentation

_S_e_p_a_r_a_b_l_e _V_i_r_t_u_a_l _P_o_p_u_l_a_t_i_o_n _A_n_a_l_y_s_i_s _i_n _F_L_R

_D_e_s_c_r_i_p_t_i_o_n:

     Performs Separable Virtual Population Analysis. Methods for
     function 'SepVPA' in package 'in Package `FLAssess''. Requires an
     object of class 'FLSepVPA.control' to be created.

_M_e_t_h_o_d_s:



     '_S_e_p_V_P_A<-_f_u_n_c_t_i_o_n(_s_t_o_c_k, _c_o_n_t_r_o_l=_F_L_S_e_p_V_P_A._c_o_n_t_r_o_l(), _r_e_f._h_a_r_v_e_s_t="_m_i_s_s_i_n_g", _f_r_a_t_i_o="_m_i_s_s_i_n_g", _f_i_t._p_l_u_s_g_r_o_u_p=_T_R_U_E, _d_e_s_c="",...)' 
          , where 'stock' is of type 'FLStock', 'control' is of type
          'FLSepVPA.control()', 'ref.harvest' and 'fratio' are numeric
          and 'fit.plusgroup' is Boolean.

_E_x_a_m_p_l_e_s:

     # Example based on ple4 dataset
     data(ple4)
     # Set up stock with correct dimensions
     my.stock <- FLStock(catch.n(ple4))
     my.stock@range["plusgroup"] <- 15
     #load catch data and mortality
     my.stock@catch.n <- ple4@catch.n
     my.stock@catch.n[my.stock@catch.n==0] <- 1
     my.stock@m <- ple4@m
     my.control <- FLSepVPA.control(sep.age = 5)
     # Set up in final year
     my.stock@stock.n[,"2001"] <- ple4@stock.n[,"2001"]
     # Run SepVPA
     my.stock.SepVPA <- SepVPA(my.stock, my.control, fit.plusgroup=TRUE)

