import numpy as np

coef_ml_drf_0 = np.array([-0.9887517])

vcov_ml_drf_0 = np.array([0.001317148]).reshape(1, 1, order='F')

cov_re_ml_drf_0 = np.array([0.2522485]).reshape(1, 1, order='F')

scale_ml_drf_0 = np.array([0.2718486])

loglike_ml_drf_0 = np.array([-240.1254])

ranef_mean_ml_drf_0 = np.array([0.04977167])

ranef_condvar_ml_drf_0 = np.array([0.130841])

coef_reml_drf_0 = np.array([-0.9887533])

vcov_reml_drf_0 = np.array([0.001323559]).reshape(1, 1, order='F')

cov_re_reml_drf_0 = np.array([0.2524129]).reshape(1, 1, order='F')

scale_reml_drf_0 = np.array([0.2733467])

loglike_reml_drf_0 = np.array([-242.5214])

ranef_mean_reml_drf_0 = np.array([0.04964696])

ranef_condvar_reml_drf_0 = np.array([0.1312315])

coef_ml_drf_1 = np.array([-0.9115929])

vcov_ml_drf_1 = np.array([0.01340632]).reshape(1, 1, order='F')

cov_re_ml_drf_1 = np.array([0]).reshape(1, 1, order='F')

scale_ml_drf_1 = np.array([4.050921])

loglike_ml_drf_1 = np.array([-538.0763])

ranef_mean_ml_drf_1 = np.array([0])

ranef_condvar_ml_drf_1 = np.array([0])

coef_reml_drf_1 = np.array([-0.9115929])

vcov_reml_drf_1 = np.array([0.01345931]).reshape(1, 1, order='F')

cov_re_reml_drf_1 = np.array([2.839777e-14]).reshape(1, 1, order='F')

scale_reml_drf_1 = np.array([4.066932])

loglike_reml_drf_1 = np.array([-539.3124])

ranef_mean_reml_drf_1 = np.array([2.424384e-14])

ranef_condvar_reml_drf_1 = np.array([2.839777e-14])

coef_ml_drf_2 = np.array([-1.012044, 0.9789052])

vcov_ml_drf_2 = np.array([
    0.00117849, 1.458744e-05,
    1.458744e-05, 0.001054926]).reshape(2, 2, order='F')

cov_re_ml_drf_2 = np.array([0.1596058]).reshape(1, 1, order='F')

scale_ml_drf_2 = np.array([0.2129146])

loglike_ml_drf_2 = np.array([-200.319])

ranef_mean_ml_drf_2 = np.array([0.3197174])

ranef_condvar_ml_drf_2 = np.array([0.09122291])

coef_reml_drf_2 = np.array([-1.012137, 0.9790792])

vcov_reml_drf_2 = np.array([
    0.001190455, 1.482666e-05,
    1.482666e-05, 0.001066002]).reshape(2, 2, order='F')

cov_re_reml_drf_2 = np.array([0.1595015]).reshape(1, 1, order='F')

scale_reml_drf_2 = np.array([0.2154276])

loglike_reml_drf_2 = np.array([-205.275])

ranef_mean_reml_drf_2 = np.array([0.3172978])

ranef_condvar_reml_drf_2 = np.array([0.09164674])

coef_ml_drf_3 = np.array([-1.028053, 0.8602685])

vcov_ml_drf_3 = np.array([
    0.01398831, 0.001592619, 0.001592619, 0.01602274]).reshape(2, 2, order='F')

cov_re_ml_drf_3 = np.array([0.8130996]).reshape(1, 1, order='F')

scale_ml_drf_3 = np.array([3.100447])

loglike_ml_drf_3 = np.array([-477.1707])

ranef_mean_ml_drf_3 = np.array([-0.2641747])

ranef_condvar_ml_drf_3 = np.array([0.6441656])

coef_reml_drf_3 = np.array([-1.027583, 0.8605714])

vcov_reml_drf_3 = np.array([
    0.01411922, 0.001607343, 0.001607343, 0.01617574]).reshape(2, 2, order='F')

cov_re_reml_drf_3 = np.array([0.8117898]).reshape(1, 1, order='F')

scale_reml_drf_3 = np.array([3.13369])

loglike_reml_drf_3 = np.array([-479.5354])

ranef_mean_reml_drf_3 = np.array([-0.2614875])

ranef_condvar_reml_drf_3 = np.array([0.6447625])

