mudata

import anndata as ad import mudata as md import numpy as np import pandas as pd # Create dummy AnnData objects for two modalities (e.g., RNA and ATAC) adata_rna = ad.AnnData( X=np.random.rand(10, 20), # 10 cells, 20 genes obs=pd.DataFrame(index=[f"cell_{i}" for i in range(10)], data={'donor': np.random.choice(['D1', 'D2'], 10)}), var=pd.DataFrame(index=[f"gene_{j}" for j in range(20)]) ) adata_atac = ad.AnnData( X=np.random.randint(0, 5, size=(10, 15)), # 10 cells, 15 peaks obs=pd.DataFrame(index=[f"cell_{i}" for i in range(10)], data={'donor': np.random.choice(['D1', 'D2'], 10)}), var=pd.DataFrame(index=[f"peak_{j}" for j in range(15)]) ) # Create a MuData object from a dictionary of AnnData objects mdata = md.MuData({'rna': adata_rna, 'atac': adata_atac}) print("MuData object created:") print(mdata) # Access a specific modality print("\nAccessing the 'rna' modality:") print(mdata['rna']) # Access shared observation metadata print("\nShared observation metadata:") print(mdata.obs.head())