Nemo  2.3.46
ttquanti.h
Go to the documentation of this file.
1 
30 #ifndef TTQUANTI_H
31 #define TTQUANTI_H
32 
33 #include <cmath>
34 #include "ttrait_with_map.h"
35 #include "filehandler.h"
36 #include "stathandler.h"
37 #include "metapop.h"
38 #include "datatable.h"
39 #ifdef HAS_GSL
40 #include <gsl/gsl_vector.h>
41 #include <gsl/gsl_matrix.h>
42 #include <gsl/gsl_eigen.h>
43 #endif
44 
45 class TTQuantiSH;
46 class TTQuantiFH;
47 class TTQFreqExtractor;
48 class TProtoQuanti;
49 
50 // ------------------------------------------------------------------------------
54 // ------------------------------------------------------------------------------
55 class TTQuanti: public TTraitWithMap {
56 
57 public:
58 
62  { }
63 
69  {
70  _phenotypes = new double [_nb_traits];
71  _init_value = new double [_nb_traits];
72  for (unsigned int i = 0; i < _nb_traits; ++i) _init_value[i] = T._init_value[i];
73  }
74 
75  ~TTQuanti () {reset();}
76 
77  //implements TTrait:
78  virtual void init ();
79  virtual void init_sequence ();
80  virtual void reset ();
81  virtual void inherit (TTrait* mother, TTrait* father);
82  virtual void mutate () {(this->*_mutationFuncPtr)();}
83  virtual void* set_trait (void* value) {return value;}
84  virtual void set_sequence (void** seq) {reset(); _sequence = (double**)seq;}
85  virtual void set_value ();
86  virtual void* getValue () const {return _phenotypes;}
87  virtual trait_t get_type () const {return QUANT;}
88  virtual void** get_sequence () const {return (void**)_sequence;}
89  virtual double get_allele_value (int loc, int all) {return (loc < (int)_seq_length && all < 2 ? _sequence[all][loc] : 0);}
90  virtual void set_allele_value (unsigned int locus, unsigned int allele, double value)
91  {assert(locus < _nb_locus && allele < 2); _sequence[allele][locus] = value;}
92  virtual void show_up ();
93  virtual TTQuanti* clone () {return new TTQuanti(*this);}
94  virtual TTQuanti& operator= (const TTrait& T);
95  virtual bool operator== (const TTrait& T);
96  virtual bool operator!= (const TTrait& T);
97 
98  double get_genotype (unsigned int trait);
99 
100  //implements StorableComponent:
101  virtual void store_data ( BinaryStorageBuffer* saver )
102  {saver->store(_sequence[0],_seq_length*sizeof(double));
103  saver->store(_sequence[1],_seq_length*sizeof(double));}
104 
105  virtual bool retrieve_data ( BinaryStorageBuffer* reader )
106  {reader->read(_sequence[0],_seq_length*sizeof(double));
107  reader->read(_sequence[1],_seq_length*sizeof(double));return true;}
108 
109  void mutate_noHC ();
110  void mutate_HC ();
111  //accessors:
112  void set_proto (TProtoQuanti* proto) {_myProto = proto;}
113  void set_nb_locus (unsigned int val) {_nb_locus = val;}
114  void set_nb_traits (unsigned int val) {_nb_traits = val;}
115  void set_seq_length (unsigned int val) {_seq_length = val;}
117  void set_init_value (double* val, unsigned int doInit);
118  void set_init_value (double* val);
119  void set_mutation_fptr (double* (TProtoQuanti::* val) (unsigned int) , bool _isHC)
121  void set_inherit_fptr (void (TProtoQuanti::* val) (sex_t, double*, double**)) {_inherit = val;}
122  void set_eVariance (double var) {_eVariance = var;}
123 
124  void set_allele (int locus, int allele, double value) {_sequence[allele][locus] = value;}
125 
126 private:
127 
128  double **_sequence;
129  double *_phenotypes;
130 
133  unsigned int _doInitMutation;
134  double _eVariance;
135 
137  double* (TProtoQuanti::* _getMutationValues) (unsigned int);
138  void (TProtoQuanti::* _inherit) (sex_t, double*, double**);
139  void (TTQuanti::* _mutationFuncPtr) (void);
140 };
141 
142 // ------------------------------------------------------------------------------
146 // ------------------------------------------------------------------------------
147 class TProtoQuanti : public TTProtoWithMap {
148 
149 public:
150 
151  TProtoQuanti ();
152  TProtoQuanti (const TProtoQuanti& T);
153  virtual ~TProtoQuanti ();
154 
155  unsigned int get_nb_traits() {return _nb_traits;}
156  unsigned int get_nb_locus() {return _nb_locus;}
157  unsigned int get_seq_length () {return _seq_length;}
158  double get_env_var () {return _eVariance;}
159  double get_trait_var (unsigned int trait) {return _mutation_matrix->get(trait, trait);}
161  unsigned int get_allele_model () {return _allele_model;}
162  double ** get_allele_values () const {return _allele_value;}
163 
165  bool setMutationParameters ();
166  bool setDiallelicMutationModel(); //@add from ttdouble
167  bool setContinuousMutationModel(); //@add from ttdouble
169 
170  double get_init_value (unsigned int i) { return _init_value[i]; }
171 
172  double* getMutationEffectMultivariateGaussian (unsigned int loc);
173  double* getMutationEffectBivariateGaussian (unsigned int loc);
174  double* getMutationEffectUnivariateGaussian (unsigned int loc);
175  double* getMutationEffectUnivariateDiallelic (unsigned int loc); //@add from ttdouble
176  double* getMutationEffectBivariateDiallelic (unsigned int loc); //@add from ttdouble
177 
178  void inherit_free (sex_t SEX, double* seq, double** parent);
179  void inherit_low (sex_t SEX, double* seq, double** parent);
180 
181  //implements TraitPrototype:
182  virtual void reset () {TTProtoWithMap::reset();}
183  virtual TTQuanti* hatch ();
184  virtual TProtoQuanti* clone () {return new TProtoQuanti(*this);}
185  virtual trait_t get_type () const {return QUANT;}
186 
187  //implementation of SimComponent:
188  virtual bool setParameters();
189  virtual void loadFileServices ( FileServices* loader );
190  virtual void loadStatServices ( StatServices* loader );
191 
192  //implementation of StorableComponent:
193  virtual void store_data ( BinaryStorageBuffer* saver )
194  {saver->store(&_seq_length,sizeof(int));}
195 
196  virtual bool retrieve_data ( BinaryStorageBuffer* reader )
197  {reader->read(&_seq_length,sizeof(int));return true;}
198 
199 private:
200 
202  unsigned int _allele_model;
203  double** _allele_value;
204 
205  //mutations:
208  gsl_vector *_eval, *_effects_multivar;
209  gsl_vector *_ws;
212  unsigned int _doInitMutation;
213 
214  //recombination:
217 
218  double _eVariance;
219 
220  double* (TProtoQuanti::* _mutation_func_ptr) (unsigned int);
221 
222  friend class TTQuanti;
223 
227 };
228 // ------------------------------------------------------------------------------
232 // ------------------------------------------------------------------------------
233 class TTQuantiSH : public TraitStatHandler<TProtoQuanti, TTQuantiSH> {
234 
235  double *_meanP, *_meanG, *_Va, *_Vb, *_Vp, *_covar,*_eigval,**_eigvect;
236  double **_pmeanP, **_pmeanG, **_pVa, **_pVp, **_pcovar, **_peigval, **_peigvect;
237 
238  unsigned int _nb_trait, _patchNbr;
239  bool _eVar;
240 
241  gsl_matrix *_G, *_evec;
242  gsl_vector *_eval;
243  gsl_eigen_symmv_workspace *_ws;
244 
247 
248 public:
249 
252  _meanP(0), _meanG(0), _Va(0), _Vb(0), _Vp(0), _covar(0), _eigval(0), _eigvect(0),
253  _pmeanP(0), _pmeanG(0), _pVa(0), _pVp(0), _pcovar(0), _peigval(0), _peigvect(0),
254  _nb_trait(0),_patchNbr(0),_G(0),_evec(0),_eval(0),_ws(0),
255  _table_set_gen(999999), _table_set_age(999999), _table_set_repl(999999)
256  {}
257 
258  virtual ~TTQuantiSH() {resetPtrs();}
259 
260  void resetPtrs();
261 
262  virtual void init ( );
263 
264  virtual bool setStatRecorders (std::string& token);
265  void addQuanti (age_t AGE);
266  void addEigen (age_t AGE);
267  void addEigenValues (age_t AGE);
268  void addEigenVect1 (age_t AGE);
269  void addQuantiPerPatch (age_t AGE);
270  void addAvgPerPatch (age_t AGE);
271  void addVarPerPatch (age_t AGE);
272  void addCovarPerPatch (age_t AGE);
273  void addEigenPerPatch (age_t AGE);
274  void addEigenValuesPerPatch (age_t AGE);
275  void addEigenVect1PerPatch (age_t AGE);
276  void addEigenStatsPerPatcg (age_t AGE);
277  void addSkewPerPatch(age_t AGE);
278 
279  void setDataTables (age_t AGE);
282  void setStats (age_t AGE);
283  double getMeanPhenot (unsigned int i) {return _meanP[i];}
284  double getVa (unsigned int i) {return _Va[i];}
285  double getVb (unsigned int i) {return _Vb[i];}
286  double getVp (unsigned int i) {return _Vp[i];}
287  double getQst (unsigned int i) {return _Vb[i]/(_Vb[i]+2*_Va[i]);}
288  double getCovar (unsigned int i) {return _covar[i];}
289  double getEigenValue (unsigned int i) {return _eigval[i];}
290  double getEigenVectorElt (unsigned int t1, unsigned int t2) {return _eigvect[t2][t1];}//eigenvectors arranged column-wise
291 
292  double getMeanPhenotPerPatch (unsigned int i, unsigned int p) {return _pmeanP[i][p];}
293  double getVaPerPatch (unsigned int i, unsigned int p) {return _pVa[i][p];}
294  double getVpPerPatch (unsigned int i, unsigned int p) {return _pVp[i][p];}
295  double getEigenValuePerPatch (unsigned int i, unsigned int p) {return _peigval[i][p];}
296  double getCovarPerPatch (unsigned int p, unsigned int i) {return _pcovar[p][i];}
297  double getEigenVectorEltPerPatch (unsigned int p, unsigned int v) {return _peigvect[p][v];}
298  double getSkewPerPatch (unsigned int i, unsigned int p);
299 };
300 // ------------------------------------------------------------------------------
304 // ------------------------------------------------------------------------------
305 class TTQuantiFH : public TraitFileHandler<TProtoQuanti> {
306 
308 
309 public:
311  virtual ~TTQuantiFH(){}
312 
313  void setOutputOption (string opt) {_output_option = opt;}
314 
315  virtual void FHwrite ();
316  void print(ofstream& FH, age_idx Ax, bool print_gene, bool print_genotype);
317  virtual void FHread (string& filename) {}
318 };
319 // ------------------------------------------------------------------------------
323 // ------------------------------------------------------------------------------
324 class TTQFreqExtractor : public TraitFileHandler<TProtoQuanti> {
325 
326  vector< string > _records;
327 
328 public:
330 
331  void resetTable();
332 
333  virtual ~TTQFreqExtractor () {}
334  virtual void FHwrite ();
335  virtual void FHread (string& filename) {}
336 };
337 
338 #endif
339 
void addEigenValuesPerPatch(age_t AGE)
Definition: ttquanti.cc:1666
void addAvgPerPatch(age_t AGE)
Definition: ttquanti.cc:1500
TTQuanti.
Definition: ttquanti.h:55
Template class for the trait's FileHandler.
Definition: filehandler.h:217
double _effects_bivar[2]
Definition: ttquanti.h:211
void set_seq_length(unsigned int val)
Definition: ttquanti.h:115
double getEigenVectorElt(unsigned int t1, unsigned int t2)
Definition: ttquanti.h:290
TTQuantiFH.
Definition: ttquanti.h:305
The Service class used to manage the StatHandler objects.
Definition: statservices.h:50
virtual double get_allele_value(int loc, int all)
Definition: ttquanti.h:89
TMatrix * _mutation_matrix
Definition: ttquanti.h:206
double get_genotype(unsigned int trait)
Definition: ttquanti.cc:1105
void inherit_free(sex_t SEX, double *seq, double **parent)
Definition: ttquanti.cc:603
double ** _pVa
Definition: ttquanti.h:236
#define ADULTS
Adults age class flag (breeders).
Definition: types.h:54
std::string trait_t
Trait types.
Definition: types.h:63
void resetTable()
Definition: ttquanti.cc:2122
void(TProtoQuanti::* _inherit)(sex_t, double *, double **)
Definition: ttquanti.h:138
double _genomic_mutation_rate
Definition: ttquanti.h:210
virtual void FHwrite()
Definition: ttquanti.cc:1993
double *(TProtoQuanti::* _getMutationValues)(unsigned int)
Definition: ttquanti.h:137
unsigned int _seq_length
Definition: ttquanti.h:131
void set_mutation_matrix_decomposition()
Definition: ttquanti.cc:488
void addEigenVect1PerPatch(age_t AGE)
Definition: ttquanti.cc:1700
virtual void set_sequence(void **seq)
Definition: ttquanti.h:84
TTQFreqExtractor(TProtoQuanti *T)
Definition: ttquanti.h:329
TTQuanti(const TTQuanti &T)
Definition: ttquanti.h:64
virtual void ** get_sequence() const
Definition: ttquanti.h:88
double ** _pcovar
Definition: ttquanti.h:236
virtual void mutate()
Definition: ttquanti.h:82
double getCovar(unsigned int i)
Definition: ttquanti.h:288
double _genomic_mutation_rate
Definition: ttquanti.h:132
age_idx
Array index of the age classes in the patch sizes and containers arrays.
Definition: types.h:41
virtual bool setParameters()
Definition: ttquanti.cc:154
virtual void * getValue() const
Definition: ttquanti.h:86
double * _Va
Definition: ttquanti.h:235
void set_eVariance(double var)
Definition: ttquanti.h:122
double * _phenotypes
Definition: ttquanti.h:129
gsl_matrix * _gsl_mutation_matrix
Definition: ttquanti.h:207
void resetPtrs()
Definition: ttquanti.cc:1156
unsigned int get_seq_length()
Definition: ttquanti.h:157
unsigned int _seq_length
Definition: ttquanti.h:201
TTQuantiSH.
Definition: ttquanti.h:233
unsigned int get_nb_traits()
Definition: ttquanti.h:155
virtual TTQuanti & operator=(const TTrait &T)
Definition: ttquanti.cc:794
virtual void * set_trait(void *value)
Definition: ttquanti.h:83
virtual void show_up()
Definition: ttquanti.cc:1119
DataTable< double > _genoTable
Definition: ttquanti.h:245
void set_nb_locus(unsigned int val)
Definition: ttquanti.h:113
void setDataTables(age_t AGE)
Definition: ttquanti.cc:1790
unsigned int age_t
Age class flags.
Definition: types.h:46
TTProtoWithMap.
Definition: ttrait_with_map.h:165
virtual void init()
Definition: ttquanti.cc:1212
void set_genomic_mutation_rate(double val)
Definition: ttquanti.h:116
double getQst(unsigned int i)
Definition: ttquanti.h:287
virtual void reset()
Definition: ttquanti.cc:1012
virtual ~TTQuantiSH()
Definition: ttquanti.h:258
double getVb(unsigned int i)
Definition: ttquanti.h:285
bool setDiallelicMutationModel()
Definition: ttquanti.cc:325
unsigned int _allele_model
Definition: ttquanti.h:202
void addSkewPerPatch(age_t AGE)
Definition: ttquanti.cc:1739
TTQuanti()
Definition: ttquanti.h:59
double * _Vp
Definition: ttquanti.h:235
void addQuantiPerPatch(age_t AGE)
Definition: ttquanti.cc:1482
void mutate_HC()
Definition: ttquanti.cc:1062
void set_proto(TProtoQuanti *proto)
Definition: ttquanti.h:112
void set_mutation_fptr(double *(TProtoQuanti::*val)(unsigned int), bool _isHC)
Definition: ttquanti.h:119
gsl_matrix * _G
Definition: ttquanti.h:241
virtual void set_allele_value(unsigned int locus, unsigned int allele, double value)
Definition: ttquanti.h:90
bool setContinuousMutationModel()
Definition: ttquanti.cc:389
virtual bool retrieve_data(BinaryStorageBuffer *reader)
Definition: ttquanti.h:105
double get(unsigned int i, unsigned int j)
Accessor to element at row i and column j.
Definition: tmatrix.h:120
Interface for all trait types, declares all basic trait operations.
Definition: ttrait.h:46
gsl_matrix * _evec
Definition: ttquanti.h:241
virtual bool retrieve_data(BinaryStorageBuffer *reader)
Definition: ttquanti.h:196
void addVarPerPatch(age_t AGE)
Definition: ttquanti.cc:1535
double * _meanP
Definition: ttquanti.h:235
virtual void FHread(string &filename)
Definition: ttquanti.h:317
double _mutation_correlation
Definition: ttquanti.h:210
void set_nb_traits(unsigned int val)
Definition: ttquanti.h:114
TProtoQuanti * _myProto
Definition: ttquanti.h:136
double getEigenVectorEltPerPatch(unsigned int p, unsigned int v)
Definition: ttquanti.h:297
TProtoQuanti.
Definition: ttquanti.h:147
double * getMutationEffectUnivariateGaussian(unsigned int loc)
Definition: ttquanti.cc:530
double * getMutationEffectBivariateDiallelic(unsigned int loc)
Definition: ttquanti.cc:546
void print(ofstream &FH, age_idx Ax, bool print_gene, bool print_genotype)
Definition: ttquanti.cc:2034
virtual void reset()
Definition: ttquanti.h:182
TProtoQuanti()
Definition: ttquanti.cc:59
void store(void *stream, unsigned int nb_bytes)
Definition: binarystoragebuffer.h:130
double ** _sequence
Definition: ttquanti.h:128
virtual ~TTQFreqExtractor()
Definition: ttquanti.h:333
~TTQuanti()
Definition: ttquanti.h:75
double get_env_var()
Definition: ttquanti.h:158
virtual void reset()
Definition: ttrait_with_map.cc:505
double getVa(unsigned int i)
Definition: ttquanti.h:284
size_t _locusByteSize
Definition: ttquanti.h:216
TTQuantiSH * _stats
Definition: ttquanti.h:224
bool * _all_chooser
Definition: ttquanti.h:215
double getVaPerPatch(unsigned int i, unsigned int p)
Definition: ttquanti.h:293
unsigned int _nb_traits
Definition: ttquanti.h:131
void set_init_value(double *val, unsigned int doInit)
Definition: ttquanti.cc:840
virtual void init_sequence()
Definition: ttquanti.cc:879
double * getMutationEffectMultivariateGaussian(unsigned int loc)
Definition: ttquanti.cc:513
Template class for the trait's StatHandler.
Definition: stathandler.h:168
double get_init_value(unsigned int i)
Definition: ttquanti.h:170
double ** _pVp
Definition: ttquanti.h:236
TTQuantiSH(TProtoQuanti *TP)
Definition: ttquanti.h:250
virtual void loadStatServices(StatServices *loader)
Definition: ttquanti.cc:777
sex_t
Sex types, males are always 0 and females 1!!
Definition: types.h:36
gsl_matrix * _evect
Definition: ttquanti.h:207
virtual bool setStatRecorders(std::string &token)
Definition: ttquanti.cc:1277
void setOffsprgStats()
Definition: ttquanti.h:281
double ** _pmeanG
Definition: ttquanti.h:236
double getMeanPhenot(unsigned int i)
Definition: ttquanti.h:283
virtual void inherit(TTrait *mother, TTrait *father)
Definition: ttquanti.cc:1029
virtual ~TProtoQuanti()
Definition: ttquanti.cc:141
double ** _eigvect
Definition: ttquanti.h:235
double * _meanG
Definition: ttquanti.h:235
TTraitWithMap.
Definition: ttrait_with_map.h:235
double * _eigval
Definition: ttquanti.h:235
gsl_vector * _eval
Definition: ttquanti.h:208
virtual trait_t get_type() const
Definition: ttquanti.h:185
double * getMutationEffectBivariateGaussian(unsigned int loc)
Definition: ttquanti.cc:521
void set_inherit_fptr(void(TProtoQuanti::*val)(sex_t, double *, double **))
Definition: ttquanti.h:121
double ** get_allele_values() const
Definition: ttquanti.h:162
double * _init_value
Definition: ttquanti.h:132
A class to store any kind of data in a char buffer before unloading it in a binary data file...
Definition: binarystoragebuffer.h:40
unsigned int _doInitMutation
Definition: ttquanti.h:133
void addEigenPerPatch(age_t AGE)
Definition: ttquanti.cc:1622
double getVpPerPatch(unsigned int i, unsigned int p)
Definition: ttquanti.h:294
void set_allele(int locus, int allele, double value)
Definition: ttquanti.h:124
size_t _sizeofLocusType
Definition: ttquanti.h:216
TTQuantiFH(TProtoQuanti *T)
Definition: ttquanti.h:310
unsigned int _nb_trait
Definition: ttquanti.h:238
unsigned int get_allele_model()
Definition: ttquanti.h:161
A class to handle matrix in params, coerces matrix into a vector of same total size.
Definition: tmatrix.h:49
virtual void set_value()
Definition: ttquanti.cc:1082
unsigned int _nb_traits
Definition: ttquanti.h:201
void addEigen(age_t AGE)
Definition: ttquanti.cc:1396
unsigned int _nb_locus
Definition: ttquanti.h:201
unsigned int _patchNbr
Definition: ttquanti.h:238
string _output_option
Definition: ttquanti.h:307
void read(void *out, unsigned int nb_bytes)
Definition: binarystoragebuffer.h:147
virtual trait_t get_type() const
Definition: ttquanti.h:87
TTQFreqExtractor * _freqExtractor
Definition: ttquanti.h:226
gsl_vector * _effects_multivar
Definition: ttquanti.h:208
double _eVariance
Definition: ttquanti.h:134
void addQuanti(age_t AGE)
Definition: ttquanti.cc:1336
double *(TProtoQuanti::* _mutation_func_ptr)(unsigned int)
Definition: ttquanti.h:220
double get_mutation_correlation()
Definition: ttquanti.h:160
void reset_mutation_pointers()
Definition: ttquanti.cc:557
unsigned int _doInitMutation
Definition: ttquanti.h:212
void setAdultStats()
Definition: ttquanti.h:280
virtual TProtoQuanti * clone()
Definition: ttquanti.h:184
unsigned int _table_set_repl
Definition: ttquanti.h:246
virtual bool operator==(const TTrait &T)
Definition: ttquanti.cc:815
#define OFFSPRG
Offspring age class flag.
Definition: types.h:50
void addEigenVect1(age_t AGE)
Definition: ttquanti.cc:1455
DataTable< double > _phenoTable
Definition: ttquanti.h:245
double * _init_value
Definition: ttquanti.h:211
virtual void loadFileServices(FileServices *loader)
Definition: ttquanti.cc:715
double getCovarPerPatch(unsigned int p, unsigned int i)
Definition: ttquanti.h:296
double * _mutation_sigma
Definition: ttquanti.h:211
unsigned int _table_set_age
Definition: ttquanti.h:246
vector< string > _records
Definition: ttquanti.h:326
virtual void store_data(BinaryStorageBuffer *saver)
Definition: ttquanti.h:101
void addEigenStatsPerPatcg(age_t AGE)
void setOutputOption(string opt)
Definition: ttquanti.h:313
double * getMutationEffectUnivariateDiallelic(unsigned int loc)
Definition: ttquanti.cc:538
virtual void init()
Definition: ttquanti.cc:863
gsl_vector * _ws
Definition: ttquanti.h:209
double getEigenValue(unsigned int i)
Definition: ttquanti.h:289
double * _Vb
Definition: ttquanti.h:235
virtual void FHwrite()
Definition: ttquanti.cc:2129
void(TTQuanti::* _mutationFuncPtr)(void)
Definition: ttquanti.h:139
bool setMutationParameters()
Definition: ttquanti.cc:269
void addEigenValues(age_t AGE)
Definition: ttquanti.cc:1427
double getVp(unsigned int i)
Definition: ttquanti.h:286
virtual TTQuanti * clone()
Definition: ttquanti.h:93
virtual void store_data(BinaryStorageBuffer *saver)
Definition: ttquanti.h:193
A class to manage the files associated with each components of the simulation.
Definition: fileservices.h:51
gsl_vector * _eval
Definition: ttquanti.h:242
#define QUANT
Definition: types.h:72
unsigned int _table_set_gen
Definition: ttquanti.h:246
unsigned int get_nb_locus()
Definition: ttquanti.h:156
double ** _peigval
Definition: ttquanti.h:236
double ** _pmeanP
Definition: ttquanti.h:236
bool _eVar
Definition: ttquanti.h:239
double ** _peigvect
Definition: ttquanti.h:236
virtual bool operator!=(const TTrait &T)
Definition: ttquanti.cc:830
double ** _allele_value
Definition: ttquanti.h:203
double getSkewPerPatch(unsigned int i, unsigned int p)
Definition: ttquanti.cc:1775
TTQuantiFH * _writer
Definition: ttquanti.h:225
void addCovarPerPatch(age_t AGE)
Definition: ttquanti.cc:1581
unsigned int _nb_locus
Definition: ttquanti.h:131
double getEigenValuePerPatch(unsigned int i, unsigned int p)
Definition: ttquanti.h:295
gsl_eigen_symmv_workspace * _ws
Definition: ttquanti.h:243
void inherit_low(sex_t SEX, double *seq, double **parent)
Definition: ttquanti.cc:621
virtual void FHread(string &filename)
Definition: ttquanti.h:335
double getMeanPhenotPerPatch(unsigned int i, unsigned int p)
Definition: ttquanti.h:292
void setStats(age_t AGE)
Definition: ttquanti.cc:1860
void mutate_noHC()
Definition: ttquanti.cc:1041
TTQFreqExtractor.
Definition: ttquanti.h:324
virtual TTQuanti * hatch()
Definition: ttquanti.cc:692
virtual ~TTQuantiFH()
Definition: ttquanti.h:311
double get_trait_var(unsigned int trait)
Definition: ttquanti.h:159
double _eVariance
Definition: ttquanti.h:218
double * _covar
Definition: ttquanti.h:235

Generated for Nemo v2.3.0 by  doxygen 1.8.8 --
Catalogued on GSR