TTQuantiFH Class Reference

TTQuantiFH. More...

#include <ttquanti.h>

Inheritance diagram for TTQuantiFH:

Collaboration diagram for TTQuantiFH:

Public Member Functions
	TTQuantiFH (TProtoQuanti *T)
virtual	~TTQuantiFH ()
void	setOutputOption (string opt)
virtual void	FHwrite ()
void	write_TABLE ()
void	print (ofstream &FH, sex_t SEX, age_idx Ax, unsigned int print_gene, bool print_genotype, bool print_additive_genotype)
void	write_PLINK ()
void	print_PLINK_PED (ofstream &FH, age_idx Ax, Patch *patch)
void	print_PLINK_FAM (ofstream &FH, age_idx Ax, Patch *patch)
virtual void	FHread (string &filename)
Public Member Functions inherited from TraitFileHandler< TProtoQuanti >
	TraitFileHandler (TProtoQuanti *trait_proto, const char *ext)
virtual	~TraitFileHandler ()
virtual void	FHwrite ()=0
virtual void	FHread (string &filename)=0
virtual void	set (bool rpl_per, bool gen_per, int rpl_occ, int gen_occ, int rank, string path, TProtoQuanti *trait_proto)
Public Member Functions inherited from FileHandler
	FileHandler (const char *ext)
virtual	~FileHandler ()
virtual void	init ()
	Called by notifier during simulation setup, performs file checking. More...
virtual vector< string >	ifExist ()
	Checks if any file associated with the current file name already exists on disk. More...
virtual void	set (bool rpl_per, bool gen_per, int rpl_occ, int gen_occ, int rank, string path)
	Sets the hanlder parameters. More...
virtual void	set_multi (bool rpl_per, bool gen_per, int rpl_occ, TMatrix *Occ, string path)
virtual void	FHwrite ()=0
	Default behavior of the class, called by Handler::update(). More...
virtual void	FHread (string &filename)=0
	Default input function. More...
virtual void	update ()
	Updates the inner replicate and generation counters and calls FHwrite if needed by the the periodicity of the file. More...
Metapop *	get_pop_ptr ()
	Returns the pointer to the current metapop through the FileServices interface. More...
void	set_pop_ptr (Metapop *pop_ptr)
FileServices *	get_service ()
	Returns pointer to the FileServices. More...
void	set_service (FileServices *srv)
std::string &	get_path ()
void	set_path ()
std::string &	get_extension ()
void	set_extension (const char *ext)
std::string &	get_filename ()
	Builds and returns the current file name depending on the periodicity of the file. More...
bool	get_isInputHandler ()
void	set_isInputHandler (bool val)
bool	get_isReplicatePeriodic ()
void	set_isReplicatePeriodic (bool val)
unsigned int	get_ReplicateOccurrence ()
void	set_ReplicateOccurrence (unsigned int val)
bool	get_isGenerationPeriodic ()
void	set_isGenerationPeriodic (bool val)
unsigned int	get_GenerationOccurrence ()
void	set_GenerationOccurrence (unsigned int val)
unsigned int	get_ExecRank ()
	unused yet... More...
void	set_ExecRank (int val)
TMatrix *	get_OccMatrix ()
void	set_OccMatrix (TMatrix *occ)
bool	get_isMasterExec ()
void	set_isMasterExec (bool is)
Public Member Functions inherited from Handler
virtual void	init ()=0
	Inits state. More...
virtual void	update ()=0
	Updates the handler state. More...
virtual	~Handler ()

Private Attributes
string	_output_option

Additional Inherited Members
Protected Attributes inherited from TraitFileHandler< TProtoQuanti >
TProtoQuanti *	_FHLinkedTrait
int	_FHLinkedTraitIndex
Protected Attributes inherited from FileHandler
Metapop *	_pop
	Pointer to the current metapop, set during initialization within the init function. More...

Detailed Description

TTQuantiFH.

Constructor & Destructor Documentation

TTQuantiFH()

TTQuantiFH::TTQuantiFH ( TProtoQuanti *  T )

inline

: TraitFileHandler<TProtoQuanti>(T,".quanti") {}

~TTQuantiFH()

virtual TTQuantiFH::~TTQuantiFH ( )

inlinevirtual

{}

