Skip to content
Snippets Groups Projects
Commit 42e6fc2e authored by Andreas Kusalananda Kähäri's avatar Andreas Kusalananda Kähäri
Browse files

Update some penalties and comments.

parent ab515d0e
No related branches found
No related tags found
No related merge requests found
...@@ -432,10 +432,9 @@ sub biotype_gene_rescore { ...@@ -432,10 +432,9 @@ sub biotype_gene_rescore {
$entry->target ); $entry->target );
if ( $source_gene->biotype() ne $target_gene->biotype() ) { if ( $source_gene->biotype() ne $target_gene->biotype() ) {
# PENALTY: The gene stable ID is now on a gene with a different # PENALTY: Lower the score for mappings that differ in biotype.
# biotype.
$matrix->set_score( $entry->source(), $entry->target(), $matrix->set_score( $entry->source(), $entry->target(),
0.75*$entry->score() ); 0.9*$entry->score() );
$i++; $i++;
} }
} }
......
...@@ -381,10 +381,10 @@ sub biotype_transcript_rescore { ...@@ -381,10 +381,10 @@ sub biotype_transcript_rescore {
if ($source_transcript->biotype() ne $target_transcript->biotype() ) if ($source_transcript->biotype() ne $target_transcript->biotype() )
{ {
# PENALTY: The transcript stable ID is now on a transcript with a # PENALTY: Lower the score for mappings to transcripts of
# different biotype. # different biotype.
$matrix->set_score( $entry->source(), $entry->target(), $matrix->set_score( $entry->source(), $entry->target(),
0.75*$entry->score() ); 0.9*$entry->score() );
$i++; $i++;
} }
} }
...@@ -430,7 +430,7 @@ sub different_translation_rescore { ...@@ -430,7 +430,7 @@ sub different_translation_rescore {
# PENALTY: The transcript stable ID is now on a transcript with a # PENALTY: The transcript stable ID is now on a transcript with a
# different translation. # different translation.
$matrix->set_score( $entry->source(), $entry->target(), $matrix->set_score( $entry->source(), $entry->target(),
0.98*$entry->score() ); 0.95*$entry->score() );
$i++; $i++;
} }
...@@ -481,7 +481,8 @@ sub non_mapped_gene_rescore { ...@@ -481,7 +481,8 @@ sub non_mapped_gene_rescore {
my $mapped_target = $gene_lookup{ $source_gene->id }; my $mapped_target = $gene_lookup{ $source_gene->id };
if ( !$mapped_target or ( $mapped_target != $target_gene->id ) ) { if ( !$mapped_target or ( $mapped_target != $target_gene->id ) ) {
# PENALTY: The transcript stable ID is now on a different gene. # PENALTY: The transcript stable ID has been mapped to an
# un-mapped gene.
$matrix->set_score( $entry->source(), $entry->target(), $matrix->set_score( $entry->source(), $entry->target(),
0.8*$entry->score() ); 0.8*$entry->score() );
$i++; $i++;
......
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