I'd like to use terastructure to check if there are multiple populations in a supposedly homogeneous dataset. For this reason I ran it first with k=2 and it converged nicely. But I need to compare the held-out likelihood for various values of k, including k=1. When I try to run it with k=1 to get a baseline likelihood value, it never stops running, although likelihood stays the same:
I'd like to use terastructure to check if there are multiple populations in a supposedly homogeneous dataset. For this reason I ran it first with
k=2and it converged nicely. But I need to compare the held-out likelihood for various values ofk, includingk=1. When I try to run it withk=1to get a baseline likelihood value, it never stops running, although likelihood stays the same: