Skip to content

Commit fd37f59

Browse files
committed
update header file for the error matrix computation functions
1 parent 9f1af38 commit fd37f59

File tree

1 file changed

+7
-0
lines changed

1 file changed

+7
-0
lines changed

model/modelgenotype.h

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -180,10 +180,17 @@ class ModelGenotype : virtual public ModelMarkov
180180

181181
// Error model parameters
182182
bool use_error_model; // whether the error model is enabled
183+
double **error_matrix; // [num_states][num_states] probability matrix
183184
double ado_rate; // allelic dropout rate
184185
double error_rate; // amplification/sequencing error rate
185186
bool fix_error_params; // whether to fix error parameters (not optimise)
186187

188+
void computeErrorMatrix();
189+
void free_error_matrix();
190+
double computeErrorProbabilities(int true_state, int obs_state, double ado, double err);
191+
bool is_heterozygote(int state);
192+
void getAlleles(int state, int &allele1, int &allele2);
193+
187194
/*
188195
* Bounds for the optimisation
189196
* lower_ado: lower bound for ADO - 0.0

0 commit comments

Comments
 (0)