just don't include profiles of missing clusters (in which case they won't align to logliks, causing flightpath trouble)
-- then should we output 2 loglik matrices, e.g. "logliks" and "logliks_including_lost_clusters"?
include the original fixedprofiles values for missing clusters
-- potentially with NAs for missing genes.
-- or potentially removing genes that are NA in the original fixedprofiles
... is there any reason to think we want to keep the complete logliks? Is anyone going to do anything with a lost cluster? Maybe people will want to hack with it, in which case outputting an extra logliks matrix is fine.
This happens at the end of insitutypeML, and maybe elsewhere, causing an error.
These lines:
Options for handling this:
... is there any reason to think we want to keep the complete logliks? Is anyone going to do anything with a lost cluster? Maybe people will want to hack with it, in which case outputting an extra logliks matrix is fine.