14 files changed, 273 insertions, 38 deletions

diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/Domain.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/Domain.xtend
new file mode 100644
index 00000000..8351e96b
--- /dev/null
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/Domain.xtend
@@ -0,0 +1,5 @@
+package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.app
+
+enum Domain{
+        Yakinduum
+}
\ No newline at end of file
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/Main.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/Main.xtend
index 15df0dde..cf871ead 100644
--- a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/Main.xtend
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/Main.xtend
@@ -1,8 +1,8 @@
 package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.app
 
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.graph.EMFGraph
-import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.input.GraphReader
-import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.output.CsvFileWriter
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.io.CsvFileWriter
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.io.GraphReader
 import hu.bme.mit.inf.dslreasoner.domains.yakindu.sgraph.yakindumm.impl.YakindummPackageImpl
 import java.util.ArrayList
 
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/PartialInterpretationMetric.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/PartialInterpretationMetric.xtend
index 0e505d30..cdd06027 100644
--- a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/PartialInterpretationMetric.xtend
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/app/PartialInterpretationMetric.xtend
@@ -1,11 +1,12 @@
 package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.app
 
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.distance.KSDistance
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.graph.PartialInterpretationGraph
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.io.CsvFileWriter
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.Metric
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.MultiplexParticipationCoefficientMetric
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.NodeActivityMetric
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.OutDegreeMetric
-import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.output.CsvFileWriter
 import hu.bme.mit.inf.dslreasoner.viatrasolver.partialinterpretationlanguage.partialinterpretation.PartialInterpretation
 import java.io.File
 import java.io.FileNotFoundException
@@ -17,10 +18,28 @@ import org.eclipse.viatra.dse.api.Solution
 
 class PartialInterpretationMetric {
         var static state = 0;
+        var static KSDistance ks;
         
         def static void initPaths(){
                 new File("debug/metric/").mkdir();
                 new File("debug/metric/trajectories/").mkdir();
+                ks = new KSDistance(Domain.Yakinduum);
+        }
+        
+        def static MetricDistanceGroup calculateMetricDistance(PartialInterpretation partial){
+                val metrics = new ArrayList<Metric>();
+                metrics.add(new OutDegreeMetric());
+                metrics.add(new NodeActivityMetric());
+                metrics.add(new MultiplexParticipationCoefficientMetric());
+                
+                val metricCalculator = new PartialInterpretationGraph(partial, metrics, null);          
+                var metricSamples = metricCalculator.evaluateAllMetricsToSamples();
+                
+                var mpc = ks.mpcDistance(metricSamples.mpcSamples);
+                var na = ks.naDistance(metricSamples.naSamples);
+                var outDegree = ks.outDegreeDistance(metricSamples.outDegreeSamples);
+                
+                return new MetricDistanceGroup(mpc, na, outDegree);
         }
         
         // calculate the metrics for a state
@@ -78,4 +97,28 @@ class PartialInterpretationMetric {
                         }
                 }                       
         }
+}
+
+class MetricDistanceGroup{
+        var double mpcDistance;
+        var double naDistance;
+        var double outDegreeDistance;
+        
+        new(double mpcDistance, double naDistance, double outDegreeDistance){
+                this.mpcDistance = mpcDistance;
+                this.naDistance = naDistance;
+                this.outDegreeDistance = outDegreeDistance;
+        }
+        
+        def double getMPCDistance(){
+                return this.mpcDistance
+        }
+        
+        def double getNADistance(){
+                return this.naDistance
+        }
+        
+        def double getOutDegreeDistance(){
+                return this.outDegreeDistance
+        }
 }
