Commit c0d8ada8 authored by Damien Leroux's avatar Damien Leroux
Browse files

WIP (attempt at fixing qtl_threshold_all_chrom...)

parent fb84c33b
......@@ -196,7 +196,7 @@ struct settings_t {
, epistasis_qtl_selection_filename()
, qtl_threshold(0)
, n_permutations(10000)
, qtl_threshold_quantile(0.025320565519104) /* FIXME : .05 by default! */
, qtl_threshold_quantile(0.05)
, cofactor_threshold(0)
, detection_method(detection_method_t::Undef)
, detection_window(10.)
......@@ -281,8 +281,8 @@ struct settings_t {
{
if (pool && msg_handler_t::color()) {
pool->set_title(t);
} else {
MSG_INFO("Current task: " << t);
/*} else {*/
/*MSG_INFO("Current task: " << t);*/
}
}
......
......@@ -584,6 +584,7 @@ double qtl_threshold_all_chromosomes(const std::string& trait_name, double quant
for (int i = 0; i < v->innerSize(); ++i) {
maxima(sz + i) = (*v)(i);
}
sz += v->innerSize();
}
return get_quantiles(maxima, {quantile})[0];
}
......
......@@ -108,6 +108,11 @@ int main(int argc, const char** argv)
value<double>(active_settings->qtl_threshold_quantile),
value<int>(active_settings->n_permutations));
if (*thresholds[0] == 0) {
std::cout << "FOIRURE" << std::endl;
return -1;
}
active_settings->set_title("Species Perscrutandis Enixe Locis Locabuntur");
for (auto& qc: active_settings->working_set) {
......@@ -120,10 +125,6 @@ int main(int argc, const char** argv)
/*all_traits(), qtl_chr,*/
/*value<double>(active_settings->qtl_threshold_quantile),*/
/*value<int>(active_settings->n_permutations));*/
if (*thresholds[0] == 0) {
std::cout << "FOIRURE" << std::endl;
return -1;
}
collection<model_block_type> popl = compute_parental_origins_multipop(all_pop, *qtl_chr, *testpos);
computation_along_chromosome cac;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment