32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
# File 'lib/rust/stats/correlation.rb', line 32
def self.test(d1, d2)
raise TypeError, "Expecting Array of numerics" if !d1.is_a?(Array) || !d1.all? { |e| e.is_a?(Numeric) }
raise TypeError, "Expecting Array of numerics" if !d2.is_a?(Array) || !d2.all? { |e| e.is_a?(Numeric) }
Rust.exclusive do
Rust['correlation.a'] = d1
Rust['correlation.b'] = d2
_, warnings = Rust._eval("correlation.result <- cor.test(correlation.a, correlation.b, method='s')", true)
result = Result.new
result.name = "Spearman's rank correlation rho"
result.statistics['S'] = Rust._pull('correlation.result$statistic')
result.pvalue = Rust._pull('correlation.result$p.value')
result.correlation = Rust._pull('correlation.result$estimate')
result.exact = !warnings.include?("Cannot compute exact p-value with ties")
return result
end
end
|