The cell type vector is generated from the following code
celltype_toy
A vector of length 5 contains 3 cell types
celltype = as.factor(sample(c(1,2,3),5,replace = TRUE))