[x] pairwise model boxplots on both apps do not display colors. It's because the test column of sig_genes. Have to remove "BayesSpace" unique(sig_genes$test)
[x] pairwise results crash at k = 16? make them smaller? filter for fdr < 0.05. Do it in a separate script and save sig_genes object.
[x] some genes are missing in the k = 9 modeling results. check for these genes in the pseudobulk object, and in the pre and post filtered spe objects. https://neuralplasticity.slack.com/archives/DU3UP2A9W/p1655841909513549. CTGF doesn't exist in the original spe or any of the downstream spe's. PDYN exists in the original spe but was filtered out in the pseudobulking process.
unique(sig_genes$test)