forester/java/src/org/forester/surfacing/CountsBasedPairwiseDomainSimilarity.java

   1 // $Id:
   2 // cmzmasek Exp $
   3 //
   4 // FORESTER -- software libraries and applications
   5 // for evolutionary biology research and applications.
   6 //
   7 // Copyright (C) 2008-2009 Christian M. Zmasek
   8 // Copyright (C) 2008-2009 Burnham Institute for Medical Research
   9 // All rights reserved
  10 //
  11 // This library is free software; you can redistribute it and/or
  12 // modify it under the terms of the GNU Lesser General Public
  13 // License as published by the Free Software Foundation; either
  14 // version 2.1 of the License, or (at your option) any later version.
  15 //
  16 // This library is distributed in the hope that it will be useful,
  17 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  18 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  19 // Lesser General Public License for more details.
  20 //
  21 // You should have received a copy of the GNU Lesser General Public
  22 // License along with this library; if not, write to the Free Software
  23 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA
  24 //
  25 // Contact: phylosoft @ gmail . com
  26 // WWW: www.phylosoft.org/forester
  27
  28 package org.forester.surfacing;
  29
  30 public class CountsBasedPairwiseDomainSimilarity implements PairwiseDomainSimilarity {
  31
  32     private final double _score;
  33     private final int    _copy_number_difference;
  34
  35     /**
  36      * counts_difference: (counts for domain 1) minus (counts for domain 2).
  37      *
  38      *
  39      * @param counts_difference value of domain_1 minus value of domain_2
  40      * @param counts_sum
  41      */
  42     public CountsBasedPairwiseDomainSimilarity( final int counts_difference, final int counts_sum ) {
  43         if ( counts_sum <= 0 ) {
  44             throw new IllegalArgumentException( "attempt to use copy sum of less than or equal to 0: " + counts_sum );
  45         }
  46         _copy_number_difference = counts_difference;
  47         final int abs_copy_number_difference = Math.abs( counts_difference );
  48         if ( abs_copy_number_difference > counts_sum ) {
  49             throw new IllegalArgumentException( "attempt to use absolute copy number difference larger than copy number sum" );
  50         }
  51         _score = 1.0 - ( ( double ) abs_copy_number_difference / counts_sum );
  52     }
  53
  54     /**
  55      * Returns (counts for domain 1) minus (counts for domain 2).
  56      *
  57      */
  58     @Override
  59     public int getDifferenceInCounts() {
  60         return _copy_number_difference;
  61     }
  62
  63     @Override
  64     public double getSimilarityScore() {
  65         return _score;
  66     }
  67 }