Update smart_pca.R

Fixed issue #7
ChristianHuber · Jan 24, 2024 · 01e53da · 01e53da
1 parent 37b1395
commit 01e53da
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/smart_pca.R b/R/smart_pca.R
@@ -375,7 +375,7 @@ smart_pca <- function(snp_data, packed_data = FALSE,
     genoMean <- genoMean[keepSNPs] # compute SNP means
     genoVar <- genoVar[keepSNPs] # compute SNP variances
     if (!isFALSE(sample_project)) {
-      snp_dat2 <- snp_dat2[keepSNPs, ] # remove invariant SNPs from ancient samples
+      snp_dat2 <- snp_dat2[keepSNPs, , drop = FALSE] # remove invariant SNPs from ancient samples
     }
   }
   rm(keepSNPs) # remove vector with indices for variant SNPs