k_factors Class Reference

Class for calculating k-factors. More...

#include <Partridge_All.h>

Public Member Functions
void	calcreal_1 ()
void	calcreal_2 ()
void	calcreal_3 ()
void	calcDickPottsk ()
void	Output_kfactors ()
void	CreateInfoDatabaseFile ()
	k_factors ()
	~k_factors ()

Public Attributes
population_attributes *	m_ThisYear
population_attributes *	m_LastYear

Protected Member Functions
void	reset ()
void	DumpInfoDatabase1 (int, int, double)
void	DumpInfoDatabase ()

Protected Attributes
FILE *	kfactorsfile
double	m_rk1
double	m_rk2
double	m_rk3
double	m_rk4
double	m_rk5
double	m_rk6
double	m_dpk1
double	m_dpk2
double	m_dpk3
double	m_dpk4
double	m_dpk5

Detailed Description

Class for calculating k-factors.

Definition at line 235 of file Partridge_All.h.

Constructor & Destructor Documentation

k_factors()

k_factors::k_factors

(

)

Definition at line 2858 of file Partridge_All.cpp.

{
  kfactorsfile = fopen("k_factors_output.txt", "w" );
  if ( !kfactorsfile ) {
      g_msg->Warn( "k_factors::k_factors()","Cannot open k_factors_output.txt" );
      exit( 1 );
  }
  m_ThisYear = new population_attributes;
  m_LastYear = new population_attributes;
  reset();
}

References g_msg, kfactorsfile, m_LastYear, m_ThisYear, reset(), and MapErrorMsg::Warn().

~k_factors()

k_factors::~k_factors

(

)

Definition at line 2871 of file Partridge_All.cpp.

{
  fclose( kfactorsfile );
  delete m_ThisYear;
  delete m_LastYear;
}

References kfactorsfile, m_LastYear, and m_ThisYear.

Member Function Documentation

calcDickPottsk()

void k_factors::calcDickPottsk

(

)

Definition at line 3026 of file Partridge_All.cpp.

{
  // Dick makes some assumptions such it is assumed old males seen on 1st sept
  // is the same number as the number of pairs in Spring.
  // He also uses a calculation for the number of chicks hatched based on
  // the geometric mean of brood size on 1st September.
  // k2 assumes the differences between old males and old females is female on nest mortality
 
  // Requires the following values to be set before entry:
  // 1) GeoMean brood size, 1st sept.
  // 2) No. old males
  // 3) No. old females
  // 4) No. Males Sept
  // 5) No. Females Sept
  // 6) No. Females Sept of last year
  // 7) No. chicks at 1st Sept.
 
  // Calculate HatchedChicks
  double HatchedChicks;
  double CSR;
  if ( m_LastYear->m_geomeanBroodSize < 10.0 ) CSR = 0.03665 * ( pow( ( double )m_LastYear->m_geomeanBroodSize, ( double )1.293 ) );
  else
    CSR = m_LastYear->m_geomeanBroodSize / 13.84;
  if ( CSR > 0 ) HatchedChicks = m_LastYear->m_NoChicksSept / CSR;
  else
    HatchedChicks = 0;
  //
  double Broods = HatchedChicks / 13.8; // 13.8 is the mean brood size used by Dick, but we use 14
  // k2 is assuming old males is the same as pairs
  m_dpk2 = log10( ( double )( 1 + m_LastYear->m_NoOldMales ) / ( double )( 1 + m_LastYear->m_NoOldFemales ) );
  // k1 assumes old males in sept is the same as no. pairs in spring
  m_dpk1 = log10( ( double )( 1 + m_LastYear->m_NoOldMales ) / ( 1 + Broods ) ) - m_dpk2;
  // k3 contains a calculated number of chicks at six weeks based on the number
  // of chicks in Sept and brood size.
  m_dpk3 = log10( ( double )( 1 + HatchedChicks ) / ( double )( 1 + m_LastYear->m_NoChicksSept ) );
  // k4 Shooting numbers recorded so should be solid
  double NoBirdsAfterShoot = ( ( m_LastYear->m_NoMalesSept + m_LastYear->m_NoFemalesSept ) - m_LastYear->m_NoShotBirds );
  m_dpk4 = log10( ( double )( ( 1 + m_LastYear->m_NoMalesSept + m_LastYear->m_NoFemalesSept ) ) / ( double )( 1 + NoBirdsAfterShoot ) );
  // k5 - tricky one. Assumes 1:1 sex ratio. Assumes that the number of females
  // at the end of year 1 is equal to the number of breeding males and this is
  // equal to the number of old males in year 2. Hence the difference is the number
  // of birds that die of everything else except hen mortality.
  m_dpk5 = log10( ( double )( 1 + (m_LastYear->m_NoFemalesSept - m_LastYear->m_NoShotFemales)) / ( double )( 1 + m_ThisYear->m_NoOldMales ) );
  DumpInfoDatabase1( ( int )Broods, ( int )HatchedChicks, CSR );
 
  m_LastYear->copyself(m_ThisYear);
  //m_ThisYear->reset();  // ***CJT*** 6/7/2010 removed because otherwise our counts are destroyed - why was this here?
}

