(counts: &[usize])
| 139 | } |
| 140 | |
| 141 | fn entropy(counts: &[usize]) -> CodependenceResult<f64> { |
| 142 | let total: usize = counts.iter().sum(); |
| 143 | if total == 0 { |
| 144 | return Err(CodependenceError::InputTooShort); |
| 145 | } |
| 146 | let total_f = total as f64; |
| 147 | |
| 148 | let mut value = 0.0; |
| 149 | for count in counts { |
| 150 | if *count == 0 { |
| 151 | continue; |
| 152 | } |
| 153 | let p = *count as f64 / total_f; |
| 154 | value -= p * p.ln(); |
| 155 | } |
| 156 | |
| 157 | Ok(value) |
| 158 | } |
| 159 | |
| 160 | pub fn angular_distance(x: &[f64], y: &[f64]) -> CodependenceResult<f64> { |
| 161 | let corr_coef = corrcoef(x, y)?; |
no outgoing calls
no test coverage detected