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')

coef_ml_irf_8 = np.array([-1.079481,0.898216])

vcov_ml_irf_8 = np.array([0.002511536,0.0001812511,0.0001812511,0.002573405]).reshape(2,2, order='F')

cov_re_ml_irf_8 = np.array([0.1498568,0,0,2.403849]).reshape(2,2, order='F')

scale_ml_irf_8 = np.array([0.2605245])

loglike_ml_irf_8 = np.array([-364.4824])

coef_reml_irf_8 = np.array([-1.07952,0.8983678])

vcov_reml_irf_8 = np.array([0.002549354,0.0001833386,0.0001833386,0.002614672]).reshape(2,2, order='F')

cov_re_reml_irf_8 = np.array([0.1497193,0,0,2.403076]).reshape(2,2, order='F')

scale_reml_irf_8 = np.array([0.2655558])

loglike_reml_irf_8 = np.array([-368.6141])

coef_ml_drf_9 = np.array([-1.272698,0.8617471])

vcov_ml_drf_9 = np.array([0.02208544,0.001527479,0.001527479,0.02597528]).reshape(2,2, order='F')

cov_re_ml_drf_9 = np.array([0.510175,0.08826114,0.08826114,3.342888]).reshape(2,2, order='F')

scale_ml_drf_9 = np.array([3.722112])

loglike_ml_drf_9 = np.array([-589.8274])

ranef_mean_ml_drf_9 = np.array([0.03253644,0.224043])

ranef_condvar_ml_drf_9 = np.array([0.3994872,0.02478884,0.02478884,1.195077]).reshape(2,2, order='F')

coef_reml_drf_9 = np.array([-1.272483,0.861814])

vcov_reml_drf_9 = np.array([0.02228589,0.001535598,0.001535598,0.0262125]).reshape(2,2, order='F')

cov_re_reml_drf_9 = np.array([0.5123204,0.08897376,0.08897376,3.341722]).reshape(2,2, order='F')

scale_reml_drf_9 = np.array([3.764058])

loglike_reml_drf_9 = np.array([-591.7188])

ranef_mean_reml_drf_9 = np.array([0.03239688,0.2230525])

ranef_condvar_reml_drf_9 = np.array([0.401762,0.02521271,0.02521271,1.203536]).reshape(2,2, order='F')

coef_ml_irf_9 = np.array([-1.277018,0.86395])

vcov_ml_irf_9 = np.array([0.02205706,0.001509887,0.001509887,0.02599941]).reshape(2,2, order='F')

cov_re_ml_irf_9 = np.array([0.5086816,0,0,3.312757]).reshape(2,2, order='F')

scale_ml_irf_9 = np.array([3.72105])

loglike_ml_irf_9 = np.array([-589.8628])

coef_reml_irf_9 = np.array([-1.276822,0.8640243])

vcov_reml_irf_9 = np.array([0.02225705,0.001517774,0.001517774,0.02623682]).reshape(2,2, order='F')

cov_re_reml_irf_9 = np.array([0.5107725,0,0,3.31152]).reshape(2,2, order='F')

scale_reml_irf_9 = np.array([3.762967])

loglike_reml_irf_9 = np.array([-591.7543])

coef_ml_drf_10 = np.array([-0.9419566,-0.02359824,1.085796])

vcov_ml_drf_10 = np.array([0.001963536,-0.0003221793,0.0001950186,-0.0003221793,0.002534251,0.0004107718,0.0001950186,0.0004107718,0.002580736]).reshape(3,3, order='F')

cov_re_ml_drf_10 = np.array([0.2040541,0.09038325,0.09038325,2.218903]).reshape(2,2, order='F')

scale_ml_drf_10 = np.array([0.2558286])

loglike_ml_drf_10 = np.array([-379.6591])

ranef_mean_ml_drf_10 = np.array([0.03876325,-0.725853])

ranef_condvar_ml_drf_10 = np.array([0.1988816,0.1872403,0.1872403,0.4052274]).reshape(2,2, order='F')

coef_reml_drf_10 = np.array([-0.9426367,-0.02336203,1.085733])

vcov_reml_drf_10 = np.array([0.002011348,-0.0003300612,0.0002002948,-0.0003300612,0.002589149,0.000418987,0.0002002948,0.000418987,0.002637433]).reshape(3,3, order='F')