References population_attributes::copyself(), DumpInfoDatabase1(), m_dpk1, m_dpk2, m_dpk3, m_dpk4, m_dpk5, population_attributes::m_geomeanBroodSize, m_LastYear, population_attributes::m_NoChicksSept, population_attributes::m_NoFemalesSept, population_attributes::m_NoMalesSept, population_attributes::m_NoOldFemales, population_attributes::m_NoOldMales, population_attributes::m_NoShotBirds, population_attributes::m_NoShotFemales, and m_ThisYear.

Referenced by Partridge_Population_Manager::DoFirst().

calcreal_1()

void k_factors::calcreal_1

(

)

Definition at line 2994 of file Partridge_All.cpp.

{
  // These can be calculated on Sept 1st.
  double anum = 1 + ( m_ThisYear->m_NoClutches - m_ThisYear->m_NoHensPredated );
  m_rk2 = log10( ( 1 + m_ThisYear->m_NoClutches ) / anum );
  anum = 1 + m_ThisYear->m_NoClutchesHatched;
  m_rk1 = log10( ( 1 + m_ThisYear->m_NoClutches ) / anum ) - m_rk2;
  anum = 1 + m_ThisYear->m_NoChicksSept;
  m_rk3 = log10( ( 1 + m_ThisYear->m_NoChicksHatched ) / anum );
}

References population_attributes::m_NoChicksHatched, population_attributes::m_NoChicksSept, population_attributes::m_NoClutches, population_attributes::m_NoClutchesHatched, population_attributes::m_NoHensPredated, m_rk1, m_rk2, m_rk3, and m_ThisYear.

Referenced by Partridge_Population_Manager::DoFirst().

calcreal_2()

void k_factors::calcreal_2

(

)

Definition at line 3006 of file Partridge_All.cpp.

{
  // This can only be calculated after shooting
  double NoBirdsAfterShoot = ( ( m_ThisYear->m_NoMalesSept + m_ThisYear->m_NoFemalesSept ) - m_ThisYear->m_NoShotBirds );
  double anum = 1 + NoBirdsAfterShoot;
  m_rk4 = log10( ( 1 + m_ThisYear->m_NoMalesSept + m_ThisYear->m_NoFemalesSept ) / anum );
}

References population_attributes::m_NoFemalesSept, population_attributes::m_NoMalesSept, population_attributes::m_NoShotBirds, m_rk4, and m_ThisYear.

Referenced by Partridge_Population_Manager::DoFirst().

calcreal_3()

void k_factors::calcreal_3

(

)

Definition at line 3015 of file Partridge_All.cpp.

