num4anova 0.0.11-java → 0.0.13-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 2a278a4c32d7e1ef21ccbb1200573d5c1b4e326c59602327f9968afb29f134b6
-  data.tar.gz: e8520453d70e0276e6bb0587aefcd31b90b84d703cf670ae25297af73172bcdd
+  metadata.gz: 84f7f49327e10a3da88a363daa4d6d5dec1dd4da00827d6f8e5c60ee6db4bc08
+  data.tar.gz: 495495a96d218034014971770b654431bb0df61aec47f8f3171fe351285022c4
 SHA512:
-  metadata.gz: 42836d15f6460c2650a7c6cb17c15de206bbabb52c40d988714fe26564d019974a416db4c22f67ed0c535694022575d380603da7c029f0b04476d4994ffca7a5
-  data.tar.gz: ed1698ed0b05da8d2f0d7417fd1271bcb41754d76fa68b19fb300245449a634c830636ea790a49b89d87066c320c55e2260464bb791ba795f8193664a0babb28
+  metadata.gz: 1b46cc7670940eb88343959c80bad55209a981388acdd823c3d9eeafa72a629d2319cfb49707f3f998d72f0d080480785129c948c8ebedb388693476be59add1
+  data.tar.gz: 83585b2c13c1f964b8bf92e507a005deb4d358a612e23b939cb9d9911175006a1378e16719c1b73073834aa510f245b4432f0aba9f5dfcdceb78ee2cf5638ba3

data/CHANGELOG.md

@@ -2,6 +2,16 @@
 
 ## Unreleased
 
+## [0.0.13] - 2024-04-04
+
+### chg
+- chg fuction of interval_estim.
+
+## [0.0.12] - 2024-03-27
+
+### add
+- add fuction of interval_estim.
+
 ## [0.0.11] - 2024-03-21
 
 ### add

data/ext/num4anova/Ancova.java

@@ -1,4 +1,5 @@
 import org.apache.commons.math3.distribution.FDistribution;
