num4anova 0.0.3-java → 0.0.4-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 89f91d9c210bb998f14a726f1c7147470cba19fba721724625b7ec87e4c88ad9
-  data.tar.gz: a79f6e3ca4adc650193065ae62037cfed0580e5f7b9c5e1d06109d74f17a05d6
+  metadata.gz: ca0164a7092861bfca67f257129332375e90bfae762365194ced78ef75a293a7
+  data.tar.gz: 16f17fb717e3d92c27c4c651e034704cd81609a3adc293555eac54f9209884be
 SHA512:
-  metadata.gz: 97a1409dfb791766292e8d6c0a49502555911529632d45b72aa6c68bc6fda7e82ed813fff4c005c9c76745a2ad05fe79d1131e23879da18afcc88d0652617c29
-  data.tar.gz: 4cd28879971daf44d20ff8122120327e8a74519f658a22e052724af316387c057cfb315a6d1cb23324bbd08634c066c1045594d6c8804f1f87e12d31b6acd0fb
+  metadata.gz: e14efb373cde03fc5872f3dfb1d4186e5d2263b0b7597970d1629c562cfb6696ecf2cdbfc3325aa2d732e252241de093c156b7d9da21b51a8ea305f3e24e0e0f
+  data.tar.gz: c9b7f310d00d060d1c6fee3e8e9f5c43fce812a65d51b192315c72a378de7e2bc775e1118346373a93047c1cd8880240d03d52309ec0d8f567481923f14e31b4

data/CHANGELOG.md

@@ -2,10 +2,15 @@
 
 ## Unreleased
 
+## [0.0.4] - 2024-01-30
+### add
+- add function of replicate_plot
+- add function of replicate_test
+
 ## [0.0.3] - 2024-01-26
 
 ### add
-- add function of dunnet_test.
+- add function of dunnet_test
 
 ## [0.0.2] - 2024-01-23
 

data/ext/num4anova/MultiComp.java