Member Function Documentation

FHread()

void TTQuantiFH::FHread ( string &  filename )

virtual

Implements FileHandler.

{
  unsigned int nb_locus = _FHLinkedTrait->get_nb_locus();
  unsigned int nb_trait = _FHLinkedTrait->get_nb_traits();
  bool has_genotype = (!_FHLinkedTrait->get_env_var().empty());
 
  unsigned int patchNbr = _pop->getPatchNbr();
 
  TTQuanti* trait;
 
  ifstream FILE(filename.c_str(),ios::in);
 
  if(!FILE) fatal("could not open QUANTI input file \"%s\"\n",filename.c_str());
 
  unsigned int age, sex, ped, origin, dummy;
  unsigned int loc;
  age_idx stagex;
  unsigned int pop;
  age_idx agex;
  Individual *ind;
  double all0, all1, pheno, geno;
  double *seq[2];
  seq[0] = new double [nb_locus*nb_trait];
  seq[1] = new double [nb_locus*nb_trait];
 
  int lnbr = 1;
 
  // swallow the header
  FILE.ignore(std::numeric_limits<std::streamsize>::max(), '\n');
 
  string str;
 
  while(FILE>>pop) {
 
    if(FILE.bad()) {
      error("Reading input genotypes from \"%s\" failed at line %i.\n",filename.c_str(), lnbr);
      FILE.clear();
      FILE >> str;
      fatal("Expecting population number as first element on the line, but received this: %s \n", str.c_str());
    }
 
    if(pop > patchNbr)
      fatal("Patch number found in file exceeds number of patches in the population.\n");
 
    for(unsigned int k = 0; k < nb_trait; ++k) {
      for(unsigned int i = 0; i < nb_locus; ++i) {
 
        loc = i * nb_trait + k;
 
        FILE>>all0>>all1;
 
        seq[0][loc] = all0;
        seq[1][loc] = all1;
      }
    }
 
 
    for(unsigned int k = 0; k < nb_trait; k++) {
      FILE >> pheno;                   //read phenotypic value
      if(has_genotype) FILE >> geno;   // genotypic values, only present if trait is influenced by env. var.
    }
 
    // read extra fields
    FILE >> age >> sex >> origin >> ped >> dummy >> dummy >> dummy >> dummy;
 
    stagex = age_idx(age);
 
    ind = _pop->makeNewIndividual(0, 0, sex_t(sex), origin - 1);
    ind->setPedigreeClass((unsigned char)ped);
    ind->setAge(age);
    trait = dynamic_cast<TTQuanti*> (ind->getTrait(_FHLinkedTraitIndex));
    trait->set_sequence((void**)seq);
    trait->set_value();
 
    //      ind->show_up();
 
    _pop->getPatch(pop-1)->add(sex_t(sex), stagex, ind);
 
    lnbr++;
  }
 
 
  FILE.close();
 
  delete seq[0];
  delete seq[1];
}

References TraitFileHandler< TProtoQuanti >::_FHLinkedTrait, TraitFileHandler< TProtoQuanti >::_FHLinkedTraitIndex, FileHandler::_pop, Patch::add(), error(), fatal(), TProtoQuanti::get_env_var(), TProtoQuanti::get_nb_locus(), TProtoQuanti::get_nb_traits(), Metapop::getPatch(), Metapop::getPatchNbr(), Individual::getTrait(), IndFactory::makeNewIndividual(), TTQuanti::set_sequence(), TTQuanti::set_value(), Individual::setAge(), and Individual::setPedigreeClass().

FHwrite()

void TTQuantiFH::FHwrite ( )

virtual

Implements TraitFileHandler< TProtoQuanti >.

{
  if (!_pop->isAlive()) return;
 
  if(!get_service()) fatal("link to file services not set!!\n");
 
  // sets the pop ptr to a sub sampled pop, if sub sampling happened
  _pop = get_service()->getSampledPop();
 
  if(_output_option == "PLINK" || _output_option == "plink")
    write_PLINK();
  else
    write_TABLE();
 
 
  // reset to pop ptr to main pop:
  _pop = get_service()->get_pop_ptr();
 
}

