Finding the variable genes between cells

[fereshtehizadi]

Sorry, I have 1563 cells in 8 time points (2h to 16 h).

If I am not wrong, I should change

var.by.stage <- lapply(seq(3,12,3), function(n) {
  findVariableGenes(axial, cells.fit=cellsInCluster(axial, "stage", stages[(n-2):n]), set.object.var.genes=F, diffCV.cutoff=0.3, mean.min=.005, mean.max=100, main.use=paste0("Stages ", stages[n-2], " to ", stages[n]), do.plot=TRUE)
})

to

var.by.stage <- lapply(seq(3,12,4), function(n) {
  findVariableGenes(axial, cells.fit=cellsInCluster(axial, "stage", stages[(n-2):n]), set.object.var.genes=F, diffCV.cutoff=0.3, mean.min=.005, mean.max=100, main.use=paste0("Stages ", stages[n-2], " to ", stages[n]), do.plot=TRUE)
})

As I 8 stages. seq(3,12,3) to seq(3,12,4)