num4normality 0.0.10-java → 0.0.11-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: bc5dd36e009ab6ce0b3a044652184f723ce70b416bab8780d01a56952bcd7b8b
-  data.tar.gz: 46fa3aefe151fb4e6dd8a47a3c4a4e38670704dd74b21a9d6fb19f6a13349a82
+  metadata.gz: e05aa008cb138a008afea68e92fd7df22fa5155f709ac2649e11db122bc28e44
+  data.tar.gz: 011e3a5473bf128c2946cc23126b59ff977b3b4d937b24405baa655a8d7b6b3f
 SHA512:
-  metadata.gz: f9933cc4c22ac645b17ae8e65b835d72fc80453e05c478903d1c748d3ff66615b9a0c631ced2ef4a2c8051297b2da16cdcd73b47d654032ed79599b0ae19f2bf
-  data.tar.gz: c74c186cc3caa4d511f8187ec9b7024d85ae8249d40692942800655c45fda60589be35745cb94659fd7d40038e65d4a0e3cb4ae91ef837d06302a6e180736e6d
+  metadata.gz: 8fdb7fe9d57d14a3bd0a7584cbd0f19803f651799ab4d3cbdf14d1ed5f057d480cbf274ca086037adfd28402a991e7b07cdb4d13ec29e6b02fde5f6985dd0906
+  data.tar.gz: b311f196404a1bd3e37b084d3dcca33100e36ff1e8ea5df3fc9a986e56e7a008693e9c01993a74bbe26ad9afbdc5b5e5d9110900b78c7be3d2335e3c25f486bb

data/CHANGELOG.md

@@ -2,6 +2,10 @@
 
 ## Unreleased
 
+## [0.0.11] - 2024-06-05
+### add
+- add fiducial interval in qqplot.
+
 ## [0.0.10] - 2024-05-18
 ### add
 - add function of jbtest.

data/ext/num4normality/Normality.java

@@ -149,6 +149,7 @@ public class Normality {
             return plotImpl.createPlot(dname, xi);
         }
         public static class QQPlot implements CreatePlot {
+            private double xbar = 0.0;
             private double[][] createData(double[] xi) {
                 DescriptiveStatistics stat = new DescriptiveStatistics();
                 NormalDistribution ndist = new NormalDistribution(0, 1);
@@ -159,15 +160,18 @@ public class Normality {
                 double sum = stat.getSum();
                 double[][] data = new double[n][2];
                 double p = 0.0;
+                double sumx = 0.0;
 
                 for (int i = 0; i < n; i++) {
                     p += xi[i] / sum;
                     double x = 
                         ndist.inverseCumulativeProbability(p * (i + 1.0) / (n + 1.0));
 
+                    sumx += x;
                     data[i][0] = x;
                     data[i][1] = xi[i];
                 }
+                xbar = sumx / n;
                 return data;
             }
             public XYPlot createPlot(String dname, double[] xi) {
@@ -224,16 +228,35 @@ public class Normality {
                 XYSeries cu     = new XYSeries("累積");
 
                 simpleReg.addData(data);
-                double a = simpleReg.getSlope();
-                double b = simpleReg.getIntercept();
-     
                 for (double x = ChartPlot.CLASS_MIN; x < ChartPlot.CLASS_MAX; x += 0.01) {
-                    double y = a * x + b;
-
-                    cu.add(x, y);
+                    cu.add(x, simpleReg.predict(x));
                 }
+                // 信頼区間
+                long n = simpleReg.getN();
+                double ve = simpleReg.getMeanSquareError();
+                double sxx = simpleReg.getXSumSquares();
+                TDistribution tDist = new TDistribution(n - 2);
+                double t = tDist.inverseCumulativeProbability(0.975);
+                XYSeries p05m     = new XYSeries("P05-");
+                XYSeries p05p     = new XYSeries("P05+");
+                
+                for(int i = 0; i < n; i++) {
+                    double x =  data[i][0];
+                    double residual = (x - xbar);
+                    double interval = t * Math.sqrt(
+                        (1 + 1.0 / n + residual * residual / sxx) * ve
+                    );
+                    double y01 = simpleReg.predict(x) - interval;
+                    double y02 = simpleReg.predict(x) + interval;
+
+                    p05m.add(x, y01);
+                    p05p.add(x, y02); 
+                }
+
                 XYSeriesCollection series = new XYSeriesCollection();
                 series.addSeries(cu);
+                series.addSeries(p05m);
+                series.addSeries(p05p);
                 return series;
             }
         }

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4normality
 version: !ruby/object:Gem::Version
-  version: 0.0.10
+  version: 0.0.11
 platform: java
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-05-18 00:00:00.000000000 Z
+date: 2024-06-05 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake