statsample 1.0.1 → 1.1.0

data/lib/statsample/vector.rb

@@ -520,6 +520,117 @@ module Statsample
       }
     end
     
+    # == Bootstrap
+    # Generate +nr+ resamples (with replacement) of size  +s+
+    # from vector, computing each estimate from +estimators+
+    # over each resample. 
+    # +estimators+ could be
+    # a) Hash with variable names as keys and lambdas as  values
+    #   a.bootstrap(:log_s2=>lambda {|v| Math.log(v.variance)},1000)
+    # b) Array with names of method to bootstrap
+    #   a.bootstrap([:mean, :sd],1000)
+    # c) A single method to bootstrap
+    #   a.jacknife(:mean, 1000)
+    # If s is nil, is set to vector size by default.
+    # 
+    # Returns a dataset where each vector is an vector
+    # of length +nr+ containing the computed resample estimates.
+    def bootstrap(estimators, nr, s=nil)
+      s||=n
+      
+      h_est, es, bss= prepare_bootstrap(estimators)
+      
+     
+      nr.times do |i|
+        bs=sample_with_replacement(s)
+        es.each do |estimator|          
+          # Add bootstrap
+          bss[estimator].push(h_est[estimator].call(bs))
+        end
+      end
+      
+      es.each do |est|
+        bss[est]=bss[est].to_scale
+        bss[est].type=:scale
+      end
+      bss.to_dataset
+      
+    end
+    
+    # == Jacknife
+    # Returns a dataset with jacknife delete-+k+ +estimators+ 
+    # +estimators+ could be:
+    # a) Hash with variable names as keys and lambdas as values
+    #   a.jacknife(:log_s2=>lambda {|v| Math.log(v.variance)})
+    # b) Array with method names to jacknife
+    #   a.jacknife([:mean, :sd])
+    # c) A single method to jacknife
+    #   a.jacknife(:mean)
+    # +k+ represent the block size for block jacknife. By default
+    # is set to 1, for classic delete-one jacknife.
+    #
+    # Returns a dataset where each vector is an vector
+    # of length +cases+/+k+ containing the computed jacknife estimates.
+    # 
+    # == Reference:
+    # * Sawyer, S. (2005). Resampling Data: Using a Statistical Jacknife.
+    def jacknife(estimators, k=1)
+      raise "n should be divisible by k:#{k}" unless n%k==0
+      
+      nb=(n / k).to_i
+      
+      
+      h_est, es, ps= prepare_bootstrap(estimators)
+
+      est_n=es.inject({}) {|h,v|
+        h[v]=h_est[v].call(self)
+        h
+      }
+      
+      
+      nb.times do |i|
+        other=@data_with_nils.dup
+        other.slice!(i*k,k)
+        other=other.to_scale
+        es.each do |estimator|
+          # Add pseudovalue
+          ps[estimator].push( nb * est_n[estimator] - (nb-1) * h_est[estimator].call(other))
+        end
+      end
+      
+      
+      es.each do |est|
+        ps[est]=ps[est].to_scale
+        ps[est].type=:scale
+      end
+      ps.to_dataset
+    end
+    
+    
+    # For an array or hash of estimators methods, returns
+    # an array with three elements
+    # 1.- A hash with estimators names as keys and lambdas as values
+    # 2.- An array with estimators names
+    # 3.- A Hash with estimators names as keys and empty arrays as values
+    def prepare_bootstrap(estimators)
+      h_est=estimators
+      
+      h_est=[h_est] unless h_est.is_a? Array or h_est.is_a? Hash
+      
+      if h_est.is_a? Array
+        h_est=h_est.inject({}) {|h,est|
+          h[est]=lambda {|v| v.send(est)}
+          h
+        }
+      end
+      
+      bss=h_est.keys.inject({}) {|h,v| h[v]=[];h}
+      
+      [h_est,h_est.keys, bss]
+      
+    end
+    private :prepare_bootstrap
+    
     # Returns an random sample of size n, with replacement,
     # only with valid data.
     #
@@ -596,7 +707,7 @@ module Statsample
     end
     
     def to_s
