Skip to content

Commit feb6601

Browse files
committed
Add debugging to seurat pipeline
1 parent 41e83f4 commit feb6601

File tree

2 files changed

+10
-1
lines changed

2 files changed

+10
-1
lines changed

singlecell/resources/chunks/Functions.R

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -62,9 +62,12 @@ printName <- function(datasetId) {
6262

6363
savedFiles <- data.frame(datasetId = character(), datasetName = character(), filename = character(), outputFileId = character(), readsetId = character())
6464
if (file.exists('/work/savedSeuratObjects.txt')) {
65+
print('Deleting pre-existing savedSeuratObjects.txt file')
6566
unlink('/work/savedSeuratObjects.txt')
6667
}
68+
6769
file.create('/work/savedSeuratObjects.txt')
70+
print(paste0('Total lines in savedSeuratObjects.txt on job start:', length(readLines('savedSeuratObjects.txt'))))
6871

6972
saveData <- function(seuratObj, datasetId) {
7073
print(paste0('Saving dataset: ', datasetId))
@@ -88,6 +91,12 @@ saveData <- function(seuratObj, datasetId) {
8891
print(paste0('readsetId: ', readsetId))
8992

9093
toAppend <- data.frame(datasetId = datasetId, datasetName = datasetName, filename = fn, outputFileId = outputFileId, readsetId = readsetId)
94+
if (nrow(toAppend) != 1) {
95+
warning(paste0('Error saving seurat objects, more than one row:'))
96+
print(toAppend)
97+
stop('Error saving seurat objects, more than one row!')
98+
}
99+
91100
write.table(toAppend, file = 'savedSeuratObjects.txt', quote = FALSE, sep = '\t', row.names = FALSE, col.names = FALSE, append = TRUE)
92101

93102
# Write cell barcodes and metadata:

singlecell/src/org/labkey/singlecell/analysis/AbstractSingleCellHandler.java

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -827,7 +827,7 @@ public List<SingleCellStep.SeuratObjectWrapper> getSeuratFromStep(int stepIdx)
827827

828828
public void setStepComplete(Logger log, SingleCellStep step, int stepIdx, RecordedAction action, List<SingleCellStep.SeuratObjectWrapper> seurat, File markdown, File html) throws PipelineJobException
829829
{
830-
log.info("Marking step complete: " + step.getProvider().getName() + ", " + stepIdx);
830+
log.info("Marking step complete: " + step.getProvider().getName() + ", " + stepIdx + ", total seurat objects: " + (seurat == null ? 0 : seurat.size()));
831831

832832
if (seurat != null)
833833
{

0 commit comments

Comments
 (0)