health_stats 0.0.1

data/LICENSE

@@ -0,0 +1,20 @@
+Copyright (c) 2009 Mark Martin
+
+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/Rakefile

@@ -0,0 +1,20 @@
+require 'rubygems'
+require 'rake'
+
+begin
+  require 'jeweler'
+  Jeweler::Tasks.new do |gemspec|
+    gemspec.name = "health_stats"
+    gemspec.summary = "Convenience methods for age, bmi, bmi percentile"
+    gemspec.email = "assplecake@gmail.com"
+    gemspec.homepage = "http://github.com/assplecake/healthdata"
+    gemspec.authors = ["Mark Martin"]
+    
+    gemspec.add_development_dependency('activesupport', '>= 2.3')
+  end
+  
+  Jeweler::GemcutterTasks.new
+  
+rescue LoadError
+  puts "Jeweler not available. Install it with: sudo gem install jeweler"
+end

data/VERSION

@@ -0,0 +1 @@
+0.0.1

data/health_stats.gemspec

@@ -0,0 +1,55 @@
+# Generated by jeweler
+# DO NOT EDIT THIS FILE DIRECTLY
+# Instead, edit Jeweler::Tasks in Rakefile, and run the gemspec command
+# -*- encoding: utf-8 -*-
+
+Gem::Specification.new do |s|
+  s.name = %q{health_stats}
+  s.version = "0.0.1"
+
+  s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
+  s.authors = ["Mark Martin"]
+  s.date = %q{2009-12-03}
+  s.email = %q{assplecake@gmail.com}
+  s.extra_rdoc_files = [
+    "LICENSE",
+     "README.md"
+  ]
+  s.files = [
+    "LICENSE",
+     "README.md",
+     "Rakefile",
+     "VERSION",
+     "health_stats.gemspec",
+     "lib/health_stats.rb",
+     "lib/health_stats/age.rb",
+     "lib/health_stats/bmi.rb",
+     "lib/health_stats/cdc_data.rb",
+     "lib/health_stats/cdc_data/bmi.rb",
+     "lib/health_stats/statistics.rb",
+     "pkg/health_stats-0.0.1.gem",
+     "spec/spec_health_stats.rb"
+  ]
+  s.homepage = %q{http://github.com/assplecake/healthdata}
+  s.rdoc_options = ["--charset=UTF-8"]
+  s.require_paths = ["lib"]
+  s.rubygems_version = %q{1.3.5}
+  s.summary = %q{Convenience methods for age, bmi, bmi percentile}
+  s.test_files = [
+    "spec/spec_health_stats.rb"
+  ]
+
+  if s.respond_to? :specification_version then
+    current_version = Gem::Specification::CURRENT_SPECIFICATION_VERSION
+    s.specification_version = 3
+
+    if Gem::Version.new(Gem::RubyGemsVersion) >= Gem::Version.new('1.2.0') then
+      s.add_development_dependency(%q<activesupport>, [">= 2.3"])
+    else
+      s.add_dependency(%q<activesupport>, [">= 2.3"])
+    end
+  else
+    s.add_dependency(%q<activesupport>, [">= 2.3"])
+  end
+end
+

data/lib/health_stats.rb

@@ -0,0 +1,15 @@
+module HealthStats
+  AttributeError = Class.new(StandardError)
+  RequiredAttributes = [:dob, :weight, :height, :gender]
+  
+  RequiredAttributes.each do |attribute|
+    define_method(attribute) do
+      message = "#{attribute} is not implemented"
+      raise HealthStats::AttributeError, message
+    end
+  end
+end
+
+require 'health_stats/age'
+require 'health_stats/statistics'
+require 'health_stats/bmi'

data/lib/health_stats/age.rb

@@ -0,0 +1,35 @@
+module HealthStats
+  module Age
+    def age_in_months
+      return if dob.nil?
+      today = Date.today
+      months = age_in_years * 12
+
+      if today.month > dob.month
+        months += (today.month - dob.month)
+      else
+        months += today.month
+      end
+
+      if today.day >= dob.day
+        months += 1 + (today.day - dob.day > 15 ? 0.5 : 0.0)
+      else
+        months += (today.day - dob.day + 30) > 15 ? -0.5 : -1.0
+      end
+    end
+
+    def age_in_years
+      return if dob.nil?
+      today = Date.today
+      
+      if today.month > dob.month || (today.month == dob.month && today.day >= dob.day)
+        today.year - dob.year
+      else
+        today.year - dob.year - 1
+      end
+    end
+    alias_method :age, :age_in_years
+  end
+  
+  include Age
+end

data/lib/health_stats/bmi.rb

@@ -0,0 +1,34 @@
+module HealthStats
+  autoload :CDCData, File.dirname(__FILE__) + '/cdc_data'
+  
+  module BMI
+    include Statistics
+    
+    def bmi    
+      if weight && height
+        (((weight * 703) * 100) / (height ** 2)) / 100.0
+      end
+    end
+
+    def bmi_percentile
+      bmi_score = bmi
+      month_age = age_in_months
+      lms_value = nil
+      
+      return unless bmi_score && month_age && gender
+      return if month_age < 24.0 || month_age > 241.0
+      
+      chart = gender == 'f' ? CDCData::BMI::Female : CDCData::BMI::Male
+      
+      while lms_value.nil?
+        lms_value = chart[month_age]
+        month_age -= 0.5
+      end
+      
+      # I feel dirty.
+      ((poz(lms(bmi_score, *lms_value)) * 1000000).to_i / 100.0).round / 100.0 
+    end
+  end
+  
+  include BMI
+end

data/lib/health_stats/cdc_data.rb

@@ -0,0 +1,8 @@
+module HealthStats
+  module CDCData
+    module BMI
+      autoload :Male, File.dirname(__FILE__) + '/cdc_data/bmi'
+      autoload :Female, File.dirname(__FILE__) + '/cdc_data/bmi'
+    end
+  end
+end

data/lib/health_stats/cdc_data/bmi.rb

@@ -0,0 +1,449 @@
+module HealthStats
+  module CDCData
+    module BMI
+      Male = {
+        24.0 => [-2.01118107, 16.575027675, 0.080592465],
+        24.5 => [-1.982373595, 16.547774867, 0.0801274288],
+        25.5 => [-1.924100169, 16.494427632, 0.0792339937],
+        26.5 => [-1.86549793, 16.442595522, 0.0783893561],
+        27.5 => [-1.807261899, 16.392243398, 0.0775935012],
+        28.5 => [-1.750118905, 16.343336543, 0.0768464622],
+        29.5 => [-1.69481584, 16.295840971, 0.0761483079],
+        30.5 => [-1.642106779, 16.249723714, 0.0754991255],
+        31.5 => [-1.592744414, 16.204952678, 0.0748989935],
+        32.5 => [-1.547442391, 16.161498714, 0.0743479966],
+        33.5 => [-1.506902601, 16.119332582, 0.0738461386],
+        34.5 => [-1.471770047, 16.078427579, 0.0733933698],
+        35.5 => [-1.442628957, 16.038758958, 0.0729895508],
+        36.5 => [-1.419991255, 16.000304012, 0.0726344323],
+        37.5 => [-1.404277619, 15.963042771, 0.072327649],
+        38.5 => [-1.39586317, 15.926954175, 0.0720686401],
+        39.5 => [-1.394935252, 15.892025816, 0.0718568052],
+        40.5 => [-1.401671596, 15.858240929, 0.071691278],
+        41.5 => [-1.416100312, 15.825588223, 0.0715710933],
+        42.5 => [-1.438164899, 15.794057282, 0.0714951131],
+        43.5 => [-1.467669032, 15.763642549, 0.0714621062],
+        44.5 => [-1.504376347, 15.734336684, 0.0714706462],
+        45.5 => [-1.547942838, 15.706135657, 0.0715192177],
+        46.5 => [-1.597896397, 15.679040623, 0.0716062769],
+        47.5 => [-1.653732283, 15.653051916, 0.071730167],
+        48.5 => [-1.714869347, 15.628172692, 0.0718892136],
+        49.5 => [-1.780673181, 15.604407997, 0.0720817373],
+        50.5 => [-1.850468473, 15.58176458, 0.0723060813],
+        51.5 => [-1.923551865, 15.560250666, 0.0725606369],
+        52.5 => [-1.999220429, 15.539874597, 0.0728438397],
+        53.5 => [-2.076707178, 15.520649927, 0.0731543235],
+        54.5 => [-2.155348017, 15.502584267, 0.0734906668],
+        55.5 => [-2.234438552, 15.485689732, 0.0738516716],
+        56.5 => [-2.313321723, 15.469977177, 0.0742362346],
+        57.5 => [-2.391381273, 15.455456916, 0.0746433738],
+        58.5 => [-2.468032491, 15.442139608, 0.0750722636],
+        59.5 => [-2.542781541, 15.430032072, 0.0755221037],
+        60.5 => [-2.61516595, 15.419141631, 0.0759922501],
+        61.5 => [-2.684789516, 15.409473561, 0.0764821284],
+        62.5 => [-2.751316949, 15.401031388, 0.0769912324],
+        63.5 => [-2.81445945, 15.393817852, 0.0775191487],
+        64.5 => [-2.87402476, 15.387830936, 0.0780653898],
+        65.5 => [-2.92984048, 15.383069448, 0.0786295919],
+        66.5 => [-2.981796828, 15.37952958, 0.0792113694],
+        67.5 => [-3.029831343, 15.37720582, 0.0798103341],
+        68.5 => [-3.073924224, 15.376091068, 0.0804260861],
+        69.5 => [-3.114093476, 15.376176765, 0.0810582059],
+        70.5 => [-3.15039004, 15.37745304, 0.0817062489],
+        71.5 => [-3.182893018, 15.379908859, 0.0823697413],
+        72.5 => [-3.21170511, 15.383532169, 0.0830481778],
+        73.5 => [-3.23694834, 15.388310046, 0.0837410207],
+        74.5 => [-3.25876011, 15.39422883, 0.0844476998],
+        75.5 => [-3.277281546, 15.40127496, 0.0851676514],
+        76.5 => [-3.292683774, 15.40943252, 0.0859001836],
+        77.5 => [-3.305124073, 15.418686911, 0.0866446671],
+        78.5 => [-3.314768951, 15.429022732, 0.0874004212],
+        79.5 => [-3.321785992, 15.440424387, 0.0881667444],
+        80.5 => [-3.326345795, 15.452875806, 0.088942897],
+        81.5 => [-3.328602731, 15.466362178, 0.0897282019],
+        82.5 => [-3.328725277, 15.480867041, 0.0905218748],
+        83.5 => [-3.32687018, 15.496374654, 0.0913231621],
+        84.5 => [-3.323188896, 15.512869363, 0.0921313054],
+        85.5 => [-3.317827016, 15.530335632, 0.0929455443],
+        86.5 => [-3.310923871, 15.548758065, 0.0937651184],
+        87.5 => [-3.302612272, 15.568121426, 0.0945892702],
+        88.5 => [-3.293018361, 15.588410651, 0.0954172465],
+        89.5 => [-3.282260813, 15.60961101, 0.0962483005],
+        90.5 => [-3.270454609, 15.63170735, 0.0970816941],
+        91.5 => [-3.257703616, 15.654685628, 0.0979166976],
+        92.5 => [-3.244108214, 15.678531387, 0.0987525931],
+        93.5 => [-3.229761713, 15.703230518, 0.0995886749],
+        94.5 => [-3.214751287, 15.728769113, 0.1004242507],
+        95.5 => [-3.199158184, 15.755133465, 0.1012586429],
+        96.5 => [-3.18305795, 15.782310065, 0.1020911894],
+        97.5 => [-3.166520664, 15.810285603, 0.1029212448],
+        98.5 => [-3.1496103, 15.839047084, 0.1037481885],
+        99.5 => [-3.132389637, 15.868581229, 0.1045713862],
+        100.5 => [-3.114911153, 15.898875618, 0.1053902685],
+        101.5 => [-3.097226399, 15.929917651, 0.1062042575],
+        102.5 => [-3.079383079, 15.961694805, 0.1070127883],
+        103.5 => [-3.061423765, 15.994194894, 0.1078153274],
+        104.5 => [-3.043386071, 16.027406071, 0.1086113736],
+        105.5 => [-3.025310003, 16.061315897, 0.1094003876],
+        106.5 => [-3.007225737, 16.095912922, 0.1101819146],
+        107.5 => [-2.989164598, 16.131185315, 0.1109554781],
+        108.5 => [-2.971148225, 16.167122344, 0.1117206908],
+        109.5 => [-2.953208047, 16.203711677, 0.1124770587],
+        110.5 => [-2.935363951, 16.240942388, 0.1132241995],
+        111.5 => [-2.917635157, 16.278803458, 0.1139617339],
+        112.5 => [-2.900039803, 16.317283847, 0.1146892914],
+        113.5 => [-2.882593796, 16.356372672, 0.1154065227],
+        114.5 => [-2.865311266, 16.396059161, 0.1161130971],
+        115.5 => [-2.848204697, 16.436332645, 0.1168087018],
+        116.5 => [-2.831285052, 16.477182556, 0.1174930418],
+        117.5 => [-2.81456189, 16.518598425, 0.1181658396],
+        118.5 => [-2.79804347, 16.560569873, 0.1188268351],
+        119.5 => [-2.781736856, 16.603086612, 0.1194757852],
+        120.5 => [-2.765648008, 16.646138439, 0.1201124636],
+        121.5 => [-2.749782197, 16.689715178, 0.1207366562],
+        122.5 => [-2.734142443, 16.733806953, 0.1213481813],
+        123.5 => [-2.718732873, 16.778403632, 0.121946849],
+        124.5 => [-2.703555506, 16.82349538, 0.1225325012],
+        125.5 => [-2.688611957, 16.869072375, 0.1231049908],
+        126.5 => [-2.673903164, 16.915124866, 0.1236641859],
+        127.5 => [-2.659429443, 16.961643168, 0.1242099694],
+        128.5 => [-2.645190534, 17.00861766, 0.1247422387],
+        129.5 => [-2.631185649, 17.056038787, 0.1252609054],
+        130.5 => [-2.617413511, 17.103897052, 0.125765895],
+        131.5 => [-2.603872392, 17.152183022, 0.1262571467],
+        132.5 => [-2.590560148, 17.200887318, 0.1267346133],
+        133.5 => [-2.577474253, 17.250000623, 0.1271982604],
+        134.5 => [-2.564611831, 17.299513673, 0.1276480666],
+        135.5 => [-2.551969684, 17.349417258, 0.128084023],
+        136.5 => [-2.539539972, 17.399703081, 0.1285061919],
+        137.5 => [-2.527325681, 17.450360715, 0.1289144974],
+        138.5 => [-2.515320235, 17.501381606, 0.1293090012],
+        139.5 => [-2.503519447, 17.552756739, 0.1296897408],
+        140.5 => [-2.491918934, 17.604477144, 0.1300567649],
+        141.5 => [-2.480514136, 17.656533895, 0.1304101325],
+        142.5 => [-2.469300331, 17.708918107, 0.1307499132],
+        143.5 => [-2.458272656, 17.761620938, 0.1310761867],
+        144.5 => [-2.447426113, 17.814633586, 0.1313890423],
+        145.5 => [-2.436755595, 17.867947289, 0.1316885791],
+        146.5 => [-2.426255887, 17.92155332, 0.1319749052],
+        147.5 => [-2.415921689, 17.975442992, 0.1322481377],
+        148.5 => [-2.405747619, 18.029607653, 0.1325084026],
+        149.5 => [-2.395728233, 18.084038684, 0.1327558343],
+        150.5 => [-2.385858029, 18.138727501, 0.1329905752],
+        151.5 => [-2.376131459, 18.193665552, 0.133212776],
+        152.5 => [-2.366542942, 18.248844314, 0.1334225948],
+        153.5 => [-2.357086871, 18.304255296, 0.1336201973],
+        154.5 => [-2.347757625, 18.359890034, 0.1338057563],
+        155.5 => [-2.338549576, 18.415740092, 0.1339794518],
+        156.5 => [-2.3294571, 18.471797059, 0.1341414703],
+        157.5 => [-2.320474586, 18.528052549, 0.1342920051],
+        158.5 => [-2.311596446, 18.584498199, 0.1344312555],
+        159.5 => [-2.302817124, 18.641125665, 0.134559427],
+        160.5 => [-2.294131107, 18.697926627, 0.1346767311],
+        161.5 => [-2.285532933, 18.754892781, 0.1347833849],
+        162.5 => [-2.277017201, 18.812015839, 0.1348796107],
+        163.5 => [-2.268578584, 18.869287528, 0.1349656365],
+        164.5 => [-2.260211837, 18.92669959, 0.1350416951],
+        165.5 => [-2.251911809, 18.984243775, 0.1351080243],
+        166.5 => [-2.243673453, 19.041911845, 0.1351648666],
+        167.5 => [-2.235491842, 19.099695568, 0.1352124688],
+        168.5 => [-2.227362173, 19.157586716, 0.1352510825],
+        169.5 => [-2.21927979, 19.215577065, 0.1352809633],
+        170.5 => [-2.211240187, 19.27365839, 0.1353023707],
+        171.5 => [-2.203239029, 19.331822466, 0.1353155684],
+        172.5 => [-2.195272161, 19.390061061, 0.1353208237],
+        173.5 => [-2.187335625, 19.448365938, 0.1353184074],
+        174.5 => [-2.179425674, 19.506728848, 0.1353085942],
+        175.5 => [-2.171538789, 19.565141532, 0.1352916617],
+        176.5 => [-2.163671689, 19.623595714, 0.1352678911],
+        177.5 => [-2.155821357, 19.682083101, 0.1352375667],
+        178.5 => [-2.147985046, 19.740595376, 0.1352009759],
+        179.5 => [-2.140160305, 19.799124201, 0.1351584088],
+        180.5 => [-2.132344989, 19.857661209, 0.1351101588],
+        181.5 => [-2.124537282, 19.916198004, 0.1350565219],
+        182.5 => [-2.116735712, 19.974726154, 0.1349977968],
+        183.5 => [-2.108939167, 20.033237194, 0.134934285],
+        184.5 => [-2.10114692, 20.091722615, 0.1348662908],
+        185.5 => [-2.093358637, 20.15017387, 0.1347941208],
+        186.5 => [-2.085574403, 20.208582361, 0.1347180845],
+        187.5 => [-2.077794735, 20.266939444, 0.1346384938],
+        188.5 => [-2.070020599, 20.325236424, 0.1345556632],
+        189.5 => [-2.062253431, 20.383464548, 0.1344699098],
+        190.5 => [-2.054495145, 20.441615008, 0.1343815533],
+        191.5 => [-2.046748156, 20.499678935, 0.1342909159],
+        192.5 => [-2.039015385, 20.557647399, 0.1341983225],
+        193.5 => [-2.031300282, 20.615511404, 0.1341041006],
+        194.5 => [-2.023606828, 20.673261889, 0.1340085806],
+        195.5 => [-2.015942013, 20.730889051, 0.1339120657],
+        196.5 => [-2.008305745, 20.788385102, 0.1338149542],
+        197.5 => [-2.000706389, 20.845740029, 0.1337175518],
+        198.5 => [-1.993150137, 20.902944494, 0.1336202002],
+        199.5 => [-1.985643741, 20.959989088, 0.1335232441],
+        200.5 => [-1.97819451, 21.01686433, 0.1334270316],
+        201.5 => [-1.970810308, 21.073560674, 0.1333319137],
+        202.5 => [-1.96349954, 21.130068501, 0.1332382449],
+        203.5 => [-1.956271141, 21.186378131, 0.1331463832],
+        204.5 => [-1.949134561, 21.242479819, 0.1330566901],
+        205.5 => [-1.942099744, 21.298363759, 0.1329695305],
+        206.5 => [-1.935177101, 21.354020094, 0.1328852735],
+        207.5 => [-1.92837748, 21.409438911, 0.1328042916],
+        208.5 => [-1.921712136, 21.464610257, 0.1327269615],
+        209.5 => [-1.915192685, 21.519524136, 0.1326536641],
+        210.5 => [-1.908831065, 21.574170525, 0.1325847841],
+        211.5 => [-1.902639482, 21.628539373, 0.1325207109],
+        212.5 => [-1.896630358, 21.682620618, 0.1324618378],
+        213.5 => [-1.890816268, 21.736404191, 0.1324085629],
+        214.5 => [-1.885209876, 21.789880029, 0.1323612888],
+        215.5 => [-1.879823505, 21.843038191, 0.1323204265],
+        216.5 => [-1.874670324, 21.895868501, 0.1322863818],
+        217.5 => [-1.869760299, 21.948361684, 0.1322595999],
+        218.5 => [-1.865113245, 22.00050569, 0.1322404176],
+        219.5 => [-1.860734944, 22.052292423, 0.1322293301],
+        220.5 => [-1.85663384, 22.103713048, 0.132226801],
+        221.5 => [-1.852827186, 22.154756029, 0.1322332005],
+        222.5 => [-1.849323204, 22.205412485, 0.1322489931],
+        223.5 => [-1.846131607, 22.255673, 0.1322746254],
+        224.5 => [-1.843261294, 22.305528306, 0.132310549],
+        225.5 => [-1.840720248, 22.354969299, 0.1323572208],
+        226.5 => [-1.83851544, 22.403987057, 0.132415103],
+        227.5 => [-1.83665586, 22.452571818, 0.132484631],
+        228.5 => [-1.835138046, 22.500717781, 0.1325663592],
+        229.5 => [-1.833972004, 22.548414372, 0.132660699],
+        230.5 => [-1.833157751, 22.595654215, 0.1327681527],
+        231.5 => [-1.83269562, 22.642429557, 0.1328892105],
+        232.5 => [-1.832584342, 22.688732921, 0.1330243684],
+        233.5 => [-1.832820974, 22.734557126, 0.1331741285],
+        234.5 => [-1.833400825, 22.779895295, 0.1333389994],
+        235.5 => [-1.834317405, 22.824740868, 0.1335194959],
+        236.5 => [-1.83555752, 22.869089116, 0.1337161923],
+        237.5 => [-1.837119466, 22.912931508, 0.1339295249],
+        238.5 => [-1.838987063, 22.956263733, 0.1341600729],
+        239.5 => [-1.841146139, 22.999080616, 0.1344083809],
+        240.0 => [-1.84233016, 23.020294238, 0.134539365],
+        240.5 => [-1.843580575, 23.041377338, 0.1346750014]
+      }
+
+      Female = {
+      	24.0 => [-0.98660853, 16.423396643, 0.085451785],
+        24.5 => [-1.024496827, 16.388040561, 0.085025838],
+        25.5 => [-1.102698353, 16.318971901, 0.0842140522],
+        26.5 => [-1.18396635, 16.252079845, 0.083455124],
+        27.5 => [-1.268071036, 16.187346686, 0.0827482842],
+        28.5 => [-1.354751525, 16.124754481, 0.0820927371],
+        29.5 => [-1.443689692, 16.064287623, 0.0814877172],
+        30.5 => [-1.53454192, 16.005930007, 0.0809324482],
+        31.5 => [-1.626928093, 15.94966631, 0.0804261754],
+        32.5 => [-1.720434829, 15.895481969, 0.0799681758],
+        33.5 => [-1.814635262, 15.843361791, 0.0795577348],
+        34.5 => [-1.909076262, 15.793291456, 0.0791941867],
+        35.5 => [-2.003296102, 15.7452564, 0.0788768946],
+        36.5 => [-2.096828937, 15.699241878, 0.0786052551],
+        37.5 => [-2.189211877, 15.655232823, 0.0783786964],
+        38.5 => [-2.279991982, 15.613213709, 0.0781966743],
+        39.5 => [-2.368732949, 15.573168427, 0.078058667],
+        40.5 => [-2.455021314, 15.53508019, 0.077964169],
+        41.5 => [-2.538471972, 15.498931449, 0.0779126837],
+        42.5 => [-2.618732901, 15.464703844, 0.0779037156],
+        43.5 => [-2.695488973, 15.432378168, 0.0779367628],
+        44.5 => [-2.768464816, 15.401934364, 0.078011309],
+        45.5 => [-2.837426693, 15.373351541, 0.0781268172],
+        46.5 => [-2.902178205, 15.346608415, 0.0782827393],
+        47.5 => [-2.962580386, 15.321681814, 0.0784784485],
+        48.5 => [-3.018521987, 15.298548972, 0.0787133246],
+        49.5 => [-3.069936555, 15.277186179, 0.0789866938],
+        50.5 => [-3.116795864, 15.257569204, 0.0792978405],
+        51.5 => [-3.159107331, 15.239673384, 0.079646006],
+        52.5 => [-3.196911083, 15.22347371, 0.0800303887],
+        53.5 => [-3.230276759, 15.208944907, 0.0804501449],
+        54.5 => [-3.259300182, 15.19606152, 0.0809043905],
+        55.5 => [-3.284099963, 15.184797987, 0.0813922027],
+        56.5 => [-3.30481415, 15.175128708, 0.0819126232],
+        57.5 => [-3.321596954, 15.167028107, 0.0824646608],
+        58.5 => [-3.334615646, 15.160470684, 0.0830472946],
+        59.5 => [-3.344047622, 15.155431067, 0.0836594775],
+        60.5 => [-3.35007771, 15.15188405, 0.0843001394],
+        61.5 => [-3.352893805, 15.149804788, 0.0849681996],
+        62.5 => [-3.352691376, 15.14916825, 0.085662539],
+        63.5 => [-3.34966438, 15.149949835, 0.086382035],
+        64.5 => [-3.343998803, 15.152125852, 0.0871255909],
+        65.5 => [-3.335889574, 15.155671862, 0.0878920466],
+        66.5 => [-3.325522491, 15.160564192, 0.0886802643],
+        67.5 => [-3.31307846, 15.166779473, 0.0894891056],
+        68.5 => [-3.298732648, 15.174294641, 0.090317434],
+        69.5 => [-3.282653831, 15.183086936, 0.0911641168],
+        70.5 => [-3.265003896, 15.193133896, 0.0920280276],
+        71.5 => [-3.245937506, 15.204413348, 0.0929080476],
+        72.5 => [-3.225606516, 15.216902957, 0.0938030328],
+        73.5 => [-3.204146115, 15.230581504, 0.0947119161],
+        74.5 => [-3.181690237, 15.245427448, 0.0956335947],
+        75.5 => [-3.158363475, 15.261419664, 0.096566992],
+        76.5 => [-3.134282833, 15.278537278, 0.0975110459],
+        77.5 => [-3.109557879, 15.296759667, 0.0984647101],
+        78.5 => [-3.084290931, 15.316066442, 0.0994269552],
+        79.5 => [-3.058577292, 15.336437447, 0.1003967693],
+        80.5 => [-3.032505499, 15.357852744, 0.1013731591],
+        81.5 => [-3.0061576, 15.380292613, 0.1023551503],
+        82.5 => [-2.979609448, 15.403737535, 0.1033417884],
+        83.5 => [-2.952930993, 15.428168191, 0.1043321392],
+        84.5 => [-2.926186592, 15.453565452, 0.1053252892],
+        85.5 => [-2.899435307, 15.479910374, 0.1063203463],
+        86.5 => [-2.872731211, 15.507184187, 0.1073164399],
+        87.5 => [-2.846123683, 15.535368293, 0.1083127212],
+        88.5 => [-2.819657704, 15.564444257, 0.1093083637],
+        89.5 => [-2.793374145, 15.594393802, 0.1103025629],
+        90.5 => [-2.767310047, 15.625198798, 0.111294537],
+        91.5 => [-2.741498897, 15.656841259, 0.1122835261],
+        92.5 => [-2.715970894, 15.689303334, 0.113268793],
+        93.5 => [-2.690753197, 15.722567299, 0.1142496222],
+        94.5 => [-2.665870146, 15.756615553, 0.1152253207],
+        95.5 => [-2.641343436, 15.791430622, 0.1161952181],
+        96.5 => [-2.617192204, 15.826995169, 0.1171586674],
+        97.5 => [-2.593430614, 15.863292407, 0.1181150731],
+        98.5 => [-2.570076037, 15.900304841, 0.1190638073],
+        99.5 => [-2.547141473, 15.938015446, 0.1200042898],
+        100.5 => [-2.524635245, 15.976407871, 0.1209359936],
+        101.5 => [-2.502569666, 16.015464834, 0.1218583548],
+        102.5 => [-2.48095189, 16.055169844, 0.1227708703],
+        103.5 => [-2.459785573, 16.09550688, 0.1236730846],
+        104.5 => [-2.439080117, 16.136458809, 0.1245644841],
+        105.5 => [-2.418838304, 16.178009551, 0.125444639],
+        106.5 => [-2.399063683, 16.220142813, 0.1263131206],
+        107.5 => [-2.379756861, 16.262842771, 0.1271695453],
+        108.5 => [-2.360920527, 16.306093162, 0.1280135154],
+        109.5 => [-2.342557728, 16.349877586, 0.1288446388],
+        110.5 => [-2.324663326, 16.39418118, 0.1296626372],
+        111.5 => [-2.307240716, 16.438987413, 0.1304671382],
+        112.5 => [-2.290287663, 16.484280823, 0.1312578524],
+        113.5 => [-2.273803847, 16.530045538, 0.1320344789],
+        114.5 => [-2.257782149, 16.57626713, 0.132796819],
+        115.5 => [-2.242227723, 16.62292864, 0.1335445247],
+        116.5 => [-2.227132805, 16.670015716, 0.1342774356],
+        117.5 => [-2.212495585, 16.717512877, 0.1349953236],
+        118.5 => [-2.19831275, 16.765404961, 0.1356979956],
+        119.5 => [-2.184580762, 16.813676886, 0.1363852755],
+        120.5 => [-2.171295888, 16.862313656, 0.1370570042],
+        121.5 => [-2.158454232, 16.911300357, 0.1377130391],
+        122.5 => [-2.146051754, 16.960622156, 0.1383532537],
+        123.5 => [-2.134084303, 17.010264304, 0.1389775374],
+        124.5 => [-2.122547629, 17.060212133, 0.1395857952],
+        125.5 => [-2.111437411, 17.110451055, 0.1401779469],
+        126.5 => [-2.100749266, 17.160966564, 0.1407539274],
+        127.5 => [-2.090478774, 17.211744236, 0.1413136859],
+        128.5 => [-2.080621484, 17.262769728, 0.1418571858],
+        129.5 => [-2.071172932, 17.314028776, 0.1423844043],
+        130.5 => [-2.062128649, 17.365507199, 0.1428953318],
+        131.5 => [-2.053484173, 17.417190895, 0.143389972],
+        132.5 => [-2.045235058, 17.469065845, 0.1438683412],
+        133.5 => [-2.03737688, 17.52111811, 0.1443304685],
+        134.5 => [-2.029906684, 17.573333469, 0.1447763715],
+        135.5 => [-2.022817914, 17.625698688, 0.1452061381],
+        136.5 => [-2.016107084, 17.678199868, 0.1456198193],
+        137.5 => [-2.009769905, 17.730823397, 0.1460174906],
+        138.5 => [-2.003802134, 17.783555746, 0.1463992386],
+        139.5 => [-1.998199572, 17.836383471, 0.1467651605],
+        140.5 => [-1.992958064, 17.889293209, 0.1471153639],
+        141.5 => [-1.988073505, 17.942271684, 0.1474499668],
+        142.5 => [-1.983541835, 17.995305704, 0.1477690965],
+        143.5 => [-1.979359041, 18.048382162, 0.1480728906],
+        144.5 => [-1.975521156, 18.101488036, 0.1483614954],
+        145.5 => [-1.972024258, 18.154610394, 0.1486350668],
+        146.5 => [-1.968864465, 18.207736386, 0.1488937694],
+        147.5 => [-1.966037938, 18.260853253, 0.1491377764],
+        148.5 => [-1.963540872, 18.313948324, 0.1493672695],
+        149.5 => [-1.961369499, 18.367009017, 0.1495824386],
+        150.5 => [-1.959520079, 18.420022839, 0.1497834816],
+        151.5 => [-1.9579889, 18.472977388, 0.1499706043],
+        152.5 => [-1.956772271, 18.525860352, 0.1501440201],
+        153.5 => [-1.95586652, 18.578659513, 0.1503039498],
+        154.5 => [-1.955267984, 18.631362745, 0.1504506214],
+        155.5 => [-1.954973011, 18.683958013, 0.1505842702],
+        156.5 => [-1.954977947, 18.736433381, 0.1507051384],
+        157.5 => [-1.955279136, 18.788777004, 0.1508134748],
+        158.5 => [-1.955872909, 18.840977134, 0.1509095352],
+        159.5 => [-1.956755579, 18.893022121, 0.1509935818],
+        160.5 => [-1.957923436, 18.944900411, 0.1510658829],
+        161.5 => [-1.959372737, 18.996600549, 0.1511267136],
+        162.5 => [-1.9610997, 19.048111179, 0.1511763547],
+        163.5 => [-1.963100496, 19.099421046, 0.1512150935],
+        164.5 => [-1.96537124, 19.150518994, 0.1512432229],
+        165.5 => [-1.967907983, 19.201393971, 0.1512610419],
+        166.5 => [-1.970706706, 19.252035026, 0.1512688553],
+        167.5 => [-1.973763307, 19.302431312, 0.1512669735],
+        168.5 => [-1.977073595, 19.352572085, 0.1512557127],
+        169.5 => [-1.980633277, 19.402446707, 0.1512353947],
+        170.5 => [-1.984437954, 19.452044646, 0.1512063468],
+        171.5 => [-1.988483106, 19.501355476, 0.1511689019],
+        172.5 => [-1.992764085, 19.550368876, 0.1511233983],
+        173.5 => [-1.997276103, 19.599074637, 0.1510701797],
+        174.5 => [-2.002014224, 19.647462655, 0.1510095954],
+        175.5 => [-2.00697335, 19.695522937, 0.1509419999],
+        176.5 => [-2.012148213, 19.743245597, 0.1508677534],
+        177.5 => [-2.017533363, 19.790620862, 0.1507872211],
+        178.5 => [-2.023123159, 19.837639068, 0.1507007738],
+        179.5 => [-2.028911755, 19.884290662, 0.1506087878],
+        180.5 => [-2.034893091, 19.930566203, 0.1505116446],
+        181.5 => [-2.041060881, 19.976456361, 0.1504097312],
+        182.5 => [-2.047408604, 20.021951917, 0.1503034402],
+        183.5 => [-2.05392949, 20.067043765, 0.1501931693],
+        184.5 => [-2.060616513, 20.11172291, 0.1500793222],
+        185.5 => [-2.067462375, 20.155980469, 0.1499623077],
+        186.5 => [-2.074459502, 20.199807672, 0.1498425404],
+        187.5 => [-2.081600029, 20.243195857, 0.1497204407],
+        188.5 => [-2.088875793, 20.286136477, 0.1495964342],
+        189.5 => [-2.096278323, 20.328621093, 0.1494709526],
+        190.5 => [-2.103798828, 20.370641376, 0.1493444333],
+        191.5 => [-2.111428194, 20.412189106, 0.1492173194],
+        192.5 => [-2.119156972, 20.453256172, 0.14909006],
+        193.5 => [-2.126975375, 20.49383457, 0.1489631101],
+        194.5 => [-2.134873266, 20.5339164, 0.1488369306],
+        195.5 => [-2.142840157, 20.573493869, 0.1487119885],
+        196.5 => [-2.150865204, 20.612559285, 0.1485887569],
+        197.5 => [-2.158937201, 20.651105058, 0.1484677151],
+        198.5 => [-2.167044578, 20.689123698, 0.1483493484],
+        199.5 => [-2.175176987, 20.726607282, 0.1482341202],
+        200.5 => [-2.183317362, 20.763550105, 0.148122614],
+        201.5 => [-2.191457792, 20.799943374, 0.1480152488],
+        202.5 => [-2.199583649, 20.83578051, 0.1479125643],
+        203.5 => [-2.207681525, 20.871054493, 0.1478150781],
+        204.5 => [-2.215737645, 20.905758394, 0.1477233147],
+        205.5 => [-2.223739902, 20.939884769, 0.1476377678],
+        206.5 => [-2.231667995, 20.973428578, 0.1475590832],
+        207.5 => [-2.239511942, 21.006381705, 0.1474877162],
+        208.5 => [-2.247257081, 21.038737402, 0.1474242097],
+        209.5 => [-2.254885145, 21.070489959, 0.1473691743],
+        210.5 => [-2.26238209, 21.101632407, 0.147323144],
+        211.5 => [-2.269731517, 21.132158447, 0.1472866982],
+        212.5 => [-2.276917229, 21.162061708, 0.1472604146],
+        213.5 => [-2.283925442, 21.191335097, 0.1472448281],
+        214.5 => [-2.290731442, 21.219974716, 0.1472406828],
+        215.5 => [-2.29732427, 21.247972623, 0.147248467],
+        216.5 => [-2.303687802, 21.275322389, 0.1472687698],
+        217.5 => [-2.309799971, 21.302019325, 0.1473022986],
+        218.5 => [-2.315651874, 21.328054894, 0.1473495144],
+        219.5 => [-2.32121731, 21.353425629, 0.1474112153],
+        220.5 => [-2.326481911, 21.378124616, 0.1474879793],
+        221.5 => [-2.331428139, 21.402145892, 0.1475804525],
+        222.5 => [-2.336038473, 21.425483514, 0.1476892889],
+        223.5 => [-2.34029545, 21.448131558, 0.1478151501],
+        224.5 => [-2.344181703, 21.470084116, 0.1479587057],
+        225.5 => [-2.34768, 21.491335286, 0.1481206332],
+        226.5 => [-2.350773286, 21.511879176, 0.1483016185],
+        227.5 => [-2.353444725, 21.531709894, 0.1485023554],
+        228.5 => [-2.355677743, 21.550821547, 0.1487235462],
+        229.5 => [-2.35745607, 21.569208237, 0.1489659018],
+        230.5 => [-2.358763788, 21.586864057, 0.1492301415],
+        231.5 => [-2.359585369, 21.603783087, 0.1495169936],
+        232.5 => [-2.359905726, 21.619959388, 0.1498271951],
+        233.5 => [-2.359710258, 21.635387002, 0.1501614923],
+        234.5 => [-2.358980464, 21.650061262, 0.150520734],
+        235.5 => [-2.357714508, 21.663972695, 0.1509054394],
+        236.5 => [-2.355892424, 21.677117355, 0.1513165313],
+        237.5 => [-2.353501353, 21.689489352, 0.1517548077],
+        238.5 => [-2.350528726, 21.701082884, 0.1522210861],
+        239.5 => [-2.346962247, 21.711892252, 0.1527162055],
+        240.0 => [-2.34495843, 21.716999342, 0.152974718],
+        240.5 => [-2.342796948, 21.721909734, 0.1532408716]
+      }
+    end
+  end
+end

data/lib/health_stats/statistics.rb

@@ -0,0 +1,97 @@
+# The following were adapted from 
+# http://www.fourmilab.ch/rpkp/experiments/analysis/zCalc.html?
+# Original comments below:
+#
+#     The following functions for calculating normal and
+#     chi-square probabilities and critical values were adapted by
+#     John Walker from C implementations
+#     written by Gary Perlman of Wang Institute, Tyngsboro, MA
+#     01879.  Both the original C code and this edition
+#     are in the public domain.  */
+# 
+#     POZ  --  probability of normal z value
+# 
+#     Adapted from a polynomial approximation in:
+#             Ibbetson D, Algorithm 209
+#             Collected Algorithms of the CACM 1963 p. 616
+#     Note:
+#             This routine has six digit accuracy, so it is only useful for absolute
+#             z values <= 6.  For z values > to 6.0, poz() returns 0.0.
+#
+module HealthStats
+  module Statistics
+    extend self
+  
+    Z_MAX = 6
+    Z_EPSILON = 0.000001
+
+    def poz(z)
+      if z == 0.0
+        x = 0.0
+      else
+        y = 0.5 * z.abs
+    
+        if y > (Z_MAX * 0.5)
+          x = 1.0
+        elsif y < 1.0
+          w = y ** 2
+          x = ((((((((0.000124818987 * w -
+                      0.001075204047) * w + 0.005198775019) * w -
+                      0.019198292004) * w + 0.059054035642) * w -
+                      0.151968751364) * w + 0.319152932694) * w -
+                      0.531923007300) * w + 0.797884560593) * y * 2.0
+        else
+          y -= 2.0
+          x = (((((((((((((-0.000045255659 * y +
+                            0.000152529290) * y - 0.000019538132) * y -
+                            0.000676904986) * y + 0.001390604284) * y -
+                            0.000794620820) * y - 0.002034254874) * y +
+                            0.006549791214) * y - 0.010557625006) * y +
+                            0.011630447319) * y - 0.009279453341) * y +
+                            0.005353579108) * y - 0.002141268741) * y +
+                            0.000535310849) * y + 0.999936657524
+        end
+      end
+  
+      z > 0.0 ? ((x + 1.0) * 0.5) : ((1.0 - x) * 0.5)
+    end
+
+    def critz(p)
+      minz = -Z_MAX
+      maxz = Z_MAX
+      zval = 0.0
+  
+      return -1 if p < 0.0 || p > 1.0
+  
+      while maxz - minz > Z_EPSILON
+        pval = poz(zval)
+    
+        if pval > p
+          maxz = zval
+        else
+          minz = zval
+        end
+    
+        zval = (maxz + minz) * 0.5
+      end
+  
+      zval
+    end
+    
+     # LMS method
+     # Box GE, Cox DR. An analysis of transformations. J Roy Stat Soc, Series 
+     #      B. 26:211-252, 1964.
+     # Cole TJ. The LMS method for constructing normalized growth standards. 
+     #      Eur J Clin Nutr. 44: 45-60, 1990
+     # Cole TJ. Fitting smoothed centile curves to reference data. 
+     #      Royal Stat Soc. 151:385-418, 1988
+     
+    def lms(x, l, m, s)
+      (((x / m) ** l) - 1) / (l * s)
+    end
+    
+    # puts poz(lms(17.08, -1.924100169, 16.494427632, 0.0792339937)) * 100
+  end
+end
+  # puts critz(0.85)
+  # 

data/spec/spec_health_stats.rb

@@ -0,0 +1,119 @@
+$:.unshift File.expand_path(File.dirname(__FILE__) + '/../lib')
+require 'health_stats'
+require 'activesupport'
+
+REQUIRED_ATTRIBUTES = [:height, :weight, :gender, :dob]
+
+describe 'HealthStats' do
+  describe 'required attributes' do
+    before do
+      @klass = Class.new do
+        include HealthStats
+      end
+    end
+    
+    REQUIRED_ATTRIBUTES.each do |attribute|
+      it "should include #{attribute}" do
+        person = @klass.new
+        lambda { person.send(attribute) }.should.raise(HealthStats::AttributeError)
+      end
+    end
+  end
+  
+  describe 'included methods' do
+    before do
+      @klass = Class.new do
+        include HealthStats
+        attr_accessor *REQUIRED_ATTRIBUTES
+      end
+      
+      @person = @klass.new
+      @person.height = 53
+      @person.weight = 100
+      @person.gender = 'f'
+      @person.dob = 14.years.ago - 5.days
+    end
+    
+    describe '#age_methods' do
+      it "should return nil if dob is nil" do
+        @person.dob = nil
+        @person.age.should.be.nil
+        @person.age_in_years.should.be.nil
+        @person.age_in_months.should.be.nil
+      end
+      
+      it "should return years as integer" do
+        @person.age_in_years.should == 14
+        @person.age.should == 14
+        
+        @person.dob = 1.month.ago
+        @person.age_in_years.should == 0
+        @person.age.should == 0
+        
+        @person.dob = 14.months.ago
+        @person.age_in_years.should == 1
+        @person.age.should == 1
+      end
+      
+      it "should include age in months as a float in .5 increments" do
+        @person.age_in_months.should == 168.0
+        
+        @person.dob = 1.month.ago + 2.days
+        @person.age_in_months.should == 0.5
+        
+        @person.dob = 1.month.ago - 3.days
+        @person.age_in_months.should == 1.0
+        
+        @person.dob = 1.year.ago - 5.days
+        @person.age_in_months.should == 12.0
+        
+        @person.dob = 1.year.ago + 5.days
+        @person.age_in_months.should == 11.5
+        
+        @person.dob = 2.years.ago - 1.month - 18.days
+        @person.age_in_months.should. == 25.5
+      end
+      
+    end
+    
+    describe '#bmi' do
+      [:height, :weight].each do |attribute|
+        it "should return nil if #{attribute} is nil" do
+          @person.send("#{attribute}=", nil)
+          @person.bmi.should.be.nil
+        end
+      end
+    
+      it 'should return the correct bmi (precision 2)' do
+        @person.bmi.should == 25.02
+      end
+    end
+  
+    describe '#bmi_percentile' do    
+      [:dob, :gender].each do |attribute|
+        it "should return nil if #{attribute} is nil" do
+          @person.send("#{attribute}=", nil)
+          @person.bmi_percentile.should.be.nil
+        end
+      end
+      
+      it "should return nil if outside cdc date range" do
+        @person.dob = 1.month.ago
+        @person.bmi_percentile.should.be.nil
+        
+        @person.dob == 20.years.ago - 10.days - 2.months
+        @person.bmi_percentile.should.be.nil
+      end
+      
+      it "should return percentile (precision 2)" do
+        @person.bmi_percentile.should == 91.02
+        
+        @person.gender = 'm'
+        @person.dob = 2.years.ago - 1.month - 18.days
+        @person.weight = 14
+        @person.height = 24
+        @person.bmi_percentile.should == 66.49
+      end
+    end
+  end
+end

metadata

@@ -0,0 +1,77 @@
+--- !ruby/object:Gem::Specification 
+name: health_stats
+version: !ruby/object:Gem::Version 
+  version: 0.0.1
+platform: ruby
+authors: 
+- Mark Martin
+autorequire: 
+bindir: bin
+cert_chain: []
+
+date: 2009-12-03 00:00:00 -08:00
+default_executable: 
+dependencies: 
+- !ruby/object:Gem::Dependency 
+  name: activesupport
+  type: :development
+  version_requirement: 
+  version_requirements: !ruby/object:Gem::Requirement 
+    requirements: 
+    - - ">="
+      - !ruby/object:Gem::Version 
+        version: "2.3"
+    version: 
+description: 
+email: assplecake@gmail.com
+executables: []
+
+extensions: []
+
+extra_rdoc_files: 
+- LICENSE
+- README.md
+files: 
+- LICENSE
+- README.md
+- Rakefile
+- VERSION
+- health_stats.gemspec
+- lib/health_stats.rb
+- lib/health_stats/age.rb
+- lib/health_stats/bmi.rb
+- lib/health_stats/cdc_data.rb
+- lib/health_stats/cdc_data/bmi.rb
+- lib/health_stats/statistics.rb
+- pkg/health_stats-0.0.1.gem
+- spec/spec_health_stats.rb
+has_rdoc: true
+homepage: http://github.com/assplecake/healthdata
+licenses: []
+
+post_install_message: 
+rdoc_options: 
+- --charset=UTF-8
+require_paths: 
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement 
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+  version: 
+required_rubygems_version: !ruby/object:Gem::Requirement 
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+  version: 
+requirements: []
+
+rubyforge_project: 
+rubygems_version: 1.3.5
+signing_key: 
+specification_version: 3
+summary: Convenience methods for age, bmi, bmi percentile
+test_files: 
+- spec/spec_health_stats.rb