\ No newline at end of file
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/distance/KSDistance.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/distance/KSDistance.xtend
new file mode 100644
index 00000000..1fb21529
--- /dev/null
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/distance/KSDistance.xtend
@@ -0,0 +1,47 @@
+package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.distance
+
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.app.Domain
+import java.util.List
+import org.apache.commons.math3.stat.inference.KolmogorovSmirnovTest
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.io.RepMetricsReader
+
+class KSDistance {
+        var static ksTester = new KolmogorovSmirnovTest();
+        var double[] mpcSamples;
+        var double[] naSamples;
+        var double[] outDegreeSamples;
+        
+        new(Domain d){
+                var metrics = RepMetricsReader.read(d);
+                mpcSamples = metrics.mpcSamples;
+                naSamples = metrics.naSamples.stream.mapToDouble([it]).toArray();
+                outDegreeSamples = metrics.outDegreeSamples.stream.mapToDouble([it]).toArray();
+        }
+        
+        def double mpcDistance(List<Double> samples){
+                // map list to array
+                var arr = samples.stream.mapToDouble([it]).toArray();
+                
+                //if the size of array is smaller than 2, ks distance cannot be performed, simply return 1
+                if(arr.size < 2) return 1;
+                return ksTester.kolmogorovSmirnovStatistic(mpcSamples, arr);
+        }
+        
+        def double naDistance(List<Double> samples){
+                // map list to array
+                var arr = samples.stream.mapToDouble([it]).toArray();
+                
+                //if the size of array is smaller than 2, ks distance cannot be performed, simply return 1
+                if(arr.size < 2) return 1;
+                return ksTester.kolmogorovSmirnovStatistic(naSamples as double[], arr);
+        }
+        
+        def double outDegreeDistance(List<Double> samples){
+                // map list to array
+                var arr = samples.stream.mapToDouble([it]).toArray();
+                
+                //if the size of array is smaller than 2, ks distance cannot be performed, simply return 1
+                if(arr.size < 2) return 1;
+                return ksTester.kolmogorovSmirnovStatistic(outDegreeSamples, arr);
+        }
+}
\ No newline at end of file
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/graph/Graph.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/graph/Graph.xtend
index 6d65367f..cf4aedba 100644
--- a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/graph/Graph.xtend
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/graph/Graph.xtend
@@ -1,8 +1,12 @@
 package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.graph
 
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.Metric
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.MetricSampleGroup
 import java.util.ArrayList
 import java.util.List
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.MultiplexParticipationCoefficientMetric
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.NodeActivityMetric
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.OutDegreeMetric
 
 abstract class Graph {
         
@@ -33,6 +37,22 @@ abstract class Graph {
                 return result;
         }       
         
+        def MetricSampleGroup evaluateAllMetricsToSamples(){
+                var sample = new MetricSampleGroup();
+                
+                for(metric : this.metrics){
+                        if(metric instanceof MultiplexParticipationCoefficientMetric){
+                                sample.mpcSamples = metric.evaluateSamples(this.statistic);                             
+                        }else if(metric instanceof NodeActivityMetric){
+                                sample.naSamples = metric.evaluateSamples(this.statistic);
+                        }else if(metric instanceof OutDegreeMetric){
+                                sample.outDegreeSamples = metric.evaluateSamples(this.statistic);
+                        }
+                }
+                
+                return sample;
+        }
+        
         def void setBasicInformation(ArrayList<ArrayList<String>> result);
         
         def GraphStatistic getStatistic();
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/input/PartialInterpretationReader.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/input/PartialInterpretationReader.xtend
deleted file mode 100644
index dbaf36b2..00000000
--- a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/input/PartialInterpretationReader.xtend
+++ /dev/null
@@ -1,7 +0,0 @@
-package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.input
-
-import hu.bme.mit.inf.dslreasoner.viatrasolver.partialinterpretationlanguage.partialinterpretation.PartialInterpretation
-
-class PartialInterpretationReader {
-        
-}
\ No newline at end of file
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/output/CsvFileWriter.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/io/CsvFileWriter.xtend
index 1dd204a0..bed356e9 100644
--- a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/output/CsvFileWriter.xtend
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/io/CsvFileWriter.xtend
@@ -1,4 +1,4 @@
-package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.output;
+package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.io;
 
 import java.io.File
 import java.io.FileNotFoundException
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/input/GraphReader.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/io/GraphReader.xtend
index 56bf95a3..fc56e142 100644
--- a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/input/GraphReader.xtend
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/io/GraphReader.xtend
@@ -1,4 +1,4 @@
-package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.input
+package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.io;
 
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.graph.EMFGraph
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.Metric
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/io/RepMetricsReader.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/io/RepMetricsReader.xtend
new file mode 100644
index 00000000..867ddd1a
--- /dev/null
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/io/RepMetricsReader.xtend
@@ -0,0 +1,79 @@
+package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.io
+
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.app.Domain
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.MultiplexParticipationCoefficientMetric
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.NodeActivityMetric
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.OutDegreeMetric
+import java.io.File
+import java.util.List
+import java.util.Scanner
+import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics.MetricSampleGroup
+
+/**
+ * Read the sample of the distribution of a metric provided the csv file of the metric
+ */
+class RepMetricsReader {
+        static def read(Domain d){
+                var domainRepPath = CsvDataName.REP_PATH + d.name + '/';
+                var rep = new MetricSampleGroup()
+                rep.mpcSamples = readFile(domainRepPath + CsvDataName.MPC_REP, MultiplexParticipationCoefficientMetric.valueName,
+                                                          MultiplexParticipationCoefficientMetric.countName).map[Double.parseDouble(it)];
+                rep.naSamples = readFile(domainRepPath+CsvDataName.NA_REP, NodeActivityMetric.valueName, NodeActivityMetric.countName
+                                                        ).map[Double.parseDouble(it)];
+                rep.outDegreeSamples = readFile(domainRepPath+CsvDataName.OUT_D_REP, OutDegreeMetric.valueName, OutDegreeMetric.countName
+                                                        ).map[Double.parseDouble(it)];
+                return rep;
+        }
+        
+        /**
+         * read metric data and parse it to samples
+         */
+        private static def List<String> readFile(String filename, String valueDataName, String countDataName){
+                var s = new Scanner(new File(filename));
+                val counts = newArrayList();
+                val values = newArrayList();
+                //read data from csv
+                while(s.hasNext()){
+                        var data = s.nextLine().split(',');     
+                        
+                        if(data.size >= 1){
+                                if(data.get(0).equals(countDataName)){
+                                        //add all data with parsing them as integers
+                                        counts.addAll(data.subList(1, data.size()).map[Integer.parseInt(it)]);
+                                }else if(data.get(0).equals(valueDataName)){
+                                        //add all data without parsing (there can be either double or string, to be parsed later)
+                                        values.addAll(data.subList(1, data.size()));
+                                }
+                        }
+                }
+                
+                return createSamples(counts, values);
+        }
+        
+        // create samples from values and counts
+        private static def List<String> createSamples(List<Integer> counts, List<String> values){
+                val samples = newArrayList();
+                
+                if(counts.size() != values.size()){
+                        throw new RuntimeException("counts and values should have the same size!");
+                }
+                
+                for(var i = 0; i < counts.size(); i++){
+                        for(var j = 0; j < counts.get(i); j++){
+                                samples.add(values.get(i));
+                        }
+                }
+                
+                return samples;
+        }
+        
+        
+}
+
+class CsvDataName{
+        public static val REP_PATH = 'data/';
+        public static val MPC_REP = 'mpc_rep.csv';
+        public static val NA_REP = 'na_rep.csv';
+        public static val OUT_D_REP = 'out_d_rep.csv';
+}
+
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/Metric.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/Metric.xtend
index ea52009a..38ef72f2 100644
--- a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/Metric.xtend
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/Metric.xtend
@@ -4,4 +4,5 @@ import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.graph.GraphStatis
 
 abstract class Metric {
         abstract def String[][] evaluate(GraphStatistic g);
+        abstract def double[] evaluateSamples(GraphStatistic g);
 }
\ No newline at end of file
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/MetricSampleGroup.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/MetricSampleGroup.xtend
new file mode 100644
index 00000000..8cd3daee
--- /dev/null
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/MetricSampleGroup.xtend
@@ -0,0 +1,9 @@
+package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics
+
+import java.util.List
+
+class MetricSampleGroup{
+        public var List<Double> mpcSamples;
+        public var List<Double> naSamples;
+        public var List<Double> outDegreeSamples;
+}
\ No newline at end of file
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/MultiplexParticipationCoefficientMetric.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/MultiplexParticipationCoefficientMetric.xtend
index 7bd48754..d9c88bb4 100644
--- a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/MultiplexParticipationCoefficientMetric.xtend
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/MultiplexParticipationCoefficientMetric.xtend
@@ -2,11 +2,13 @@ package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics
 
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.graph.GraphStatistic
 import java.text.DecimalFormat
+import java.util.ArrayList
 import java.util.HashMap
+import org.eclipse.emf.ecore.EObject
 
 class MultiplexParticipationCoefficientMetric extends Metric {
-        static val countName = "MPCCount";
-        static val valueName = "MPCValue";
+        public static val countName = "MPCCount";
+        public static val valueName = "MPCValue";
         
         
         override evaluate(GraphStatistic g) {
@@ -18,26 +20,7 @@ class MultiplexParticipationCoefficientMetric extends Metric {
                 val map = new HashMap<String, Integer>();
                 //calculate the metric distribution
                 g.allNodes.forEach[n|
-                        val edgeCounts = g.outDegree(n) + g.inDegree(n);
-                        
-                        var coef = 0.0;
-                        
-                        for(type : g.allTypes){
-                                val degree = g.dimentionalDegree(n, type) as double;
-                                coef += Math.pow(degree / edgeCounts, 2);
-                        }               
-                        coef = 1 - coef;                
-                        coef = coef * typeCounts / (typeCounts-1);
-                        
-                        //Consider the case that either typeCounts-1 or the edgeCounts could be 0 in some situation
-                        //in this case the metric should be evaluated to 0
-                        if(typeCounts == 1){
-                                println("bad");
-                        }
-                        
-                        if(Double.isNaN(coef)){
-                                coef = 0;
-                        }
+                        var coef = calculateMPC(n, g, typeCounts);
                         
                         //format number to String
                         val value = formatter.format(coef);
@@ -62,4 +45,36 @@ class MultiplexParticipationCoefficientMetric extends Metric {
                 
                 return datas;
         }
+        
+        override evaluateSamples(GraphStatistic g){
+                val samples = new ArrayList<Double>();
+                val typeCounts = g.allTypes.size;
+                //calculate the metric distribution
+                g.allNodes.forEach[
+                        samples.add(calculateMPC(it, g, typeCounts));
+                ]
+                
+                return samples;
+        }
+        
+        def double calculateMPC(EObject n, GraphStatistic g, int typeCounts){
+                val edgeCounts = g.outDegree(n) + g.inDegree(n);
+                        
+                var coef = 0.0;
+                
+                for(type : g.allTypes){
+                        val degree = g.dimentionalDegree(n, type) as double;
+                        coef += Math.pow(degree / edgeCounts, 2);
+                }               
+                coef = 1 - coef;                
+                coef = coef * typeCounts / (typeCounts-1);
+                
+                //Consider the case that either typeCounts-1 or the edgeCounts could be 0 in some situation
+                //in this case the metric should be evaluated to 0
+                if(Double.isNaN(coef)){
+                        coef = 0;
+                }
+                
+                return coef;
+        }
 }
\ No newline at end of file
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/NodeActivityMetric.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/NodeActivityMetric.xtend
index 297bdd14..fbf06c47 100644
--- a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/NodeActivityMetric.xtend
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/NodeActivityMetric.xtend
@@ -1,11 +1,12 @@
 package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics
 
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.graph.GraphStatistic
+import java.util.ArrayList
 import java.util.HashMap
 
 class NodeActivityMetric extends Metric {
-        static val countName = "NACount";
-        static val valueName = "NAValue";
+        public static val countName = "NACount";
+        public static val valueName = "NAValue";
         
         override evaluate(GraphStatistic g) {
                 val map = new HashMap<Integer, Integer>();
@@ -34,4 +35,15 @@ class NodeActivityMetric extends Metric {
                 
                 return datas;
         }
+        
+        override evaluateSamples(GraphStatistic g){
+                val samples = new ArrayList<Double>();
+                
+                //calculate the metric distribution
+                g.allNodes.forEach[
+                        samples.add(g.numOfEdgeTypes(it) as double);
+                ]
+                
+                return samples;
+        }
 }
\ No newline at end of file
diff --git a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/OutDegreeMetric.xtend b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/OutDegreeMetric.xtend
index 1ba6c8c7..55046b14 100644
--- a/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/OutDegreeMetric.xtend
+++ b/Metrics/Metrics-Calculation/ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator/src/ca/mcgill/ecse/dslreasoner/realistic/metrics/calculator/metrics/OutDegreeMetric.xtend
@@ -1,11 +1,12 @@
 package ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.metrics
 
 import ca.mcgill.ecse.dslreasoner.realistic.metrics.calculator.graph.GraphStatistic
+import java.util.ArrayList
 import java.util.HashMap
 
 class OutDegreeMetric extends Metric {
-        static val countName = "OutDegreeCount";
-        static val valueName = "OutDegreeValue";
+        public static val countName = "OutDegreeCount";
+        public static val valueName = "OutDegreeValue";
         
         override evaluate(GraphStatistic g) {
                 val map = new HashMap<Integer, Integer>();
@@ -34,5 +35,15 @@ class OutDegreeMetric extends Metric {
                 return datas;
         }
         
+        override evaluateSamples(GraphStatistic g){
+                val samples = new ArrayList<Double>();
+                
+                //calculate the metric distribution
+                g.allNodes.forEach[
+                        samples.add(g.outDegree(it) as double);
+                ]
+                
+                return samples;
+        }
         
 }
\ No newline at end of file