Update lawlor data preperation

code/prep-lawlor.R

@@ -44,9 +44,14 @@ lawlor <- computeSumFactors(lawlor, clusters = clusters)
 lawlor <- logNormCounts(lawlor)
 
 dec_lawlor <- modelGeneVarByPoisson(lawlor)
-top_lawlor <- getTopHVGs(dec_lawlor, n = 5000)
+top_lawlor <- getTopHVGs(dec_lawlor, n = 10000)
 lawlor <- lawlor[top_lawlor, ]
 
+old_rownames <- rownames(lawlor)
+rownames(lawlor) <- rowData(lawlor)$SYMBOL
+is_na_rownames <- which(is.na(rownames(lawlor)))
+rownames(lawlor)[is_na_rownames] <- old_rownames[is_na_rownames]
+
 lawlor <- runPCA(lawlor)
 lawlor <- runTSNE(lawlor, dimred = "PCA")
 lawlor <- runUMAP(lawlor, dimred = "PCA")