docs: def fit_df(self, dfs, pstate_col=PSTATE_COL):
          """
          Convenience function to fit a model from a list of dataframes
          """
          obs_cols = list(self.emission_name)
          obs = [df[df.columns.difference([pstate_col])][obs_cols].values for df in dfs]
          pstates = [df[pstate_col].values for df in dfs]
          return self.fit(obs, pstates)