my $self = shift;
my $num = shift;
my $gsb = shift;
my $mappings = shift;
my $gene_scores = shift;
unless ($gene_scores->loaded) {
$self->logger->info("Synteny Framework building...\n", 0, 'stamped');
my $dump_path = path_append($self->conf->param('basedir'), 'mapping');
-DUMP_PATH => $dump_path,
-CACHE_FILE => 'synteny_framework.ser',
-LOGGER => $self->logger,
-CONF => $self->conf,
-CACHE => $self->cache,
);
# use it to rescore the genes
$self->logger->info("\nSynteny assisted mapping...\n", 0, 'stamped');
$gene_scores = $sf->rescore_gene_matrix_lsf($gene_scores);
# checkpoint
$gene_scores->write_to_file;
}
my $new_mappings = $self->basic_mapping($gene_scores, "gene_mappings$num");
$num++;
my $new_scores = $gsb->create_shrinked_matrix($gene_scores, $new_mappings,
"gene_matrix$num");
return ($new_scores, $new_mappings);
}