@cache_readonly
    def bse_re(self):
        p = self.model.exog.shape[1]
        return np.sqrt(self.scale * np.diag(self.cov_params())[p:])
Returns the standard errors of the variance parameters.  Note
dataset/ETHPy150Open statsmodels/statsmodels/statsmodels/regression/mixed_linear_model.py/MixedLMResults.bse_re