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

RunPCA(data, features = VariableFeatures(object = data))

RunUMAP(data, dims = 1:20, min.dist = 0.75)

RunTSNE(data, dims = 1:10, nthreads = 4, max_iter = 1000,
k.seed = 12345)