Skip to content
Snippets Groups Projects
Commit ae9e4069 authored by Veronique Legrand's avatar Veronique Legrand
Browse files

limit to 2 digits after dot the number of digits displayed for flase positive probability

parent cc11d121
No related branches found
No related tags found
No related merge requests found
...@@ -83,7 +83,7 @@ const void printVerboseInfo(const ROCKparams& Rparms,const Filter& filter,const ...@@ -83,7 +83,7 @@ const void printVerboseInfo(const ROCKparams& Rparms,const Filter& filter,const
cout<<setw(40)<<left<<"expected no. distinct k-mers"<<Rparms.get_expected_nbKmers()<<endl; cout<<setw(40)<<left<<"expected no. distinct k-mers"<<Rparms.get_expected_nbKmers()<<endl;
cout<<setw(40)<<left<<"no. hash. function(s)"<<parms.lambda<<endl; cout<<setw(40)<<left<<"no. hash. function(s)"<<parms.lambda<<endl;
cout<<setw(40)<<left<<"expected false positive proba."<<Rparms.get_expected_collision_proba()<<" (cov. depth > 1)"<<endl; cout<<setw(40)<<left<<"expected false positive proba."<<setprecision (2)<<fixed<<Rparms.get_expected_collision_proba()<<" (cov. depth > 1)"<<endl;
for (int i=0;i<parms.lambda;i++) { for (int i=0;i<parms.lambda;i++) {
string tmp="no. buckets for hash. "; string tmp="no. buckets for hash. ";
...@@ -105,13 +105,13 @@ const void printVerboseInfo(const ROCKparams& Rparms,const Filter& filter,const ...@@ -105,13 +105,13 @@ const void printVerboseInfo(const ROCKparams& Rparms,const Filter& filter,const
int smallest_kappa=parms.kappa; int smallest_kappa=parms.kappa;
if (parms.kappa_prime>0) smallest_kappa=parms.kappa_prime; if (parms.kappa_prime>0) smallest_kappa=parms.kappa_prime;
float p =getCollisionProba(smallest_kappa,approx_nbDiffKm,Pi_js[0],parms.lambda); float p =getCollisionProba(smallest_kappa,approx_nbDiffKm,Pi_js[0],parms.lambda);
cout<<setw(40)<<left<<"approx. false positive proba."<<p<<endl; cout<<setw(40)<<left<<"approx. false positive proba."<<setprecision (2)<<p<<endl;
cout<<setw(40)<<left<<"no. selected reads"<<rc.nb_output_reads<<endl; cout<<setw(40)<<left<<"no. selected reads"<<rc.nb_output_reads<<endl;
for (it_pe=v_PE_files.begin();it_pe!=v_PE_files.end();it_pe++) { for (it_pe=v_PE_files.begin();it_pe!=v_PE_files.end();it_pe++) {
cout<<setw(40)<<left<<"PE output files"<<it_pe->PE1.out_fq_file<<" "<<it_pe->PE2.out_fq_file<<" "<<it_pe->PE1.undef_fq_file<<" "<<it_pe->PE2.undef_fq_file<<endl; cout<<setw(40)<<left<<"PE output files"<<it_pe->PE1.out_fq_file<<" "<<it_pe->PE2.out_fq_file<<" "<<it_pe->PE1.undef_fq_file<<" "<<it_pe->PE2.undef_fq_file<<endl;
} }
for (it_s=single_files.begin();it_s!=single_files.end();it_s++) { for (it_s=single_files.begin();it_s!=single_files.end();it_s++) {
cout<<setw(40)<<left<<"SE output file(s)"<<it_s->out_fq_file<<" "<<it_s->undef_fq_file<<endl; cout<<setw(40)<<left<<"SE output files"<<it_s->out_fq_file<<" "<<it_s->undef_fq_file<<endl;
} }
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment