Class: Statsample::Regression::Multiple::RubyEngine

Inherits:
MatrixEngine show all
Defined in:
lib/statsample/regression/multiple/rubyengine.rb

Overview

lr=Statsample::Regression::Multiple::RubyEngine.new(ds,‘y’)

Instance Attribute Summary

Attributes inherited from MatrixEngine

#cases, #digits, #x_mean, #x_sd, #y_mean, #y_sd

Attributes inherited from BaseEngine

#cases, #digits, #name, #total_cases, #valid_cases

Instance Method Summary collapse

Methods inherited from MatrixEngine

#coeffs, #coeffs_se, #constant, #constant_t, #df_e, #df_r, #r, #r2, #sst, #standarized_coeffs, #tolerance

Methods inherited from BaseEngine

#anova, #assign_names, #coeffs_se, #coeffs_t, #coeffs_tolerances, #constant_t, #df_e, #df_r, #estimated_variance_covariance_matrix, #f, #mse, #msr, #predicted, #probability, #process, #r, #r2_adjusted, #report_building, #residuals, #se_estimate, #se_r2, #sse, #sse_direct, #ssr, #ssr_direct, #sst, #standarized_predicted, #tolerance, univariate?

Methods included from Summarizable

#summary

Constructor Details

#initialize(ds, y_var, opts = Hash.new) ⇒ RubyEngine 



19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
 
# File 'lib/statsample/regression/multiple/rubyengine.rb', line 19

def initialize(ds,y_var, opts=Hash.new)
  matrix=ds.correlation_matrix
  fields_indep=ds.fields-[y_var]
  default={
    :y_mean=>ds[y_var].mean,
    :x_mean=>fields_indep.inject({}) {|ac,f|  ac[f]=ds[f].mean; ac},
    :y_sd=>ds[y_var].sd,
    :x_sd=>fields_indep.inject({}) {|ac,f|  ac[f]=ds[f].sd; ac},
    :cases=>Statsample::Bivariate.min_n_valid(ds)
  }
  opts=opts.merge(default)
  super(matrix, y_var, opts)
  @ds=ds
  @dy=ds[@y_var]
  @ds_valid=ds.dup_only_valid
  @total_cases=@ds.cases
  @valid_cases=@ds_valid.cases
  @ds_indep = ds.dup(ds.fields-[y_var])
  set_dep_columns
end

Instance Method Details

#constant_se ⇒ Object

Standard error for constant



84
85
86
 
# File 'lib/statsample/regression/multiple/rubyengine.rb', line 84

def constant_se
  estimated_variance_covariance_matrix[0,0]
end

#fix_with_mean ⇒ Object 



47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
 
# File 'lib/statsample/regression/multiple/rubyengine.rb', line 47

def fix_with_mean
  i=0
  @ds_indep.each do |row|
    empty=[]
    row.each do |k,v|
      empty.push(k) if v.nil?
    end
    if empty.size==1
      @ds_indep[empty[0]][i]=@ds[empty[0]].mean
    end
    i+=1
  end
  @ds_indep.update_valid_data
  set_dep_columns
end

#fix_with_regression ⇒ Object 



62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
 
# File 'lib/statsample/regression/multiple/rubyengine.rb', line 62

def fix_with_regression
  i=0
  @ds_indep.each{|row|
    empty=[]
    row.each{|k,v|
        empty.push(k) if v.nil?
    }
    if empty.size==1
      field=empty[0]
      lr=MultipleRegression.new(@ds_indep,field)
      fields=[]
      @ds_indep.fields.each{|f|
          fields.push(row[f]) unless f==field
      }
      @ds_indep[field][i]=lr.process(fields)
    end
    i+=1
  }
  @ds_indep.update_valid_data
  set_dep_columns
end

#set_dep_columns ⇒ Object 



40
41
42
43
44
45
 
# File 'lib/statsample/regression/multiple/rubyengine.rb', line 40

def set_dep_columns
  @dep_columns=[]
  @ds_indep.each_vector{|k,v|
    @dep_columns.push(v.data_with_nils)
  }                
end