Getting Started With SEraster (original) (raw)
Rasterized cell-type labels can be used to analyze pair-wise cell-type co-enrichment To do so, we binarize the rasterized cell-type labels using a relative enrichment metric and a previously developed tool called CooccurrenceAffinity
. Please refer to our paper for more details about the methodology and CooccurrenceAffinityfor more details about the package.
# extract cell-type labels
ct_labels <- as.factor(colData(merfish_mousePOA)$celltype)
# compute relative enrichment (RE) metric
mat <- assay(rastCt, "pixelval")
mat_re <- do.call(rbind, lapply(rownames(rastCt), function(ct_label) {
mat[ct_label,] / (sum(mat[ct_label,]) / sum(mat) * colSums(mat))
}))
rownames(mat_re) <- rownames(mat)
# binarize
mat_bin <- ifelse(mat_re >= 1, 1, 0)
# add RE and binarized layers to SpatialExperiment object
assays(rastCt) <- list(pixelval = assay(rastCt, "pixelval"), re = mat_re, bin = mat_bin)
ct_interest <- "Ependymal"
# plot pixel value for a cell-type of interest
plotRaster(rastCt, assay_name = "pixelval", feature_name = ct_interest, name = "cell-type counts", option = "inferno")
# plot RE value for a cell-type of interest
plotRaster(rastCt, assay_name = "re", feature_name = ct_interest, name = "RE", option = "inferno")
# plot binarized value for a cell-type of interest
plotRaster(rastCt, assay_name = "bin", feature_name = ct_interest, factor_levels = c(0,1), name = "binarized", option = "inferno")
# run CooccurrenceAffinity
ct_coocc <- CooccurrenceAffinity::affinity(data = mat_bin, row.or.col = "row", squarematrix = c("all"))
#>
#> ------ as expected, the data ready for analysis has only 1 and 0... 1 = present, 0 = absent used for the interpretation ------
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AcceptAffCI(x, marg, lev, Int2): NAs introduced by coercion
#> Warning in AcceptAffCI(x, marg, lev, Int2): NAs introduced by coercion
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = Infty is capped, along with upper confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = Infty is capped, along with upper confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> Warning in AlphInts(x, marg, lev = lev, scal = scal, pvalType = pvalType): MLE = -Infty is capped, along with lower confidence limits
#> ~~~~~~~~~~ printing head of all elements of the output list ~~~~~~~~~~
#> $all
#> entity_1 entity_2 entity_1_count_mA entity_2_count_mB obs_cooccur_X
#> 1 Ambiguous Astrocyte 561 557 232
#> 2 Ambiguous Endothelial 1 561 308 136
#> 3 Ambiguous Endothelial 2 561 20 9
#> 4 Ambiguous Endothelial 3 561 102 49
#> 5 Ambiguous Ependymal 561 47 4
#> 6 Ambiguous Excitatory 561 504 192
#> total_N p_value exp_cooccur alpha_mle alpha_medianInt conf_level
#> 1 1301 0.3657 240.182 -0.105 [-0.111, -0.098] 0.95
#> 2 1301 0.693 132.812 0.055 [0.047, 0.064] 0.95
#> 3 1301 1 8.624 0.077 [-0.021, 0.181] 0.95
#> 4 1301 0.2998 43.983 0.215 [0.194, 0.237] 0.95
#> 5 1301 2.1435e-07 20.267 -2.150 [-2.252, -1.986] 0.95
#> 6 1301 0.0041 217.328 -0.337 [-0.343, -0.33] 0.95
#> ci_blaker ci_cp ci_midQ ci_midP jaccard
#> 1 [-0.331, 0.119] [-0.333, 0.123] [-0.327, 0.117] [-0.327, 0.117] 0.262
#> 2 [-0.208, 0.313] [-0.212, 0.321] [-0.203, 0.312] [-0.203, 0.313] 0.186
#> 3 [-0.894, 1.04] [-0.934, 1.061] [-0.823, 0.959] [-0.846, 0.98] 0.016
#> 4 [-0.201, 0.625] [-0.212, 0.641] [-0.19, 0.62] [-0.192, 0.622] 0.080
#> 5 [-3.27, -1.164] [-3.5, -1.128] [-3.291, -1.223] [-3.341, -1.2] 0.007
#> 6 [-0.569, -0.104] [-0.571, -0.104] [-0.565, -0.11] [-0.565, -0.11] 0.220
#> sorensen simpson errornote
#> 1 0.415 0.417 NA
#> 2 0.313 0.442 NA
#> 3 0.031 0.450 NA
#> 4 0.148 0.480 NA
#> 5 0.013 0.085 NA
#> 6 0.361 0.381 NA
#>
#> $occur_mat
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3 Ependymal
#> pixel23 1 0 0 0 0 0
#> pixel24 0 0 1 0 0 0
#> pixel25 0 1 0 0 0 0
#> pixel26 0 0 0 0 0 0
#> pixel27 1 0 1 0 0 0
#> pixel28 0 1 0 0 0 0
#> Excitatory Inhibitory Microglia OD Immature 1 OD Immature 2 OD Mature 1
#> pixel23 0 1 0 0 0 0
#> pixel24 0 1 0 0 0 0
#> pixel25 0 0 0 1 0 0
#> pixel26 1 0 0 0 0 0
#> pixel27 0 0 0 0 0 0
#> pixel28 1 0 0 0 0 0
#> OD Mature 2 OD Mature 3 OD Mature 4 Pericytes
#> pixel23 0 0 0 0
#> pixel24 0 0 0 0
#> pixel25 0 0 0 0
#> pixel26 0 0 0 0
#> pixel27 0 0 0 0
#> pixel28 0 0 0 0
#>
#> $alpha_mle
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous NA NA NA NA NA
#> Astrocyte -0.105 NA NA NA NA
#> Endothelial 1 0.055 -0.015 NA NA NA
#> Endothelial 2 0.077 -0.118 0.779 NA NA
#> Endothelial 3 0.215 -0.029 0.322 1.902 NA
#> Ependymal -2.150 -0.807 -0.143 0.345 -0.667
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#>
#> $alpha_mle_sig
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA 1.902 NA
#> Ependymal -2.15 -0.807 NA NA NA
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#>
#> $p_value
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous <NA> <NA> <NA> <NA> NA
#> Astrocyte 0.3657 <NA> <NA> <NA> NA
#> Endothelial 1 0.693 0.9475 <NA> <NA> NA
#> Endothelial 2 1 0.825 0.1079 <NA> NA
#> Endothelial 3 0.2998 0.9173 0.1813 5e-04 NA
#> Ependymal 2.1435e-07 0.0159 0.7323 1 0.4274
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous <NA> <NA> <NA> NA NA
#> Astrocyte <NA> <NA> <NA> NA NA
#> Endothelial 1 <NA> <NA> <NA> NA NA
#> Endothelial 2 <NA> <NA> <NA> NA NA
#> Endothelial 3 <NA> <NA> <NA> NA NA
#> Ependymal <NA> <NA> <NA> NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#>
#> $cooccur.null
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous NA NA NA NA NA
#> Astrocyte 240.182 NA NA NA NA
#> Endothelial 1 132.812 131.865 NA NA NA
#> Endothelial 2 8.624 8.563 4.735 NA NA
#> Endothelial 3 43.983 43.669 24.148 1.568 NA
#> Ependymal 20.267 20.122 11.127 0.723 3.685
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#>
#> $cooccur.obs
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous NA NA NA NA NA
#> Astrocyte 232 NA NA NA NA
#> Endothelial 1 136 131 NA NA NA
#> Endothelial 2 9 8 8 NA NA
#> Endothelial 3 49 43 30 7 NA
#> Ependymal 4 12 10 1 2
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#>
#> $jaccard
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous NA NA NA NA NA
#> Astrocyte 0.262 NA NA NA NA
#> Endothelial 1 0.186 0.178 NA NA NA
#> Endothelial 2 0.016 0.014 0.025 NA NA
#> Endothelial 3 0.080 0.070 0.079 0.061 NA
#> Ependymal 0.007 0.020 0.029 0.015 0.014
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#>
#> $jaccard_sig
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA 0.061 NA
#> Ependymal 0.007 0.02 NA NA NA
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#>
#> $sorensen
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous NA NA NA NA NA
#> Astrocyte 0.415 NA NA NA NA
#> Endothelial 1 0.313 0.303 NA NA NA
#> Endothelial 2 0.031 0.028 0.049 NA NA
#> Endothelial 3 0.148 0.131 0.146 0.115 NA
#> Ependymal 0.013 0.040 0.056 0.030 0.027
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#>
#> $sorensen_sig
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA 0.115 NA
#> Ependymal 0.013 0.04 NA NA NA
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#>
#> $simpson
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous NA NA NA NA NA
#> Astrocyte 0.417 NA NA NA NA
#> Endothelial 1 0.442 0.425 NA NA NA
#> Endothelial 2 0.450 0.400 0.400 NA NA
#> Endothelial 3 0.480 0.422 0.294 0.35 NA
#> Ependymal 0.085 0.255 0.213 0.05 0.043
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#>
#> $simpson_sig
#> Ambiguous Astrocyte Endothelial 1 Endothelial 2 Endothelial 3
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA 0.35 NA
#> Ependymal 0.085 0.255 NA NA NA
#> Ependymal Excitatory Inhibitory Microglia OD Immature 1
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> OD Immature 2 OD Mature 1 OD Mature 2 OD Mature 3 OD Mature 4
#> Ambiguous NA NA NA NA NA
#> Astrocyte NA NA NA NA NA
#> Endothelial 1 NA NA NA NA NA
#> Endothelial 2 NA NA NA NA NA
#> Endothelial 3 NA NA NA NA NA
#> Ependymal NA NA NA NA NA
#> Pericytes
#> Ambiguous NA
#> Astrocyte NA
#> Endothelial 1 NA
#> Endothelial 2 NA
#> Endothelial 3 NA
#> Ependymal NA
#> ~~~~~~~~~~ COMPLETED: printing head of all elements of the output list ~~~~~~~~~~
# plot maximum likelihood estimates of affinity metric (alpha MLE)
CooccurrenceAffinity::plotgg(data = ct_coocc, variable = "alpha_mle", legendlimit = "datarange")
#> you can hide the printed values with show.value=F
#> use the argument value.digit to change number of digits and text.size to adjust the text size