docs: def get_statsmodels_summary(self,
                                  title=None,
                                  alpha=.05):
          """
          Parameters
          ----------
          title : str, or None, optional.
              Will be the title of the returned summary. If None, the default
              title is used.
          alpha : float, optional.
              Should be between 0.0 and 1.0. Determines the width of the
              displayed, (1 - alpha)% confidence interval.
  
          Returns
          -------
          statsmodels.summary object or None.
          """
          try:
              # Get the statsmodels Summary class
              from statsmodels.iolib.summary import Summary
          except ImportError:
              print("statsmodels not installed. Resorting to standard summary")
              return self.print_summaries()
  
          if not hasattr(self, "estimation_success"):
              msg = "Must estimate a model before a summary can be returned."
              raise NotImplementedError(msg)
  
          # Get an instantiation of the Summary class.
          smry = Summary()
  
          # Get the yname and yname_list.
          # Note I'm not really sure what the yname_list is.
          new_yname, new_yname_list = self.choice_col, None
  
          # Get the model name
          model_name = self.model_type
  
          ##########
          # Note the following commands are basically directly from
          # statsmodels.discrete.discrete_model
          ##########
          top_left = [('Dep. Variable:', None),
                      ('Model:', [model_name]),
                      ('Method:', ['MLE']),
                      ('Date:', None),
                      ('Time:', None),
                      ('AIC:', ["{:,.3f}".format(self.aic)]),
                      ('BIC:', ["{:,.3f}".format(self.bic)])
                      ]
  
          top_right = [('No. Observations:', ["{:,}".format(self.nobs)]),
                       ('Df Residuals:', ["{:,}".format(self.df_resid)]),
                       ('Df Model:', ["{:,}".format(self.df_model)]),
                       ('Pseudo R-squ.:',
                        ["{:.3f}".format(self.rho_squared)]),
                       ('Pseudo R-bar-squ.:',
                        ["{:.3f}".format(self.rho_bar_squared)]),
                       ('Log-Likelihood:', ["{:,.3f}".format(self.llf)]),
                       ('LL-Null:',
                        ["{:,.3f}".format(self.null_log_likelihood)]),
                       ]
  
          if title is None:
              title = model_name + ' ' + "Regression Results"
  
          xnames = self.params.index.tolist()
  
          # for top of table
          smry.add_table_2cols(self,
                               gleft=top_left,
                               gright=top_right,  # [],
                               yname=new_yname,
                               xname=xnames,
                               title=title)
          # for parameters, etc
          smry.add_table_params(self,
                                yname=[new_yname_list],
                                xname=xnames,
                                alpha=alpha,
                                use_t=False)
          return smry