Good afternoon,

I have a question about seurat analysis problem.

My study combined three dataset. After I integrate the data from three database together, I run following steps:

DefaultAssay(pancreatic) <- "integrated"
pancreatic<-ScaleData(pancreatic)

pancreatic <- RunPCA(pancreatic, npcs = 20, verbose = FALSE)

pdf("elbow_combine.pdf", width = 6, height = 6)
ElbowPlot(pancreatic) #this is to help you select how many pcs to use
dev.off()

saveRDS(pancreatic,'pancreatic.rds'). 

Then I meet an error, please give me some suggestions to solve the problem.

Error in idx[i, ] <- res[[i]][[1]] : 
  number of items to replace is not a multiple of replacement length