SciRuby/statsample

#!/usr/bin/ruby
$:.unshift(File.dirname(__FILE__)+'/../lib')

# == Description
#
# Reliability Scale Analysis with statsample
require 'statsample'
Statsample::Analysis.store(Statsample::Reliability) do 
  samples=100
  a=rnorm(samples)
  
  ds = Daru::DataFrame.new({})
  
  20.times do |i|
    ds["v#{i}".to_sym]= a + rnorm(samples,0,0.2)
  end
  
  rel=Statsample::Reliability::ScaleAnalysis.new(ds)
  summary rel
  
  ms=Statsample::Reliability::MultiScaleAnalysis.new(:name=>"Multi Scale analyss") do |m|
    m.scale "Scale 1", ds.clone([:v1, :v2, :v3, :v4, :v5, :v6, :v7, :v8, :v9, :v10])
    m.scale "Scale 2", ds.clone([:v11, :v12, :v13, :v14, :v15, :v16, :v17, :v18, :v19])
  end
  
  summary ms
end

if __FILE__==$0
   Statsample::Analysis.run_batch
end