-    sprintf("Vector(type:%s, n:%d)[%s]",@type.to_s,@data.size, @data.collect{|d| d.nil? ? "nil":d}.join(","))
+      sprintf("Vector(type:%s, n:%d)[%s]",@type.to_s,@data.size, @data.collect{|d| d.nil? ? "nil":d}.join(","))
     end
     # Ugly name. Really, create a Vector for standard 'matrix' package.
     # <tt>dir</tt> could be :horizontal or :vertical
@@ -666,18 +777,27 @@ module Statsample
       b.section(:name=>name) do |s|
         s.text _("n :%d") % n        
         s.text _("n valid:%d") % n_valid
-        s.text  _("factors:%s") % factors.join(",")
-        s.text   _("mode: %s") % mode
-        s.table(:name=>_("Distribution")) do |t|
-          frequencies.sort.each do |k,v|
-            key=labels.has_key?(k) ? labels[k]:k
-            t.row [key, v , ("%0.2f%%" % (v.quo(n_valid)*100))]
+        if @type==:nominal
+          s.text  _("factors:%s") % factors.join(",") 
+          s.text   _("mode: %s") % mode 
+          
+          s.table(:name=>_("Distribution")) do |t|
+            frequencies.sort.each do |k,v|
+              key=labels.has_key?(k) ? labels[k]:k
+              t.row [key, v , ("%0.2f%%" % (v.quo(n_valid)*100))]
+            end
           end
         end
-        s.text _("median: %s") % median.to_s if(@type==:ordinal)
+        
+        s.text _("median: %s") % median.to_s if(@type==:ordinal or @type==:scale)
         if(@type==:scale)
           s.text _("mean: %0.4f") % mean
-          s.text _("sd: %0.4f") % sd.to_s
+          if sd
+            s.text _("std.dev.: %0.4f") % sd
+            s.text _("std.err.: %0.4f") % se
+            s.text _("skew: %0.4f") % skew
+            s.text _("kurtosis: %0.4f") % kurtosis
+          end
         end
       end
     end
@@ -822,12 +942,18 @@ module Statsample
     end
 
     # Population average deviation (denominator N)
+    # author: Al Chou
+   
     def average_deviation_population( m = nil )
       check_type :scale
       m ||= mean
       ( @scale_data.inject( 0 ) { |a, x| ( x - m ).abs + a } ).quo( n_valid )
     end
-
+    def median_absolute_deviation
+      med=median
+      recode {|x| (x-med).abs}.median
+    end
+    alias  :mad :median_absolute_deviation
     # Sample Variance (denominator n-1)
     def variance_sample(m=nil)
       check_type :scale
@@ -894,6 +1020,13 @@ module Statsample
         check_type :scale
         standard_deviation_sample.quo(mean)
     end
+    # Standard error of the distribution mean
+    # Calculated using sd/sqrt(n)
+    def standard_error
+      standard_deviation_sample.quo(Math.sqrt(valid_data.size))
+    end
+    alias :se :standard_error
+    
     alias_method :sdp, :standard_deviation_population
     alias_method :sds, :standard_deviation_sample
     alias_method :adp, :average_deviation_population

data/po/es/statsample.po

@@ -1,8 +1,8 @@
 msgid ""
 msgstr ""
-"Project-Id-Version: statsample 1.0.0\n"
-"POT-Creation-Date: 2011-01-28 18:40-0300\n"
-"PO-Revision-Date: 2011-01-28 18:41-0300\n"
+"Project-Id-Version: statsample 1.0.1\n"
+"POT-Creation-Date: 2011-03-03 12:03-0300\n"
+"PO-Revision-Date: 2011-03-03 12:05-0300\n"
 "Last-Translator: Claudio Bustos <clbustos@gmail.com>\n"
 "Language-Team: Desarrollador\n"
 "MIME-Version: 1.0\n"
@@ -15,49 +15,81 @@ msgstr ""
 msgid "F Test"
 msgstr "Prueba F"
 
