stat 0.0.1

data/LICENSE

@@ -0,0 +1,22 @@
+Copyright (c) 2012 arvicco
+
+MIT License
+
+Permission is hereby granted, free of charge, to any person obtaining
+a copy of this software and associated documentation files (the
+"Software"), to deal in the Software without restriction, including
+without limitation the rights to use, copy, modify, merge, publish,
+distribute, sublicense, and/or sell copies of the Software, and to
+permit persons to whom the Software is furnished to do so, subject to
+the following conditions:
+
+The above copyright notice and this permission notice shall be
+included in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

data/README.md

@@ -0,0 +1,29 @@
+# Stat
+
+Support for doing stats in Ruby with a little help from R and Octave.
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+    gem 'stat'
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install stat
+
+## Usage
+
+No documentation, this gem is pre-alpa as of now.
+
+## Contributing
+
+1. Fork it
+2. Create your feature branch (`git checkout -b my-new-feature`)
+3. Commit your changes (`git commit -am 'Added some feature'`)
+4. Push to the branch (`git push origin my-new-feature`)
+5. Create new Pull Request

data/Rakefile

@@ -0,0 +1,2 @@
+#!/usr/bin/env rake
+require "bundler/gem_tasks"

data/bin/preg/convert.rb

@@ -0,0 +1,11 @@
+#!/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

@@ -0,0 +1,32 @@
+#!/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

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

data/bin/stat.rb

@@ -0,0 +1,8 @@
+#!/usr/bin/env ruby
+# This loads bundler and sets loadpath so that you don't have to.
+
+require 'rubygems'
+require 'bundler/setup'
+$LOAD_PATH.unshift File.expand_path(File.dirname(__FILE__) + '/../lib')
+
+require 'stat'

data/lib/preg.rb

@@ -0,0 +1,54 @@
+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

data/lib/stat/record.rb

@@ -0,0 +1,18 @@
+require 'ostruct'
+
+module Stat
+
+  ### Record is a single row in a data table
+  class Record < OpenStruct
+
+    def initialize fields, line
+      super()
+      fields.each do |(name, start, finish, type)|
+        val = line[start-1...finish].strip.send("to_#{type}") rescue nil
+        send "#{name}=", val
+      end
+    end
+
+  end # class Record
+end # module Stat
+

data/lib/stat/table.rb

@@ -0,0 +1,49 @@
+require 'pathname'
+require 'stat/record'
+
+module Stat
+  ### Table represents a table data file (actual file format may differ)
+  class Table
+    attr_reader :records
+
+    def initialize
+      @records = []
+    end
+
+    def length
+      @records.length
+    end
+
+    # Use it to define fields in subclasses
+    def fields
+      []
+    end
+
+    # Reads in data file
+    def read filename, constructor, n=nil
+      file = Pathname.new(filename)
+      file.open('r') do |file|
+        file.readlines.each_with_index do |line, i|
+          break if i == n
+          @records << constructor.new(fields, line)
+        end
+      end
+      self
+    end
+
+    # Child classes can override this to recode values
+    def recode
+      self
+    end
+
+    # Turn Table into CSV format
+    def to_csv
+      fields.map(&:first).join(',') + "\n" +
+          records.map do |rec|
+            fields.map { |(name, _, _, _)| rec.send(name) || '' }.join(',')
+          end.join("\n")
+    end
+
+  end # class Table
+end # module Stat
+

data/lib/stat/version.rb

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

data/lib/stat.rb

@@ -0,0 +1,8 @@
+# This file contains code library for use with "Think Stats"
+
+module Stat
+end
+
+require "stat/version"
+require "stat/table"
+

metadata

@@ -0,0 +1,88 @@
+--- !ruby/object:Gem::Specification
+name: stat
+version: !ruby/object:Gem::Version
+  prerelease:
+  version: 0.0.1
+platform: ruby
+authors:
+- arvicco
+autorequire:
+bindir: bin
+cert_chain: []
+date: 2012-04-30 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: 1.1.3
+    none: false
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: 1.1.3
+    none: false
+  prerelease: false
+  type: :runtime
+description: Support for doing stats in Ruby
+email:
+- arvicco@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- !binary |-
+  YmluL3N0YXQucmI=
+- !binary |-
+  YmluL3ByZWcvY29udmVydC5yYg==
+- !binary |-
+  YmluL3ByZWcvZmlyc3QucmI=
+- !binary |-
+  YmluL3ByZWcvaGVscGVyLnJi
+- !binary |-
+  bGliL3ByZWcucmI=
+- !binary |-
+  bGliL3N0YXQucmI=
+- !binary |-
+  bGliL3N0YXQvcmVjb3JkLnJi
+- !binary |-
+  bGliL3N0YXQvdGFibGUucmI=
+- !binary |-
+  bGliL3N0YXQvdmVyc2lvbi5yYg==
+- !binary |-
+  UmFrZWZpbGU=
+- !binary |-
+  UkVBRE1FLm1k
+- !binary |-
+  TElDRU5TRQ==
+homepage: ''
+licenses: []
+post_install_message:
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: !binary |-
+        MA==
+  none: false
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: !binary |-
+        MA==
+  none: false
+requirements: []
+rubyforge_project:
+rubygems_version: 1.8.23
+signing_key:
+specification_version: 3
+summary: Support for doing stats in Ruby with a little help from R and Octave
+test_files: []
+...