Improve logging in PerformTcrClustering

Author: bbimber

singlecell/resources/chunks/PerformTcrClustering.R

@@ -2,6 +2,7 @@ for (datasetId in names(seuratObjects)) {
     printName(datasetId)
     seuratObj <- readSeuratRDS(seuratObjects[[datasetId]])
 
+    print(paste0('Calculating distances for: ', datasetId))
     seuratObj <- tcrClustR::CalculateTcrDistances(
       inputData = seuratObj,
       chains = c('TRA', 'TRB', 'TRG', 'TRD'),
@@ -10,6 +11,7 @@ for (datasetId in names(seuratObjects)) {
       calculateChainPairs = TRUE
     )
 
+    print('Performing TCR Clustering')
     seuratObj <- tcrClustR::RunTcrClustering(
       seuratObj_TCR = seuratObj,
       dianaHeight = 20,
@@ -22,7 +24,8 @@ for (datasetId in names(seuratObjects)) {
     } else {
       for (an in names(seuratObj@misc$TCR_Distances)) {
         ad <- seuratObj@misc$TCR_Distances[[an]]
-        print(paste0('Assay: ', an, ', total clones: ', nrow(ad)))
+        fn <- length(unique(seuratObj[[paste0(an, '_ClusterIdx')]]))
+        print(paste0('Assay: ', an, ', total clones: ', nrow(ad), '. Distinct families: ', fn))
       }
     }