5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
# File 'lib/rust/stats/correlation.rb', line 5
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='p')", true)
result = Result.new
result.name = "Pearson's product-moment correlation"
result.statistics['t'] = 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
|