While running code "run_cicero", I keep getting "The previous R session was abnormally terminated due to an unexpected crash"
I don't know what to do,I also test the other data, meet the same issue!
Here is my code
library(cicero)
library(Signac)
library(Seurat)
library(SeuratWrappers)
library(ggplot2)
library(patchwork)
library(monocle3)
bone.cds <- as.cell_data_set(x = bone, assay = "ATAC")
bone.cds <- monocle3::detect_genes(bone.cds)
#Ensure there are no peaks included with zero reads
bone.cds <- bone.cds[Matrix::rowSums(exprs(bone.cds)) != 0,]
bone.cds <- detect_genes(bone.cds)
bone.cds <- estimate_size_factors(bone.cds)
bone.cds <- preprocess_cds(bone.cds, method = "LSI")
bone.cds <- reduce_dimension(bone.cds, reduction_method = "UMAP", preprocess_method = "LSI")
bone.cicero <- make_cicero_cds(bone.cds,reduced_coordinates = reducedDims(bone.cds)$UMAP)
genome <- seqlengths(bone@assays$ATAC@annotation@seqinfo)
use chromosome 1 to save some time
omit this step to run on the whole genome
genome <- genome[30]
convert chromosome sizes to a dataframe
genome.df <- data.frame("chr" = names(genome), "length" = genome)
run cicero
conns <- run_cicero(bone.cicero, genomic_coords = genome.df, sample_num = 100)
And I don't want to update my seurat!!
other attached packages:
[1] patchwork_1.1.3 ggplot2_3.4.4 SeuratWrappers_0.2.0
[4] SeuratObject_4.1.3 Seurat_4.3.0 Signac_1.11.0
[7] cicero_1.3.9 Gviz_1.46.1 monocle3_1.3.4
[10] SingleCellExperiment_1.24.0 SummarizedExperiment_1.32.0 GenomicRanges_1.54.1
[13] GenomeInfoDb_1.38.4 IRanges_2.36.0 S4Vectors_0.40.2
[16] MatrixGenerics_1.14.0 matrixStats_1.2.0 Biobase_2.62.0
[19] BiocGenerics_0.48.1