Regulon_specificity_scores question

[dilkamur]

Dear Team:
I am interested to compare regulons in a specific cell type at specific condition: Seems that I can calculate only rss by cell type or by condition. My question is how to achive the goal to calculate rss for cell type and condition. I created the following code but not sure if this approach is correct:

filtered_data = {}
filtered_annotations = {}
rss_results = {}

 #Loop through each unique cell type from the DataFrame
for cell_type in cellAnnot['Cell_Type'].unique():
#Filter the cellAnnot DataFrame for the current cell type_
    cell_mask = cellAnnot['Cell_Type'] == cell_type
    filtered_data[cell_type] = auc_mtx.loc[cell_mask]
    
    #Filter the cellAnnot DataFrame for the same cells
    filtered_annotations[cell_type] = cellAnnot.loc[cell_mask]
    
     #Calculate RSS and store the result_
    rss_results[cell_type] = regulon_specificity_scores(filtered_data[cell_type],
                                                        filtered_annotations[cell_type]['Condition'])

#Example: Access the RSS results for a specific cell type
rss_results['CellType1']  

So this will give me rss for each regulon for each condition for a specific cell type. And then I would like to compare regulons by condition as the following

(rss_results['CellType1'].loc["condition1"]-rss_results['CellType2'].loc['condition2']).

Could you please provide your feedback/suggestions? Thank you in advance!!