@@ -249,6 +249,7 @@ public class MultiComp {
             private double[] n = null;
             protected int getK() { return k;}
             protected int getV() { return v;}
+
             public double[][] calcTestStatistic(double[][] xi) {
                 k = xi.length;
                 mean = new double[k];
@@ -267,8 +268,10 @@ public class MultiComp {
             private double calcVe(double[][] xi) {
                 double sumSq = 0.0;
                 int sumN = 0;
+
                 for(int i = 0; i < k; i++) {
                     DescriptiveStatistics stat = new DescriptiveStatistics();
+
                     Arrays.stream(xi[i]).forEach(stat::addValue);
                     mean[i] = stat.getMean();
                     n[i] = stat.getN();
@@ -286,7 +289,7 @@ public class MultiComp {
                 int v = super.getV();
                 int k = super.getK();
                 double den = k - 1;
-                double p = 1.0 - a / den;
+
                 TDistribution tDist = new TDistribution(v);
                 double l_val = tDist.inverseCumulativeProbability(a  / den);
                 double r_val = tDist.inverseCumulativeProbability(1.0 - a / den);

data/ext/num4anova/OneWayLayout.java

@@ -28,6 +28,7 @@ import java.util.ArrayList;
 import org.apache.commons.math3.stat.inference.OneWayAnova;
 import org.apache.commons.math3.stat.descriptive.DescriptiveStatistics;
 import org.apache.commons.math3.distribution.ChiSquaredDistribution;
+import org.apache.commons.math3.distribution.FDistribution;
 import java.util.Map;
 public class OneWayLayout {
     private static OneWayLayout oneWay = new OneWayLayout();
@@ -59,6 +60,18 @@ public class OneWayLayout {
         double statistic = oneway.calcTestStatistic(xi);
         return oneway.test(statistic, a);
     }
+    public void replicatePlot(String dname, Map<String, double[]> vals) {
+        ChartPlot plot = new ReplicateChartPlot();
+
+        JFreeChart chart = plot.createChart("反復測定", dname, vals);
+        plot.writeJPEG("replicate.jpeg", chart, 800, 500);
+    }
+    public boolean replicateTest(double[][] xi, double a) {
+        OneWayAnovaTest oneway = new ReplicateTest();
+
+        double statistic = oneway.calcTestStatistic(xi);
+        return oneway.test(statistic, a);
+    }
     /*********************************/
     /* interface define              */
     /*********************************/
@@ -196,6 +209,7 @@ public class OneWayLayout {
             }
         }
     }
+    // バートレット検定
     private class BartletTest implements OneWayAnovaTest {
         private int n = 0;
         public double calcTestStatistic(double[][] xi) {
@@ -226,6 +240,7 @@ public class OneWayLayout {
             double invSumN = 0.0;
             int sumN = 0;
             DescriptiveStatistics stat = new DescriptiveStatistics();
+
             for(int i = 0; i < n; i++) {
                 Arrays.stream(xi[i]).forEach(stat::addValue);
                 invSumN += 1.0 / (stat.getN() - 1.0);
@@ -243,5 +258,127 @@ public class OneWayLayout {
             return (r_val < statistic) ? true : false;
         }
     }
+    // 反復測定Plot
+    private class ReplicateChartPlot implements ChartPlot {
+        public JFreeChart createChart(String title, String dname, Map<String, double[]> vals) {
+            CategoryPlot plot = createPlot(dname, vals);
+            ChartFactory.setChartTheme(StandardChartTheme.createLegacyTheme());
+            JFreeChart chart = new JFreeChart(title, plot);
+
+            ChartUtils.applyCurrentTheme(chart);
+            return chart;
+        }
+        private CategoryPlot createPlot( String dname, Map<String, double[]> vals) {
+            CreatePlot plotImpl = new ReplicatePlot();
+
+            return plotImpl.createPlot(dname, vals);
+        }
+        private class ReplicatePlot implements CreatePlot {
+            public CategoryPlot createPlot(String dname, Map<String, double[]> vals) {
+                LineAndShapeRenderer renderer = new LineAndShapeRenderer(true, true);
+
+                renderer.setDefaultToolTipGenerator(
+                    new StandardCategoryToolTipGenerator()
+                );
+                CategoryPlot plot = new CategoryPlot();
+
+                plot.setOrientation(PlotOrientation.VERTICAL);
+                plot.mapDatasetToRangeAxis(0,0);
+	        plot.setDatasetRenderingOrder(DatasetRenderingOrder.FORWARD);
+                
+                /*--- 横軸 ---*/
+                CategoryAxis categoryAxis = new CategoryAxis("因子");
+                plot.setDomainAxis(categoryAxis);
+
+                /*--- 縦軸 ---*/
+                NumberAxis valueAxis0 = new NumberAxis("推定周辺平均");
+                plot.setRangeAxis(valueAxis0);
+
+                plot.setRenderer(0, renderer);
+                plot.setDataset(0, createDataset(dname, vals));
+                return plot;
+            }
+            private CategoryDataset createDataset(String dname, Map<String, double[]> vals) {
+                DefaultCategoryDataset data = new DefaultCategoryDataset();
+                DescriptiveStatistics stat = new DescriptiveStatistics();
+
+                for(Map.Entry<String, double[]> entry : vals.entrySet()) {
+                    double[] v = entry.getValue();
+
+                    Arrays.stream(v).forEach(stat::addValue);
+                    data.addValue(stat.getMean(), dname, entry.getKey());
+                    stat.clear();
+                }
+                return data;
+            }
+        }
+    }
+    private class ReplicateTest implements OneWayAnovaTest {
+        private int a1 = 0;
+        private int b1 = 0;
+        public double calcTestStatistic(double[][] xi) {
+            b1 = xi[0].length;
+            a1 = xi.length;
+            double st = calcSt(xi, a1, b1);
+            double sa = calcSa(xi, a1, b1);
+            double sb = calcSb(xi, a1, b1);       
+            double se = st - sa - sb;
+            double meanSa = sa / (a1 - 1);
+            double meanSe = se / ((a1 - 1) * (b1 - 1));
+
+            return meanSa / meanSe ;
+        }
+        private double calcSt(double[][] xi, int a, int b) {
+            DescriptiveStatistics stat = new DescriptiveStatistics();
+            double sumSt1 = 0.0;
+            double sumSt2 = 0.0;
+
+            for(int i = 0; i < a; i++) {
+                Arrays.stream(xi[i]).forEach(stat::addValue);
+                sumSt1 += stat.getSumsq();
+                sumSt2 += stat.getSum();
+                stat.clear();
+            }
+            return sumSt1 - sumSt2 * sumSt2 / (a * b);
+        }
+        private double calcSa(double[][] xi, int a, int b) {
+            double sumSa1 = 0.0;
+            double sumSa2 = 0.0;
+            double[] an = new double[a];
+
+            for(int i = 0; i < a; i++) {
+                for(int j = 0; j < b; j++) {
+                    an[i] += xi[i][j];
+                    sumSa2 += xi[i][j];
+                }
+            }
+            DescriptiveStatistics stat = new DescriptiveStatistics();
+            Arrays.stream(an).forEach(stat::addValue);
+            sumSa1 = stat.getSumsq() / b;
+            return sumSa1 - sumSa2 * sumSa2 / (a * b);
+        }
+        private double calcSb(double[][] xi, int a, int b) {
+            double[] bn = new double[b];
+            double sumSb1 = 0.0;
+            double sumSb2 = 0.0;
+
+            for(int i = 0; i < a; i++) {
+                for(int j = 0; j < b; j++) {
+                    bn[j] += xi[i][j];
+                    sumSb2 += xi[i][j];
+                }
+            }
+            DescriptiveStatistics stat = new DescriptiveStatistics();
+            Arrays.stream(bn).forEach(stat::addValue);
+            sumSb1 = stat.getSumsq() / a;
+            return sumSb1 - sumSb2 * sumSb2 / (a * b);
+        }
+        public boolean test(double statistic, double a) {
+            FDistribution fDist = new FDistribution(a1 - 1, (a1 - 1) * (b1 - 1));
+            double f = fDist.inverseCumulativeProbability(1.0 - a);
+            
+            return (statistic >= f) ? true : false;
+        }
+    }
 }
 

data/lib/num4anova.rb

@@ -72,13 +72,13 @@ module Num4AnovaLib
         #   @param [double] a         有意水準
         #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
         # @example
-        #    xi = [
-        #        [12.2, 18.8, 18.2],
-        #        [22.2, 20.5, 14.6],
-        #        [20.8, 19.5, 26.3],
-        #        [26.4, 32.5, 31.3],
-        #        [24.5, 21.2, 22.4],
-        #    ]
+        #   xi = [
+        #       [12.2, 18.8, 18.2],
+        #       [22.2, 20.5, 14.6],
+        #       [20.8, 19.5, 26.3],
+        #       [26.4, 32.5, 31.3],
+        #       [24.5, 21.2, 22.4],
+        #   ]
         #   oneWay = Num4AnovaLib::OneWayLayoutLib.new 
         #   oneWay.oneWay.oneway_anova(xi, 0.05)
         #   => true
@@ -92,19 +92,63 @@ module Num4AnovaLib
         #   @param [double] a         有意水準
         #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
         # @example
-        #    xi = [
-        #        [12.2, 18.8, 18.2],
-        #        [22.2, 20.5, 14.6],
-        #        [20.8, 19.5, 26.3],
-        #        [26.4, 32.5, 31.3],
-        #        [24.5, 21.2, 22.4],
-        #    ]
+        #   xi = [
+        #       [12.2, 18.8, 18.2],
+        #       [22.2, 20.5, 14.6],
+        #       [20.8, 19.5, 26.3],
+        #       [26.4, 32.5, 31.3],
+        #       [24.5, 21.2, 22.4],
+        #   ]
         #   oneWay = Num4AnovaLib::OneWayLayoutLib.new 
         #   oneWay.bartlet(xi, 0.05)
         #   => true
         def bartlet(xi, a)
             return @oneWay.bartletTest(xi.to_java(Java::double[]), a)
         end
+        # 反復測定Plot
+        #
+        # @overload replicate_plot(dname, vals)
+        #   @param [String] dname データ名
+        #   @param [Hash] vals Hash(String, double[])
+        #   @return [void]  replicate.jpegファイルを出力
+        # @example
+        #   vals = {
+        #        "stageA1" => [27, 52, 18, 21, 32],
+        #        "stageA2" => [52, 72, 31, 50, 45],
+        #        "stageA3" => [47, 54, 29, 43, 32],
+        #        "stageA4" => [28, 50, 22, 26, 29],
+        #       }
+        #   oneWay = Num4AnovaLib::OneWayLayoutLib.new 
+        #   oneWay.replicate_plot("LDH", vals)
+        #   => replicate.jpeg
+        # @note
+        #   グラフは、jfreechartを使用
+        def replicate_plot(dname, vals)
+            o = HashMap.new
+            vals.each{|k, v|
+                o[k] = v.to_java(Java::double)
+            }
+            return @oneWay.replicatePlot(dname, o)
+        end
+        # 反復測定検定
+        #
+        # @overload replicate_test(xi, a)
+        #   @param [array]  xi データ(double[][])
+        #   @param [double] a         有意水準
+        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
+        # @example
+        #   xi = [
+        #       [27, 52, 18, 21, 32],
+        #       [52, 72, 31, 50, 45],
+        #       [47, 54, 29, 43, 32],
+        #       [28, 50, 22, 26, 29],
+        #   ]
+        #   oneWay = Num4AnovaLib::OneWayLayoutLib.new 
+        #   oneWay.replicate_test("LDH", vals)
+        #   => true
+        def replicate_test(xi, a)
+            return @oneWay.replicateTest(xi.to_java(Java::double[]), a)
+        end
     end
 end
 

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4anova
 version: !ruby/object:Gem::Version
-  version: 0.0.3
+  version: 0.0.4
 platform: java
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-01-26 00:00:00.000000000 Z
+date: 2024-01-30 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake