By the way, note that at 08_layer_differential_expression the color legend includes colors not used. It's a small detail that we don't need to deal with right now, but well, just take in into account when updating 09_region_differential_expression.
This makes me think that some other parts of the DE analysis at the region level are incomplete or not up to date. I know this was by design since we wanted to figure out the code well with the layer DE analysis before doing the region DE one. So we might need to chat about filterByExpr() + the statistical model used for the region analysis.
thanks @lcolladotor. as you noted the region analysis is not finished yet. The model is not correct and there are probably some other issues still. That's why the plots haven't been polished yet.
The PCA plots (currently not version controlled) from https://github.com/LieberInstitute/spatialDLPFC/tree/main/plots/09_region_differential_expression don't match the style from those at https://github.com/LieberInstitute/spatialDLPFC/tree/main/plots/08_layer_differential_expression, so I know that they are older plots.
Some things these don't have are:
Compare against https://github.com/LieberInstitute/spatialDLPFC/blob/main/plots/08_layer_differential_expression/pca_explanatory_variables_k2.pdf.
By the way, note that at
08_layer_differential_expression
the color legend includes colors not used. It's a small detail that we don't need to deal with right now, but well, just take in into account when updating09_region_differential_expression
.This makes me think that some other parts of the DE analysis at the region level are incomplete or not up to date. I know this was by design since we wanted to figure out the code well with the layer DE analysis before doing the region DE one. So we might need to chat about
filterByExpr()
+ the statistical model used for the region analysis.