References _output_option, FileHandler::_pop, fatal(), FileServices::get_pop_ptr(), FileHandler::get_service(), FileServices::getSampledPop(), Metapop::isAlive(), write_PLINK(), and write_TABLE().

print()

void TTQuantiFH::print	(	ofstream &	FH,
		sex_t	SEX,
		age_idx	Ax,
		unsigned int	print_gene,
		bool	print_genotype,
		bool	print_additive_genotype
	)

{
  int patchNbr = _pop->getPatchNbr();
  Patch* current_patch;
  Individual* ind;
  TTQuanti* trait;
  double* Tval;
  double **genes;
  double **allele_values = _FHLinkedTrait->get_allele_values(); //
  unsigned int loc, nb_trait=_FHLinkedTrait->get_nb_traits(), nb_locus = _FHLinkedTrait->get_nb_locus();
 
  for(int i = 0; i < patchNbr; i++) {
 
    current_patch = _pop->getPatch(i);
 
    for(unsigned int j = 0, size = current_patch->size(SEX, Ax); j < size; j++) {
 
      ind = current_patch->get(SEX, Ax, j);
      trait = dynamic_cast<TTQuanti*> (ind->getTrait(_FHLinkedTraitIndex));
 
      FH<<current_patch->getID() +1<<" ";
      Tval = (double*)trait->getValue();
 
      if(print_gene){
        genes = (double**)trait->get_sequence();
 
        FH.precision(6);
 
        for(unsigned int k = 0; k < nb_trait; k++) {
          for(unsigned int l = 0; l < nb_locus; l++) {
            loc = l * nb_trait + k;
 
            if(print_gene == 1) {
 
              FH<<genes[FEM][loc]<<" "<<genes[MAL][loc]<<" ";
 
            } else {
 
              FH << int(genes[0][loc] !=  allele_values[l][0]) + int(genes[1][loc] !=  allele_values[l][0]) <<" ";
 
            }
          }
        }
      }
 
      FH.precision(4);
      for(unsigned int k = 0; k < nb_trait; k++) {
        FH<<Tval[k]<<" ";
        if(print_genotype) FH << trait->get_full_genotype(k) << " ";
        if(print_additive_genotype) FH << trait->get_additive_genotype(k) << " ";
      }
 
      FH<<Ax<<" "<<ind->getSex()<<" "<<ind->getHome()+1<<" "<<ind->getPedigreeClass()<<" "
          << (ind->getFather() && ind->getMother() ?
              (ind->getFather()->getHome()!= current_patch->getID() ) + (ind->getMother()->getHome()!= current_patch->getID() ) : 0)
              <<" "<<ind->getFatherID()<<" "<<ind->getMotherID()<<" "<<ind->getID()<<std::endl;
    }
  }
}

References TraitFileHandler< TProtoQuanti >::_FHLinkedTrait, TraitFileHandler< TProtoQuanti >::_FHLinkedTraitIndex, FileHandler::_pop, FEM, Patch::get(), TTQuanti::get_additive_genotype(), TProtoQuanti::get_allele_values(), TTQuanti::get_full_genotype(), TProtoQuanti::get_nb_locus(), TProtoQuanti::get_nb_traits(), TTQuanti::get_sequence(), Individual::getFather(), Individual::getFatherID(), Individual::getHome(), Individual::getID(), Patch::getID(), Individual::getMother(), Individual::getMotherID(), Metapop::getPatch(), Metapop::getPatchNbr(), Individual::getPedigreeClass(), Individual::getSex(), Individual::getTrait(), TTQuanti::getValue(), MAL, and Patch::size().

Referenced by write_TABLE().

Here is the caller graph for this function:

print_PLINK_FAM()

void TTQuantiFH::print_PLINK_FAM	(	ofstream &	FH,
		age_idx	Ax,
		Patch *	patch
	)

