alk.cpp

Go to the documentation of this file.

#include "statsLib.h"

dvar_matrix ageLengthKey(const dvar_vector& mu, const dvar_vector& sig, const dvector& x)
{
  RETURN_ARRAYS_INCREMENT();
  int i, j;
  dvariable z1;
  dvariable z2;
  int si,ni; si=mu.indexmin(); ni=mu.indexmax();
  int sj,nj; sj=x.indexmin(); nj=x.indexmax()-1;
  dvar_matrix pdf(si,ni,sj,nj);
  pdf.initialize();
  //double xs=0.5*(x[sj+1]-x[sj]);
  for(i=si;i<=ni;i++) //loop over ages
  {
     for(j=sj;j<=nj;j++) //loop over length bins
    {
      z1=(x(j)-mu(i))/sig(i);
      z2=(x(j+1)-mu(i))/sig(i);
      pdf(i,j)=cumd_norm(z2)-cumd_norm(z1);
    }//end nbins
    pdf(i)/=sum(pdf(i));
  }//end nage
  
  RETURN_ARRAYS_DECREMENT();
  return(pdf);
}


dmatrix ageLengthKey(const dvector& mu, const dvector& sig, const dvector& x)
{
  RETURN_ARRAYS_INCREMENT();
  int i, j;
  double z1;
  double z2;
  int si,ni; si=mu.indexmin(); ni=mu.indexmax();
  int sj,nj; sj=x.indexmin(); nj=x.indexmax()-1;
  dmatrix pdf(si,ni,sj,nj);
  pdf.initialize();
  //double xs=0.5*(x[sj+1]-x[sj]);
  for(i=si;i<=ni;i++) //loop over ages
  {
     for(j=sj;j<=nj;j++) //loop over length bins
    {
      z1=(x(j)-mu(i))/sig(i);
      z2=(x(j+1)-mu(i))/sig(i);
      pdf(i,j)=cumd_norm(z2)-cumd_norm(z1);
    }//end nbins
    pdf(i)/=sum(pdf(i));
  }//end nage
  
  RETURN_ARRAYS_DECREMENT();
  return(pdf);
}

dvar_matrix ALK(dvar_vector mu, dvar_vector sig, dvector x)
{
  RETURN_ARRAYS_INCREMENT();
  int i, j;
  dvariable z1;
  dvariable z2;
  int si,ni; si=mu.indexmin(); ni=mu.indexmax();
  int sj,nj; sj=x.indexmin(); nj=x.indexmax();
  dvar_matrix pdf(si,ni,sj,nj);
  pdf.initialize();
  double xs=0.5*(x[sj+1]-x[sj]);
  for(i=si;i<=ni;i++) //loop over ages
  {
     for(j=sj;j<=nj;j++) //loop over length bins
    {
      z1=((x(j)-xs)-mu(i))/sig(i);
      z2=((x(j)+xs)-mu(i))/sig(i);
      pdf(i,j)=cumd_norm(z2)-cumd_norm(z1);
    }//end nbins
    pdf(i)/=sum(pdf(i));
  }//end nage
  //pdf/=sum(pdf);
  RETURN_ARRAYS_DECREMENT();
  return(pdf);
}