{
  // K5 is the number of males and females in Sept/those that survive to
  // breeding time. This means k5 can only be calculated in April.
  double anum = 1 + m_ThisYear->m_NoFemalesSept - m_ThisYear->m_NoShotBirds / 2;
  m_rk5 = log10( anum / ( double )( 1 + m_ThisYear->m_NoFemalesApril ) );
  anum = 1 + m_ThisYear->m_NoMalesSept + m_ThisYear->m_NoFemalesSept - m_ThisYear->m_NoShotBirds;
  m_rk6 = log10( anum / ( double )( 1 + m_ThisYear->m_NoBirdsApril ) );
}

References population_attributes::m_NoBirdsApril, population_attributes::m_NoFemalesApril, population_attributes::m_NoFemalesSept, population_attributes::m_NoMalesSept, population_attributes::m_NoShotBirds, m_rk5, m_rk6, and m_ThisYear.

Referenced by Partridge_Population_Manager::DoFirst().

CreateInfoDatabaseFile()

void k_factors::CreateInfoDatabaseFile

(

)

Open a k-factor output file.

Definition at line 3145 of file Partridge_All.cpp.

{
  FILE * idb = fopen("k_facInfoBase.txt", "w" );
  if (!idb) {
      g_msg->Warn("k_factors::CreateInfoDatabaseFile","Cannot open k_facInfoBase.txt" );
      exit( 1 );
  }
  fclose( idb );
}

References g_msg, and MapErrorMsg::Warn().

Referenced by Partridge_Population_Manager::Init().

DumpInfoDatabase()

void k_factors::DumpInfoDatabase ( )

protected

Definition at line 3088 of file Partridge_All.cpp.

{
  FILE * idb = fopen("k_facInfoBase.txt", "a" );
  if (!idb) {
      g_msg->Warn( "k_factors::DumpInfoDatabase","Cannot open k_facInfoBase.txt" );
      exit( 1 );
  }
  fprintf( idb, "%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%g\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\n",
     m_ThisYear->m_NoClutches,m_ThisYear->m_NoClutchesHatched, m_ThisYear->m_NoHensPredated, m_ThisYear->m_NoChicksSixWeeks, m_ThisYear->m_NoChicksSept, m_ThisYear->m_NoMalesSept, m_ThisYear->m_NoFemalesSept,  m_ThisYear->m_NoDeadFemales, m_ThisYear->m_NoDeadMales, m_ThisYear->m_NoShotBirds,
     m_ThisYear->m_NoShotFemales, m_ThisYear->m_NoChicksHatched, m_ThisYear->m_NoOldMales, m_ThisYear->m_NoOldFemales, m_ThisYear->m_geomeanBroodSize,
     m_ThisYear->m_NoClutchesPredated, m_ThisYear->m_NoFemalesApril, m_ThisYear->m_NoTerritorialFemalesApril, m_ThisYear->m_NoDeadChicks, m_ThisYear->m_NoDeadClutches,
     m_ThisYear->m_NoClutchesReplacements, m_ThisYear->m_NoAgDeadClutches, m_ThisYear->m_NoPairsApril, m_ThisYear->m_NoMalesApril,m_ThisYear->m_NoTerritorialFemalesMay, m_ThisYear->m_NoNonTerritorialFemalesMay,
     m_ThisYear->m_NoMalesDec, m_ThisYear->m_NoFemalesDec, m_ThisYear->m_NoChick1sAug, m_ThisYear->m_NoChick2sAug,
     m_ThisYear->m_NoMalesAug, m_ThisYear->m_NoFemsAug, m_ThisYear->m_NoStarvedChicks );
 
  /*
       1m_NoClutches,
       2m_NoClutchesHatched,
       3m_NoHensPredated,
       4m_NoChicksSixWeeks,
       5m_NoChicksSept,
       6m_NoMalesSept,
       7m_NoFemalesSept,
       8m_NoDeadFemales,
       9m_NoDeadMales,
       10m_NoShotBirds,
       11m_NoChicksHatched,
       12m_NoOldMales,
       13m_NoOldFemales,
       14m_geomeanBroodSize,
       15m_NoClutchesPredated,
       16m_NoFemalesApril,
       17m_NoTerritorialFemalesApril,
       18m_NoDeadChicks,
       19m_NoDeadClutches,
       20m_NoClutchesReplacements,
       21m_NoAgDeadClutches,
       22m_NoPairsApril,
       23m_NoMalesApril,
       24m_NoTerritorialFemalesMay,
       25m_NoNonTerritorialFemalesMay
       26m_NoMalesDec
       27m_NoFemalesDec
       28m_NoChick1sAug,
       29m_NoChick2sAug,
       30m_NoMalesAug,
       31m_NoFemsAug
       32m_NoStarvedChicks
       */
 
  fclose( idb );
}

