RunBanksy(codex.obj, lambda = 0.8, verbose = TRUE, assay = "Spatial",
slot = "data", features = "variable", k_geom = 5)

RunPCA(object = codex.obj, npcs = max_pc, verbose = FALSE,
reduction.name = "pca")

RunPCA(codex.obj, assay = "BANKSY", reduction.name = "pca.banksy",
features = rownames(codex.obj), npcs = 30)

RunUMAP(object = codex.obj, dims = 1:npc, verbose = FALSE,
reduction = "pca", reduction.name = paste0("umap_pc", npc))

FindNeighbors(object = codex.obj, dims = 1:npc, verbose = FALSE,
reduction = "pca")

FindNeighbors(codex.obj, reduction = "pca.banksy", dims = 1:npc)