coef_ml_drf_4 = np.array([-1.005151, -0.003657404, 1.054786])

vcov_ml_drf_4 = np.array([
    0.001190639, -5.327162e-05, 5.992985e-05, -5.327162e-05,
    0.001460303, -2.662532e-05, 5.992985e-05, -2.662532e-05,
    0.00148609]).reshape(3, 3, order='F')

cov_re_ml_drf_4 = np.array([0.1703249]).reshape(1, 1, order='F')

scale_ml_drf_4 = np.array([0.251763])

loglike_ml_drf_4 = np.array([-231.6389])

ranef_mean_ml_drf_4 = np.array([-0.2063164])

ranef_condvar_ml_drf_4 = np.array([0.0459578])

coef_reml_drf_4 = np.array([-1.005067, -0.003496032, 1.054666])

vcov_reml_drf_4 = np.array([
    0.001206925, -5.4182e-05, 6.073475e-05, -5.4182e-05,
    0.001479871, -2.723494e-05, 6.073475e-05, -2.723494e-05,
    0.001506198]).reshape(3, 3, order='F')

cov_re_reml_drf_4 = np.array([0.1705659]).reshape(1, 1, order='F')

scale_reml_drf_4 = np.array([0.2556394])

loglike_reml_drf_4 = np.array([-238.761])

ranef_mean_reml_drf_4 = np.array([-0.2055303])

ranef_condvar_reml_drf_4 = np.array([0.04649027])

coef_ml_drf_5 = np.array([-0.8949725, 0.08141558, 1.052529])

vcov_ml_drf_5 = np.array([
    0.01677563, 0.0008077524, -0.001255011, 0.0008077524,
    0.01719346, 0.0009266736, -0.001255011, 0.0009266736,
    0.01608435]).reshape(3, 3, order='F')

cov_re_ml_drf_5 = np.array([0.3444677]).reshape(1, 1, order='F')

scale_ml_drf_5 = np.array([4.103944])

loglike_ml_drf_5 = np.array([-579.4568])

ranef_mean_ml_drf_5 = np.array([0.08254713])

ranef_condvar_ml_drf_5 = np.array([0.3177935])

coef_reml_drf_5 = np.array([-0.8946164, 0.08134261, 1.052486])

vcov_reml_drf_5 = np.array([
    0.0169698, 0.0008162714, -0.001268635, 0.0008162714,
    0.01739219, 0.0009345538, -0.001268635, 0.0009345538,
    0.01627074]).reshape(3, 3, order='F')

cov_re_reml_drf_5 = np.array([0.3420993]).reshape(1, 1, order='F')

scale_reml_drf_5 = np.array([4.155737])

loglike_reml_drf_5 = np.array([-582.8377])

ranef_mean_reml_drf_5 = np.array([0.08111449])

ranef_condvar_reml_drf_5 = np.array([0.3160797])

coef_ml_drf_6 = np.array([-0.8885425])

vcov_ml_drf_6 = np.array([0.002443738]).reshape(1, 1, order='F')

cov_re_ml_drf_6 = np.array([
    0.2595201, 0.04591071,
    0.04591071, 2.204612]).reshape(2, 2, order='F')

scale_ml_drf_6 = np.array([0.243133])

loglike_ml_drf_6 = np.array([-382.551])

ranef_mean_ml_drf_6 = np.array([-0.0597406, 0.6037288])

ranef_condvar_ml_drf_6 = np.array([
    0.2420741, 0.2222169,
    0.2222169, 0.4228908]).reshape(2, 2, order='F')

coef_reml_drf_6 = np.array([-0.8883881])

vcov_reml_drf_6 = np.array([0.002461581]).reshape(1, 1, order='F')

cov_re_reml_drf_6 = np.array([
    0.2595767, 0.04590012,
    0.04590012, 2.204822]).reshape(2, 2, order='F')

scale_reml_drf_6 = np.array([0.2453537])

loglike_reml_drf_6 = np.array([-384.6373])

ranef_mean_reml_drf_6 = np.array([-0.05969892, 0.6031793])

ranef_condvar_reml_drf_6 = np.array([
    0.2421365, 0.2221108,
    0.2221108, 0.4244443]).reshape(2, 2, order='F')

coef_ml_irf_6 = np.array([-0.8874992])

vcov_ml_irf_6 = np.array([0.002445505]).reshape(1, 1, order='F')