References g_msg, population_attributes::m_geomeanBroodSize, population_attributes::m_NoAgDeadClutches, population_attributes::m_NoChick1sAug, population_attributes::m_NoChick2sAug, population_attributes::m_NoChicksHatched, population_attributes::m_NoChicksSept, population_attributes::m_NoChicksSixWeeks, population_attributes::m_NoClutches, population_attributes::m_NoClutchesHatched, population_attributes::m_NoClutchesPredated, population_attributes::m_NoClutchesReplacements, population_attributes::m_NoDeadChicks, population_attributes::m_NoDeadClutches, population_attributes::m_NoDeadFemales, population_attributes::m_NoDeadMales, population_attributes::m_NoFemalesApril, population_attributes::m_NoFemalesDec, population_attributes::m_NoFemalesSept, population_attributes::m_NoFemsAug, population_attributes::m_NoHensPredated, population_attributes::m_NoMalesApril, population_attributes::m_NoMalesAug, population_attributes::m_NoMalesDec, population_attributes::m_NoMalesSept, population_attributes::m_NoNonTerritorialFemalesMay, population_attributes::m_NoOldFemales, population_attributes::m_NoOldMales, population_attributes::m_NoPairsApril, population_attributes::m_NoShotBirds, population_attributes::m_NoShotFemales, population_attributes::m_NoStarvedChicks, population_attributes::m_NoTerritorialFemalesApril, population_attributes::m_NoTerritorialFemalesMay, m_ThisYear, and MapErrorMsg::Warn().

Referenced by Output_kfactors().

DumpInfoDatabase1()

void k_factors::DumpInfoDatabase1 ( int  Broods, int  HC, double  CSR  )

protected

Definition at line 3076 of file Partridge_All.cpp.

{
  FILE * idb = fopen("k_facInfoBase.txt", "a" );
  if ( !idb ) {
      g_msg->Warn( "k_factors::DumpInfoDatabase1","Cannot open k_facInfoBase.txt" );
      exit( 1 );
  }
  fprintf( idb, "%d\t%d\t%d\t%g\t", HC, Broods, m_ThisYear->m_NoOldMales, CSR );
  fclose( idb );
}

References g_msg, population_attributes::m_NoOldMales, m_ThisYear, and MapErrorMsg::Warn().

Referenced by calcDickPottsk().

Output_kfactors()

void k_factors::Output_kfactors

(

)

Definition at line 2893 of file Partridge_All.cpp.

{
  //calcDickPottsk();
  fprintf( kfactorsfile, "%g\t%g\t%g\t%g\t%g\t%g\t%g\t%g\t%g\t%g\t%g\n", m_rk1, m_rk2, m_rk3, m_rk4, m_rk5, m_rk6, m_dpk1,
       m_dpk2, m_dpk3, m_dpk4, m_dpk5 );
  fflush( kfactorsfile );
  DumpInfoDatabase();
  reset();
}

References DumpInfoDatabase(), kfactorsfile, m_dpk1, m_dpk2, m_dpk3, m_dpk4, m_dpk5, m_rk1, m_rk2, m_rk3, m_rk4, m_rk5, m_rk6, and reset().

Referenced by Partridge_Population_Manager::DoFirst().

reset()

void k_factors::reset ( )

protected

Definition at line 2978 of file Partridge_All.cpp.

{
  m_ThisYear-> reset();
  m_rk1 = 0;
  m_rk2 = 0;
  m_rk3 = 0;
  m_rk4 = 0;
  m_rk5 = 0;
  m_dpk1 = 0;
  m_dpk2 = 0;
  m_dpk3 = 0;
  m_dpk4 = 0;
  m_dpk5 = 0;
 
}

