scRNAseq
scRNAseq copied to clipboard
Missing G2 phase in LengESCData
Column metadata for this dataset supposedly contains experimentally determined cell cycle phase for each cell, however:
sce.ref <- scRNAseq::LengESCData()
table(sce.ref$Phase, useNA="ifany")
G1 S <NA>
91 80 289
So, I looked back at the data and the code used to create this dataset, and found that G2 in fact is present:
phase <- sub("_Exp.*", "", colnames(sce.ref))
table(phase)
phase
G1 G2 H1 S
91 76 213 80
So, this line of code:
phase[!phase %in% c("G1", "S", "G2M")] <- NA
should be changed into
phase[!phase %in% c("G1", "S", "G2")] <- NA
Oops. How embarrassing; fixed in 2.17.10. (This will release as 2.18.0 next week with the new BioC release.)