WIP

Author: choroba

t/01-deltas.t

@@ -4,8 +4,10 @@ use strict;
 
 use Statistics::Krippendorff;
 
+use List::Util qw{ sum };
+
 use Test2::V0;
-plan 16;
+plan 20;
 
 my $sk = 'Statistics::Krippendorff'->new(units => []);
 
@@ -29,3 +31,25 @@ is $sk->delta_jaccard('b,a', 'c,a'), 2/3, 'ba ca';
 is $sk->delta_jaccard('b,a', 'a,c'), 2/3, 'ba ac';
 
 is $sk->delta_jaccard('a,b,c', 'b,c,d,e'), 3/5, 'abc bcde';
+
+my @figures = (
+    [['x,y', 'x,y,z'], ['x,y', 'x,y,z'], ['x', 'x,y,z']],
+    [['x,y', 'x'], ['x,y', 'y,z'], ['z', 'y,z']]
+);
+
+my %means = (jaccard => [4 / 9, 5 / 9],
+             masi    => [10 / 27, 6 / 27]);
+for my $figure_index (0 .. $#figures) {
+    my $figure = $figures[$figure_index];
+    my @jaccard_deltas = map $sk->delta_jaccard(@$_), @$figure;
+    my $mean = sum(@jaccard_deltas) / @jaccard_deltas;
+    is $mean, shift @{ $means{jaccard} }, "fig$figure_index jaccard mean";
+
+    if (1 == $figure_index) {
+        $figure->[1] = ['x', 'z'];
+    }
+
+    my @masi_deltas = map $sk->delta_masi(@$_), @$figure;
+    my $mean = 1 - sum(@masi_deltas) / @masi_deltas;
+    is $mean, shift @{ $means{masi} }, "fig$figure_index masi mean";
+}