+import org.apache.commons.math3.distribution.TDistribution;
 
 public class Ancova {
     private static Ancova ancova = new Ancova();
@@ -23,6 +24,11 @@ public class Ancova {
         double statistic = hypoth.calcTestStatistic(xi);
         return hypoth.executeTest(statistic, a);
     }
+    public Interval[] intervalEstim(double[][][] xi, double a) {
+        Estim estim = new IntervalEstim();
+
+        return estim.calcInterval(xi, a);
+    }
     /*********************************/
     /* interface define              */
     /*********************************/
@@ -30,9 +36,22 @@ public class Ancova {
         double calcTestStatistic(double[][][] xi);
         boolean executeTest(double statistic, double a);
     }
+    private interface Estim {
+        Interval[] calcInterval(double[][][] xi, double a);
+    }
     /*********************************/
     /* class define                  */
     /*********************************/
+    public class Interval {
+        private double min;
+        private double max;
+        public Interval(double min, double max) {
+            this.min = min;
+            this.max = max;
+        }
+        public double getMin() { return this.min; }
+        public double getMax() { return this.max; }
+    }
     private class RegressionLine {
         protected int calcSumn(double[][][] xi) {
             int sum = 0;
@@ -321,5 +340,93 @@ public class Ancova {
             return (sumey * sumex - sumeyx * sumeyx) / (m * sumex);
         }
     }
+    //
+    private class IntervalEstim extends RegressionLine 
+                                implements Estim {
+        private int n = 0;
+        private int[] ni = null;
+        private double sumex = 0.0;
+        public Interval[] calcInterval(double[][][] xi, double a) {
+            Interval[] ret = new Interval[xi.length];
+            ni = calcNi(xi);
+            int sumn = calcSumn(xi);
+            n = sumn - xi.length - 1;
+            sumex = calcSex(xi, sumn);
+            double ve = calcVe(xi, sumn);
+            double b = calcB(xi, sumn);
+
+            double[] meanyi = calcMeanyi(xi);
+            double[] meanxi = calcMeanxi(xi);
+            double meanx = calcMeanx(xi);
+
+            TDistribution tDist = new TDistribution(n);
+            double t = tDist.inverseCumulativeProbability(1.0 - a / 2.0);
+            for(int i = 0; i < ret.length; i++) {
+                double wk = (meanxi[i] - meanx);
+                double wk2 = t * Math.sqrt((1/ni[i] + wk * wk / sumex) * ve);
+                double min = meanyi[i] - b * wk - wk2;
+                double max = meanyi[i] - b * wk + wk2;
+
+                ret[i] = new Interval(min, max);
+            }
+            
+            return ret;
+        }
+        private int[] calcNi(double[][][] xi) {
+            int[] ni = new int[xi.length];
+
+            for(int i = 0; i < xi.length; i++) {
+                ni[i] = xi[i].length;
+            }
+            return ni;
+        }
+        private double calcVe(double[][][] xi, int sumn) {
+            double sumey = calcSey(xi, sumn);
+            double sumeyx = calcSeyx(xi, sumn);
+
+            return (sumey * sumex - sumeyx * sumeyx) / (n * sumex);
+        }
+        private double calcB(double[][][] xi, int sumn) {
+            double sumeyx = calcSeyx(xi, sumn);
+            double sex = calcSex(xi, sumn);
+            
+            return sumeyx / sex;
+        }
+        private double[] calcMeanyi(double[][][] xi) {
+            double[] meanyi = new double[xi.length];
+
+            for(int i = 0; i < xi.length; i++) {
+                double sum = 0.0;
+                for(int j = 0; j < xi[i].length; j++) {
+                    sum += xi[i][j][0];
+                }
+                meanyi[i] = sum / xi[i].length;
+            }
+            return meanyi;
+        }
+        private double[] calcMeanxi(double[][][] xi) {
+            double[] meanxi = new double[xi.length];
+
+            for(int i = 0; i < xi.length; i++) {
+                double sum = 0.0;
+                for(int j = 0; j < xi[i].length; j++) {
+                    sum += xi[i][j][1];
+                }
+                meanxi[i] = sum / xi[i].length;
+            }
+            return meanxi;
+        }
+        private double calcMeanx(double[][][] xi) {
+            double sum = 0.0;
+            double n = 0;
+            for(int i = 0; i < xi.length; i++) {
+                for(int j = 0; j < xi[i].length; j++) {
+                    sum += xi[i][j][1];
+                    n++;
+                }
+            }
+            return sum / n;
+        }
+    }
 }
 

data/lib/num4anova.rb

@@ -314,6 +314,52 @@ module Num4AnovaLib
         def difference_test(xi, a)
             @ancova.differenceTest(xi.to_java(Java::double[][]), a)
         end
+        # 区間推定
+        #
+        # @overload interval_estim(xi, a)
+        #   @param [array]  xi データ(double[][][])
+        #   @param [double] a  有意水準
+        #   @return [Hash]  (min:下限信頼区間 max:上限信頼区間)
+        # @example
+        #   xi = [
+        #      [
+        #          [3,35], [5,38], [3,39],
+        #      ],
+        #      [
+        #          [3,36], [3,39], [8,54],
+        #      ],
+        #      [
+        #          [2,40], [2,45], [2,39],
+        #      ],
+        #      [
+        #          [3,47], [4,52], [2,48],
+        #      ],
+        #      [
+        #          [1,64], [2,80], [0,70],
+        #      ],
+        #    ]
+        #    ancova = Num4AnovaLib::Num4AncovaLib.new
+        #    ancova.interval_estim(xi, 0.05)
+        #    =>
+        #       {:min=>4.466605469341916,  :max=>7.1909253948556096}
+        #       {:min=>5.05699825110459,   :max=>6.386335082228742}
+        #       {:min=>2.510804295684195,  :max=>4.250430272217034}
+        #       {:min=>2.8089257316042135, :max=>2.9566298239513418}
+        #       {:min=>-6.303283147572267, :max=>-0.6577045067487104}
+        def interval_estim(xi, a)
+            retRb = []
+            retJava = @ancova.intervalEstim(xi.to_java(Java::double[][]), a)
+            sz = retJava.size
+            sz.times do |i|
+                retRb.push(
+                  {
+                    "min": retJava[i].getMin(),
+                    "max": retJava[i].getMax()
+                  }
+                )
+            end
+            return retRb
+        end
     end
 end
 

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4anova
 version: !ruby/object:Gem::Version
-  version: 0.0.11
+  version: 0.0.13
 platform: java
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-03-21 00:00:00.000000000 Z
+date: 2024-04-04 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake