ValueError: could not broadcast input array from shape

[wariobrega]

Hello,

I have a dataset made of unpaired RNA/ATAC seq data.

I performed PCA / nearest neighbors on each modality using a standard approach and I was going for the joint neighbors embedding.

I obtain however the following error:

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
[/home/dcapocef/Projects/EPINet/scripts/screen_orig_data.ipynb](https://file+.vscode-resource.vscode-cdn.net/home/dcapocef/Projects/EPINet/scripts/screen_orig_data.ipynb) Cell 31 line 1
----> [1](vscode-notebook-cell:/home/dcapocef/Projects/EPINet/scripts/screen_orig_data.ipynb#X54sZmlsZQ%3D%3D?line=0) mu.pp.neighbors(ctrl)

File [~/miniconda3/envs/muon/lib/python3.10/site-packages/muon/_core/preproc.py:465](https://file+.vscode-resource.vscode-cdn.net/home/dcapocef/Projects/EPINet/scripts/~/miniconda3/envs/muon/lib/python3.10/site-packages/muon/_core/preproc.py:465), in neighbors(mdata, n_neighbors, n_bandwidth_neighbors, n_multineighbors, neighbor_keys, metric, low_memory, key_added, weight_key, add_weights_to_modalities, eps, copy, random_state)
    [463](https://file+.vscode-resource.vscode-cdn.net/home/dcapocef/Projects/EPINet/scripts/~/miniconda3/envs/muon/lib/python3.10/site-packages/muon/_core/preproc.py:463)         currtheta = theta
    [464](https://file+.vscode-resource.vscode-cdn.net/home/dcapocef/Projects/EPINet/scripts/~/miniconda3/envs/muon/lib/python3.10/site-packages/muon/_core/preproc.py:464)     else:
--> [465](https://file+.vscode-resource.vscode-cdn.net/home/dcapocef/Projects/EPINet/scripts/~/miniconda3/envs/muon/lib/python3.10/site-packages/muon/_core/preproc.py:465)         thetas[:, lasti] = theta
    [466](https://file+.vscode-resource.vscode-cdn.net/home/dcapocef/Projects/EPINet/scripts/~/miniconda3/envs/muon/lib/python3.10/site-packages/muon/_core/preproc.py:466)         lasti += 1
    [467](https://file+.vscode-resource.vscode-cdn.net/home/dcapocef/Projects/EPINet/scripts/~/miniconda3/envs/muon/lib/python3.10/site-packages/muon/_core/preproc.py:467) ratios[ratioidx, i1] = currtheta [/](https://file+.vscode-resource.vscode-cdn.net/) (np.max(thetas, axis=1) + eps)

ValueError: could not broadcast input array from shape (0,) into shape (5106,)

This usually happens after a couple of minutes after running the command, so way along the line

The data is made of 5106 Cells from ATAC and 4215 Cells from RNA.

I am using the latest stable version of muon on a conda environment running python 3.10. The conda env contains only muon and its minimal dependencies + ipython kernel/jupyter specs.

Can you help me understand this issue and how to fix it?

Thanks,

Daniele