References m_dpk1, m_dpk2, m_dpk3, m_dpk4, m_dpk5, m_rk1, m_rk2, m_rk3, m_rk4, m_rk5, and m_ThisYear.

Referenced by k_factors(), and Output_kfactors().

Member Data Documentation

kfactorsfile

FILE* k_factors::kfactorsfile

protected

This class contains the functionality required to make two types of k-factors the first type is the 'real' k-factor, the second is the way Dick Potts calculated his.
Not yet functional/tested - 5thOct2008

Definition at line 243 of file Partridge_All.h.

Referenced by k_factors(), Output_kfactors(), and ~k_factors().

m_dpk1

double k_factors::m_dpk1

protected

Definition at line 245 of file Partridge_All.h.

Referenced by calcDickPottsk(), Output_kfactors(), and reset().

m_dpk2

double k_factors::m_dpk2

protected

Definition at line 245 of file Partridge_All.h.

Referenced by calcDickPottsk(), Output_kfactors(), and reset().

m_dpk3

double k_factors::m_dpk3

protected

Definition at line 245 of file Partridge_All.h.

Referenced by calcDickPottsk(), Output_kfactors(), and reset().

m_dpk4

double k_factors::m_dpk4

protected

Definition at line 245 of file Partridge_All.h.

Referenced by calcDickPottsk(), Output_kfactors(), and reset().

m_dpk5

double k_factors::m_dpk5

protected

Definition at line 245 of file Partridge_All.h.

Referenced by calcDickPottsk(), Output_kfactors(), and reset().

m_LastYear

population_attributes* k_factors::m_LastYear

Definition at line 252 of file Partridge_All.h.

Referenced by calcDickPottsk(), k_factors(), and ~k_factors().

m_rk1

double k_factors::m_rk1

protected

Definition at line 244 of file Partridge_All.h.

Referenced by calcreal_1(), Output_kfactors(), and reset().

m_rk2

double k_factors::m_rk2

protected

Definition at line 244 of file Partridge_All.h.

Referenced by calcreal_1(), Output_kfactors(), and reset().

m_rk3

double k_factors::m_rk3

protected

Definition at line 244 of file Partridge_All.h.

Referenced by calcreal_1(), Output_kfactors(), and reset().

m_rk4

double k_factors::m_rk4

protected

Definition at line 244 of file Partridge_All.h.

Referenced by calcreal_2(), Output_kfactors(), and reset().

m_rk5

double k_factors::m_rk5

protected

Definition at line 244 of file Partridge_All.h.

Referenced by calcreal_3(), Output_kfactors(), and reset().

m_rk6

double k_factors::m_rk6

protected

Definition at line 244 of file Partridge_All.h.

Referenced by calcreal_3(), and Output_kfactors().

m_ThisYear

population_attributes* k_factors::m_ThisYear

Definition at line 251 of file Partridge_All.h.

Referenced by Partridge_Clutch::AgDying(), calcDickPottsk(), calcreal_1(), calcreal_2(), calcreal_3(), Partridge_Chick::ChDeveloping(), Partridge_Chick::ChDying(), Partridge_Clutch::ClDying(), Partridge_Population_Manager::CreateObjects(), Partridge_Population_Manager::DoFirst(), DumpInfoDatabase(), DumpInfoDatabase1(), Partridge_Female::FDying(), Partridge_Female::FIncubating(), Partridge_Female::FLaying(), Partridge_Female::FMakingNest(), Partridge_Population_Manager::Hunting(), Partridge_Population_Manager::HuntingDifferentiatedBeetleBankArea(), Partridge_Population_Manager::HuntingGrid(), Partridge_Population_Manager::Init(), k_factors(), Partridge_Male::MDying(), reset(), and ~k_factors().

---

The documentation for this class was generated from the following files:

• C:/Repo/ALMaSS_All/Partridge/Partridge_All.h
• C:/Repo/ALMaSS_All/Partridge/Partridge_All.cpp