-#: lib/statsample/test/t.rb:209
+#: lib/statsample/test/t.rb:82
+msgid "T Test"
+msgstr "Prueba T"
+
+#: lib/statsample/test/t.rb:83
+msgid "Estimate"
+msgstr "Estimado"
+
+#: lib/statsample/test/t.rb:84
+msgid "Std.Err.of Estimate"
+msgstr "Err.Est. del Estimado"
+
+#: lib/statsample/test/t.rb:114
+msgid "%s: %0.4f | %s: %0.4f"
+msgstr "%s: %0.4f | %s: %0.4f"
+
+#: lib/statsample/test/t.rb:120
+msgid "t(%d) = %0.4f, p=%0.4f (%s tails)"
+msgstr "t(%d) = %0.4f, p=%0.4f (%s colas)"
+
+#: lib/statsample/test/t.rb:121
+msgid "CI(%d%%): %0.4f - %0.4f"
+msgstr "IC(%d%%): %0.4f - %0.4f"
+
+#: lib/statsample/test/t.rb:190
+msgid "Sample mean: %0.4f | Sample sd: %0.4f | se : %0.4f"
+msgstr "Media de la muestra: %0.4f | DE de la muestra: %0.4f | EE : %0.4f"
+
+#: lib/statsample/test/t.rb:191
+msgid "Population mean: %0.4f"
+msgstr "Promedio población: %0.4f"
+
+#: lib/statsample/test/t.rb:292
 msgid "Mean and standard deviation"
 msgstr "Promedio y desviación estándar"
 
-#: lib/statsample/test/t.rb:209
+#: lib/statsample/test/t.rb:292
 #: lib/statsample/regression/simple.rb:109
 #: lib/statsample/factor/pca.rb:216
 #: lib/statsample/factor/principalaxis.rb:202
 msgid "Variable"
 msgstr "Variable"
 
-#: lib/statsample/test/t.rb:209
+#: lib/statsample/test/t.rb:292
 #: lib/statsample/dominanceanalysis/bootstrap.rb:208
 msgid "mean"
 msgstr "promedio"
 
-#: lib/statsample/test/t.rb:209
+#: lib/statsample/test/t.rb:292
 msgid "sd"
 msgstr "de"
 
-#: lib/statsample/test/t.rb:209
+#: lib/statsample/test/t.rb:292
 #: lib/statsample/factor/parallelanalysis.rb:103
 #: lib/statsample/factor/parallelanalysis.rb:111
 msgid "n"
 msgstr "n"
 
-#: lib/statsample/test/t.rb:213
+#: lib/statsample/test/t.rb:296
 msgid "Levene test for equality of variances"
 msgstr "Test de Levene para igualdad de variancas"
 
-#: lib/statsample/test/t.rb:215
+#: lib/statsample/test/t.rb:298
 msgid "T statistics"
 msgstr "Estadístico T"
 
-#: lib/statsample/test/t.rb:216
+#: lib/statsample/test/t.rb:299
 msgid "Equal variance"
 msgstr "Varianza Igual"
 
-#: lib/statsample/test/t.rb:217
+#: lib/statsample/test/t.rb:300
 msgid "Non equal variance"
 msgstr "Varianza Desigual"
 
-#: lib/statsample/test/t.rb:219
+#: lib/statsample/test/t.rb:302
 msgid "Effect size"
 msgstr "Tamaño del efecto"
 
@@ -272,7 +304,7 @@ msgid "Anova Two-Way on %s"
 msgstr "Anova de dos vías en %s"
 
 #: lib/statsample/anova/twoway.rb:184
-#: lib/statsample/anova/oneway.rb:101
+#: lib/statsample/anova/oneway.rb:127
 msgid "Test of Homogeneity of variances (Levene)"
 msgstr "Test de homogeneidad de varianza (Levene)"
 
@@ -289,22 +321,38 @@ msgstr "Varianza explicada"
 msgid "Unexplained variance"
 msgstr "Varianza sin explicar"
 
-#: lib/statsample/anova/oneway.rb:91
+#: lib/statsample/anova/oneway.rb:97
 msgid "Anova One-Way"
 msgstr "Anova de una vía"
 