{
  Individual *ind;
  unsigned int nb_trait = _FHLinkedTrait->get_nb_traits();
 
  // SPECIFICATION FOR THE .fam FILE = 6 first values in .ped files:
  //  .fam (PLINK sample information file)
  //
  //  Sample information file accompanying a .bed binary genotype table.
  //  Also generated by "--recode lgen" and "--recode rlist".
  //
  //  A text file with no header line, and one line per sample with the following six fields:
  //
  //      1. Family ID ('FID')
  //      2. Within-family ID ('IID'; cannot be '0')
  //      3. Within-family ID of father ('0' if father isn't in dataset)
  //      4. Within-family ID of mother ('0' if mother isn't in dataset)
  //      5. Sex code ('1' = male, '2' = female, '0' = unknown)
  //      6. Phenotype value ('1' = control, '2' = case, '-9'/'0'/non-numeric = missing data if case/control)
  //
  //  If there are any numeric phenotype values other than {-9, 0, 1, 2}, the phenotype is interpreted as a
  //  quantitative trait instead of case/control status. In this case, -9 normally still designates a missing
  //  phenotype; use --missing-phenotype if this is problematic.
 
 
  // here, we save the phenotypic value of the first trait by default
  // the complete set of phenotypic values for all traits is saved in the .fam file
 
  for (unsigned int j = 0; j < patch->size(FEM, Ax); ++j) {
 
    ind = patch->get(FEM, Ax, j);
 
    FH<<"fam"<< ind->getHome()+1
        <<" "<< ind->getID()+1
        <<" 0 0" //<<ind->getFatherID()<<" "<<ind->getMotherID() //not in file for adults
        <<" 2";
 
    for(unsigned int k = 0; k < nb_trait; ++k) {
      FH << " " << ((double*)ind->getTraitValue(_FHLinkedTraitIndex))[k];
    }
 
    FH <<std::endl;
 
  }
 
  for (unsigned int j = 0; j < patch->size(MAL, Ax); ++j) {
 
    ind = patch->get(MAL, Ax, j);
 
    FH<<"fam"<< ind->getHome()+1
        <<" "<< ind->getID()+1
        <<" 0 0"  //<<ind->getFatherID()<<" "<<ind->getMotherID()
        <<" 1";
 
    for(unsigned int k = 0; k < nb_trait; ++k) {
      FH << " " << ((double*)ind->getTraitValue(_FHLinkedTraitIndex))[k];
    }
 
    FH<<std::endl;
 
  }
}

References TraitFileHandler< TProtoQuanti >::_FHLinkedTrait, TraitFileHandler< TProtoQuanti >::_FHLinkedTraitIndex, FEM, Patch::get(), TProtoQuanti::get_nb_traits(), Individual::getHome(), Individual::getID(), Individual::getTraitValue(), MAL, and Patch::size().

Referenced by write_PLINK().

Here is the caller graph for this function:

print_PLINK_PED()

void TTQuantiFH::print_PLINK_PED	(	ofstream &	FH,
		age_idx	Ax,
		Patch *	patch
	)

