checkpoint

BEAST2-Dev · May 4, 2015 · cc059e6 · cc059e6
1 parent 841f4cb
commit cc059e6
Showing 1 changed file with 16 additions and 0 deletions.
diff --git a/src/snap/Data.java b/src/snap/Data.java
@@ -500,4 +500,20 @@ protected boolean isEqual(int iSite1, int iSite2) {
     public int[] getPatternLineagCounts(int id) {
         return m_nPatternLineageCounts[id];
     }
+
+    public double getProportionZeros() {
+    	int zeroCount = 0;
+    	int oneCount  = 0;
+    	for (int i = 0; i < sitePatterns.length; i++) {
+    		int [] p = getPattern(i);
+    		int [] n = getPatternLineagCounts(i);
+    		int w = getPatternWeight(i);
+    		for (int j = 0; j < p.length; j++) {
+    			oneCount += w * p[j];
+    			zeroCount += w * (n[j] - p[j]);
+    		}
+    	}
+    	return (double) zeroCount / ((double) zeroCount + oneCount);
+    }
+
 }