-#: lib/statsample/anova/oneway.rb:92
+#: lib/statsample/anova/oneway.rb:98
 msgid "Between Groups"
 msgstr "Entre grupos"
 
-#: lib/statsample/anova/oneway.rb:93
+#: lib/statsample/anova/oneway.rb:99
 msgid "Within Groups"
 msgstr "Dentro de grupos"
 
-#: lib/statsample/anova/oneway.rb:137
+#: lib/statsample/anova/oneway.rb:119
+msgid "Contrast for %s"
+msgstr "Contraste para %s"
+
+#: lib/statsample/anova/oneway.rb:163
 msgid "Descriptives"
 msgstr "Descriptivos"
 
+#: lib/statsample/anova/contrast.rb:13
+msgid "Psi estimate"
+msgstr "Psi Estimado"
+
+#: lib/statsample/anova/contrast.rb:14
+msgid "Contrast"
+msgstr "Contraste"
+
+#: lib/statsample/anova/contrast.rb:73
+msgid "Contrast:%s"
+msgstr "Contraste: %s"
+
 #: lib/statsample/graph/scatterplot.rb:72
 msgid "Scatterplot (%s - %s)"
 msgstr "Diagrama de dispersión (%s - %s)"
@@ -540,27 +588,27 @@ msgstr "Análisis de múltiples escalas"
 msgid "Scale %s"
 msgstr "Escala %s"
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:143
+#: lib/statsample/reliability/multiscaleanalysis.rb:145
 msgid "Reliability analysis of scales"
 msgstr "Análisis de confiabilidad de escalas"
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:149
+#: lib/statsample/reliability/multiscaleanalysis.rb:151
 msgid "Correlation matrix for %s"
 msgstr "Matriz de correlaciones para %s"
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:154
+#: lib/statsample/reliability/multiscaleanalysis.rb:156
 msgid "PCA for %s"
 msgstr "ACP para %s"
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:159
+#: lib/statsample/reliability/multiscaleanalysis.rb:161
 msgid "Principal Axis for %s"
 msgstr "Ejes principales para %s"
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:165
+#: lib/statsample/reliability/multiscaleanalysis.rb:167
 msgid "Parallel Analysis for %s"
 msgstr "Análisis Paralelo para %s"
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:170
+#: lib/statsample/reliability/multiscaleanalysis.rb:172
 msgid "MAP for %s"
 msgstr "MAP para %s"
 
@@ -600,96 +648,96 @@ msgstr "p: %0.3f"
 msgid "No problematic items"
 msgstr "Sin ítems problemáticos"
 
-#: lib/statsample/reliability/scaleanalysis.rb:45
+#: lib/statsample/reliability/scaleanalysis.rb:44
 msgid "Reliability Analisis"
 msgstr "Análisis de confiabilidad"
 
-#: lib/statsample/reliability/scaleanalysis.rb:158
+#: lib/statsample/reliability/scaleanalysis.rb:157
 msgid "Summary for %s with all items"
 msgstr "Sumario para %s con todos los ítems"
 
-#: lib/statsample/reliability/scaleanalysis.rb:159
+#: lib/statsample/reliability/scaleanalysis.rb:158
 msgid "Items"
 msgstr "Ítems"
 
-#: lib/statsample/reliability/scaleanalysis.rb:160
-#: lib/statsample/reliability/scaleanalysis.rb:177
+#: lib/statsample/reliability/scaleanalysis.rb:159
+#: lib/statsample/reliability/scaleanalysis.rb:176
 msgid "Sum mean"
 msgstr "Promedio de suma"
 
-#: lib/statsample/reliability/scaleanalysis.rb:161
+#: lib/statsample/reliability/scaleanalysis.rb:160
 msgid "S.d. mean"
 msgstr "Promedio de d.e."
 
-#: lib/statsample/reliability/scaleanalysis.rb:163
+#: lib/statsample/reliability/scaleanalysis.rb:162
 msgid "Deleted items"
 msgstr "Ítems eliminados"
 
-#: lib/statsample/reliability/scaleanalysis.rb:173
+#: lib/statsample/reliability/scaleanalysis.rb:172
 msgid "Summary for %s"
 msgstr "Sumario para %s"
 
