analysis/scripts/run_spase_cere3_overall_bias.R

library(spacexr)
library(spASE)

res <- readRDS('results/rctd_cere_3.rds')
maternal_counts_matrix <- res@originalSpatialRNA@maternalCounts
paternal_counts_matrix <- res@originalSpatialRNA@paternalCounts
myfit <- scase(
  maternal_counts_matrix,
  paternal_counts_matrix,
  min.cells = 2^7, cores=8, verbose=F
)
saveRDS(myfit, 'results/results_overall_bias_cere_3.rds')
lulizou/spASE documentation built on May 22, 2024, 5:24 a.m.