From f7870ad65f6506d9332e8bdb62d823485e4e8bf8 Mon Sep 17 00:00:00 2001
From: Sahil Sahni <86493572+sahil-sahni@users.noreply.github.com>
Date: Thu, 5 Mar 2026 15:28:18 -0500
Subject: [PATCH] Correct receptor expression references from ct1 to ct2

Updated receptor expression calculations to use ct2 instead of ct1.
---
 R/SOCIAL.R | 15 ++++++++-------
 1 file changed, 8 insertions(+), 7 deletions(-)

diff --git a/R/SOCIAL.R b/R/SOCIAL.R
index 42f57de..ac3d361 100644
--- a/R/SOCIAL.R
+++ b/R/SOCIAL.R
@@ -1,5 +1,6 @@
 #IRIS package
-#Sahil Sahni & Kun Wang 07/31/24
+#Sahil Sahni & Kun Wang 07/31/24;
+#Corrected by Sahil Sahni 03/05/26;
 
 ## SOCIAL
 
@@ -81,9 +82,9 @@ calculate_interaction_score <- function(expr,ct_map, pairs, n_iterations = 100){
       else
         ligand_exp = as.numeric(summary_expr[g1,ct1])
       if(length(g2) > 1)
-        receptor_exp = min(as.numeric(summary_expr[g2,ct1]))
+        receptor_exp = min(as.numeric(summary_expr[g2,ct2])) # 3/5/26: correct ct1 to ct2
       else
-        receptor_exp = as.numeric(summary_expr[g2,ct1])
+        receptor_exp = as.numeric(summary_expr[g2,ct2]) # 3/5/26: correct ct1 to ct2
 
       #interaction score
       sc = ligand_exp*receptor_exp
@@ -119,9 +120,9 @@ calculate_interaction_score <- function(expr,ct_map, pairs, n_iterations = 100){
       score = c()
       receptor_exp = NULL
       if(length(g2) > 1)
-        receptor_exp = min(as.numeric(summary_expr[g2,ct1]))
+        receptor_exp = min(as.numeric(summary_expr[g2,ct2])) # 3/5/26: correct ct1 to ct2
       else
-        receptor_exp = as.numeric(summary_expr[g2,ct1])
+        receptor_exp = as.numeric(summary_expr[g2,ct2]) # 3/5/26: correct ct1 to ct2
 
       #receptor expression
       sc = receptor_exp
@@ -171,9 +172,9 @@ calculate_interaction_score <- function(expr,ct_map, pairs, n_iterations = 100){
         else
           ligand_exp = as.numeric(summary_expr[g1,ct1])
         if(length(g2) > 1)
-          receptor_exp = min(as.numeric(summary_expr[g2,ct1]))
+          receptor_exp = min(as.numeric(summary_expr[g2,ct2])) # 3/5/26: correct ct1 to ct2
         else
-          receptor_exp = as.numeric(summary_expr[g2,ct1])
+          receptor_exp = as.numeric(summary_expr[g2,ct2]) # 3/5/26: correct ct1 to ct2
 
         #null score
         sc = ligand_exp*receptor_exp