-#: lib/statsample/reliability/scaleanalysis.rb:174
+#: lib/statsample/reliability/scaleanalysis.rb:173
 msgid "Valid Items"
 msgstr "Ítems Válidos"
 
-#: lib/statsample/reliability/scaleanalysis.rb:176
+#: lib/statsample/reliability/scaleanalysis.rb:175
 msgid "Valid cases"
 msgstr "casos válidos"
 
-#: lib/statsample/reliability/scaleanalysis.rb:178
+#: lib/statsample/reliability/scaleanalysis.rb:177
 msgid "Sum sd"
 msgstr "d.e. de suma"
 
-#: lib/statsample/reliability/scaleanalysis.rb:180
+#: lib/statsample/reliability/scaleanalysis.rb:179
 msgid "Sum median"
 msgstr "Mediana de suma"
 
-#: lib/statsample/reliability/scaleanalysis.rb:182
+#: lib/statsample/reliability/scaleanalysis.rb:181
 msgid "Item mean"
 msgstr "Promedio de los ítemes"
 
-#: lib/statsample/reliability/scaleanalysis.rb:183
+#: lib/statsample/reliability/scaleanalysis.rb:182
 msgid "Item sd"
 msgstr "DE de Items"
 
-#: lib/statsample/reliability/scaleanalysis.rb:185
+#: lib/statsample/reliability/scaleanalysis.rb:184
 msgid "Skewness"
 msgstr "Sesgo"
 
-#: lib/statsample/reliability/scaleanalysis.rb:186
+#: lib/statsample/reliability/scaleanalysis.rb:185
 msgid "Kurtosis"
 msgstr "Curtosis"
 
-#: lib/statsample/reliability/scaleanalysis.rb:188
+#: lib/statsample/reliability/scaleanalysis.rb:187
 msgid "Cronbach's alpha"
 msgstr "Alfa de Cronbach"
 
-#: lib/statsample/reliability/scaleanalysis.rb:189
+#: lib/statsample/reliability/scaleanalysis.rb:188
 msgid "Standarized Cronbach's alpha"
 msgstr "Alfa de Cronbach estandarizado"
 
-#: lib/statsample/reliability/scaleanalysis.rb:190
+#: lib/statsample/reliability/scaleanalysis.rb:189
 msgid "Mean rpb"
 msgstr "rbp medio"
 
-#: lib/statsample/reliability/scaleanalysis.rb:192
+#: lib/statsample/reliability/scaleanalysis.rb:191
 msgid "Variances mean"
 msgstr "Promedio de las varianzas"
 
-#: lib/statsample/reliability/scaleanalysis.rb:193
+#: lib/statsample/reliability/scaleanalysis.rb:192
 msgid "Covariances mean"
 msgstr "Promedio de las covarianzas"
 
-#: lib/statsample/reliability/scaleanalysis.rb:197
+#: lib/statsample/reliability/scaleanalysis.rb:196
 msgid "Items for obtain alpha(0.8) : %d"
 msgstr "Ítems para obtener alfa(0,8): %d"
 
-#: lib/statsample/reliability/scaleanalysis.rb:198
+#: lib/statsample/reliability/scaleanalysis.rb:197
 msgid "Items for obtain alpha(0.9) : %d"
 msgstr "Ítems para obtener alfa(0,9): %d"
 
-#: lib/statsample/reliability/scaleanalysis.rb:206
+#: lib/statsample/reliability/scaleanalysis.rb:205
 msgid "Items report for %s"
 msgstr "Reporte de ítems para %s"
 
@@ -857,104 +905,55 @@ msgstr "%s(centrado)"
 msgid "%s(percentil)"
 msgstr "%s(percentil)"
 
-#: lib/statsample/vector.rb:667
+#: lib/statsample/vector.rb:778
 msgid "n :%d"
 msgstr "n: %s"
 
-#: lib/statsample/vector.rb:668
+#: lib/statsample/vector.rb:779
 msgid "n valid:%d"
 msgstr "n válido: %d"
 
