stat 0.0.1 → 0.0.2

data/bin/octave_test

@@ -0,0 +1,18 @@
+#!/usr/bin/env octave -qf
+
+printf ("%s", program_name());
+
+for i = 1:nargin
+  printf (" %s", argv(){i});
+end
+
+printf ("\n");
+printf ("Hello, world!\n");
+
+r = rand(50) * 100;
+whos
+
+setenv("GNUTERM", "x11");
+hist(r);
+sleep(8);
+

data/bin/octave_test.rb

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require_relative 'stat'
+require 'octave'
+
+engine = Octave::Engine.new
+engine.eval "123.456 * 789.101112"
+engine.rand(10)
+
+matrix = Octave::Matrix.new(20, 400)
+20.times { |m| 400.times { |n| matrix[m, n] = rand } }
+engine.put_variable("m", matrix)
+
+engine.save "/tmp/20_x_400_matrix"

data/bin/r_test.rb

@@ -0,0 +1,68 @@
+#!/usr/bin/env ruby
+
+require_relative 'stat'
+
+p MY_DIR = `pwd`.strip
+p R >> 'getwd()'
+
+# converse: casts the result with 'to_ruby'
+p R.converse "mean(c(1,2,3))" # -> 2.0
+
+# convert: returns the raw Rserve::REXP result (like eval)
+rexp = R.convert "mean(c(1,2,3))" # -> #<Rserve::REXP::Double:0x00000002296498 @payload=[2.0], @attr=nil>
+rexp.as_doubles # -> [2.0]
+p rexp.to_ruby # -> 2.0
+
+# command: tell R what to do and only expect boolean reply
+R.command "z <- mean(c(1,2,3))" # -> true
+# all variables are persistent in the session and can be retrieved later
+p R >> "z" # -> 2.0
+
+p 'Correlation:'
+p R.>> "cor(a,b)", a: [1, 2, 3], b: [1, 2, 3]
+
+p hash = {
+    var1: [1, 2, 3, 4],
+    fac1: [3, 4, 5, 6],
+    res1: [4, 5, 6, 7]
+}
+
+## convert with hash.to_dataframe or Rserve::DataFrame.new(hash)
+
+R.command(dfr: hash.to_dataframe) do
+  %Q[
+    pdf("#{MY_DIR}/out/r.pdf")
+    hist(dfr$var1)
+    dev.off()
+    pdf("#{MY_DIR}/out/s.pdf")
+    plot(dfr)
+    dev.off()
+   ]
+end
+
+R.>> "pdf('#{MY_DIR}/out/t.pdf')",
+     "plot(c(1,2,3),c(4,5,6))",
+     "dev.off()"
+
+DataRow = Struct.new(:var1, :fac1, :res1)
+
+structs = [
+    DataRow.new(1, 3, 4),
+    DataRow.new(2, 4, 5),
+    DataRow.new(3, 5, 6),
+    DataRow.new(4, 6, 7)
+]
+
+datafr = Rserve::DataFrame.from_structs(structs)
+reply = R.converse('summary(df)', df: datafr).each_cons(6).to_a
+datafr.colnames.zip(reply) { |name, data| puts [name].+(data).join("\n ") }
+
+R >> 'plot(c(1,2,3), c(4,5,6))'
+R >> 'dev.hold'
+sleep(2)
+
+R >> 'r = rand(50) * 100'
+R >> 'hist(r)'
+sleep(6)
+
+#R.pause

data/bin/stat.rb

@@ -6,3 +6,4 @@ require 'bundler/setup'
 $LOAD_PATH.unshift File.expand_path(File.dirname(__FILE__) + '/../lib')
 
 require 'stat'
+require 'rserve/simpler/R'

data/lib/stat/version.rb

@@ -1,3 +1,3 @@
 module Stat
-  VERSION = "0.0.1"
+  VERSION = "0.0.2"
 end

metadata

@@ -2,14 +2,14 @@
 name: stat
 version: !ruby/object:Gem::Version
   prerelease:
-  version: 0.0.1
+  version: 0.0.2
 platform: ruby
 authors:
 - arvicco
 autorequire:
 bindir: bin
 cert_chain: []
-date: 2012-04-30 00:00:00.000000000 Z
+date: 2012-05-01 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -27,6 +27,22 @@ dependencies:
     none: false
   prerelease: false
   type: :runtime
+- !ruby/object:Gem::Dependency
+  name: rserve-simpler
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: 0.0.6
+    none: false
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: 0.0.6
+    none: false
+  prerelease: false
+  type: :runtime
 description: Support for doing stats in Ruby
 email:
 - arvicco@gmail.com
