FindVariableFeatures(data, selection.method = "vst", nfeatures = 2000,
verbose = FALSE)

ScaleData(data, features = VariableFeatures(data), verbose = FALSE)

RunPCA(data, npcs = pc, verbose = FALSE)

RunUMAP(data, features = VariableFeatures(data), verbose = FALSE)

FindNeighbors(data, dims = 1:pc, verbose = FALSE)

FindClusters(data, verbose = FALSE, resolution = 0.5)

FindVariableFeatures(data, selection.method = "vst", nfeatures = 2000,
verbose = FALSE)

FindVariableFeatures(data, selection.method = "vst", nfeatures = 2000,
verbose = FALSE)

FindVariableFeatures(data, selection.method = "vst", nfeatures = 2000,
verbose = FALSE)

FindVariableFeatures(data, selection.method = "vst", nfeatures = 2000,
verbose = FALSE)