-#: lib/statsample/vector.rb:669
+#: lib/statsample/vector.rb:780
 msgid "factors:%s"
 msgstr "factores:%s"
 
-#: lib/statsample/vector.rb:670
+#: lib/statsample/vector.rb:781
 msgid "mode: %s"
 msgstr "modo: %s"
 
-#: lib/statsample/vector.rb:671
+#: lib/statsample/vector.rb:782
 msgid "Distribution"
 msgstr "Distribución"
 
-#: lib/statsample/vector.rb:677
+#: lib/statsample/vector.rb:788
 msgid "median: %s"
 msgstr "Mediana: %s"
 
-#: lib/statsample/vector.rb:679
+#: lib/statsample/vector.rb:790
 msgid "mean: %0.4f"
 msgstr "promedio: %0.3f"
 
-#: lib/statsample/vector.rb:680
+#: lib/statsample/vector.rb:791
 msgid "sd: %0.4f"
 msgstr "d.e.: %0.3f"
 
-#: lib/statsample/dataset.rb:134
+#: lib/statsample/dataset.rb:161
 msgid "Dataset %d"
 msgstr "Dataset %d"
 
-#: lib/statsample/dataset.rb:427
+#: lib/statsample/dataset.rb:457
 msgid "Sum from %s"
 msgstr "Suma para %s"
 
-#: lib/statsample/dataset.rb:480
+#: lib/statsample/dataset.rb:510
 msgid "Means from %s"
 msgstr "Media desde %s"
 
-#: lib/statsample/dataset.rb:704
+#: lib/statsample/dataset.rb:734
 msgid "%s(filtered)"
 msgstr "%s(filtrado)"
 
-#: lib/statsample/dataset.rb:926
+#: lib/statsample/dataset.rb:956
 msgid "Cases: %d"
 msgstr "Casos: %s"
 
-#~ msgid "R=%0.3f"
-#~ msgstr "R: %0.3f"
-#~ msgid "Sum variance"
-#~ msgstr "Varianza de suma"
-#~ msgid "actual"
-#~ msgstr "real"
-#~ msgid "Polychoric correlation"
-#~ msgstr "Correlación policórica"
-#~ msgid "Minimizing using GSL Brent method\n"
-#~ msgstr "Minimizando usando método GSL Brent\n"
-#~ msgid "Two step minimization using %s method\n"
-#~ msgstr "Minimización en dos etapas usando método %s\n"
-#~ msgid "Contingence Table"
-#~ msgstr "Tabla de Contingencia"
-#~ msgid "Thresholds"
-#~ msgstr "Umbrales"
-
-#, fuzzy
-#~ msgid "Threshold X %d"
-#~ msgstr "Umbral X: %0.3f"
-
-#, fuzzy
-#~ msgid "Threshold Y %d"
-#~ msgstr "Umbral Y:%0.3f"
-#~ msgid "Test of bivariate normality: X2 = %0.3f, df = %d, p= %0.5f"
-#~ msgstr "Prueba de normalidad bivariada: X2 = %0.3f, g.l. = %d, p= %0.5f"
-#~ msgid "Threshold X: %0.3f "
-#~ msgstr "Umbral X: %0.3f"
-#~ msgid "Threshold Y: %0.3f "
-#~ msgstr "Umbral Y:%0.3f"
-#~ msgid "Factor Analysis: "
-#~ msgstr "Análisis de Factores:"
-#~ msgid "DAB: "
-#~ msgstr "RAD:"
-#~ msgid "Total mean"
-#~ msgstr "Promedio total"
-#~ msgid "Total sd"
-#~ msgstr "d.e. total"
-#~ msgid "Median"
-#~ msgstr "Mediana"
-#~ msgid "\"Anova Two-Way on #{@ds[dep_var].name}\""
-#~ msgstr "\"Anova de dos vías en #{@ds[dep_var].name}\""
-#~ msgid "Crosstab"
-#~ msgstr "Tabulación cruzada"
-#~ msgid "Crosstab: "
-#~ msgstr "Tabulación cruzada:"
-#~ msgid ")}"
-#~ msgstr ")}"
-