TTQuanti_continuous_full_pleio Class Reference

TTQuanti_continuous_full_pleio : universal pleiotropy. More...

#include <ttquanti.h>

Inheritance diagram for TTQuanti_continuous_full_pleio:

Collaboration diagram for TTQuanti_continuous_full_pleio:

Public Member Functions
	TTQuanti_continuous_full_pleio ()
	TTQuanti_continuous_full_pleio (const TTQuanti_continuous_full_pleio &TT)
virtual	~TTQuanti_continuous_full_pleio ()
virtual void	init_sequence ()
virtual TTQuanti_continuous_full_pleio *	clone ()
virtual void	show_up ()
virtual double	get_additive_genotype (const unsigned int trait) const
virtual double	get_dominant_genotype (const unsigned int trait) const
virtual void	copy_sequence_block (sex_t SEX, unsigned int strand, unsigned int from_locus, unsigned int to_locus, const TTQuanti *parent)
virtual void	copy_sequence_1locus (sex_t SEX, unsigned int strand, unsigned int at, const TTQuanti *parent)
Public Member Functions inherited from TTQuanti_continuous
	TTQuanti_continuous ()
	TTQuanti_continuous (const TTQuanti &T)
virtual	~TTQuanti_continuous ()
virtual void	reset ()
virtual void	init ()
virtual void	set_sequence (void **seq)
virtual void **	get_sequence () const
virtual unsigned int	get_allele (int loc, int all) const
virtual double	get_allele_value (int loc, int all) const
virtual void	set_allele_value (unsigned int locus, unsigned int allele, double value)
virtual TTQuanti_continuous &	operator= (const TTrait &T)
virtual bool	operator== (const TTrait &T)
virtual bool	operator!= (const TTrait &T)
virtual void	store_data (BinaryStorageBuffer *saver)
virtual bool	retrieve_data (BinaryStorageBuffer *reader)
virtual double	get_full_genotype (unsigned int trait)
virtual void	set_allele (int locus, int allele, double value)
virtual void	mutate_add (unsigned int position, unsigned int allele, double value)
virtual void	mutate_inplace (unsigned int position, unsigned int allele, double value)
virtual bool	get_allele_bit (unsigned int position, unsigned int allele) const
virtual void	set_allele_bit (unsigned int position, unsigned int allele, bool value)
Public Member Functions inherited from TTQuanti
	TTQuanti ()
	TTQuanti (const TTQuanti &T)
virtual	~TTQuanti ()
virtual trait_t	get_type () const
virtual void	mutate ()
virtual void	inherit (const TTrait *mother, const TTrait *father)
virtual void *	set_trait (void *value)
virtual void	set_value ()
virtual void *	getValue () const
void	set_proto (TProtoQuanti *proto)
TProtoQuanti *	get_proto ()
double	get_phenotype (unsigned int trai)
void	set_phenotype (unsigned int trait, double value)
void	reset_phenotype_to_genotypic_value ()
Public Member Functions inherited from TTrait
virtual	~TTrait ()
Public Member Functions inherited from StorableComponent
virtual	~StorableComponent ()

Additional Inherited Members
Protected Attributes inherited from TTQuanti_continuous
double **	_sequence
Protected Attributes inherited from TTQuanti
double *	_phenotypes
double *	_genotypic_values
TProtoQuanti *	_myProto

Detailed Description

TTQuanti_continuous_full_pleio : universal pleiotropy.

Constructor & Destructor Documentation

TTQuanti_continuous_full_pleio() [1/2]

TTQuanti_continuous_full_pleio::TTQuanti_continuous_full_pleio ( )

inline

: TTQuanti_continuous() {}

Referenced by clone().

TTQuanti_continuous_full_pleio() [2/2]

TTQuanti_continuous_full_pleio::TTQuanti_continuous_full_pleio ( const TTQuanti_continuous_full_pleio &  TT )

inline

: TTQuanti_continuous(TT) {}

~TTQuanti_continuous_full_pleio()

virtual TTQuanti_continuous_full_pleio::~TTQuanti_continuous_full_pleio ( )

inlinevirtual

{}

Member Function Documentation

clone()

virtual TTQuanti_continuous_full_pleio* TTQuanti_continuous_full_pleio::clone ( )

