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)