cheap-stats 0.1.0 → 0.1.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 0e5a30b7e09b21f08359019f9a4a33bb8bc8b4e1e103bea1eed739f12de818f0
-  data.tar.gz: dc3146cd65e34779373ffe61b5f32c90ac1cfb89218e5fefd81804bf6d40ac59
+  metadata.gz: f5a99f35a58f32ff2cfb41964a1ec75a82cf249200ad6b2b2392a9ffd2571976
+  data.tar.gz: 6b262daad757204e7d6842de7a43b682f7143d22697b5526d58126e3e9c036c9
 SHA512:
-  metadata.gz: 5860b61e54bc57a8ade4724d4f5e5ac62811079b75869177e774c0eb281341fa0a2c7ee7f41f59fbb409355b4a2e5757ccf14a7179d6912d0919d240928759fe
-  data.tar.gz: a9b091cfcb393b1236dd5eca02eaf2d5742c51ce0ea94da38827171c6752c0a6140856c50679b97d363eceabffbdd2876ffb5795aa99c0e365ffe81ad7f85ffd
+  metadata.gz: 54ca8675dcef936be225e9ecc7747aa9e709796fc45b1449b1df81162895427b50349f5497ac1d575f01d13984ac6bfbe89a84045228c0a6222f652fa6eaa1ef
+  data.tar.gz: 2488dc611af4b88933755af4b583a712401f8be45d93e7b304249a997491d2acc96af46412ee355985a207245aafb6fc47cb20a49d7e96f3ad35251cb878ab8f

data/ext/cheap_stats/cheap_stats.c

@@ -603,3 +603,37 @@ cheap_stats_pearson_skewness(cheap_stats_t* ptr, double* dst)
   return ret;
 }
 
+int
+cheap_stats_z_score(cheap_stats_t* ptr, double v, double* dst)
+{
+  int ret;
+
+  /*
+   * initialize
+   */
+  ret = 0;
+
+  /*
+   * argument check
+   */
+  do {
+    if (ptr == NULL) {
+      ret = DEFAULT_ERROR;
+      break;
+    }
+
+    if (dst == NULL) {
+      ret = DEFAULT_ERROR;
+      break;
+    }
+  } while (0);
+
+  /*
+   * do test (by Smirnov-Grubbs test)
+   */
+  if (!ret) {
+    *dst = (v - ptr->mean) / ptr->std;
+  }
+
+  return ret;
+}

data/ext/cheap_stats/cheap_stats.h

@@ -33,5 +33,6 @@ int cheap_stats_central_moment(cheap_stats_t* obj, double k, double* dst);
 int cheap_stats_std_moment(cheap_stats_t* obj, double k, double* dst);
 int cheap_stats_skewness(cheap_stats_t* obj, double* dst);
 int cheap_stats_pearson_skewness(cheap_stats_t* obj, double* dst);
+int cheap_stats_z_score(cheap_stats_t* obj, double v, double* res);
 
 #endif /* !defined(__SMALL_STATS_H__) */

data/ext/cheap_stats/rb_cheap_stats.c

@@ -484,6 +484,37 @@ rb_cheap_stats_pearson_skewness(VALUE self)
   return DBL2NUM(ret);
 }
 
+/**
+ * calc Z-score
+ *
+ * @param [Numeric] v   target value
+ *
+ * @return [Float] skewness value
+ */
+static VALUE
+rb_cheap_stats_z_score(VALUE self, VALUE v)
+{
+  rb_cheap_stats_t* ptr;
+  int err;
+  double ret;
+
+  /*
+   * strip context data
+   */
+  TypedData_Get_Struct(self, rb_cheap_stats_t, &rb_cheap_stats_data_type, ptr);
+
+  /*
+   * check argument
+   */
+  err = cheap_stats_z_score(ptr->stats, rb_num2dbl(v), &ret);
+  if (err) {
+    RUNTIME_ERROR("cheap_stats_z_score() failed [err=%d]", err);
+  }
+
+  return DBL2NUM(ret);
+}
+
+
 void
 Init_cheap_stats()
 {
@@ -507,6 +538,7 @@ Init_cheap_stats()
   rb_define_method(klass, "std_moment", rb_cheap_stats_std_moment, 1);
   rb_define_method(klass, "skewness", rb_cheap_stats_skewness, 0);
   rb_define_method(klass, "pearson_skewness",rb_cheap_stats_pearson_skewness,0);
+  rb_define_method(klass, "z_score",rb_cheap_stats_z_score, 1);
 
   rb_alias(klass, rb_intern("average"), rb_intern("mean"));
   rb_alias(klass, rb_intern("sigma"), rb_intern("std"));

data/lib/cheap_stats/version.rb

@@ -1,3 +1,3 @@
 class CheapStats
-  VERSION = "0.1.0"
+  VERSION = "0.1.1"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: cheap-stats
 version: !ruby/object:Gem::Version
-  version: 0.1.0
+  version: 0.1.1
 platform: ruby
 authors:
 - Hiroshi Kuwagata
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2019-04-17 00:00:00.000000000 Z
+date: 2019-04-24 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler