Just came across this cool preprint from the Pachter lab about improving the ordering of variables in alluvial plots and the R implementation is still a work in progress at https://github.com/pachterlab/wompwomp so just putting this here for the future.