analysis/scripts/run_spase_hippo1_overall_spatial.R

In lulizou/spASE: Detecting Allele-Specific Expression in Spatial Transcriptomics

library(spacexr)
library(spASE)

res <- readRDS('results/rctd_hippo_1.rds')
maternal_counts_matrix <- res@originalSpatialRNA@maternalCounts
paternal_counts_matrix <- res@originalSpatialRNA@paternalCounts
coords <- res@originalSpatialRNA@coords
coords$bead <- rownames(coords)
coords <- coords[c('bead','x','y')]

myfit <- spase(maternal_counts_matrix, paternal_counts_matrix, coords,min.pixels=2^7,cores=8,verbose=T)
saveRDS(myfit, file = 'results/results_overall_spatial_hippo_1.rds')