59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
# File 'lib/rust/stats/correlation.rb', line 59
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='k')", true)
result = Result.new
result.name = "Kendall's rank correlation tau"
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
|