Updated to use read_stanfit()