{
  Individual *ind;
  double** seq;
  char BASE[2] = {'A','G'};
  double** ref_all = _FHLinkedTrait->get_allele_values();
  unsigned int loc;
  unsigned int nb_trait = _FHLinkedTrait->get_nb_traits();
  unsigned int nb_locus = _FHLinkedTrait->get_nb_locus();
 
  // SPECIFICATION FOR THE .fam FILE = 6 first values in .ped files:
  //  .fam (PLINK sample information file)
  //
  //  Sample information file accompanying a .bed binary genotype table.
  //  Also generated by "--recode lgen" and "--recode rlist".
  //
  //  A text file with no header line, and one line per sample with the following six fields:
  //
  //      1. Family ID ('FID')
  //      2. Within-family ID ('IID'; cannot be '0')
  //      3. Within-family ID of father ('0' if father isn't in dataset)
  //      4. Within-family ID of mother ('0' if mother isn't in dataset)
  //      5. Sex code ('1' = male, '2' = female, '0' = unknown)
  //      6. Phenotype value ('1' = control, '2' = case, '-9'/'0'/non-numeric = missing data if case/control)
  //
  //  If there are any numeric phenotype values other than {-9, 0, 1, 2}, the phenotype is interpreted as a quantitative trait instead of case/control status. In this case, -9 normally still designates a missing phenotype; use --missing-phenotype if this is problematic.
 
 
  // here, we save the phenotypic value of the first trait by default
  // the complete set of phenotypic values for all traits is saved in the .fam file
 
  for (unsigned int j = 0; j < patch->size(FEM, Ax); ++j) {
 
    ind = patch->get(FEM, Ax, j);
 
    FH<<"fam"<<ind->getHome()+1<<" "<<ind->getID()+1<<" 0 0" //<<ind->getFatherID()<<" "<<ind->getMotherID() //not in file for adults
        <<" 2 "<<((double*)ind->getTraitValue(_FHLinkedTraitIndex))[0]<<" ";
 
    seq = (double**)ind->getTrait(_FHLinkedTraitIndex)->get_sequence();
 
    for(unsigned int k = 0; k < nb_trait; ++k) {
      for (unsigned int l = 0; l < nb_locus; ++l) {
 
        loc = l * nb_trait + k;
 
        FH<< BASE[ (seq[FEM][loc] == ref_all[l][1]) ]<<" "<< BASE[ (seq[MAL][loc] == ref_all[l][1]) ]<<" ";
      }
    }
 
    FH <<std::endl;
 
  }
 
  for (unsigned int j = 0; j < patch->size(MAL, Ax); ++j) {
 
    ind = patch->get(MAL, Ax, j);
 
    FH<<"fam"<<ind->getHome()+1<<" "<<ind->getID()+1<<" 0 0"  //<<ind->getFatherID()<<" "<<ind->getMotherID()
        <<" 1 "<<((double*)ind->getTraitValue(_FHLinkedTraitIndex))[0]<<" ";
 
    seq = (double**)ind->getTrait(_FHLinkedTraitIndex)->get_sequence();
 
    for(unsigned int k = 0; k < nb_trait; ++k) {
      for (unsigned int l = 0; l < nb_locus; ++l) {
 
        loc = l * nb_trait + k;
 
        FH<< BASE[ (seq[FEM][loc] == ref_all[l][1]) ]<<" "<< BASE[ (seq[MAL][loc] == ref_all[l][1]) ]<<" ";
      }
    }
    FH<<std::endl;
 
  }
}

References TraitFileHandler< TProtoQuanti >::_FHLinkedTrait, TraitFileHandler< TProtoQuanti >::_FHLinkedTraitIndex, FEM, Patch::get(), TProtoQuanti::get_allele_values(), TProtoQuanti::get_nb_locus(), TProtoQuanti::get_nb_traits(), TTrait::get_sequence(), Individual::getHome(), Individual::getID(), Individual::getTrait(), Individual::getTraitValue(), MAL, and Patch::size().

Referenced by write_PLINK().

Here is the caller graph for this function:

setOutputOption()

void TTQuantiFH::setOutputOption

(

string

opt

)

{
  _output_option = opt;
 
  if( !(_output_option == "genotypes" ||
      _output_option == "genotype" ||
      _output_option == "snp_genotypes" ||
      _output_option == "snp_genotype" || _output_option == "1" || _output_option == "0"))
    fatal("option \"%s\" to parameter \"quanti_output\" is not recognized\n", opt.c_str());
}

References _output_option, and fatal().

Referenced by TProtoQuanti::loadFileServices().

Here is the caller graph for this function:

write_PLINK()

void TTQuantiFH::write_PLINK

(

)

{
  if(_FHLinkedTrait->get_allele_model() > 2)
    fatal("PLINK file output for the quanti trait is only possible for di-allelic loci (for now)\n");
 
  unsigned int patchNbr = _pop->getPatchNbr();
  Patch* current_patch;
 
  std::string filename = get_path() + this->get_service()->getGenerationReplicateFileName() + ".ped";
 
  // the PED file -------------------------------------------------------------------------
  ofstream PED (filename.c_str(), ios::out);
 
  if(!PED) fatal("could not open plink .ped output file!!\n");
 
#ifdef _DEBUG_
  message("TTQuantiFH::write_PLINK (%s)\n",filename.c_str());
#endif
 
  for (unsigned int i = 0; i < patchNbr; ++i) {
 
    current_patch = _pop->getPatch(i);
 
    if( _pop->getCurrentAge() & ADULTS)
      print_PLINK_PED(PED, ADLTx, current_patch);
 
    if( _pop->getCurrentAge() & OFFSPRG)
      print_PLINK_PED(PED, OFFSx, current_patch);
 
  }
 
 
  PED.close();
  // the FAM file -------------------------------------------------------------------------
  // we write the .fam file only when more than one trait; phenotype columns are added here
 
  if(_FHLinkedTrait->get_nb_traits() > 1) {
 
    filename = get_path() + this->get_service()->getGenerationReplicateFileName() + ".fam";
 
    ofstream FAM (filename.c_str(), ios::out);
 
    if(!FAM) fatal("could not open plink .fam output file!!\n");
 
#ifdef _DEBUG_
    message("TTQuantiFH::write_PLINK (%s)\n",filename.c_str());
#endif
 
    for (unsigned int i = 0; i < patchNbr; ++i) {
 
      current_patch = _pop->getPatch(i);
 
      if( _pop->getCurrentAge() & ADULTS)
        print_PLINK_FAM(FAM, ADLTx, current_patch);
 
      if( _pop->getCurrentAge() & OFFSPRG)
        print_PLINK_FAM(FAM, OFFSx, current_patch);
    }
 
    FAM.close();
  }
 
  // the MAP file -------------------------------------------------------------------------
  filename = get_path() + this->get_service()->getGenerationReplicateFileName() + ".map";
 
  ofstream MAP (filename.c_str(), ios::out);
 
  if(!MAP) fatal("could not open plink .map output file!!\n");
 
#ifdef _DEBUG_
  message("TTQuantiFH::write_PLINK (%s)\n",filename.c_str());
#endif
 
  unsigned int nb_locus = _FHLinkedTrait->get_nb_locus();
 
  double **map;
  map = new double* [nb_locus];
  for(unsigned int k = 0; k < nb_locus; ++k) 
    map[k] = new double[2]; // [0]: chr; [1]: map pos in cM
 
  bool found = _FHLinkedTrait->_map.getGeneticMap(_FHLinkedTrait->get_type(), map, nb_locus);
 
  if( found ) {
 
    for(unsigned int k = 0; k < _FHLinkedTrait->get_nb_traits() ; ++k) {  //_FHLinkedTrait->get_nb_traits()
      for (unsigned int l = 0; l < nb_locus; ++l) { // this will give the locus ID
 
        // MAP FORMAT (PLINK1.9): chrmsm ID; Loc ID; position (cM); bp ID
        MAP<<map[l][0]+1<<" T"<<k+1<<"loc"<<l+1<<" "<<map[l][1]<<" "<<l+1<<endl;
 
      }
    }
  } else { // trait didn't register a genetic map, loci are unlinked (free recombination)
 
    warning("PLINK .map file: we assume QTL are unlinked, separated by 50M and on a single chromosome\n");
 
    // we're gonna set all loci on a single chrmsme, but 50M apart
    for(unsigned int k = 0; k < _FHLinkedTrait->get_nb_traits() ; ++k) {  //_FHLinkedTrait->get_nb_traits()
      for (unsigned int l = 0; l < nb_locus; ++l) { // this will give the locus ID
 
        // write: chromosome, locus name, map position, locus number
        MAP<< 1 <<" T"<<k+1<<"loc"<<l+1<<" "<< l*5000.0 + 1.0 <<" "<<l+1<<endl;
 
      }
    }
 
  }
  MAP.close();
 
  for(unsigned int k = 0; k < nb_locus; ++k) delete [] map[k];
  delete [] map;
}

References TraitFileHandler< TProtoQuanti >::_FHLinkedTrait, TTProtoWithMap::_map, FileHandler::_pop, ADLTx, ADULTS, fatal(), TProtoQuanti::get_allele_model(), TProtoQuanti::get_nb_locus(), TProtoQuanti::get_nb_traits(), FileHandler::get_path(), FileHandler::get_service(), TProtoQuanti::get_type(), Metapop::getCurrentAge(), FileServices::getGenerationReplicateFileName(), GeneticMap::getGeneticMap(), Metapop::getPatch(), Metapop::getPatchNbr(), message(), OFFSPRG, OFFSx, print_PLINK_FAM(), print_PLINK_PED(), and warning().

Referenced by FHwrite().

Here is the caller graph for this function:

write_TABLE()

void TTQuantiFH::write_TABLE

(

)

