avoiding wrong cbinds

FertigLab · May 20, 2024 · 5868616 · 5868616
1 parent 86729d5
commit 5868616
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/modules/local/spacemarkers.nf b/modules/local/spacemarkers.nf
@@ -37,8 +37,10 @@ process SPACEMARKERS {
       dataMatrix <- dataMatrix[keepGenes,]
       
       coords <- load10XCoords("$data");
+      rownames(coords) <- coords$barcode;
       features <- getSpatialFeatures("$cogapsResult");
-      spPatterns <- cbind(coords, features);
+      barcodes <- intersect(rownames(features), rownames(coords))
+      spPatterns <- cbind(coords[barcodes,], features[barcodes,]);
       
       saveRDS(spPatterns, file = "${prefix}/spPatterns.rds");