statsmodels.regression.mixed_linear_model.MixedLMResults.bootstrap

MixedLMResults.bootstrap(nrep=100, method='nm', disp=0, store=1)

simple bootstrap to get mean and variance of estimator

see notes

Parameters:

nrep : int

number of bootstrap replications

method : str

optimization method to use

disp : bool

If true, then optimization prints results

store : bool

If true, then parameter estimates for all bootstrap iterations are attached in self.bootstrap_results

Returns:

mean : array

mean of parameter estimates over bootstrap replications

std : array

standard deviation of parameter estimates over bootstrap replications

Notes

This was mainly written to compare estimators of the standard errors of the parameter estimates. It uses independent random sampling from the original endog and exog, and therefore is only correct if observations are independently distributed.

This will be moved to apply only to models with independently distributed observations.