{
  std::string filename = get_filename();
 
  std::ofstream FILE (filename.c_str(), ios::out);
 
  if(!FILE) fatal("could not open \"%s\" output file!!\n",filename.c_str());
 
#ifdef _DEBUG_
  message("TTQuantiFH::write_TABLE (%s)\n", filename.c_str());
#endif
 
  //print the genotype at each locus
  unsigned int print_gene = 0;
 
  if(_output_option == "genotypes" || _output_option == "genotype")
    print_gene = 1;
  else if(_output_option == "snp_genotypes" || _output_option == "snp_genotype")
    print_gene = 2;
 
 
  //print the genotypic value of the individual in sus of the phenotype
  bool print_genotype = (!_FHLinkedTrait->get_env_var().empty());
 
  bool print_additive_genotype = (_FHLinkedTrait->get_dominance_model() != 0);
 
  // First column is the patch ID:
  FILE<<"pop ";
 
  // print the allelic values?
  if(print_gene) {
 
    // we print allelic values ordered trait:locus, alleles at loci affecting one trait are grouped
    for(unsigned int k = 0; k < _FHLinkedTrait->get_nb_traits(); k++)
      for(unsigned int l = 0; l < _FHLinkedTrait->get_nb_locus(); l++)
 
        if(print_gene == 1) {
          // "x" maternally inherited, "y" paternally inherited
          FILE<<"t"<<k+1<<"l"<<l+1<<"x "<<"t"<<k+1<<"l"<<l+1<<"y ";
 
        } else {
          // snp-encoded genotypes: 0, 1, 2
          FILE<<"t"<<k+1<<"l"<<l+1<<" ";
        }
  }
 
  for(unsigned int k = 0; k < _FHLinkedTrait->get_nb_traits(); k++) {
    FILE<<"P"<<k+1<< " "; 
    if(print_genotype) FILE<<"G"<<k+1<< " ";
    if(print_additive_genotype) FILE<<"A"<<k+1<< " ";
  }
 
  FILE<<"age sex home ped isMigrant father mother ID\n";
 
  // add the dominance effects if any:
  if(_FHLinkedTrait->get_dominance_model() != 0 && print_gene) {
    FILE<< "NA ";
    for(unsigned int k = 0; k < _FHLinkedTrait->get_nb_traits(); k++)
      for(unsigned int l = 0; l < _FHLinkedTrait->get_nb_locus(); l++)
 
        if(print_gene == 1) {
          FILE << _FHLinkedTrait->get_dominance(l, k) << " NA ";
        } else {
          FILE << _FHLinkedTrait->get_dominance(l, k) << " ";
        }
    FILE << "NA "; //phenotype
    if(print_genotype) FILE<<"NA ";
    if(print_additive_genotype) FILE<<"NA ";
    FILE << "NA NA NA NA NA NA NA NA\n";
  }
 
  age_t pop_age = _pop->getCurrentAge(); //flag telling which age class should contain individuals
 
  //we print anything that is present in the pop:
  if( (pop_age & OFFSPRG) != 0) {
    print(FILE, FEM, OFFSx, print_gene, print_genotype, print_additive_genotype);
    print(FILE, MAL, OFFSx, print_gene, print_genotype, print_additive_genotype);
  }
  if( (pop_age & ADULTS) != 0) {
    print(FILE, FEM, ADLTx, print_gene, print_genotype, print_additive_genotype);
    print(FILE, MAL, ADLTx, print_gene, print_genotype, print_additive_genotype);
  }
 
  FILE.close();  
 
  // reset to pop ptr to main pop:
  _pop = get_service()->get_pop_ptr();
}

References TraitFileHandler< TProtoQuanti >::_FHLinkedTrait, _output_option, FileHandler::_pop, ADLTx, ADULTS, fatal(), FEM, TProtoQuanti::get_dominance(), TProtoQuanti::get_dominance_model(), TProtoQuanti::get_env_var(), FileHandler::get_filename(), TProtoQuanti::get_nb_locus(), TProtoQuanti::get_nb_traits(), FileServices::get_pop_ptr(), FileHandler::get_service(), Metapop::getCurrentAge(), MAL, message(), OFFSPRG, OFFSx, and print().

Referenced by FHwrite().

Here is the caller graph for this function:

Member Data Documentation

_output_option

string TTQuantiFH::_output_option

private

Referenced by FHwrite(), setOutputOption(), and write_TABLE().

---

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

• ttquanti.h
• ttquanti.cc