cov_re_ml_irf_6 = np.array([
    0.2587624, 0,
    0, 2.188653]).reshape(2, 2, order='F')

scale_ml_irf_6 = np.array([0.2432694])

loglike_ml_irf_6 = np.array([-382.6581])

coef_reml_irf_6 = np.array([-0.8873394])

vcov_reml_irf_6 = np.array([0.002463375]).reshape(1, 1, order='F')

cov_re_reml_irf_6 = np.array([
    0.2588157, 0,
    0, 2.188876]).reshape(2, 2, order='F')

scale_reml_irf_6 = np.array([0.2454935])

loglike_reml_irf_6 = np.array([-384.7441])

coef_ml_drf_7 = np.array([-0.9645281])

vcov_ml_drf_7 = np.array([0.01994]).reshape(1, 1, order='F')

cov_re_ml_drf_7 = np.array([
    0.2051329, 0.0734377,
    0.0734377, 3.285381]).reshape(2, 2, order='F')

scale_ml_drf_7 = np.array([3.423247])

loglike_ml_drf_7 = np.array([-587.7101])

ranef_mean_ml_drf_7 = np.array([0.07007965, -0.2920284])

ranef_condvar_ml_drf_7 = np.array([
    0.1823183, 0.02247519,
    0.02247519, 1.125011]).reshape(2, 2, order='F')

coef_reml_drf_7 = np.array([-0.9647862])

vcov_reml_drf_7 = np.array([0.02002546]).reshape(1, 1, order='F')

cov_re_reml_drf_7 = np.array([
    0.2056226, 0.0726139,
    0.0726139, 3.2876]).reshape(2, 2, order='F')

scale_reml_drf_7 = np.array([3.440244])

loglike_reml_drf_7 = np.array([-588.7476])

ranef_mean_reml_drf_7 = np.array([0.07000628, -0.2916737])

ranef_condvar_reml_drf_7 = np.array([
    0.1828266, 0.02229138,
    0.02229138, 1.128947]).reshape(2, 2, order='F')

coef_ml_irf_7 = np.array([-0.9665524])

vcov_ml_irf_7 = np.array([0.01998144]).reshape(1, 1, order='F')

cov_re_ml_irf_7 = np.array([
    0.2021561, 0, 0, 3.270735]).reshape(2, 2, order='F')

scale_ml_irf_7 = np.array([3.423186])

loglike_ml_irf_7 = np.array([-587.7456])

coef_reml_irf_7 = np.array([-0.9667854])

vcov_reml_irf_7 = np.array([0.02006657]).reshape(1, 1, order='F')

cov_re_reml_irf_7 = np.array([
    0.2026938, 0, 0, 3.273129]).reshape(2, 2, order='F')

scale_reml_irf_7 = np.array([3.440197])

loglike_reml_irf_7 = np.array([-588.782])

coef_ml_drf_8 = np.array([-1.083882, 0.8955623])

vcov_ml_drf_8 = np.array([
    0.002491643, 0.0001693531,
    0.0001693531, 0.00253309]).reshape(2, 2, order='F')

cov_re_ml_drf_8 = np.array([
    0.1506188, 0.126091, 0.126091, 2.485462]).reshape(2, 2, order='F')

scale_ml_drf_8 = np.array([0.2586519])

loglike_ml_drf_8 = np.array([-363.6234])

ranef_mean_ml_drf_8 = np.array([0.2852326, -0.5047804])

ranef_condvar_ml_drf_8 = np.array([
    0.05400391, 0.002330104,
    0.002330104, 0.122761]).reshape(2, 2, order='F')

coef_reml_drf_8 = np.array([-1.083938, 0.8956893])

vcov_reml_drf_8 = np.array([
    0.002528969, 0.0001712206,
    0.0001712206, 0.002573335]).reshape(2, 2, order='F')

cov_re_reml_drf_8 = np.array([
    0.1505098, 0.1256311,
    0.1256311, 2.484219]).reshape(2, 2, order='F')

scale_reml_drf_8 = np.array([0.2635901])

loglike_reml_drf_8 = np.array([-367.7667])

ranef_mean_reml_drf_8 = np.array([0.2829798, -0.5042857])

ranef_condvar_reml_drf_8 = np.array([
    0.05463632, 0.002393538,
    0.002393538, 0.1249828]).reshape(2, 2, order='F')