Equivalent to gstat::plot.gstatVariogram
, but using ggplot2 to be more
customizable.
Arguments
- res
Cross variogram results for one sample, from
calculateBivariate
. Global bivariate results are not stored in the SFE object.- plot_np
Logical, whether to plot the number of pairs in each distance bin instead of the variance.
Examples
library(SFEData)
library(scater)
sfe <- McKellarMuscleData()
#> see ?SFEData and browseVignettes('SFEData') for documentation
#> loading from cache
sfe <- sfe[,sfe$in_tissue]
sfe <- logNormCounts(sfe)
res <- calculateBivariate(sfe, type = "cross_variogram_map",
feature1 = c("Myh1", "Myh2", "Csrp3"), swap_rownames = "symbol",
width = 500, cutoff = 2000)
plotCrossVariogramMap(res)