write new file now that status error is fixed
add a few more complex models

species*strata, pguild, psi*strata
need to determine best way to set breakpoint
need to determine best way to set p(1-->2) == p(2-->2)