@@ -35,15 +51,13 @@ extensions: []
 extra_rdoc_files: []
 files:
 - !binary |-
-  YmluL3N0YXQucmI=
-- !binary |-
-  YmluL3ByZWcvY29udmVydC5yYg==
+  YmluL29jdGF2ZV90ZXN0
 - !binary |-
-  YmluL3ByZWcvZmlyc3QucmI=
+  YmluL29jdGF2ZV90ZXN0LnJi
 - !binary |-
-  YmluL3ByZWcvaGVscGVyLnJi
+  YmluL3JfdGVzdC5yYg==
 - !binary |-
-  bGliL3ByZWcucmI=
+  YmluL3N0YXQucmI=
 - !binary |-
   bGliL3N0YXQucmI=
 - !binary |-

data/bin/preg/convert.rb

@@ -1,11 +0,0 @@
-#!/usr/bin/env ruby
-
-require_relative 'helper'
-
-output_file = ARGV[0] || "#{DATA_DIR}/preg.csv"
-
-preg = Pregnancies.new.read #(100)
-
-File.open(output_file, 'w') do |output|
-  output.puts preg.to_csv
-end

data/bin/preg/first.rb

@@ -1,32 +0,0 @@
-#!/usr/bin/env ruby
-
-require_relative 'helper'
-
-resp = Respondents.new.read.recode #(10)
-puts 'Number of respondents', resp.length
-
-preg = Pregnancies.new.read.recode #(10)
-puts 'Number of pregnancies', preg.length
-
-#p preg.records[1]
-
-live_births = preg.records.find_all { |rec| rec.outcome == 1 }
-puts "Live births: #{live_births.size}"
-
-first_births = live_births.find_all { |rec| rec.birthord == 1 }
-first_avg_length = first_births.inject(0.0) { |sum, rec| sum + rec.prglength } / first_births.size
-puts "First births: #{first_births.size} total, #{first_avg_length} average length"
-
-other_births = live_births.find_all { |rec| rec.birthord != 1 }
-other_avg_length = other_births.inject(0.0) { |sum, rec| sum + rec.prglength } / other_births.size
-puts "Other births: #{other_births.size} total, #{other_avg_length} average length"
-
-puts 'Difference first-others (weeks)'
-p (first_avg_length - other_avg_length)
-
-puts 'Difference (hours)'
-p (first_avg_length - other_avg_length) * 7 * 24
-
-puts 'Difference %'
-
-p (first_avg_length - other_avg_length) / other_avg_length * 100

data/bin/preg/helper.rb

@@ -1,6 +0,0 @@
-require_relative '../stat'
-require 'preg'
-include Preg
-
-DATA_DIR = 'data/preg'
-

data/lib/preg.rb

@@ -1,54 +0,0 @@
-require 'stat'
-
-# This module contains code for use with "Think Stats" pregnancy examples
-module Preg
-
-  class Respondent < Stat::Record
-  end
-
-  class Pregnancy < Stat::Record
-  end
-
-  class Respondents < Stat::Table
-    def read n=nil
-      super 'data/preg/2002FemResp.dat', Respondent, n
-    end
-
-    def fields
-      [['caseid', 1, 12, :i]
-      ]
-    end
-  end
-
-  class Pregnancies < Stat::Table
-    def read n=nil
-      super 'data/preg/2002FemPreg.dat', Pregnancy, n
-    end
-
-    def fields
-      [['caseid', 1, 12, :i],
-       ['nbrnaliv', 22, 22, :i],
-       ['babysex', 56, 56, :i],
-       ['birthwgt_lb', 57, 58, :i],
-       ['birthwgt_oz', 59, 60, :i],
-       ['prglength', 275, 276, :i],
-       ['outcome', 277, 277, :i],
-       ['birthord', 278, 279, :i],
-       ['agepreg', 284, 287, :i],
-       ['finalwgt', 423, 440, :f],
-      ]
-    end
-
-    def recode
-      @records.each do |rec|
-        rec.agepreg /= 100.0 if rec.agepreg
-        if rec.birthwgt_lb && rec.birthwgt_oz &&
-            rec.birthwgt_lb < 20 && rec.birthwgt_oz <= 16
-          rec.totalwgt_oz = rec.birthwgt_lb * 16 + rec.birthwgt_oz
-        end
-      end
-
-      super
-    end
-  end # class Pregnancies
-end # module Preg