inlinevirtual

Implements TTrait.

{return new TTQuanti_continuous_full_pleio(*this);}

References TTQuanti_continuous_full_pleio().

copy_sequence_1locus()

void TTQuanti_continuous_full_pleio::copy_sequence_1locus ( sex_t  SEX, unsigned int  strand, unsigned int  at, const TTQuanti *  parent  )

inlinevirtual

Implements TTQuanti.

{
  const double *orig = (const double*)parent->get_sequence()[strand];
  double *seq = _sequence[SEX];
  unsigned int from_pos = at * _myProto->get_num_traits();
 
  memcpy(&seq[from_pos], &orig[from_pos], _myProto->get_locus_byte_size());
}

References TTQuanti::_myProto, TTQuanti_continuous::_sequence, TProtoQuanti::get_locus_byte_size(), TProtoQuanti::get_num_traits(), and TTrait::get_sequence().

copy_sequence_block()

void TTQuanti_continuous_full_pleio::copy_sequence_block ( sex_t  SEX, unsigned int  strand, unsigned int  from_locus, unsigned int  to_locus, const TTQuanti *  parent  )

inlinevirtual

Implements TTQuanti.

{
  const double *orig = (const double*)parent->get_sequence()[strand];
  double *seq = _sequence[SEX];
  unsigned int from_pos = from_locus * _myProto->get_num_traits();
 
  assert(to_locus >= from_locus);
 
  size_t block_size = (to_locus - from_locus) * _myProto->get_locus_byte_size();
 
  memcpy(&seq[from_pos], &orig[from_pos], block_size);
 
}

References TTQuanti::_myProto, TTQuanti_continuous::_sequence, TProtoQuanti::get_locus_byte_size(), TProtoQuanti::get_num_traits(), and TTrait::get_sequence().

get_additive_genotype()

double TTQuanti_continuous_full_pleio::get_additive_genotype ( const unsigned int  trait ) const

inlinevirtual

Implements TTQuanti.

{
  unsigned int loc = trait;  // initial position in the sequence, in the first locus
  double genotype = 0;
 
 
  for(unsigned int j = 0; j < _myProto->get_num_locus(); ++j) {
    // alleles are ordered locus-trait: seq[0]:l0t0,l0t1,l0t2, ..., l1t0,l1t1,l1t2, ...
    genotype +=  _sequence[0][loc] + _sequence[1][loc];
    loc += _myProto->get_num_traits(); // stride size on the genome, between loci
  }
 
  return genotype;
}

References TTQuanti::_myProto, TTQuanti_continuous::_sequence, TProtoQuanti::get_num_locus(), and TProtoQuanti::get_num_traits().

get_dominant_genotype()

double TTQuanti_continuous_full_pleio::get_dominant_genotype ( const unsigned int  trait ) const

virtual

Implements TTQuanti.

{
  unsigned int loc = trait;
  double genotype = 0;
 
  for(unsigned int j = 0; j < _myProto->get_num_locus(); ++j) {
    genotype += _myProto->get_genotype_with_dominance( _sequence[0][loc] , _sequence[1][loc], j, trait);
    loc += _myProto->get_num_traits(); // stride size on the genome, between loci
  }
 
  return genotype;
}

References TTQuanti::_myProto, TTQuanti_continuous::_sequence, TProtoQuanti::get_genotype_with_dominance(), TProtoQuanti::get_num_locus(), and TProtoQuanti::get_num_traits().

init_sequence()

void TTQuanti_continuous_full_pleio::init_sequence ( )

inlinevirtual

Implements TTrait.

{
//  cout<<" TTQuanti_continuous_full_pleio::init_sequence\n";
 
  //options:
  //0: no variation, init value = (trait value)/(2*_num_locus)
  //1: init value = (trait value)/(2*_num_locus) + 1 mutation/locus
  //2: init value = (trait value)/(2*_num_locus) + random deviate ~N(0, _init_variance/2*_num_locus)
 
  //Note: the initial values may have been set individually by LCE_quanti
  //      it wouldn't make sense then to store the init values in the prototype
  //      because init values are set patch-specific by LCE_quanti
  //      Problem: without LCE_quanti, the original init values must not change from one replicate
  //      to the other, so we use a local variable
 
  double my_init[_myProto->get_num_traits()];
 
  for(unsigned int j = 0; j < _myProto->get_num_traits(); j++)
    my_init[j] = _myProto->get_init_value(j) / (2*_myProto->get_num_locus());
 
  unsigned int pos = 0;
 
  for(unsigned int j = 0; j < _myProto->get_num_traits(); ++j) {
    pos = j;
    //set the allele values from the init value, mutations may be added later
    for(unsigned int i = 0; i < _myProto->get_num_locus(); ++i) {
      _sequence[0][pos] = my_init[j];
      _sequence[1][pos] = my_init[j];
      pos += _myProto->get_num_traits();
    }
  }
 
  //add random effects to allele values if needed
  if(_myProto->get_doInitMutation() == 1) {
 
    double *mut;
 
    for(unsigned int i = 0; i < _myProto->get_num_locus(); i++) {
 
      pos = i * _myProto->get_num_traits();
 
      // change only one of the two allele at the locus, randomly
      mut = (_myProto->*_myProto->_getMutationValues)(i);
 
      for(unsigned int j = 0; j < _myProto->get_num_traits(); j++) {
        _sequence[RAND::RandBool()][pos] += mut[j];
        pos++;
      }
    } // end for_num_locus
  } // end if_doInit == 1
 
  // random deviate to initial additive effect
  if(_myProto->get_doInitMutation() == 2) {
 
    double initSD[_myProto->get_num_traits()];
 
    for(unsigned int j = 0; j < _myProto->get_num_traits(); j++) {
      initSD[j] = sqrt(_myProto->get_init_variance(j) / (2*_myProto->get_num_locus()) );
    }
 
    for(unsigned int i = 0; i < _myProto->get_num_locus(); i++) {
 
      pos = i * _myProto->get_num_traits();
 
      for(unsigned int j = 0; j < _myProto->get_num_traits(); j++) {
        _sequence[0][pos] += RAND::Gaussian(initSD[j]);
        _sequence[1][pos] += RAND::Gaussian(initSD[j]);
        pos++;
      }
    } // end for_num_locus
  } // end if_doInit == 2
 
//  cout<<" TTQuanti_continuous_full_pleio::init_sequence::end\n";
}

References TProtoQuanti::_getMutationValues, TTQuanti::_myProto, TTQuanti_continuous::_sequence, RAND::Gaussian(), TProtoQuanti::get_doInitMutation(), TProtoQuanti::get_init_value(), TProtoQuanti::get_init_variance(), TProtoQuanti::get_num_locus(), TProtoQuanti::get_num_traits(), and RAND::RandBool().

show_up()

void TTQuanti_continuous_full_pleio::show_up ( )

virtual

Implements TTrait.

{
  message("\
          Trait's type: QUANTI (continuous, full pleiotropy)\n\
          traits: %i\n\
          loci: %i\n\
          seq length: %i\n",_myProto->get_num_traits(),_myProto->get_num_locus()
          ,_myProto->get_seq_length());
 
  for(unsigned int i = 0; i < _myProto->get_num_traits(); i++)
    message("phenotype %i: %f\n",i+1,_phenotypes[i]);
 
  message("_sequence: \n0:");
  unsigned int loc;
  for(unsigned int i = 0; i < _myProto->get_num_traits(); ++i) {
    message("\nt%i:",i+1);
    for(unsigned int j = 0; (j < _myProto->get_num_locus()); ++j) {
      loc = j * _myProto->get_num_traits() + i;
      message("%.3f,",_sequence[0][loc]);
    }
  }
  message("\n1:");
  for(unsigned int i = 0; i < _myProto->get_num_traits(); ++i) {
    message("\nt%i:",i+1);
    for(unsigned int j = 0; (j < _myProto->get_num_locus()); ++j) {
      loc = j * _myProto->get_num_traits() + i;
      message("%.3f,",_sequence[1][loc]);
    }
  }
  message("\n");
 
}

References TTQuanti::_myProto, TTQuanti::_phenotypes, TTQuanti_continuous::_sequence, TProtoQuanti::get_num_locus(), TProtoQuanti::get_num_traits(), TProtoQuanti::get_seq_length(), and message().

---

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

• ttquanti.h
• ttquanti.cc