floodPseudotime

[allisoncheney]

Hi, I'm having trouble with the floodPseudotime function.

rootcells <- cellsInCluster(hffURD, "ARCcluster", "12")

hffURD.floods <- floodPseudotime(hffURD, root.cells = rootcells, n=50, minimum.cells.flooded = 2,  verbose=TRUE)

[1] "Starting flood number 1"
Error in intI(j, n = d[2], dn[[2]], give.dn = FALSE) : no 'dimnames[[.]]': cannot use character indexing

here's the traceback:
11. stop("no 'dimnames[[.]]': cannot use character indexing")
10. intI(j, n = d[2], dn[[2]], give.dn = FALSE)
9. subCsp_ij(x, i, j, drop = drop)
8. tm.flood[cells.not.visited, cells.visited, drop = F]
7. tm.flood[cells.not.visited, cells.visited, drop = F]
6. apply(tm.flood[cells.not.visited, cells.visited, drop = F], 1, combine.probs)
5. floodPseudotimeCalc(object, start.cells = root.cells, minimum.cells.flooded = minimum.cells.flooded, tm.flood = tm.flood, verbose.freq = verbose.freq)
4. FUN(X[[i]], ...)
3. lapply(1:n, function(x) { if (verbose) print(paste("Starting flood number", x)) floodPseudotimeCalc(object, start.cells = root.cells, minimum.cells.flooded = minimum.cells.flooded, ...
2. as.data.frame(lapply(1:n, function(x) { if (verbose) print(paste("Starting flood number", x)) floodPseudotimeCalc(object, start.cells = root.cells, minimum.cells.flooded = minimum.cells.flooded, ...

1. floodPseudotime(hffURD, root.cells = rootcells, n = 50, minimum.cells.flooded = 2, verbose = TRUE)

I thought maybe it was that the cell names somehow weren't present in the data, but
all(rootcells %in% colnames([email protected]))
returns TRUE.

Thank you!

[farrellja]

Hi @allisoncheney,

the floodPseudotime function looks at [email protected]. Does that slot have data with cell names as columns and genes as rows?