<s>[INST] Possible bug with statistical variant in Integral Discrete Ordination
Hello, I am performing a statistical (also known as probability-based) variant of Integral Discrete Ordination on a set of binary data. I am concerned about the existence of a bug in the algorithm and I wanted to confirm it with the developers before I open an issue on GitHub. 
I'm using the `integral_discrete_ordination()` function to perform the analysis. However, I notice that the resulting eigenvalues appear to be not entirely sorted in ascending order as they should be.
This is the minimal reproducible code to recreate the issue:
```python
import numpy as np
from integrated-ordination import integral_discrete_ordination

# Create a sample dataset
np.random.seed(42)
n_rows, n_cols = 100, 10
data = np.random.randint(2, size=(n_rows, n_cols))

# Perform integral discrete ordination
eigenvalues, coordinates = integral_discrete_ordination(data, method="statistical")

# Verify eigenvalues are sorted
if not np.all(np.diff(eigenvalues) <= 1e-6):
    print("Eigenvalues are not sorted in ascending order.")
```
Have I truly uncovered a bug in Integral Discrete Ordination, or am I misunderstanding its implementation?[/INST]