To perform clustering on a seuset
object, the function FindClusters()
from the package Seurat
can be used:
seuset <- FindClusters(
object = seuset,
print.output = FALSE
)
We now want to compare our clustering to the clustering from the published mouse epithelium paper, which used the code:
pub <- seuset@meta.data[, "publishedClusters"]
table(seuset@ident, pub)
The output is as follows:
published
7 9 15
0 41 1 6
1 0 45 1
Is our new clustering the same as the clustering from the mouse epithelium paper?