ADMB-users,<br><br>I'm having problems to fit a ssb-rec Ricker curve.<br>The problem is with the rec array indexing.<br>Recruitment strength is measured at age 3 so what I do is to create<br>two extras arrays holding the ssb-rec pairs with a 3-yrs lag.<br>
I thought this solution would make it easier but I get the following<br>error:<br><br>LRC $ ./catage6109_rec_mod_surv_ric<br>ok to here<br>Incompatible bounds in dvar_vector& dvar_vector::operator = (_CONST dvar_vector& t)<br>
<br>Snipet of the tpl file:<br><br>nyrs=49<br>.......................<br>  sdreport_vector ssb(1,nyrs)<br>  sdreport_vector ssb_m3(1,nyrs-3)      // SAME AS SSB with dimension (1,46)<br>  sdreport_vector rec(1,nyrs)<br>  sdreport_vector rec_m3(1,nyrs-3)      // SAME AS REC with dimension (1,46)<br>
  sdreport_vector pred_rec(1,nyrs)<br>  sdreport_vector pred_rec_m3(1,nyrs-3)    // SAME AS PRED_REC with dimension (1,46)<br>.........................<br><br>FUNCTION ssb_rec_ricker<br>  ssb=rowsum(elem_prod(elem_prod(N,weight_at_age),mat_at_age))/1000;<br>
  ssb_m3=ssb.sub(1,46);<br>  rec=column(N,1)/1000;<br>  cout<<"ok to here"<<endl;          // HERE IS THE<br>  rec_m3=rec.sub(3,49);              // PROBLEM <br><br>  for (int i=1;i<=nyrs-3;i++)<br>
  {<br>  pred_rec_m3(i)=exp(log_alfa)*ssb_m3(i)*exp(-exp(log_beta)*ssb_m3(i));<br>  }  <br><br>............................<br>
<br>The valid index for rec is (1,49)<br>Why can't I subset it to (3,49)?<br><br>Thanks in advance,<br>Luis<br>