Differential Expression Analysis using Harmony Integrated Clusters
Nasir Mahmood Abbasi
2025-01-18
1. Load Libraries
2. Load Seurat Object
#Load Seurat Object merged from cell lines and a control(PBMC) after filtration
load("0-imp_Robj/Harmony_integrated_CD4Tcells_harmony_integrated.Robj")
All_samples_Merged An object of class Seurat
62926 features across 49360 samples within 6 assays
Active assay: SCT (26174 features, 3000 variable features)
3 layers present: counts, data, scale.data
5 other assays present: RNA, ADT, prediction.score.celltype.l1, prediction.score.celltype.l2, prediction.score.celltype.l3
6 dimensional reductions calculated: integrated_dr, ref.umap, pca, umap, harmony, umap.harmony3. Set Up Identifiers for Clustering
# Assign cluster identities to the Seurat object
Idents(All_samples_Merged) <- "seurat_clusters"
DimPlot(All_samples_Merged, reduction = "umap.harmony", group.by = "seurat_clusters",label = T, label.box = T, repel = T) +
ggtitle("Harmony Integration - By Clusters")
4. Differential Expression Analysis
4.1 Using MAST with Latent Variables (Batch Correction)
markers_mast_batch <- FindMarkers(
All_samples_Merged,
ident.1 = c("0", "1", "2", "3", "4", "5", "7", "9"), # Cell lines
ident.2 = c("6", "8", "10", "11"), # CD4 T cells
assay = "SCT",
test.use = "MAST",
latent.vars = c("cell_line"), # Adjust for batch effects
min.pct = 0,
logfc.threshold = 0
)Avis : sparse->dense coercion: allocating vector of size 9.6 GiBAvis : Dropping illegal slot(s) norm.method for class BayesGLMlike.
This likely indicates a bug in an upstream package.
Completed [-----------------------------------------------------------------------------] 0% with 0 failuresRegistered S3 method overwritten by 'rmarkdown':
method from
print.paged_df
Completed [-----------------------------------------------------------------------------] 1% with 0 failures
Completed [>----------------------------------------------------------------------------] 1% with 0 failures
Completed [>----------------------------------------------------------------------------] 2% with 0 failures
Completed [=>---------------------------------------------------------------------------] 2% with 0 failures
Completed [=>---------------------------------------------------------------------------] 3% with 0 failures
Completed [==>--------------------------------------------------------------------------] 3% with 0 failures
Completed [==>--------------------------------------------------------------------------] 4% with 0 failures
Completed [==>--------------------------------------------------------------------------] 5% with 0 failures
Completed [===>-------------------------------------------------------------------------] 5% with 0 failures
Completed [===>-------------------------------------------------------------------------] 6% with 0 failures
Completed [====>------------------------------------------------------------------------] 6% with 0 failures
Completed [====>------------------------------------------------------------------------] 7% with 0 failures
Completed [=====>-----------------------------------------------------------------------] 7% with 0 failures
Completed [=====>-----------------------------------------------------------------------] 8% with 0 failures
Completed [======>----------------------------------------------------------------------] 8% with 0 failures
Completed [======>----------------------------------------------------------------------] 9% with 0 failures
Completed [======>----------------------------------------------------------------------] 10% with 0 failures
Completed [=======>---------------------------------------------------------------------] 10% with 0 failures
Completed [=======>---------------------------------------------------------------------] 11% with 0 failures
Completed [========>--------------------------------------------------------------------] 11% with 0 failures
Completed [========>--------------------------------------------------------------------] 12% with 0 failures
Completed [=========>-------------------------------------------------------------------] 12% with 0 failures
Completed [=========>-------------------------------------------------------------------] 13% with 0 failures
Completed [=========>-------------------------------------------------------------------] 14% with 0 failures
Completed [==========>------------------------------------------------------------------] 14% with 0 failures
Completed [==========>------------------------------------------------------------------] 15% with 0 failures
Completed [===========>-----------------------------------------------------------------] 15% with 0 failures
Completed [===========>-----------------------------------------------------------------] 16% with 0 failures
Completed [============>----------------------------------------------------------------] 16% with 0 failures
Completed [============>----------------------------------------------------------------] 17% with 0 failures
Completed [============>----------------------------------------------------------------] 18% with 0 failures
Completed [=============>---------------------------------------------------------------] 18% with 0 failures
Completed [=============>---------------------------------------------------------------] 19% with 0 failures
Completed [==============>--------------------------------------------------------------] 19% with 0 failures
Completed [==============>--------------------------------------------------------------] 20% with 0 failures
Completed [===============>-------------------------------------------------------------] 20% with 0 failures
Completed [===============>-------------------------------------------------------------] 21% with 0 failures
Completed [================>------------------------------------------------------------] 21% with 0 failures
Completed [================>------------------------------------------------------------] 22% with 0 failures
Completed [================>------------------------------------------------------------] 23% with 0 failures
Completed [=================>-----------------------------------------------------------] 23% with 0 failures
Completed [=================>-----------------------------------------------------------] 24% with 0 failures
Completed [==================>----------------------------------------------------------] 24% with 0 failures
Completed [==================>----------------------------------------------------------] 25% with 0 failures
Completed [===================>---------------------------------------------------------] 25% with 0 failures
Completed [===================>---------------------------------------------------------] 26% with 0 failures
Completed [===================>---------------------------------------------------------] 27% with 0 failures
Completed [====================>--------------------------------------------------------] 27% with 0 failures
Completed [====================>--------------------------------------------------------] 28% with 0 failures
Completed [=====================>-------------------------------------------------------] 28% with 0 failures
Completed [=====================>-------------------------------------------------------] 29% with 0 failures
Completed [======================>------------------------------------------------------] 29% with 0 failures
Completed [======================>------------------------------------------------------] 30% with 0 failures
Completed [======================>------------------------------------------------------] 31% with 0 failures
Completed [=======================>-----------------------------------------------------] 31% with 0 failures
Completed [=======================>-----------------------------------------------------] 32% with 0 failures
Completed [========================>----------------------------------------------------] 32% with 0 failures
Completed [========================>----------------------------------------------------] 33% with 0 failures
Completed [=========================>---------------------------------------------------] 33% with 0 failures
Completed [=========================>---------------------------------------------------] 34% with 0 failures
Completed [==========================>--------------------------------------------------] 34% with 0 failures
Completed [==========================>--------------------------------------------------] 35% with 0 failures
Completed [==========================>--------------------------------------------------] 36% with 0 failures
Completed [===========================>-------------------------------------------------] 36% with 0 failures
Completed [===========================>-------------------------------------------------] 37% with 0 failures
Completed [============================>------------------------------------------------] 37% with 0 failures
Completed [============================>------------------------------------------------] 38% with 0 failures
Completed [=============================>-----------------------------------------------] 38% with 0 failures
Completed [=============================>-----------------------------------------------] 39% with 0 failures
Completed [=============================>-----------------------------------------------] 40% with 0 failures
Completed [==============================>----------------------------------------------] 40% with 0 failures
Completed [==============================>----------------------------------------------] 41% with 0 failures
Completed [===============================>---------------------------------------------] 41% with 0 failures
Completed [===============================>---------------------------------------------] 42% with 0 failures
Completed [================================>--------------------------------------------] 42% with 0 failures
Completed [================================>--------------------------------------------] 43% with 0 failures
Completed [================================>--------------------------------------------] 44% with 0 failures
Completed [=================================>-------------------------------------------] 44% with 0 failures
Completed [=================================>-------------------------------------------] 45% with 0 failures
Completed [==================================>------------------------------------------] 45% with 0 failures
Completed [==================================>------------------------------------------] 46% with 0 failures
Completed [===================================>-----------------------------------------] 46% with 0 failures
Completed [===================================>-----------------------------------------] 47% with 0 failures
Completed [====================================>----------------------------------------] 47% with 0 failures
Completed [====================================>----------------------------------------] 48% with 0 failures
Completed [====================================>----------------------------------------] 49% with 0 failures
Completed [=====================================>---------------------------------------] 49% with 0 failures
Completed [=====================================>---------------------------------------] 50% with 0 failures
Completed [======================================>--------------------------------------] 50% with 0 failures
Completed [======================================>--------------------------------------] 51% with 0 failures
Completed [=======================================>-------------------------------------] 51% with 0 failures
Completed [=======================================>-------------------------------------] 52% with 0 failures
Completed [=======================================>-------------------------------------] 53% with 0 failures
Completed [========================================>------------------------------------] 53% with 0 failures
Completed [========================================>------------------------------------] 54% with 0 failures
Completed [=========================================>-----------------------------------] 54% with 0 failures
Completed [=========================================>-----------------------------------] 55% with 0 failures
Completed [==========================================>----------------------------------] 55% with 0 failures
Completed [==========================================>----------------------------------] 56% with 0 failures
Completed [===========================================>---------------------------------] 56% with 0 failures
Completed [===========================================>---------------------------------] 57% with 0 failures
Completed [===========================================>---------------------------------] 58% with 0 failures
Completed [============================================>--------------------------------] 58% with 0 failures
Completed [============================================>--------------------------------] 59% with 0 failures
Completed [=============================================>-------------------------------] 59% with 0 failures
Completed [=============================================>-------------------------------] 60% with 0 failures
Completed [==============================================>------------------------------] 60% with 0 failures
Completed [==============================================>------------------------------] 61% with 0 failures
Completed [==============================================>------------------------------] 62% with 0 failures
Completed [===============================================>-----------------------------] 62% with 0 failures
Completed [===============================================>-----------------------------] 63% with 0 failures
Completed [================================================>----------------------------] 63% with 0 failures
Completed [================================================>----------------------------] 64% with 0 failures
Completed [=================================================>---------------------------] 64% with 0 failures
Completed [=================================================>---------------------------] 65% with 0 failures
Completed [=================================================>---------------------------] 66% with 0 failures
Completed [==================================================>--------------------------] 66% with 0 failures
Completed [==================================================>--------------------------] 67% with 0 failures
Completed [===================================================>-------------------------] 67% with 0 failures
Completed [===================================================>-------------------------] 68% with 0 failures
Completed [====================================================>------------------------] 68% with 0 failures
Completed [====================================================>------------------------] 69% with 0 failures
Completed [=====================================================>-----------------------] 69% with 0 failures
Completed [=====================================================>-----------------------] 70% with 0 failures
Completed [=====================================================>-----------------------] 71% with 0 failures
Completed [======================================================>----------------------] 71% with 0 failures
Completed [======================================================>----------------------] 72% with 0 failures
Completed [=======================================================>---------------------] 72% with 0 failures
Completed [=======================================================>---------------------] 73% with 0 failures
Completed [========================================================>--------------------] 73% with 0 failures
Completed [========================================================>--------------------] 74% with 0 failures
Completed [========================================================>--------------------] 75% with 0 failures
Completed [=========================================================>-------------------] 75% with 0 failures
Completed [=========================================================>-------------------] 76% with 0 failures
Completed [==========================================================>------------------] 76% with 0 failures
Completed [==========================================================>------------------] 77% with 0 failures
Completed [===========================================================>-----------------] 77% with 0 failures
Completed [===========================================================>-----------------] 78% with 0 failures
Completed [===========================================================>-----------------] 79% with 0 failures
Completed [============================================================>----------------] 79% with 0 failures
Completed [============================================================>----------------] 80% with 0 failures
Completed [=============================================================>---------------] 80% with 0 failures
Completed [=============================================================>---------------] 81% with 0 failures
Completed [==============================================================>--------------] 81% with 0 failures
Completed [==============================================================>--------------] 82% with 0 failures
Completed [===============================================================>-------------] 82% with 0 failures
Completed [===============================================================>-------------] 83% with 0 failures
Completed [===============================================================>-------------] 84% with 0 failures
Completed [================================================================>------------] 84% with 0 failures
Completed [================================================================>------------] 85% with 0 failures
Completed [=================================================================>-----------] 85% with 0 failures
Completed [=================================================================>-----------] 86% with 0 failures
Completed [==================================================================>----------] 86% with 0 failures
Completed [==================================================================>----------] 87% with 0 failures
Completed [==================================================================>----------] 88% with 0 failures
Completed [===================================================================>---------] 88% with 0 failures
Completed [===================================================================>---------] 89% with 0 failures
Completed [====================================================================>--------] 89% with 0 failures
Completed [====================================================================>--------] 90% with 0 failures
Completed [=====================================================================>-------] 90% with 0 failures
Completed [=====================================================================>-------] 91% with 0 failures
Completed [=====================================================================>-------] 92% with 0 failures
Completed [======================================================================>------] 92% with 0 failures
Completed [======================================================================>------] 93% with 0 failures
Completed [=======================================================================>-----] 93% with 0 failures
Completed [=======================================================================>-----] 94% with 0 failures
Completed [========================================================================>----] 94% with 0 failures
Completed [========================================================================>----] 95% with 0 failures
Completed [=========================================================================>---] 95% with 0 failures
Completed [=========================================================================>---] 96% with 0 failures
Completed [=========================================================================>---] 97% with 0 failures
Completed [==========================================================================>--] 97% with 0 failures
Completed [==========================================================================>--] 98% with 0 failures
Completed [===========================================================================>-] 98% with 0 failures
Completed [===========================================================================>-] 99% with 0 failures
Completed [============================================================================>] 99% with 0 failures
Completed [============================================================================>] 100% with 0 failures
Completed [=============================================================================] 100% with 0 failures
Done!
Combining coefficients and standard errors
Calculating log-fold changes
Processed 219380 groups out of 471132. 47% done. Time elapsed: 3s. ETA: 3s.
Processed 289789 groups out of 471132. 62% done. Time elapsed: 4s. ETA: 2s.
Processed 361491 groups out of 471132. 77% done. Time elapsed: 5s. ETA: 1s.
Processed 436911 groups out of 471132. 93% done. Time elapsed: 6s. ETA: 0s.
Processed 471132 groups out of 471132. 100% done. Time elapsed: 6s. ETA: 0s.Calculating likelihood ratio tests
Refitting on reduced model...
Completed [---------------------------------------------------------------------------------------] 0% with 0 failures
Completed [---------------------------------------------------------------------------------------] 1% with 0 failures
Completed [>--------------------------------------------------------------------------------------] 1% with 0 failures
Completed [>--------------------------------------------------------------------------------------] 2% with 0 failures
Completed [=>-------------------------------------------------------------------------------------] 2% with 0 failures
Completed [=>-------------------------------------------------------------------------------------] 3% with 0 failures
Completed [==>------------------------------------------------------------------------------------] 3% with 0 failures
Completed [==>------------------------------------------------------------------------------------] 4% with 0 failures
Completed [===>-----------------------------------------------------------------------------------] 4% with 0 failures
Completed [===>-----------------------------------------------------------------------------------] 5% with 0 failures
Completed [====>----------------------------------------------------------------------------------] 5% with 0 failures
Completed [====>----------------------------------------------------------------------------------] 6% with 0 failures
Completed [=====>---------------------------------------------------------------------------------] 6% with 0 failures
Completed [=====>---------------------------------------------------------------------------------] 7% with 0 failures
Completed [======>--------------------------------------------------------------------------------] 7% with 0 failures
Completed [======>--------------------------------------------------------------------------------] 8% with 0 failures
Completed [======>--------------------------------------------------------------------------------] 9% with 0 failures
Completed [=======>-------------------------------------------------------------------------------] 9% with 0 failures
Completed [=======>-------------------------------------------------------------------------------] 10% with 0 failures
Completed [========>------------------------------------------------------------------------------] 10% with 0 failures
Completed [========>------------------------------------------------------------------------------] 11% with 0 failures
Completed [=========>-----------------------------------------------------------------------------] 11% with 0 failures
Completed [=========>-----------------------------------------------------------------------------] 12% with 0 failures
Completed [==========>----------------------------------------------------------------------------] 12% with 0 failures
Completed [==========>----------------------------------------------------------------------------] 13% with 0 failures
Completed [===========>---------------------------------------------------------------------------] 13% with 0 failures
Completed [===========>---------------------------------------------------------------------------] 14% with 0 failures
Completed [============>--------------------------------------------------------------------------] 14% with 0 failures
Completed [============>--------------------------------------------------------------------------] 15% with 0 failures
Completed [============>--------------------------------------------------------------------------] 16% with 0 failures
Completed [=============>-------------------------------------------------------------------------] 16% with 0 failures
Completed [=============>-------------------------------------------------------------------------] 17% with 0 failures
Completed [==============>------------------------------------------------------------------------] 17% with 0 failures
Completed [==============>------------------------------------------------------------------------] 18% with 0 failures
Completed [===============>-----------------------------------------------------------------------] 18% with 0 failures
Completed [===============>-----------------------------------------------------------------------] 19% with 0 failures
Completed [================>----------------------------------------------------------------------] 19% with 0 failures
Completed [================>----------------------------------------------------------------------] 20% with 0 failures
Completed [=================>---------------------------------------------------------------------] 20% with 0 failures
Completed [=================>---------------------------------------------------------------------] 21% with 0 failures
Completed [==================>--------------------------------------------------------------------] 21% with 0 failures
Completed [==================>--------------------------------------------------------------------] 22% with 0 failures
Completed [===================>-------------------------------------------------------------------] 22% with 0 failures
Completed [===================>-------------------------------------------------------------------] 23% with 0 failures
Completed [===================>-------------------------------------------------------------------] 24% with 0 failures
Completed [====================>------------------------------------------------------------------] 24% with 0 failures
Completed [====================>------------------------------------------------------------------] 25% with 0 failures
Completed [=====================>-----------------------------------------------------------------] 25% with 0 failures
Completed [=====================>-----------------------------------------------------------------] 26% with 0 failures
Completed [======================>----------------------------------------------------------------] 26% with 0 failures
Completed [======================>----------------------------------------------------------------] 27% with 0 failures
Completed [=======================>---------------------------------------------------------------] 27% with 0 failures
Completed [=======================>---------------------------------------------------------------] 28% with 0 failures
Completed [========================>--------------------------------------------------------------] 28% with 0 failures
Completed [========================>--------------------------------------------------------------] 29% with 0 failures
Completed [=========================>-------------------------------------------------------------] 29% with 0 failures
Completed [=========================>-------------------------------------------------------------] 30% with 0 failures
Completed [==========================>------------------------------------------------------------] 30% with 0 failures
Completed [==========================>------------------------------------------------------------] 31% with 0 failures
Completed [==========================>------------------------------------------------------------] 32% with 0 failures
Completed [===========================>-----------------------------------------------------------] 32% with 0 failures
Completed [===========================>-----------------------------------------------------------] 33% with 0 failures
Completed [============================>----------------------------------------------------------] 33% with 0 failures
Completed [============================>----------------------------------------------------------] 34% with 0 failures
Completed [=============================>---------------------------------------------------------] 34% with 0 failures
Completed [=============================>---------------------------------------------------------] 35% with 0 failures
Completed [==============================>--------------------------------------------------------] 35% with 0 failures
Completed [==============================>--------------------------------------------------------] 36% with 0 failures
Completed [===============================>-------------------------------------------------------] 36% with 0 failures
Completed [===============================>-------------------------------------------------------] 37% with 0 failures
Completed [================================>------------------------------------------------------] 37% with 0 failures
Completed [================================>------------------------------------------------------] 38% with 0 failures
Completed [================================>------------------------------------------------------] 39% with 0 failures
Completed [=================================>-----------------------------------------------------] 39% with 0 failures
Completed [=================================>-----------------------------------------------------] 40% with 0 failures
Completed [==================================>----------------------------------------------------] 40% with 0 failures
Completed [==================================>----------------------------------------------------] 41% with 0 failures
Completed [===================================>---------------------------------------------------] 41% with 0 failures
Completed [===================================>---------------------------------------------------] 42% with 0 failures
Completed [====================================>--------------------------------------------------] 42% with 0 failures
Completed [====================================>--------------------------------------------------] 43% with 0 failures
Completed [=====================================>-------------------------------------------------] 43% with 0 failures
Completed [=====================================>-------------------------------------------------] 44% with 0 failures
Completed [======================================>------------------------------------------------] 44% with 0 failures
Completed [======================================>------------------------------------------------] 45% with 0 failures
Completed [=======================================>-----------------------------------------------] 45% with 0 failures
Completed [=======================================>-----------------------------------------------] 46% with 0 failures
Completed [=======================================>-----------------------------------------------] 47% with 0 failures
Completed [========================================>----------------------------------------------] 47% with 0 failures
Completed [========================================>----------------------------------------------] 48% with 0 failures
Completed [=========================================>---------------------------------------------] 48% with 0 failures
Completed [=========================================>---------------------------------------------] 49% with 0 failures
Completed [==========================================>--------------------------------------------] 49% with 0 failures
Completed [==========================================>--------------------------------------------] 50% with 0 failures
Completed [===========================================>-------------------------------------------] 50% with 0 failures
Completed [===========================================>-------------------------------------------] 51% with 0 failures
Completed [============================================>------------------------------------------] 51% with 0 failures
Completed [============================================>------------------------------------------] 52% with 0 failures
Completed [=============================================>-----------------------------------------] 52% with 0 failures
Completed [=============================================>-----------------------------------------] 53% with 0 failures
Completed [==============================================>----------------------------------------] 53% with 0 failures
Completed [==============================================>----------------------------------------] 54% with 0 failures
Completed [==============================================>----------------------------------------] 55% with 0 failures
Completed [===============================================>---------------------------------------] 55% with 0 failures
Completed [===============================================>---------------------------------------] 56% with 0 failures
Completed [================================================>--------------------------------------] 56% with 0 failures
Completed [================================================>--------------------------------------] 57% with 0 failures
Completed [=================================================>-------------------------------------] 57% with 0 failures
Completed [=================================================>-------------------------------------] 58% with 0 failures
Completed [==================================================>------------------------------------] 58% with 0 failures
Completed [==================================================>------------------------------------] 59% with 0 failures
Completed [===================================================>-----------------------------------] 59% with 0 failures
Completed [===================================================>-----------------------------------] 60% with 0 failures
Completed [====================================================>----------------------------------] 60% with 0 failures
Completed [====================================================>----------------------------------] 61% with 0 failures
Completed [=====================================================>---------------------------------] 61% with 0 failures
Completed [=====================================================>---------------------------------] 62% with 0 failures
Completed [=====================================================>---------------------------------] 63% with 0 failures
Completed [======================================================>--------------------------------] 63% with 0 failures
Completed [======================================================>--------------------------------] 64% with 0 failures
Completed [=======================================================>-------------------------------] 64% with 0 failures
Completed [=======================================================>-------------------------------] 65% with 0 failures
Completed [========================================================>------------------------------] 65% with 0 failures
Completed [========================================================>------------------------------] 66% with 0 failures
Completed [=========================================================>-----------------------------] 66% with 0 failures
Completed [=========================================================>-----------------------------] 67% with 0 failures
Completed [==========================================================>----------------------------] 67% with 0 failures
Completed [==========================================================>----------------------------] 68% with 0 failures
Completed [===========================================================>---------------------------] 68% with 0 failures
Completed [===========================================================>---------------------------] 69% with 0 failures
Completed [===========================================================>---------------------------] 70% with 0 failures
Completed [============================================================>--------------------------] 70% with 0 failures
Completed [============================================================>--------------------------] 71% with 0 failures
Completed [=============================================================>-------------------------] 71% with 0 failures
Completed [=============================================================>-------------------------] 72% with 0 failures
Completed [==============================================================>------------------------] 72% with 0 failures
Completed [==============================================================>------------------------] 73% with 0 failures
Completed [===============================================================>-----------------------] 73% with 0 failures
Completed [===============================================================>-----------------------] 74% with 0 failures
Completed [================================================================>----------------------] 74% with 0 failures
Completed [================================================================>----------------------] 75% with 0 failures
Completed [=================================================================>---------------------] 75% with 0 failures
Completed [=================================================================>---------------------] 76% with 0 failures
Completed [==================================================================>--------------------] 76% with 0 failures
Completed [==================================================================>--------------------] 77% with 0 failures
Completed [==================================================================>--------------------] 78% with 0 failures
Completed [===================================================================>-------------------] 78% with 0 failures
Completed [===================================================================>-------------------] 79% with 0 failures
Completed [====================================================================>------------------] 79% with 0 failures
Completed [====================================================================>------------------] 80% with 0 failures
Completed [=====================================================================>-----------------] 80% with 0 failures
Completed [=====================================================================>-----------------] 81% with 0 failures
Completed [======================================================================>----------------] 81% with 0 failures
Completed [======================================================================>----------------] 82% with 0 failures
Completed [=======================================================================>---------------] 82% with 0 failures
Completed [=======================================================================>---------------] 83% with 0 failures
Completed [========================================================================>--------------] 83% with 0 failures
Completed [========================================================================>--------------] 84% with 0 failures
Completed [=========================================================================>-------------] 84% with 0 failures
Completed [=========================================================================>-------------] 85% with 0 failures
Completed [=========================================================================>-------------] 86% with 0 failures
Completed [==========================================================================>------------] 86% with 0 failures
Completed [==========================================================================>------------] 87% with 0 failures
Completed [===========================================================================>-----------] 87% with 0 failures
Completed [===========================================================================>-----------] 88% with 0 failures
Completed [============================================================================>----------] 88% with 0 failures
Completed [============================================================================>----------] 89% with 0 failures
Completed [=============================================================================>---------] 89% with 0 failures
Completed [=============================================================================>---------] 90% with 0 failures
Completed [==============================================================================>--------] 90% with 0 failures
Completed [==============================================================================>--------] 91% with 0 failures
Completed [===============================================================================>-------] 91% with 0 failures
Completed [===============================================================================>-------] 92% with 0 failures
Completed [===============================================================================>-------] 93% with 0 failures
Completed [================================================================================>------] 93% with 0 failures
Completed [================================================================================>------] 94% with 0 failures
Completed [=================================================================================>-----] 94% with 0 failures
Completed [=================================================================================>-----] 95% with 0 failures
Completed [==================================================================================>----] 95% with 0 failures
Completed [==================================================================================>----] 96% with 0 failures
Completed [===================================================================================>---] 96% with 0 failures
Completed [===================================================================================>---] 97% with 0 failures
Completed [====================================================================================>--] 97% with 0 failures
Completed [====================================================================================>--] 98% with 0 failures
Completed [=====================================================================================>-] 98% with 0 failures
Completed [=====================================================================================>-] 99% with 0 failures
Completed [======================================================================================>] 99% with 0 failures
Completed [======================================================================================>] 100% with 0 failures
Completed [=======================================================================================] 100% with 0 failures
Done!gc() used (Mb) gc trigger (Mb) max used (Mb)
Ncells 8388915 448.1 15661334 836.5 15661334 836.5
Vcells 1232996630 9407.1 6392433025 48770.4 7990541281 60963.0gc() used (Mb) gc trigger (Mb) max used (Mb)
Ncells 8388837 448.1 15661334 836.5 15661334 836.5
Vcells 1232961244 9406.8 5113946420 39016.4 7990541281 60963.0write.csv(markers1, file = "0-imp_Robj/1-MAST_with_batch_as_Covariate.csv", row.names = TRUE)Erreur dans is.data.frame(x) : objet 'markers1' introuvable4.2 Using Wilcoxon Test (Without Batch Correction)
4.3 Using MAST Without Latent Variables (No Batch Correction)
markers_mast_no_batch <- FindMarkers(
All_samples_Merged,
ident.1 = c("0", "1", "2", "3", "4", "5", "7", "9"), # Cell lines
ident.2 = c("6", "8", "10", "11"), # CD4 T cells
assay = "SCT",
test.use = "MAST",
min.pct = 0,
logfc.threshold = 0
)Avis : sparse->dense coercion: allocating vector of size 9.6 GiBAvis : Dropping illegal slot(s) norm.method for class BayesGLMlike.
This likely indicates a bug in an upstream package.
Completed [----------------------------------------------------------] 0% with 0 failures
Completed [----------------------------------------------------------] 1% with 0 failures
Completed [>---------------------------------------------------------] 1% with 0 failures
Completed [>---------------------------------------------------------] 2% with 0 failures
Completed [>---------------------------------------------------------] 3% with 0 failures
Completed [=>--------------------------------------------------------] 3% with 0 failures
Completed [=>--------------------------------------------------------] 4% with 0 failures
Completed [==>-------------------------------------------------------] 4% with 0 failures
Completed [==>-------------------------------------------------------] 5% with 0 failures
Completed [==>-------------------------------------------------------] 6% with 0 failures
Completed [===>------------------------------------------------------] 6% with 0 failures
Completed [===>------------------------------------------------------] 7% with 0 failures
Completed [===>------------------------------------------------------] 8% with 0 failures
Completed [====>-----------------------------------------------------] 8% with 0 failures
Completed [====>-----------------------------------------------------] 9% with 0 failures
Completed [=====>----------------------------------------------------] 9% with 0 failures
Completed [=====>----------------------------------------------------] 10% with 0 failures
Completed [=====>----------------------------------------------------] 11% with 0 failures
Completed [======>---------------------------------------------------] 11% with 0 failures
Completed [======>---------------------------------------------------] 12% with 0 failures
Completed [======>---------------------------------------------------] 13% with 0 failures
Completed [=======>--------------------------------------------------] 13% with 0 failures
Completed [=======>--------------------------------------------------] 14% with 0 failures
Completed [=======>--------------------------------------------------] 15% with 0 failures
Completed [========>-------------------------------------------------] 15% with 0 failures
Completed [========>-------------------------------------------------] 16% with 0 failures
Completed [=========>------------------------------------------------] 16% with 0 failures
Completed [=========>------------------------------------------------] 17% with 0 failures
Completed [=========>------------------------------------------------] 18% with 0 failures
Completed [==========>-----------------------------------------------] 18% with 0 failures
Completed [==========>-----------------------------------------------] 19% with 0 failures
Completed [==========>-----------------------------------------------] 20% with 0 failures
Completed [===========>----------------------------------------------] 20% with 0 failures
Completed [===========>----------------------------------------------] 21% with 0 failures
Completed [===========>----------------------------------------------] 22% with 0 failures
Completed [============>---------------------------------------------] 22% with 0 failures
Completed [============>---------------------------------------------] 23% with 0 failures
Completed [=============>--------------------------------------------] 23% with 0 failures
Completed [=============>--------------------------------------------] 24% with 0 failures
Completed [=============>--------------------------------------------] 25% with 0 failures
Completed [==============>-------------------------------------------] 25% with 0 failures
Completed [==============>-------------------------------------------] 26% with 0 failures
Completed [==============>-------------------------------------------] 27% with 0 failures
Completed [===============>------------------------------------------] 27% with 0 failures
Completed [===============>------------------------------------------] 28% with 0 failures
Completed [================>-----------------------------------------] 28% with 0 failures
Completed [================>-----------------------------------------] 29% with 0 failures
Completed [================>-----------------------------------------] 30% with 0 failures
Completed [=================>----------------------------------------] 30% with 0 failures
Completed [=================>----------------------------------------] 31% with 0 failures
Completed [=================>----------------------------------------] 32% with 0 failures
Completed [==================>---------------------------------------] 32% with 0 failures
Completed [==================>---------------------------------------] 33% with 0 failures
Completed [==================>---------------------------------------] 34% with 0 failures
Completed [===================>--------------------------------------] 34% with 0 failures
Completed [===================>--------------------------------------] 35% with 0 failures
Completed [====================>-------------------------------------] 35% with 0 failures
Completed [====================>-------------------------------------] 36% with 0 failures
Completed [====================>-------------------------------------] 37% with 0 failures
Completed [=====================>------------------------------------] 37% with 0 failures
Completed [=====================>------------------------------------] 38% with 0 failures
Completed [=====================>------------------------------------] 39% with 0 failures
Completed [======================>-----------------------------------] 39% with 0 failures
Completed [======================>-----------------------------------] 40% with 0 failures
Completed [======================>-----------------------------------] 41% with 0 failures
Completed [=======================>----------------------------------] 41% with 0 failures
Completed [=======================>----------------------------------] 42% with 0 failures
Completed [========================>---------------------------------] 42% with 0 failures
Completed [========================>---------------------------------] 43% with 0 failures
Completed [========================>---------------------------------] 44% with 0 failures
Completed [=========================>--------------------------------] 44% with 0 failures
Completed [=========================>--------------------------------] 45% with 0 failures
Completed [=========================>--------------------------------] 46% with 0 failures
Completed [==========================>-------------------------------] 46% with 0 failures
Completed [==========================>-------------------------------] 47% with 0 failures
Completed [===========================>------------------------------] 47% with 0 failures
Completed [===========================>------------------------------] 48% with 0 failures
Completed [===========================>------------------------------] 49% with 0 failures
Completed [============================>-----------------------------] 49% with 0 failures
Completed [============================>-----------------------------] 50% with 0 failures
Completed [============================>-----------------------------] 51% with 0 failures
Completed [=============================>----------------------------] 51% with 0 failures
Completed [=============================>----------------------------] 52% with 0 failures
Completed [=============================>----------------------------] 53% with 0 failures
Completed [==============================>---------------------------] 53% with 0 failures
Completed [==============================>---------------------------] 54% with 0 failures
Completed [===============================>--------------------------] 54% with 0 failures
Completed [===============================>--------------------------] 55% with 0 failures
Completed [===============================>--------------------------] 56% with 0 failures
Completed [================================>-------------------------] 56% with 0 failures
Completed [================================>-------------------------] 57% with 0 failures
Completed [================================>-------------------------] 58% with 0 failures
Completed [=================================>------------------------] 58% with 0 failures
Completed [=================================>------------------------] 59% with 0 failures
Completed [==================================>-----------------------] 59% with 0 failures
Completed [==================================>-----------------------] 60% with 0 failures
Completed [==================================>-----------------------] 61% with 0 failures
Completed [===================================>----------------------] 61% with 0 failures
Completed [===================================>----------------------] 62% with 0 failures
Completed [===================================>----------------------] 63% with 0 failures
Completed [====================================>---------------------] 63% with 0 failures
Completed [====================================>---------------------] 64% with 0 failures
Completed [====================================>---------------------] 65% with 0 failures
Completed [=====================================>--------------------] 65% with 0 failures
Completed [=====================================>--------------------] 66% with 0 failures
Completed [======================================>-------------------] 66% with 0 failures
Completed [======================================>-------------------] 67% with 0 failures
Completed [======================================>-------------------] 68% with 0 failures
Completed [=======================================>------------------] 68% with 0 failures
Completed [=======================================>------------------] 69% with 0 failures
Completed [=======================================>------------------] 70% with 0 failures
Completed [========================================>-----------------] 70% with 0 failures
Completed [========================================>-----------------] 71% with 0 failures
Completed [========================================>-----------------] 72% with 0 failures
Completed [=========================================>----------------] 72% with 0 failures
Completed [=========================================>----------------] 73% with 0 failures
Completed [==========================================>---------------] 73% with 0 failures
Completed [==========================================>---------------] 74% with 0 failures
Completed [==========================================>---------------] 75% with 0 failures
Completed [===========================================>--------------] 75% with 0 failures
Completed [===========================================>--------------] 76% with 0 failures
Completed [===========================================>--------------] 77% with 0 failures
Completed [============================================>-------------] 77% with 0 failures
Completed [============================================>-------------] 78% with 0 failures
Completed [=============================================>------------] 78% with 0 failures
Completed [=============================================>------------] 79% with 0 failures
Completed [=============================================>------------] 80% with 0 failures
Completed [==============================================>-----------] 80% with 0 failures
Completed [==============================================>-----------] 81% with 0 failures
Completed [==============================================>-----------] 82% with 0 failures
Completed [===============================================>----------] 82% with 0 failures
Completed [===============================================>----------] 83% with 0 failures
Completed [===============================================>----------] 84% with 0 failures
Completed [================================================>---------] 84% with 0 failures
Completed [================================================>---------] 85% with 0 failures
Completed [=================================================>--------] 85% with 0 failures
Completed [=================================================>--------] 86% with 0 failures
Completed [=================================================>--------] 87% with 0 failures
Completed [==================================================>-------] 87% with 0 failures
Completed [==================================================>-------] 88% with 0 failures
Completed [==================================================>-------] 89% with 0 failures
Completed [===================================================>------] 89% with 0 failures
Completed [===================================================>------] 90% with 0 failures
Completed [===================================================>------] 91% with 0 failures
Completed [====================================================>-----] 91% with 0 failures
Completed [====================================================>-----] 92% with 0 failures
Completed [=====================================================>----] 92% with 0 failures
Completed [=====================================================>----] 93% with 0 failures
Completed [=====================================================>----] 94% with 0 failures
Completed [======================================================>---] 94% with 0 failures
Completed [======================================================>---] 95% with 0 failures
Completed [======================================================>---] 96% with 0 failures
Completed [=======================================================>--] 96% with 0 failures
Completed [=======================================================>--] 97% with 0 failures
Completed [========================================================>-] 97% with 0 failures
Completed [========================================================>-] 98% with 0 failures
Completed [========================================================>-] 99% with 0 failures
Completed [=========================================================>] 99% with 0 failures
Completed [=========================================================>] 100% with 0 failures
Completed [==========================================================] 100% with 0 failures
Done!
Combining coefficients and standard errors
Calculating log-fold changes
Calculating likelihood ratio tests
Refitting on reduced model...
Completed [----------------------------------------------------------] 0% with 0 failures
Completed [----------------------------------------------------------] 1% with 0 failures
Completed [>---------------------------------------------------------] 1% with 0 failures
Completed [>---------------------------------------------------------] 2% with 0 failures
Completed [>---------------------------------------------------------] 3% with 0 failures
Completed [=>--------------------------------------------------------] 3% with 0 failures
Completed [=>--------------------------------------------------------] 4% with 0 failures
Completed [==>-------------------------------------------------------] 4% with 0 failures
Completed [==>-------------------------------------------------------] 5% with 0 failures
Completed [==>-------------------------------------------------------] 6% with 0 failures
Completed [===>------------------------------------------------------] 6% with 0 failures
Completed [===>------------------------------------------------------] 7% with 0 failures
Completed [===>------------------------------------------------------] 8% with 0 failures
Completed [====>-----------------------------------------------------] 8% with 0 failures
Completed [====>-----------------------------------------------------] 9% with 0 failures
Completed [=====>----------------------------------------------------] 9% with 0 failures
Completed [=====>----------------------------------------------------] 10% with 0 failures
Completed [=====>----------------------------------------------------] 11% with 0 failures
Completed [======>---------------------------------------------------] 11% with 0 failures
Completed [======>---------------------------------------------------] 12% with 0 failures
Completed [======>---------------------------------------------------] 13% with 0 failures
Completed [=======>--------------------------------------------------] 13% with 0 failures
Completed [=======>--------------------------------------------------] 14% with 0 failures
Completed [=======>--------------------------------------------------] 15% with 0 failures
Completed [========>-------------------------------------------------] 15% with 0 failures
Completed [========>-------------------------------------------------] 16% with 0 failures
Completed [=========>------------------------------------------------] 16% with 0 failures
Completed [=========>------------------------------------------------] 17% with 0 failures
Completed [=========>------------------------------------------------] 18% with 0 failures
Completed [==========>-----------------------------------------------] 18% with 0 failures
Completed [==========>-----------------------------------------------] 19% with 0 failures
Completed [==========>-----------------------------------------------] 20% with 0 failures
Completed [===========>----------------------------------------------] 20% with 0 failures
Completed [===========>----------------------------------------------] 21% with 0 failures
Completed [===========>----------------------------------------------] 22% with 0 failures
Completed [============>---------------------------------------------] 22% with 0 failures
Completed [============>---------------------------------------------] 23% with 0 failures
Completed [=============>--------------------------------------------] 23% with 0 failures
Completed [=============>--------------------------------------------] 24% with 0 failures
Completed [=============>--------------------------------------------] 25% with 0 failures
Completed [==============>-------------------------------------------] 25% with 0 failures
Completed [==============>-------------------------------------------] 26% with 0 failures
Completed [==============>-------------------------------------------] 27% with 0 failures
Completed [===============>------------------------------------------] 27% with 0 failures
Completed [===============>------------------------------------------] 28% with 0 failures
Completed [================>-----------------------------------------] 28% with 0 failures
Completed [================>-----------------------------------------] 29% with 0 failures
Completed [================>-----------------------------------------] 30% with 0 failures
Completed [=================>----------------------------------------] 30% with 0 failures
Completed [=================>----------------------------------------] 31% with 0 failures
Completed [=================>----------------------------------------] 32% with 0 failures
Completed [==================>---------------------------------------] 32% with 0 failures
Completed [==================>---------------------------------------] 33% with 0 failures
Completed [==================>---------------------------------------] 34% with 0 failures
Completed [===================>--------------------------------------] 34% with 0 failures
Completed [===================>--------------------------------------] 35% with 0 failures
Completed [====================>-------------------------------------] 35% with 0 failures
Completed [====================>-------------------------------------] 36% with 0 failures
Completed [====================>-------------------------------------] 37% with 0 failures
Completed [=====================>------------------------------------] 37% with 0 failures
Completed [=====================>------------------------------------] 38% with 0 failures
Completed [=====================>------------------------------------] 39% with 0 failures
Completed [======================>-----------------------------------] 39% with 0 failures
Completed [======================>-----------------------------------] 40% with 0 failures
Completed [======================>-----------------------------------] 41% with 0 failures
Completed [=======================>----------------------------------] 41% with 0 failures
Completed [=======================>----------------------------------] 42% with 0 failures
Completed [========================>---------------------------------] 42% with 0 failures
Completed [========================>---------------------------------] 43% with 0 failures
Completed [========================>---------------------------------] 44% with 0 failures
Completed [=========================>--------------------------------] 44% with 0 failures
Completed [=========================>--------------------------------] 45% with 0 failures
Completed [=========================>--------------------------------] 46% with 0 failures
Completed [==========================>-------------------------------] 46% with 0 failures
Completed [==========================>-------------------------------] 47% with 0 failures
Completed [===========================>------------------------------] 47% with 0 failures
Completed [===========================>------------------------------] 48% with 0 failures
Completed [===========================>------------------------------] 49% with 0 failures
Completed [============================>-----------------------------] 49% with 0 failures
Completed [============================>-----------------------------] 50% with 0 failures
Completed [============================>-----------------------------] 51% with 0 failures
Completed [=============================>----------------------------] 51% with 0 failures
Completed [=============================>----------------------------] 52% with 0 failures
Completed [=============================>----------------------------] 53% with 0 failures
Completed [==============================>---------------------------] 53% with 0 failures
Completed [==============================>---------------------------] 54% with 0 failures
Completed [===============================>--------------------------] 54% with 0 failures
Completed [===============================>--------------------------] 55% with 0 failures
Completed [===============================>--------------------------] 56% with 0 failures
Completed [================================>-------------------------] 56% with 0 failures
Completed [================================>-------------------------] 57% with 0 failures
Completed [================================>-------------------------] 58% with 0 failures
Completed [=================================>------------------------] 58% with 0 failures
Completed [=================================>------------------------] 59% with 0 failures
Completed [==================================>-----------------------] 59% with 0 failures
Completed [==================================>-----------------------] 60% with 0 failures
Completed [==================================>-----------------------] 61% with 0 failures
Completed [===================================>----------------------] 61% with 0 failures
Completed [===================================>----------------------] 62% with 0 failures
Completed [===================================>----------------------] 63% with 0 failures
Completed [====================================>---------------------] 63% with 0 failures
Completed [====================================>---------------------] 64% with 0 failures
Completed [====================================>---------------------] 65% with 0 failures
Completed [=====================================>--------------------] 65% with 0 failures
Completed [=====================================>--------------------] 66% with 0 failures
Completed [======================================>-------------------] 66% with 0 failures
Completed [======================================>-------------------] 67% with 0 failures
Completed [======================================>-------------------] 68% with 0 failures
Completed [=======================================>------------------] 68% with 0 failures
Completed [=======================================>------------------] 69% with 0 failures
Completed [=======================================>------------------] 70% with 0 failures
Completed [========================================>-----------------] 70% with 0 failures
Completed [========================================>-----------------] 71% with 0 failures
Completed [========================================>-----------------] 72% with 0 failures
Completed [=========================================>----------------] 72% with 0 failures
Completed [=========================================>----------------] 73% with 0 failures
Completed [==========================================>---------------] 73% with 0 failures
Completed [==========================================>---------------] 74% with 0 failures
Completed [==========================================>---------------] 75% with 0 failures
Completed [===========================================>--------------] 75% with 0 failures
Completed [===========================================>--------------] 76% with 0 failures
Completed [===========================================>--------------] 77% with 0 failures
Completed [============================================>-------------] 77% with 0 failures
Completed [============================================>-------------] 78% with 0 failures
Completed [=============================================>------------] 78% with 0 failures
Completed [=============================================>------------] 79% with 0 failures
Completed [=============================================>------------] 80% with 0 failures
Completed [==============================================>-----------] 80% with 0 failures
Completed [==============================================>-----------] 81% with 0 failures
Completed [==============================================>-----------] 82% with 0 failures
Completed [===============================================>----------] 82% with 0 failures
Completed [===============================================>----------] 83% with 0 failures
Completed [===============================================>----------] 84% with 0 failures
Completed [================================================>---------] 84% with 0 failures
Completed [================================================>---------] 85% with 0 failures
Completed [=================================================>--------] 85% with 0 failures
Completed [=================================================>--------] 86% with 0 failures
Completed [=================================================>--------] 87% with 0 failures
Completed [==================================================>-------] 87% with 0 failures
Completed [==================================================>-------] 88% with 0 failures
Completed [==================================================>-------] 89% with 0 failures
Completed [===================================================>------] 89% with 0 failures
Completed [===================================================>------] 90% with 0 failures
Completed [===================================================>------] 91% with 0 failures
Completed [====================================================>-----] 91% with 0 failures
Completed [====================================================>-----] 92% with 0 failures
Completed [=====================================================>----] 92% with 0 failures
Completed [=====================================================>----] 93% with 0 failures
Completed [=====================================================>----] 94% with 0 failures
Completed [======================================================>---] 94% with 0 failures
Completed [======================================================>---] 95% with 0 failures
Completed [======================================================>---] 96% with 0 failures
Completed [=======================================================>--] 96% with 0 failures
Completed [=======================================================>--] 97% with 0 failures
Completed [========================================================>-] 97% with 0 failures
Completed [========================================================>-] 98% with 0 failures
Completed [========================================================>-] 99% with 0 failures
Completed [=========================================================>] 99% with 0 failures
Completed [=========================================================>] 100% with 0 failures
Completed [==========================================================] 100% with 0 failures
Done!View(markers_mast_batch)
View(markers_wilcox)
View(markers_wilcox)
View(markers_wilcox)
write.csv(markers_mast_no_batch, file = "0-imp_Robj/3-MAST_without_batch_as_Covariate.csv", row.names = TRUE)5. Compute Mean Expression for Groups
# Define cell groups
group1_cells <- WhichCells(All_samples_Merged, idents = c("0", "1", "2", "3", "4", "5", "7", "9"))
group2_cells <- WhichCells(All_samples_Merged, idents = c("6", "8", "10", "11"))
# Extract normalized expression values
expression_data <- GetAssayData(All_samples_Merged, slot = "data") # Log-normalized valuesAvis : The `slot` argument of `GetAssayData()` is deprecated as of SeuratObject 5.0.0.
Please use the `layer` argument instead.# Calculate mean expression for each group
calculate_mean_expression <- function(markers, group1_cells, group2_cells, expression_data) {
group1_mean <- rowMeans(expression_data[, group1_cells])
group2_mean <- rowMeans(expression_data[, group2_cells])
markers <- markers %>%
rownames_to_column("gene") %>%
mutate(mean_expr_group1 = group1_mean[gene],
mean_expr_group2 = group2_mean[gene])
return(markers)
}
markers_mast_batch <- calculate_mean_expression(markers_mast_batch, group1_cells, group2_cells, expression_data)
markers_wilcox <- calculate_mean_expression(markers_wilcox, group1_cells, group2_cells, expression_data)
markers_mast_no_batch <- calculate_mean_expression(markers_mast_no_batch, group1_cells, group2_cells, expression_data)
# Save updated results with mean expression
write.csv(markers_mast_batch, file = "0-imp_Robj/1-MAST_with_batch_as_Covariate_with_meanExpression.csv", row.names = TRUE)
write.csv(markers_wilcox, file = "0-imp_Robj/2-Wilcox_min.pct_logfcT-0_with_meanExpression.csv", row.names = TRUE)
write.csv(markers_mast_no_batch, file = "0-imp_Robj/3-MAST_without_batch_as_Covariate_with_meanExpression.csv", row.names = TRUE)6. Summarize Results
6.1 Summary Function
summarize_markers <- function(markers) {
num_pval0 <- sum(markers$p_val_adj == 0)
num_pval1 <- sum(markers$p_val_adj == 1)
num_significant <- sum(markers$p_val_adj < 0.05)
num_upregulated <- sum(markers$avg_log2FC > 1)
num_downregulated <- sum(markers$avg_log2FC < -1)
cat("Number of genes with p_val_adj = 0:", num_pval0, "\n")
cat("Number of genes with p_val_adj = 1:", num_pval1, "\n")
cat("Number of significant genes (p_val_adj < 0.05):", num_significant, "\n")
cat("Number of upregulated genes (avg_log2FC > 1):", num_upregulated, "\n")
cat("Number of downregulated genes (avg_log2FC < -1):", num_downregulated, "\n")
}
6.2 Summarize Markers1 (MAST with Batch Correction)
cat("Markers1 Summary (MAST with Batch Correction):\n")Markers1 Summary (MAST with Batch Correction):summarize_markers(markers_mast_batch)Number of genes with p_val_adj = 0: 1
Number of genes with p_val_adj = 1: 16799
Number of significant genes (p_val_adj < 0.05): 8601
Number of upregulated genes (avg_log2FC > 1): 14796
Number of downregulated genes (avg_log2FC < -1): 2494 6.3 Summary for Markers2 (Wilcoxon Test)
cat("Markers2 Summary (Wilcoxon Test):\n")Markers2 Summary (Wilcoxon Test):summarize_markers(markers_wilcox)Number of genes with p_val_adj = 0: 4502
Number of genes with p_val_adj = 1: 10049
Number of significant genes (p_val_adj < 0.05): 15357
Number of upregulated genes (avg_log2FC > 1): 14796
Number of downregulated genes (avg_log2FC < -1): 2494 6.4 Summary for Markers3 (MAST without Batch Correction)
cat("Markers3 Summary (MAST without Batch Correction):\n")Markers3 Summary (MAST without Batch Correction):summarize_markers(markers_mast_no_batch)Number of genes with p_val_adj = 0: 4957
Number of genes with p_val_adj = 1: 9570
Number of significant genes (p_val_adj < 0.05): 15882
Number of upregulated genes (avg_log2FC > 1): 14796
Number of downregulated genes (avg_log2FC < -1): 2494 7. Visualization
7.1 Volcano Plot for MAST with Batch Correction
EnhancedVolcano(markers_mast_batch,
lab = markers_mast_batch$gene,
x = "avg_log2FC",
y = "p_val_adj",
title = "MAST with Batch Correction",
pCutoff = 0.05,
FCcutoff = 1.0)Avis : One or more p-values is 0. Converting to 10^-1 * current lowest non-zero p-value...
7.2 Volcano Plot for Wilcoxon Test
EnhancedVolcano(markers_wilcox,
lab = markers_wilcox$gene,
x = "avg_log2FC",
y = "p_val_adj",
title = "MAST with Batch Correction",
pCutoff = 0.05,
FCcutoff = 1.0)Avis : One or more p-values is 0. Converting to 10^-1 * current lowest non-zero p-value...
7.3 Volcano Plot for MAST without Batch Correction
EnhancedVolcano(markers_mast_no_batch,
lab = markers_mast_no_batch$gene,
x = "avg_log2FC",
y = "p_val_adj",
title = "MAST with Batch Correction",
pCutoff = 0.05,
FCcutoff = 1.0)Avis : One or more p-values is 0. Converting to 10^-1 * current lowest non-zero p-value...
7.4 Heatmap of Top Differentially Expressed Genes
# Load the pheatmap package
library(pheatmap)
# Define cell groups
group1_cells <- WhichCells(All_samples_Merged, idents = c("0", "1", "2", "3", "4", "5", "7", "9"))
group2_cells <- WhichCells(All_samples_Merged, idents = c("6", "8", "10", "11"))
# Define top genes based on adjusted p-value
top_genes <- head(markers_mast_batch$gene[order(markers_mast_batch$p_val_adj)], 20)
# Extract the corresponding expression data
expression_data <- GetAssayData(All_samples_Merged, slot = "data") # Log-normalized values
heatmap_data <- expression_data[top_genes, c(group1_cells, group2_cells)]
# Create annotation data
annotation_data <- data.frame(
CellLine = factor(All_samples_Merged$cell_line)[c(group1_cells, group2_cells)],
Cluster = factor(All_samples_Merged$seurat_clusters[c(group1_cells, group2_cells)])
)
# Set row names for annotation data to match column names of heatmap_data
rownames(annotation_data) <- colnames(heatmap_data)
# Generate the heatmap
pheatmap(heatmap_data,
cluster_rows = TRUE,
cluster_cols = TRUE,
show_rownames = TRUE,
show_colnames = FALSE, # Hide cell labels below the heatmap
annotation_col = annotation_data)
LS0tCnRpdGxlOiAiRGlmZmVyZW50aWFsIEV4cHJlc3Npb24gQW5hbHlzaXMgdXNpbmcgSGFybW9ueSBJbnRlZ3JhdGVkIENsdXN0ZXJzIgphdXRob3I6IE5hc2lyIE1haG1vb2QgQWJiYXNpCmRhdGU6ICJgciBTeXMuRGF0ZSgpYCIKb3V0cHV0OgogIGh0bWxfbm90ZWJvb2s6CiAgICB0b2M6IHRydWUKICAgIHRvY19mbG9hdDogdHJ1ZQogICAgdG9jX2NvbGxhcHNlZDogdHJ1ZQotLS0KCiMjIDEuIExvYWQgTGlicmFyaWVzCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFfQpsaWJyYXJ5KFNldXJhdCkKbGlicmFyeShoYXJtb255KQpsaWJyYXJ5KGRwbHlyKQpsaWJyYXJ5KEVuaGFuY2VkVm9sY2FubykKbGlicmFyeShnZ3Bsb3QyKQpsaWJyYXJ5KHBoZWF0bWFwKQpsaWJyYXJ5KHRpYmJsZSkKYGBgCgoKIyAyLiBMb2FkIFNldXJhdCBPYmplY3QgCmBgYHtyIGxvYWRfc2V1cmF0fQoKI0xvYWQgU2V1cmF0IE9iamVjdCBtZXJnZWQgZnJvbSBjZWxsIGxpbmVzIGFuZCBhIGNvbnRyb2woUEJNQykgYWZ0ZXIgZmlsdHJhdGlvbgpsb2FkKCIwLWltcF9Sb2JqL0hhcm1vbnlfaW50ZWdyYXRlZF9DRDRUY2VsbHNfaGFybW9ueV9pbnRlZ3JhdGVkLlJvYmoiKQoKQWxsX3NhbXBsZXNfTWVyZ2VkIAoKCgpgYGAKCiMgMy4gU2V0IFVwIElkZW50aWZpZXJzIGZvciBDbHVzdGVyaW5nCmBgYHtyfQojIEFzc2lnbiBjbHVzdGVyIGlkZW50aXRpZXMgdG8gdGhlIFNldXJhdCBvYmplY3QKSWRlbnRzKEFsbF9zYW1wbGVzX01lcmdlZCkgPC0gInNldXJhdF9jbHVzdGVycyIKCgpEaW1QbG90KEFsbF9zYW1wbGVzX01lcmdlZCwgcmVkdWN0aW9uID0gInVtYXAuaGFybW9ueSIsIGdyb3VwLmJ5ID0gInNldXJhdF9jbHVzdGVycyIsbGFiZWwgPSBULCBsYWJlbC5ib3ggPSBULCByZXBlbCA9IFQpICsgCiAgZ2d0aXRsZSgiSGFybW9ueSBJbnRlZ3JhdGlvbiAtIEJ5IENsdXN0ZXJzIikKYGBgCgojIDQuICBEaWZmZXJlbnRpYWwgRXhwcmVzc2lvbiBBbmFseXNpcwoKIyMgNC4xIFVzaW5nIE1BU1Qgd2l0aCBMYXRlbnQgVmFyaWFibGVzIChCYXRjaCBDb3JyZWN0aW9uKQpgYGB7cn0KbWFya2Vyc19tYXN0X2JhdGNoIDwtIEZpbmRNYXJrZXJzKAogIEFsbF9zYW1wbGVzX01lcmdlZCwKICBpZGVudC4xID0gYygiMCIsICIxIiwgIjIiLCAiMyIsICI0IiwgIjUiLCAiNyIsICI5IiksICMgQ2VsbCBsaW5lcwogIGlkZW50LjIgPSBjKCI2IiwgIjgiLCAiMTAiLCAiMTEiKSwgIyBDRDQgVCBjZWxscwogIGFzc2F5ID0gIlNDVCIsCiAgdGVzdC51c2UgPSAiTUFTVCIsCiAgbGF0ZW50LnZhcnMgPSBjKCJjZWxsX2xpbmUiKSwgIyBBZGp1c3QgZm9yIGJhdGNoIGVmZmVjdHMKICBtaW4ucGN0ID0gMCwKICBsb2dmYy50aHJlc2hvbGQgPSAwCikKd3JpdGUuY3N2KG1hcmtlcnNfbWFzdF9iYXRjaCwgZmlsZSA9ICIwLWltcF9Sb2JqLzEtTUFTVF93aXRoX2JhdGNoX2FzX0NvdmFyaWF0ZS5jc3YiLCByb3cubmFtZXMgPSBUUlVFKQpgYGAKCiMjIDQuMiBVc2luZyBXaWxjb3hvbiBUZXN0IChXaXRob3V0IEJhdGNoIENvcnJlY3Rpb24pCmBgYHtyfQptYXJrZXJzX3dpbGNveCA8LSBGaW5kTWFya2VycygKICBBbGxfc2FtcGxlc19NZXJnZWQsCiAgaWRlbnQuMSA9IGMoIjAiLCAiMSIsICIyIiwgIjMiLCAiNCIsICI1IiwgIjciLCAiOSIpLCAjIENlbGwgbGluZXMKICBpZGVudC4yID0gYygiNiIsICI4IiwgIjEwIiwgIjExIiksICMgQ0Q0IFQgY2VsbHMKICBhc3NheSA9ICJTQ1QiLAogIHRlc3QudXNlID0gIndpbGNveCIsCiAgbWluLnBjdCA9IDAsCiAgbG9nZmMudGhyZXNob2xkID0gMAopCndyaXRlLmNzdihtYXJrZXJzX3dpbGNveCwgZmlsZSA9ICIwLWltcF9Sb2JqLzItV2lsY294X21pbi5wY3RfbG9nZmNULTAuY3N2Iiwgcm93Lm5hbWVzID0gVFJVRSkKYGBgCgojIyA0LjMgVXNpbmcgTUFTVCBXaXRob3V0IExhdGVudCBWYXJpYWJsZXMgKE5vIEJhdGNoIENvcnJlY3Rpb24pCmBgYHtyfQptYXJrZXJzX21hc3Rfbm9fYmF0Y2ggPC0gRmluZE1hcmtlcnMoCiAgQWxsX3NhbXBsZXNfTWVyZ2VkLAogIGlkZW50LjEgPSBjKCIwIiwgIjEiLCAiMiIsICIzIiwgIjQiLCAiNSIsICI3IiwgIjkiKSwgIyBDZWxsIGxpbmVzCiAgaWRlbnQuMiA9IGMoIjYiLCAiOCIsICIxMCIsICIxMSIpLCAjIENENCBUIGNlbGxzCiAgYXNzYXkgPSAiU0NUIiwKICB0ZXN0LnVzZSA9ICJNQVNUIiwKICBtaW4ucGN0ID0gMCwKICBsb2dmYy50aHJlc2hvbGQgPSAwCikKd3JpdGUuY3N2KG1hcmtlcnNfbWFzdF9ub19iYXRjaCwgZmlsZSA9ICIwLWltcF9Sb2JqLzMtTUFTVF93aXRob3V0X2JhdGNoX2FzX0NvdmFyaWF0ZS5jc3YiLCByb3cubmFtZXMgPSBUUlVFKQpgYGAKCiMgNS4gQ29tcHV0ZSBNZWFuIEV4cHJlc3Npb24gZm9yIEdyb3VwcwpgYGB7cn0KIyBEZWZpbmUgY2VsbCBncm91cHMKZ3JvdXAxX2NlbGxzIDwtIFdoaWNoQ2VsbHMoQWxsX3NhbXBsZXNfTWVyZ2VkLCBpZGVudHMgPSBjKCIwIiwgIjEiLCAiMiIsICIzIiwgIjQiLCAiNSIsICI3IiwgIjkiKSkKZ3JvdXAyX2NlbGxzIDwtIFdoaWNoQ2VsbHMoQWxsX3NhbXBsZXNfTWVyZ2VkLCBpZGVudHMgPSBjKCI2IiwgIjgiLCAiMTAiLCAiMTEiKSkKCiMgRXh0cmFjdCBub3JtYWxpemVkIGV4cHJlc3Npb24gdmFsdWVzCmV4cHJlc3Npb25fZGF0YSA8LSBHZXRBc3NheURhdGEoQWxsX3NhbXBsZXNfTWVyZ2VkLCBzbG90ID0gImRhdGEiKSAgIyBMb2ctbm9ybWFsaXplZCB2YWx1ZXMKCiMgQ2FsY3VsYXRlIG1lYW4gZXhwcmVzc2lvbiBmb3IgZWFjaCBncm91cApjYWxjdWxhdGVfbWVhbl9leHByZXNzaW9uIDwtIGZ1bmN0aW9uKG1hcmtlcnMsIGdyb3VwMV9jZWxscywgZ3JvdXAyX2NlbGxzLCBleHByZXNzaW9uX2RhdGEpIHsKICBncm91cDFfbWVhbiA8LSByb3dNZWFucyhleHByZXNzaW9uX2RhdGFbLCBncm91cDFfY2VsbHNdKQogIGdyb3VwMl9tZWFuIDwtIHJvd01lYW5zKGV4cHJlc3Npb25fZGF0YVssIGdyb3VwMl9jZWxsc10pCiAgCiAgbWFya2VycyA8LSBtYXJrZXJzICU+JQogICAgcm93bmFtZXNfdG9fY29sdW1uKCJnZW5lIikgJT4lCiAgICBtdXRhdGUobWVhbl9leHByX2dyb3VwMSA9IGdyb3VwMV9tZWFuW2dlbmVdLAogICAgICAgICAgIG1lYW5fZXhwcl9ncm91cDIgPSBncm91cDJfbWVhbltnZW5lXSkKICAKICByZXR1cm4obWFya2VycykKfQoKbWFya2Vyc19tYXN0X2JhdGNoIDwtIGNhbGN1bGF0ZV9tZWFuX2V4cHJlc3Npb24obWFya2Vyc19tYXN0X2JhdGNoLCBncm91cDFfY2VsbHMsIGdyb3VwMl9jZWxscywgZXhwcmVzc2lvbl9kYXRhKQptYXJrZXJzX3dpbGNveCA8LSBjYWxjdWxhdGVfbWVhbl9leHByZXNzaW9uKG1hcmtlcnNfd2lsY294LCBncm91cDFfY2VsbHMsIGdyb3VwMl9jZWxscywgZXhwcmVzc2lvbl9kYXRhKQptYXJrZXJzX21hc3Rfbm9fYmF0Y2ggPC0gY2FsY3VsYXRlX21lYW5fZXhwcmVzc2lvbihtYXJrZXJzX21hc3Rfbm9fYmF0Y2gsIGdyb3VwMV9jZWxscywgZ3JvdXAyX2NlbGxzLCBleHByZXNzaW9uX2RhdGEpCgojIFNhdmUgdXBkYXRlZCByZXN1bHRzIHdpdGggbWVhbiBleHByZXNzaW9uCndyaXRlLmNzdihtYXJrZXJzX21hc3RfYmF0Y2gsIGZpbGUgPSAiMC1pbXBfUm9iai8xLU1BU1Rfd2l0aF9iYXRjaF9hc19Db3ZhcmlhdGVfd2l0aF9tZWFuRXhwcmVzc2lvbi5jc3YiLCByb3cubmFtZXMgPSBUUlVFKQp3cml0ZS5jc3YobWFya2Vyc193aWxjb3gsIGZpbGUgPSAiMC1pbXBfUm9iai8yLVdpbGNveF9taW4ucGN0X2xvZ2ZjVC0wX3dpdGhfbWVhbkV4cHJlc3Npb24uY3N2Iiwgcm93Lm5hbWVzID0gVFJVRSkKd3JpdGUuY3N2KG1hcmtlcnNfbWFzdF9ub19iYXRjaCwgZmlsZSA9ICIwLWltcF9Sb2JqLzMtTUFTVF93aXRob3V0X2JhdGNoX2FzX0NvdmFyaWF0ZV93aXRoX21lYW5FeHByZXNzaW9uLmNzdiIsIHJvdy5uYW1lcyA9IFRSVUUpCmBgYAoKIyA2LiBTdW1tYXJpemUgUmVzdWx0cwoKIyMgNi4xIFN1bW1hcnkgRnVuY3Rpb24KYGBge3J9CnN1bW1hcml6ZV9tYXJrZXJzIDwtIGZ1bmN0aW9uKG1hcmtlcnMpIHsKICBudW1fcHZhbDAgPC0gc3VtKG1hcmtlcnMkcF92YWxfYWRqID09IDApCiAgbnVtX3B2YWwxIDwtIHN1bShtYXJrZXJzJHBfdmFsX2FkaiA9PSAxKQogIG51bV9zaWduaWZpY2FudCA8LSBzdW0obWFya2VycyRwX3ZhbF9hZGogPCAwLjA1KQogIG51bV91cHJlZ3VsYXRlZCA8LSBzdW0obWFya2VycyRhdmdfbG9nMkZDID4gMSkKICBudW1fZG93bnJlZ3VsYXRlZCA8LSBzdW0obWFya2VycyRhdmdfbG9nMkZDIDwgLTEpCiAgCiAgY2F0KCJOdW1iZXIgb2YgZ2VuZXMgd2l0aCBwX3ZhbF9hZGogPSAwOiIsIG51bV9wdmFsMCwgIlxuIikKICBjYXQoIk51bWJlciBvZiBnZW5lcyB3aXRoIHBfdmFsX2FkaiA9IDE6IiwgbnVtX3B2YWwxLCAiXG4iKQogIGNhdCgiTnVtYmVyIG9mIHNpZ25pZmljYW50IGdlbmVzIChwX3ZhbF9hZGogPCAwLjA1KToiLCBudW1fc2lnbmlmaWNhbnQsICJcbiIpCiAgY2F0KCJOdW1iZXIgb2YgdXByZWd1bGF0ZWQgZ2VuZXMgKGF2Z19sb2cyRkMgPiAxKToiLCBudW1fdXByZWd1bGF0ZWQsICJcbiIpCiAgY2F0KCJOdW1iZXIgb2YgZG93bnJlZ3VsYXRlZCBnZW5lcyAoYXZnX2xvZzJGQyA8IC0xKToiLCBudW1fZG93bnJlZ3VsYXRlZCwgIlxuIikKfQoKYGBgCgojIyA2LjIgU3VtbWFyaXplIE1hcmtlcnMxIChNQVNUIHdpdGggQmF0Y2ggQ29ycmVjdGlvbikKYGBge3J9CmNhdCgiTWFya2VyczEgU3VtbWFyeSAoTUFTVCB3aXRoIEJhdGNoIENvcnJlY3Rpb24pOlxuIikKc3VtbWFyaXplX21hcmtlcnMobWFya2Vyc19tYXN0X2JhdGNoKQoKYGBgCgojIyA2LjMgU3VtbWFyeSBmb3IgTWFya2VyczIgKFdpbGNveG9uIFRlc3QpCmBgYHtyfQpjYXQoIk1hcmtlcnMyIFN1bW1hcnkgKFdpbGNveG9uIFRlc3QpOlxuIikKc3VtbWFyaXplX21hcmtlcnMobWFya2Vyc193aWxjb3gpCmBgYAoKIyMgNi40IFN1bW1hcnkgZm9yIE1hcmtlcnMzIChNQVNUIHdpdGhvdXQgQmF0Y2ggQ29ycmVjdGlvbikKYGBge3J9CmNhdCgiTWFya2VyczMgU3VtbWFyeSAoTUFTVCB3aXRob3V0IEJhdGNoIENvcnJlY3Rpb24pOlxuIikKc3VtbWFyaXplX21hcmtlcnMobWFya2Vyc19tYXN0X25vX2JhdGNoKQpgYGAKCgojIDcuIFZpc3VhbGl6YXRpb24KCiMjIDcuMSBWb2xjYW5vIFBsb3QgZm9yIE1BU1Qgd2l0aCBCYXRjaCBDb3JyZWN0aW9uCmBgYHtyLCBmaWcuaGVpZ2h0PTEyLCBmaWcud2lkdGg9MTJ9CkVuaGFuY2VkVm9sY2FubyhtYXJrZXJzX21hc3RfYmF0Y2gsCiAgICAgICAgICAgICAgICBsYWIgPSBtYXJrZXJzX21hc3RfYmF0Y2gkZ2VuZSwKICAgICAgICAgICAgICAgIHggPSAiYXZnX2xvZzJGQyIsCiAgICAgICAgICAgICAgICB5ID0gInBfdmFsX2FkaiIsCiAgICAgICAgICAgICAgICB0aXRsZSA9ICJNQVNUIHdpdGggQmF0Y2ggQ29ycmVjdGlvbiIsCiAgICAgICAgICAgICAgICBwQ3V0b2ZmID0gMC4wNSwKICAgICAgICAgICAgICAgIEZDY3V0b2ZmID0gMS4wKQoKYGBgCgojIyA3LjIgVm9sY2FubyBQbG90IGZvciBXaWxjb3hvbiBUZXN0CmBgYHtyLCBmaWcuaGVpZ2h0PTEyLCBmaWcud2lkdGg9MTJ9CkVuaGFuY2VkVm9sY2FubyhtYXJrZXJzX3dpbGNveCwKICAgICAgICAgICAgICAgIGxhYiA9IG1hcmtlcnNfd2lsY294JGdlbmUsCiAgICAgICAgICAgICAgICB4ID0gImF2Z19sb2cyRkMiLAogICAgICAgICAgICAgICAgeSA9ICJwX3ZhbF9hZGoiLAogICAgICAgICAgICAgICAgdGl0bGUgPSAiTUFTVCB3aXRoIEJhdGNoIENvcnJlY3Rpb24iLAogICAgICAgICAgICAgICAgcEN1dG9mZiA9IDAuMDUsCiAgICAgICAgICAgICAgICBGQ2N1dG9mZiA9IDEuMCkKCmBgYAoKIyMgNy4zIFZvbGNhbm8gUGxvdCBmb3IgTUFTVCB3aXRob3V0IEJhdGNoIENvcnJlY3Rpb24KYGBge3IsIGZpZy5oZWlnaHQ9MTIsIGZpZy53aWR0aD0xMn0KRW5oYW5jZWRWb2xjYW5vKG1hcmtlcnNfbWFzdF9ub19iYXRjaCwKICAgICAgICAgICAgICAgIGxhYiA9IG1hcmtlcnNfbWFzdF9ub19iYXRjaCRnZW5lLAogICAgICAgICAgICAgICAgeCA9ICJhdmdfbG9nMkZDIiwKICAgICAgICAgICAgICAgIHkgPSAicF92YWxfYWRqIiwKICAgICAgICAgICAgICAgIHRpdGxlID0gIk1BU1Qgd2l0aCBCYXRjaCBDb3JyZWN0aW9uIiwKICAgICAgICAgICAgICAgIHBDdXRvZmYgPSAwLjA1LAogICAgICAgICAgICAgICAgRkNjdXRvZmYgPSAxLjApCgpgYGAKCiMjIDcuNCBIZWF0bWFwIG9mIFRvcCBEaWZmZXJlbnRpYWxseSBFeHByZXNzZWQgR2VuZXMKYGBge3IsIGZpZy5oZWlnaHQ9NiwgZmlnLndpZHRoPTEwfQojIExvYWQgdGhlIHBoZWF0bWFwIHBhY2thZ2UKbGlicmFyeShwaGVhdG1hcCkKCiMgRGVmaW5lIGNlbGwgZ3JvdXBzCmdyb3VwMV9jZWxscyA8LSBXaGljaENlbGxzKEFsbF9zYW1wbGVzX01lcmdlZCwgaWRlbnRzID0gYygiMCIsICIxIiwgIjIiLCAiMyIsICI0IiwgIjUiLCAiNyIsICI5IikpCmdyb3VwMl9jZWxscyA8LSBXaGljaENlbGxzKEFsbF9zYW1wbGVzX01lcmdlZCwgaWRlbnRzID0gYygiNiIsICI4IiwgIjEwIiwgIjExIikpCgojIERlZmluZSB0b3AgZ2VuZXMgYmFzZWQgb24gYWRqdXN0ZWQgcC12YWx1ZQp0b3BfZ2VuZXMgPC0gaGVhZChtYXJrZXJzX21hc3RfYmF0Y2gkZ2VuZVtvcmRlcihtYXJrZXJzX21hc3RfYmF0Y2gkcF92YWxfYWRqKV0sIDIwKQoKIyBFeHRyYWN0IHRoZSBjb3JyZXNwb25kaW5nIGV4cHJlc3Npb24gZGF0YQpleHByZXNzaW9uX2RhdGEgPC0gR2V0QXNzYXlEYXRhKEFsbF9zYW1wbGVzX01lcmdlZCwgc2xvdCA9ICJkYXRhIikgICMgTG9nLW5vcm1hbGl6ZWQgdmFsdWVzCmhlYXRtYXBfZGF0YSA8LSBleHByZXNzaW9uX2RhdGFbdG9wX2dlbmVzLCBjKGdyb3VwMV9jZWxscywgZ3JvdXAyX2NlbGxzKV0KCiMgQ3JlYXRlIGFubm90YXRpb24gZGF0YQphbm5vdGF0aW9uX2RhdGEgPC0gZGF0YS5mcmFtZSgKICBDZWxsTGluZSA9IGZhY3RvcihBbGxfc2FtcGxlc19NZXJnZWQkY2VsbF9saW5lKVtjKGdyb3VwMV9jZWxscywgZ3JvdXAyX2NlbGxzKV0sCiAgQ2x1c3RlciA9IGZhY3RvcihBbGxfc2FtcGxlc19NZXJnZWQkc2V1cmF0X2NsdXN0ZXJzW2MoZ3JvdXAxX2NlbGxzLCBncm91cDJfY2VsbHMpXSkKKQoKIyBTZXQgcm93IG5hbWVzIGZvciBhbm5vdGF0aW9uIGRhdGEgdG8gbWF0Y2ggY29sdW1uIG5hbWVzIG9mIGhlYXRtYXBfZGF0YQpyb3duYW1lcyhhbm5vdGF0aW9uX2RhdGEpIDwtIGNvbG5hbWVzKGhlYXRtYXBfZGF0YSkKCiMgR2VuZXJhdGUgdGhlIGhlYXRtYXAKcGhlYXRtYXAoaGVhdG1hcF9kYXRhLCAKICAgICAgICAgY2x1c3Rlcl9yb3dzID0gVFJVRSwgCiAgICAgICAgIGNsdXN0ZXJfY29scyA9IFRSVUUsIAogICAgICAgICBzaG93X3Jvd25hbWVzID0gVFJVRSwKICAgICAgICAgc2hvd19jb2xuYW1lcyA9IEZBTFNFLCAjIEhpZGUgY2VsbCBsYWJlbHMgYmVsb3cgdGhlIGhlYXRtYXAKICAgICAgICAgYW5ub3RhdGlvbl9jb2wgPSBhbm5vdGF0aW9uX2RhdGEpCmBgYAoKCgoK