cov_re_reml_drf_10 = np.array([0.2034827,0.09063836,0.09063836,2.219191]).reshape(2,2, order='F')

scale_reml_drf_10 = np.array([0.2630213])

loglike_reml_drf_10 = np.array([-386.0008])

ranef_mean_reml_drf_10 = np.array([0.03838686,-0.7240812])

ranef_condvar_reml_drf_10 = np.array([0.1983981,0.1865469,0.1865469,0.4100937]).reshape(2,2, order='F')

coef_ml_irf_10 = np.array([-0.9441033,-0.01755913,1.088568])

vcov_ml_irf_10 = np.array([0.001960114,-0.0003215658,0.0001944005,-0.0003215658,0.00253441,0.0004061179,0.0001944005,0.0004061179,0.002589158]).reshape(3,3, order='F')

cov_re_ml_irf_10 = np.array([0.2032228,0,0,2.192893]).reshape(2,2, order='F')

scale_ml_irf_10 = np.array([0.2553399])

loglike_ml_irf_10 = np.array([-380.162])

coef_reml_irf_10 = np.array([-0.9448257,-0.01722993,1.088557])

vcov_reml_irf_10 = np.array([0.00200783,-0.0003294349,0.0001996613,-0.0003294349,0.00258937,0.0004141667,0.0001996613,0.0004141667,0.002646242]).reshape(3,3, order='F')

cov_re_reml_irf_10 = np.array([0.2026653,0,0,2.193124]).reshape(2,2, order='F')

scale_reml_irf_10 = np.array([0.2625147])

loglike_reml_irf_10 = np.array([-386.5024])

coef_ml_drf_11 = np.array([-1.36971,0.1596278,0.8588724])

vcov_ml_drf_11 = np.array([0.0232326,0.00172214,0.002275343,0.00172214,0.02318941,0.0004755663,0.002275343,0.0004755663,0.02123474]).reshape(3,3, order='F')

cov_re_ml_drf_11 = np.array([0.3719096,0.332198,0.332198,1.120588]).reshape(2,2, order='F')

scale_ml_drf_11 = np.array([4.849781])

loglike_ml_drf_11 = np.array([-601.6432])

ranef_mean_ml_drf_11 = np.array([-0.4256917,-0.3907759])

ranef_condvar_ml_drf_11 = np.array([0.2987928,0.1992074,0.1992074,0.7477486]).reshape(2,2, order='F')

coef_reml_drf_11 = np.array([-1.370236,0.1597671,0.8585994])

vcov_reml_drf_11 = np.array([0.02351795,0.001749756,0.002301599,0.001749756,0.02346869,0.0004785668,0.002301599,0.0004785668,0.02149093]).reshape(3,3, order='F')

cov_re_reml_drf_11 = np.array([0.3680346,0.3324419,0.3324419,1.118623]).reshape(2,2, order='F')

scale_reml_drf_11 = np.array([4.922222])

loglike_reml_drf_11 = np.array([-604.5746])

ranef_mean_reml_drf_11 = np.array([-0.4168539,-0.3879533])

ranef_condvar_reml_drf_11 = np.array([0.2965372,0.2010191,0.2010191,0.7503986]).reshape(2,2, order='F')

coef_ml_irf_11 = np.array([-1.370117,0.1414964,0.8466083])

vcov_ml_irf_11 = np.array([0.02319951,0.001705996,0.002265252,0.001705996,0.02345623,0.000514879,0.002265252,0.000514879,0.02153162]).reshape(3,3, order='F')

cov_re_ml_irf_11 = np.array([0.4004789,0,0,1.108087]).reshape(2,2, order='F')

scale_ml_irf_11 = np.array([4.78776])

loglike_ml_irf_11 = np.array([-602.308])

coef_reml_irf_11 = np.array([-1.370663,0.1417561,0.8464232])

vcov_reml_irf_11 = np.array([0.02348548,0.001734072,0.002291519,0.001734072,0.02373715,0.0005177618,0.002291519,0.0005177618,0.02178966]).reshape(3,3, order='F')

cov_re_reml_irf_11 = np.array([0.3966454,0,0,1.106551]).reshape(2,2, order='F')

scale_reml_irf_11 = np.array([4.860342])

loglike_reml_irf_11 = np.array([-605.2274])