Generalized Linear Models
=========================


.. _glm_notebook:

`Link to Notebook GitHub <https://github.com/statsmodels/statsmodels/blob/master/examples/notebooks/glm.ipynb>`_

.. raw:: html

   
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[1]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="kn">from</span> <span class="nn">__future__</span> <span class="kn">import</span> <span class="n">print_function</span>
   <span class="kn">import</span> <span class="nn">numpy</span> <span class="kn">as</span> <span class="nn">np</span>
   <span class="kn">import</span> <span class="nn">statsmodels.api</span> <span class="kn">as</span> <span class="nn">sm</span>
   <span class="kn">from</span> <span class="nn">scipy</span> <span class="kn">import</span> <span class="n">stats</span>
   <span class="kn">from</span> <span class="nn">matplotlib</span> <span class="kn">import</span> <span class="n">pyplot</span> <span class="k">as</span> <span class="n">plt</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <h2 id="glm-binomial-response-data">GLM: Binomial response data</h2>
   <h3 id="load-data">Load data</h3>
   <p> In this example, we use the Star98 dataset which was taken with permission
    from Jeff Gill (2000) Generalized linear models: A unified approach. Codebook
    information can be obtained by typing: </p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[2]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="k">print</span><span class="p">(</span><span class="n">sm</span><span class="o">.</span><span class="n">datasets</span><span class="o">.</span><span class="n">star98</span><span class="o">.</span><span class="n">NOTE</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
   ::
   
       Number of Observations - 303 (counties in California).
   
       Number of Variables - 13 and 8 interaction terms.
   
       Definition of variables names::
   
           NABOVE   - Total number of students above the national median for the
                      math section.
           NBELOW   - Total number of students below the national median for the
                      math section.
           LOWINC   - Percentage of low income students
           PERASIAN - Percentage of Asian student
           PERBLACK - Percentage of black students
           PERHISP  - Percentage of Hispanic students
           PERMINTE - Percentage of minority teachers
           AVYRSEXP - Sum of teachers&apos; years in educational service divided by the
                   number of teachers.
           AVSALK   - Total salary budget including benefits divided by the number
                      of full-time teachers (in thousands)
           PERSPENK - Per-pupil spending (in thousands)
           PTRATIO  - Pupil-teacher ratio.
           PCTAF    - Percentage of students taking UC/CSU prep courses
           PCTCHRT  - Percentage of charter schools
           PCTYRRND - Percentage of year-round schools
   
           The below variables are interaction terms of the variables defined
           above.
   
           PERMINTE_AVYRSEXP
           PEMINTE_AVSAL
           AVYRSEXP_AVSAL
           PERSPEN_PTRATIO
           PERSPEN_PCTAF
           PTRATIO_PCTAF
           PERMINTE_AVTRSEXP_AVSAL
           PERSPEN_PTRATIO_PCTAF
   
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p>Load the data and add a constant to the exogenous (independent) variables:</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[3]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">data</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">datasets</span><span class="o">.</span><span class="n">star98</span><span class="o">.</span><span class="n">load</span><span class="p">()</span>
   <span class="n">data</span><span class="o">.</span><span class="n">exog</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">add_constant</span><span class="p">(</span><span class="n">data</span><span class="o">.</span><span class="n">exog</span><span class="p">,</span> <span class="n">prepend</span><span class="o">=</span><span class="bp">False</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p> The dependent variable is N by 2 (Success: NABOVE, Failure: NBELOW): </p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[4]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="k">print</span><span class="p">(</span><span class="n">data</span><span class="o">.</span><span class="n">endog</span><span class="p">[:</span><span class="mi">5</span><span class="p">,:])</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
   [[ 452.  355.]
    [ 144.   40.]
    [ 337.  234.]
    [ 395.  178.]
    [   8.   57.]]
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p> The independent variables include all the other variables described above, as
    well as the interaction terms:</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[5]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="k">print</span><span class="p">(</span><span class="n">data</span><span class="o">.</span><span class="n">exog</span><span class="p">[:</span><span class="mi">2</span><span class="p">,:])</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
   [[  3.43973000e+01   2.32993000e+01   1.42352800e+01   1.14111200e+01
       1.59183700e+01   1.47064600e+01   5.91573200e+01   4.44520700e+00
       2.17102500e+01   5.70327600e+01   0.00000000e+00   2.22222200e+01
       2.34102872e+02   9.41688110e+02   8.69994800e+02   9.65065600e+01
       2.53522420e+02   1.23819550e+03   1.38488985e+04   5.50403520e+03
       1.00000000e+00]
    [  1.73650700e+01   2.93283800e+01   8.23489700e+00   9.31488400e+00
       1.36363600e+01   1.60832400e+01   5.95039700e+01   5.26759800e+00
       2.04427800e+01   6.46226400e+01   0.00000000e+00   0.00000000e+00
       2.19316851e+02   8.11417560e+02   9.57016600e+02   1.07684350e+02
       3.40406090e+02   1.32106640e+03   1.30502233e+04   6.95884680e+03
       1.00000000e+00]]
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <h3 id="fit-and-summary">Fit and summary</h3>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[6]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">glm_binom</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">GLM</span><span class="p">(</span><span class="n">data</span><span class="o">.</span><span class="n">endog</span><span class="p">,</span> <span class="n">data</span><span class="o">.</span><span class="n">exog</span><span class="p">,</span> <span class="n">family</span><span class="o">=</span><span class="n">sm</span><span class="o">.</span><span class="n">families</span><span class="o">.</span><span class="n">Binomial</span><span class="p">())</span>
   <span class="n">res</span> <span class="o">=</span> <span class="n">glm_binom</span><span class="o">.</span><span class="n">fit</span><span class="p">()</span>
   <span class="k">print</span><span class="p">(</span><span class="n">res</span><span class="o">.</span><span class="n">summary</span><span class="p">())</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
                    Generalized Linear Model Regression Results                  
   ==============================================================================
   Dep. Variable:           [&apos;y1&apos;, &apos;y2&apos;]   No. Observations:                  303
   Model:                            GLM   Df Residuals:                      282
   Model Family:                Binomial   Df Model:                           20
   Link Function:                  logit   Scale:                             1.0
   Method:                          IRLS   Log-Likelihood:                -2998.6
   Date:                Wed, 20 May 2015   Deviance:                       4078.8
   Time:                        21:52:38   Pearson chi2:                 4.05e+03
   No. Iterations:                     7                                         
   ==============================================================================
                    coef    std err          z      P&gt;|z|      [95.0% Conf. Int.]
   ------------------------------------------------------------------------------
   x1            -0.0168      0.000    -38.749      0.000        -0.018    -0.016
   x2             0.0099      0.001     16.505      0.000         0.009     0.011
   x3            -0.0187      0.001    -25.182      0.000        -0.020    -0.017
   x4            -0.0142      0.000    -32.818      0.000        -0.015    -0.013
   x5             0.2545      0.030      8.498      0.000         0.196     0.313
   x6             0.2407      0.057      4.212      0.000         0.129     0.353
   x7             0.0804      0.014      5.775      0.000         0.053     0.108
   x8            -1.9522      0.317     -6.162      0.000        -2.573    -1.331
   x9            -0.3341      0.061     -5.453      0.000        -0.454    -0.214
   x10           -0.1690      0.033     -5.169      0.000        -0.233    -0.105
   x11            0.0049      0.001      3.921      0.000         0.002     0.007
   x12           -0.0036      0.000    -15.878      0.000        -0.004    -0.003
   x13           -0.0141      0.002     -7.391      0.000        -0.018    -0.010
   x14           -0.0040      0.000     -8.450      0.000        -0.005    -0.003
   x15           -0.0039      0.001     -4.059      0.000        -0.006    -0.002
   x16            0.0917      0.015      6.321      0.000         0.063     0.120
   x17            0.0490      0.007      6.574      0.000         0.034     0.064
   x18            0.0080      0.001      5.362      0.000         0.005     0.011
   x19            0.0002   2.99e-05      7.428      0.000         0.000     0.000
   x20           -0.0022      0.000     -6.445      0.000        -0.003    -0.002
   const          2.9589      1.547      1.913      0.056        -0.073     5.990
   ==============================================================================
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <h3 id="quantities-of-interest">Quantities of interest</h3>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[7]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="k">print</span><span class="p">(</span><span class="s">&#39;Total number of trials:&#39;</span><span class="p">,</span>  <span class="n">data</span><span class="o">.</span><span class="n">endog</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span><span class="o">.</span><span class="n">sum</span><span class="p">())</span>
   <span class="k">print</span><span class="p">(</span><span class="s">&#39;Parameters: &#39;</span><span class="p">,</span> <span class="n">res</span><span class="o">.</span><span class="n">params</span><span class="p">)</span>
   <span class="k">print</span><span class="p">(</span><span class="s">&#39;T-values: &#39;</span><span class="p">,</span> <span class="n">res</span><span class="o">.</span><span class="n">tvalues</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
   Total number of trials: 807.0
   Parameters:  [ -1.68150366e-02   9.92547661e-03  -1.87242148e-02  -1.42385609e-02
      2.54487173e-01   2.40693664e-01   8.04086739e-02  -1.95216050e+00
     -3.34086475e-01  -1.69022168e-01   4.91670212e-03  -3.57996435e-03
     -1.40765648e-02  -4.00499176e-03  -3.90639579e-03   9.17143006e-02
      4.89898381e-02   8.04073890e-03   2.22009503e-04  -2.24924861e-03
      2.95887793e+00]
   T-values:  [-38.74908321  16.50473627 -25.1821894  -32.81791308   8.49827113
      4.21247925   5.7749976   -6.16191078  -5.45321673  -5.16865445
      3.92119964 -15.87825999  -7.39093058  -8.44963886  -4.05916246
      6.3210987    6.57434662   5.36229044   7.42806363  -6.44513698
      1.91301155]
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p>First differences: We hold all explanatory variables constant at their means and manipulate the percentage of low income households to assess its impact on the response variables: </p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[8]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">means</span> <span class="o">=</span> <span class="n">data</span><span class="o">.</span><span class="n">exog</span><span class="o">.</span><span class="n">mean</span><span class="p">(</span><span class="n">axis</span><span class="o">=</span><span class="mi">0</span><span class="p">)</span>
   <span class="n">means25</span> <span class="o">=</span> <span class="n">means</span><span class="o">.</span><span class="n">copy</span><span class="p">()</span>
   <span class="n">means25</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">=</span> <span class="n">stats</span><span class="o">.</span><span class="n">scoreatpercentile</span><span class="p">(</span><span class="n">data</span><span class="o">.</span><span class="n">exog</span><span class="p">[:,</span><span class="mi">0</span><span class="p">],</span> <span class="mi">25</span><span class="p">)</span>
   <span class="n">means75</span> <span class="o">=</span> <span class="n">means</span><span class="o">.</span><span class="n">copy</span><span class="p">()</span>
   <span class="n">means75</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">=</span> <span class="n">lowinc_75per</span> <span class="o">=</span> <span class="n">stats</span><span class="o">.</span><span class="n">scoreatpercentile</span><span class="p">(</span><span class="n">data</span><span class="o">.</span><span class="n">exog</span><span class="p">[:,</span><span class="mi">0</span><span class="p">],</span> <span class="mi">75</span><span class="p">)</span>
   <span class="n">resp_25</span> <span class="o">=</span> <span class="n">res</span><span class="o">.</span><span class="n">predict</span><span class="p">(</span><span class="n">means25</span><span class="p">)</span>
   <span class="n">resp_75</span> <span class="o">=</span> <span class="n">res</span><span class="o">.</span><span class="n">predict</span><span class="p">(</span><span class="n">means75</span><span class="p">)</span>
   <span class="n">diff</span> <span class="o">=</span> <span class="n">resp_75</span> <span class="o">-</span> <span class="n">resp_25</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p>The interquartile first difference for the percentage of low income households in a school district is:</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[9]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="k">print</span><span class="p">(</span><span class="s">&quot;</span><span class="si">%2.4f%%</span><span class="s">&quot;</span> <span class="o">%</span> <span class="p">(</span><span class="n">diff</span><span class="o">*</span><span class="mi">100</span><span class="p">))</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
   -11.8753%
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <h3 id="plots">Plots</h3>
   <p> We extract information that will be used to draw some interesting plots: </p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[10]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">nobs</span> <span class="o">=</span> <span class="n">res</span><span class="o">.</span><span class="n">nobs</span>
   <span class="n">y</span> <span class="o">=</span> <span class="n">data</span><span class="o">.</span><span class="n">endog</span><span class="p">[:,</span><span class="mi">0</span><span class="p">]</span><span class="o">/</span><span class="n">data</span><span class="o">.</span><span class="n">endog</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="mi">1</span><span class="p">)</span>
   <span class="n">yhat</span> <span class="o">=</span> <span class="n">res</span><span class="o">.</span><span class="n">mu</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p>Plot yhat vs y:</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[11]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="kn">from</span> <span class="nn">statsmodels.graphics.api</span> <span class="kn">import</span> <span class="n">abline_plot</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[12]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">fig</span><span class="p">,</span> <span class="n">ax</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">subplots</span><span class="p">()</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">scatter</span><span class="p">(</span><span class="n">yhat</span><span class="p">,</span> <span class="n">y</span><span class="p">)</span>
   <span class="n">line_fit</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">OLS</span><span class="p">(</span><span class="n">y</span><span class="p">,</span> <span class="n">sm</span><span class="o">.</span><span class="n">add_constant</span><span class="p">(</span><span class="n">yhat</span><span class="p">,</span> <span class="n">prepend</span><span class="o">=</span><span class="bp">True</span><span class="p">))</span><span class="o">.</span><span class="n">fit</span><span class="p">()</span>
   <span class="n">abline_plot</span><span class="p">(</span><span class="n">model_results</span><span class="o">=</span><span class="n">line_fit</span><span class="p">,</span> <span class="n">ax</span><span class="o">=</span><span class="n">ax</span><span class="p">)</span>
   
   
   <span class="n">ax</span><span class="o">.</span><span class="n">set_title</span><span class="p">(</span><span class="s">&#39;Model Fit Plot&#39;</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">set_ylabel</span><span class="p">(</span><span class="s">&#39;Observed values&#39;</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">set_xlabel</span><span class="p">(</span><span class="s">&#39;Fitted values&#39;</span><span class="p">);</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   
   
   <div class="output_png output_subarea ">
   <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAYYAAAEZCAYAAACTsIJzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz
   AAALEgAACxIB0t1+/AAAIABJREFUeJztnXl8VNXZ+L9PCJFUQAiJqOyiFSkKg1Zp0SYuLNqWvkpb
   f65ULWjrVhsVceXVUF6tuGAXt6pUq9bWpbgl4BLcFxQRF1oXXABrCWhViAbI8/vj3klmnzs3c2cm
   k+f7+ZxPZu6ce865VzzPeZ7nnOcRVcUwDMMwwpTkewCGYRhGYWGCwTAMw4jCBINhGIYRhQkGwzAM
   IwoTDIZhGEYUJhgMwzCMKEwwGF0SERkqIq0ikvb/ARH5mYg8lWH7+4vISv8jjGrL81gNIxvYPzSj
   4BGR90XkaxHpF3N9mTthDs7j2MKT9hcRZZmqPqWqIyLqvS8iB6Zopyainc9FZKWI/MzHeGaLyG0+
   H8cwACjN9wAMwwMKvAccCfwOQET2AMrd3wqB7VS1NcXvCkiaNtao6iAAEfkR8HcReR74KktjNAxP
   mMZgdBZuB46L+D4N+DMRk62IbCcifxaR/7gr9PNFRNzfSkTkChFZJyLvAt+PbNy9908islZEVovI
   pR0x3bgawEfu59uAwcADrkZwVrr7VfUfwKfAyARt7yQiC0VkvYi8LSI/d69PBmYBR4Q1F7/jN7o2
   JhiMzsLzQG8RGSEi3YAjcIRFJNcCvYBhQDWOIDne/W0GjjAYA+wN/JhobeNWoAUYDoSAicDPMxhf
   Um1AVY8FPgR+oKq9VPWKlA05QuwwoA+wIkGVu9z2dsR5jt+IyAGqWg/8BrjL7SeUwfgNow0TDEZn
   4jacyX4C8CawJvxDhLCYpaobVfUDYB5wrFvlp8BVqrpGVT/FmUDD2kR/4BDgTFVtVtV1wNXA/8tg
   bE0i8qlbfu3z+XYSkU+BdcCFwDGq+nZkBREZBHwXmKmqLaq6HLiJdm1KSG+yMoyUmI/B6CwojmB4
   CkcjiDIjAZVAd+CDiGsfAgPczzsCH8X8FmaIe+/HruUJnEVTZJ109Iv0MYhITQb3hlkb9jGkYCdg
   g6pujLj2IY4WZBhZwQSD0WlQ1Q9F5D2c1f0JMT83AZuBocBb7rXBwGr388fudyJ+C/MR8DUxk3uW
   yZaTfC1QISI9VfVL91rkcxaKM97oxJgpyehsnAgcqKrNkRdVdStwNzBHRHqKyBDgTNr9EHcDp4vI
   ABHpC5wbce/HwCLgShHp5dr4h4vI97I47k9w/BcdQlU/Ap4F5orINiKyJ46QDD/nv4GhEqH6GEam
   mGAwOhWq+p6qvhJ5KeLzacBGnK2tTwF/AW5xf7sRaACWA0uBe2LuPQ4ow/FdbAD+BuwQ0UeqlXiy
   3yKvzwUuSOOD8NrHkTia0VrgXuAiVX3c/e1v7t/1IrI0RXuGkRQJMlGPiNyMsxPkP6q6R5I683FM
   A5uAn6mqbbEzDMPII0FrDLcAk5P9KCKHAruo6q442wn/GPB4DMMwjDQEKhhU9SmcQzrJmAIscOu+
   APRxtw4ahmEYeSLfPoYBRG8hXA0MzNNYDMMwDPIvGCD+MI5ttzMMw8gj+T7HsAaIPNAzkIjTrGFE
   xISFYRiGD1Q1463L+dYYFuIe5ReRccBnqvpJooqqakWViy++OO9jKJRi78Lehb2L1MUvgWoMInIn
   TjCzSjfS5MU4oQdQ1etV9WEROVRE3sHZf3588tYMwzAML2zZAvPn+78/UMGgqkd6qHNqkGMwDMPo
   Srz+OpxwAvTs6b+NfJuSjAypqanJ9xAKBnsX7di7aKervovNm6GuDmpq4MQT4bHH/LcV6MnnbCEi
   2hnGaRiGkQ+WL4fjj4eqKrjxRhjshogUEbQTOp8NwzAMn7S0wMUXw8EHw6mnQn19u1DoCPnermoY
   hmH44OWXHS1h8GB49VUYMCD9PV4xjcEwDKMT8fXXcN55cOihcM458MAD2RUKYBqDYRhGp+GFFxwt
   YcQIx6+www7p7/GDCQbDMIwCp7kZLroIbrsNrrkGfvpTCDIVk5mSDMMwCphnnoExY+CDD+C11+CI
   I4IVCmAag2EYRkGycSOcfz7cfTdcey1MnZq7vk1jMAzDiKGhoYGJE6cyceJUGhoact5/YyOMHg1N
   TbBiRW6FAtgBN8MwjCgaGho47LBpNDdfBkB5+Uzuu28BkyZNCrzvL76Ac8+Ff/wD/vAHmDKlY+3Z
   ATfDMIwsMG/eDa5QmAY4AmLevBsC7/fRR2GPPWDTJkdL6KhQ6AjmYzAMw8gj//0vnHUWNDTA9dfD
   IYfke0SmMRiGYURRWzuD8vKZOOnoF1BePpPa2hmB9PXII46WIOJoCYUgFMB8DIZhGHE0NDS0mY9q
   a2dk3b/w6adw5pmwZIkT9O7gg7PafBt+fQwmGAzDMHLIwoXwi1/AYYfB3LnQq1dwffkVDOZjMAzD
   yAHr18PppzthLe64A6qr8z2i5JiPwTAMIwkNDQ2MHVtDv367MHbsfr7PNNxzD4waBdtv78Q4KmSh
   AGZKMgzDSEhDQwNTphxLS8tv3StnUVa2hYUL7/Lsc/jPf+CUU5xQFjffDOPHBzfeRNg5BsMwOh3Z
   PmGcqD2/fcybd4MrFKa55QpaWkZ4OtOgCnfe6ew4GjbMyZeQa6HQEczHYBhGXog9Yfz009M6dMI4
   UXvnn38ac+Zcm7U+vPDxx45z+e23HUfzvvsG1lVwqGrBF2eYhmEUExMmHK5wqzrra1W4VSdMODyr
   7VVUDPfdR319vZaVVbn336pQqWVlfbS+vj5h/dZW1QULVKuqVM87T/Wrr3w/StZw586M51zTGAzD
   MBIwadIkFi68jVmz5vLBB6sZMmQ35s69MKG2sWYNzJgBq1c7eZfHjs3DgLOJH2mS64JpDIZRdNTX
   12t5ef+2FXl5ef+kq3G/7dXV1WW1j1haW1Vvukm1slJ19mzVr7/OWtNZAZ8ag+1KMgwjb2T7hHGi
   9oI6xfzhhzB9uhMa+5ZbYM89s9JsVrGTz4ZhGDmgtRVuuAEuuAB+/Ws4+2zo3j3fo0qMnXw2DMMI
   mPfeg5//HL780olz9K1v5XtEwWDnGAzDMNLQ2uqk19xnHycC6rPPFq9QANMYDMMwUvL223DiibB1
   KzzzDOy2W75HFDymMRiGYSRg61a48kr4znfg8MPhySe7hlAA0xgMwzDieOstOOEEKCuD55+HXXbJ
   94hyi2kMhmF0KrIRXylZG1u2wP/9H+y/PxxzDDzxRNcTCmDbVQ3D6ETExkMqL5+ZceyjZG0MGDCJ
   44+H7bZzsqoNGxbII+QUO8dgGEbRM3HiVBYvnoIT7RRgARMmLGTRons60MafGT68O//975HMmeMc
   WpOMp9LCxM4xGIZh+OL7fPbZByxbBoMG5XsshYH5GAzDCJxUfoFMfAa1tTMoL58JLAAWUF4+k9ra
   GRmNpbZ2Bj16XAAsB5rp3v0Sbr99nQmFSPwEWPJagMnASuBtYGaC3yuBeuBV4HXgZ0nayUY8KcMw
   8kCqYHl+AunV19frhAmH64QJh/sKiPfSS6pDh36ulZUv6v77n5jVoHqFBj6D6AUpFLoB7wBDge7u
   5L97TJ3ZwFxtFxLrgdIEbQXz1gzDCJxUeReynZMhFc3NqjNnqm6/verttzuRUYsdv4IhSFPSPsA7
   qvq+qm4G7gJ+FFPnY6C3+7k3sF5VtwQ4JsMwuiDPPQehkHOKeflyOPro4nEwB0GQzucBwEcR31cD
   sUnubgQeF5G1QC/gpwGOxzCMLJBpGOva2hk8/fQ0mpud745fYEHa37LBpk1w4YXwl7/A/Pnwk5+Y
   QPBCkILBy/7S84BXVbVGRIYDi0VktKp+EVtx9uzZbZ9ramqoqanJ1jgNw/CInzzNkyZN4r77FkQI
   k/b6qX7rKE895cQ4GjsWVqyAqqqsNFvQNDY20tjY2OF2AjvHICLjgNmqOtn9PgtoVdXLIuo8DMxR
   1Wfc74/hOKmXxrSlQY3TMAzvZOMcQTKylVBn40Y47zz429/g97+Hww7r8NA6LX7PMQTpY1gK7Coi
   Q0WkDDgCWBhTZyVwMICI9Ad2A94LcEyGYeQYL9tRw5rI4sVTWLx4CocdNs1XuIsnnnAyqX36Kbz+
   etcWCh3Cj8faawEOAf6JsztplnvtJOAkbd+J9ADOhuIVwFFJ2smal94wDP9kur3Ua/2O7k76/HPV
   k09WHThQ9YEHfD1aUYLPXUmBnnxW1UeAR2KuXR/xuQn4YZBjMAwjHr9mm0x8Ag0NDRx11Ck0Nw8D
   dgAm0dwM8+bdkDU/AsCiRTBjBhx0kONL6NMna013XfxIk1wXTGMwjKzh51BZR/uA/gr1STUBP2P6
   7DPVE09UHTRI9ZFHsjr8ooFC1BgMwyg85s27wd1V5DiQk63iY7WK8L0A1dVjWbLklbbfYu+N7cNh
   NuXlqxJuR810d9LDD8NJJ8H3v+/4Enr3TlrV8IEJBsMw4ojdlrpkybHAZlpargZWsHjx5cB8wNuW
   VYCKinXccUfyepMmTUrbxoYNcOaZzlbUW291zEdGAPhRM3JdMFOSUeR0NP5Ppn2lM9skcgbDOPdz
   ekexlz6SPXOy6/ffr7rTTqqnnqr6xRcBvJgiBDMlGUbnxM+hsY6QzGwTaTpqalofSB9hkj0zEHf9
   1lvv4L77DmTpUrjzTvje9zo0NMMLfqRJrgumMRhFTLYDyfnRPmJX+GVlfbSsrCrie5WWlm7rag0j
   FHp3yHmd7Jnjrz+hZWUb9MwzVTdu9PM2ujaYxmAYhl/tI9ZZ3NICodCNVFY6Z1Krq8/gkkuuAU4G
   oLT0DPbY4xYqK/tlNYxFvKYSYsiQ87jyyt9npX3DGyYYDCPPZDOQnNcdR16orOzfFupi4sSptLT8
   tq3dLVugstJ/KIxkz3zuuZcCz+LE03wHOJSePYf46sPwjwkGw8gzQQaS80JDQwNNTZ9QUnImra3O
   tWxHOY0l0TPvueck1qzpCwwGLgbeBY6gsnJVYOMwkuDH/pTrgvkYDMMTHQtZUaslJf00FKpOuIMo
   qENxra2qt96qWlWleuSR72iPHoMCPXzXlcCnjyGw6KrZxKKrGoZ3Mgl3kUm01GxFP41k9WonnMXa
   tXDLLU4ynSD66ar4ja5qgsEwujBBhtFOhSrcdJMTHvu00+Dcc6GsLNAuuyR+BYP5GAyjCxN0BrVE
   vP8+TJ/uhMZ+/HFYu7aBH/zANIRCIsh8DIZhFDhhJ/CECQuZMGFhoAfrWlvhj3+EvfeGAw+E5593
   hEI28jAY2cVMSYZhBM577zlpNpub4eabYeRI53q+TFldhULM4GYYRg7wkiEtX7S2wvz5sM8+8IMf
   wDPPtAsFo3AxH4NhdGJyHWcpE/71LzjhBBCBZ5+Fb34zvk4+fBxGekxjMIxOTPRJZ0dAhLd6BoEX
   7WTrVrjiCvjud+GnP4UlSxILBcitj8PwjmkMhmF4wot28uabjpZQXg4vvADDh6dv10seBiO3mPPZ
   MAoAv4e6Yifr8vKZvlfd6cbgOIqHAeEQFcMIhZ6nsrI/ra0lDBo0hwce+CaXXupkVysxe0Te8et8
   znu4Cy8FC4lhFDEdDTeRjSQ/XsYQCo1XqIzI41yhIj0V7ldYpyUljXrrrY2++jeCAZ8hMTKdoLsB
   vf101JFigsEoZpLlJsjWhD9hwuEaClVrKDQ+YVv19fVaUTHczbVQnzQnRChUHTPOBQo3uJ9XKIzT
   iorhFtuogPArGNL6GETkTuAkYCvwErCdiFyjqpdnrJ4YhuGJpqZPOrzbKNbMBGcB46Laiq8zDUi8
   K6iysl/MlR8AK4ElwBHAZWzY4GRgMydyJyed5ACWu3+PBuYB3YEVfqSQ34JpDEYnxcuqP1H2tF69
   BnvK6paq/cR5mw+PaitZbudkOZqdyKfLFTZpt26navfufVxNI3sZ6IzsgU+NwYt7qFREugP/Azyg
   qpsB8wQbRhrCq/F04R4it2yGQjcC3fnii50ybn/y5CPYZZc9O3zIraJiXcIVf0XFJPr3f4Oqqq/Y
   f//TeOihH/DAA3dRUbGuQ/0ZBUg6yQGcDqwBHsE59zAUeMqPFPJbMI3B6IT4yeXcfk+9QmpncLLV
   fllZldbX18dpIo7juDaqLS9O5+Zm1XPOUd1+e9W//MXJnxBJkLkajI5BUD4GVZ0PzA9/F5EPgAOy
   LqEMw4hhF+ACevUS/vY3rzb7nWhpOZl5825g0aJ72rKkObmUd6OyclVUhrh02eOefdY5lzBqFLz2
   GvTvH99jvjPQGQGQTnIAOwB/Aurd7yOBE/1IIb8F0xiMToiflXR9fb2WlfWJ2hYa1gAS1S0p6Ruh
   EfR3NY2O2/g3blQ980zVHXZQ/dvfOtSUkUcIarsqUI+z5eA193t34HU/nfktJhiMbJONraBB9RO/
   LTT5RF9XV6cifRUqFGpTChKvLFmiussuqkceqbpuXXaeycgPQQqGpe7fZRHXXvXTmd9igsHIJoVu
   E8/UN1FfX6+h0HitqBieMF+zV774QvXUU1V32kn1/vuT91XI786IJkjB0Aj0CwsGYBywxE9nfosJ
   BiOb+HEK5xK/k2+mK/nI+nPnvqhDh6oed5zq+vXJ7yn0d2dE41cweAmiVws8AOwsIs8CVcCPfTk0
   DMNIix9nbqbht9vrXwnsxaOP9mT27Je56KK9svosRifFi/TA8SuMckt3PxKoIwXTGIwskm9zSBA2
   +kxX8k79RRH1/+Jp5Z/vd2dkBgGGxJiGc6AtHKFvrBux78+BSCrDCJh8bq8shMQ6n30Gb7zxS2Df
   iKubaWpaz8SJU4HkEV5ta2rXIG3YbRH5He0nnXsABwGvqGrOzEkWdtsoFoLKcew1/PaDD8LJJ0Mo
   9CGPPnoAX311EQBlZWcDm2lpuTrl/Ubnwm/YbS8H3E6N6agP8FePg5oMXI0TlfUmVb0sQZ0a4Coc
   c1WTqtZ4adswjHbCK/lZsy7lgw/+zZAhI6J+37ABzjjDObB2221wwAGDaWj4Q9vKv6npmyxbNp2w
   wGpudrLDmWDomvjJ4LYJGJaukoh0A34HHIwTUuMlEVmoqm9F1OkD/B6YpKqrRaTSx3gMo9MQdI7j
   lSvfobk5Osrppk2TOOUU+MlPnNPL227r1I3MnBY2IRkG4Gm76gMR5SGc9E2XebjvO7inpd3v5wLn
   xtT5JXCJh7ay6Y8xjLySzvns1zkd74C+U/v3f0p33VX1qafSj8mrU7krHHArlmckwHMMNRFlP2CQ
   p4adLa03Rnw/Brg2ps5VOFrFE8BS4NgkbQX24gyjkOjIrp94wbBJhwy5Xzdu9N53piHCi3FXUjE9
   Y2CCwW8BpnoQDL8DngXKcQ7R/QvYNUFbgbw0w+gIhbDtNHY8PXp8S+F9hU91m20OiRpXNsYbP75a
   ragY3ulX1pEU0yE+v4IhqY9BRL4ked4FVdXeKSxU4PgVBkV8HwSsjqnzEY7DuRloFpEngdHA27GN
   zZ49u+1zTU0NNTU1abo3jOAohG2nkahCU9MkevSoYYcdHmLnne/mnHPOaBvPnDlzuOiiq2ht3RUY
   n6XxNgAL2LDhChYvzv87MKCxsZHGxsaON+RHmngpOI7td3HyN5QBrwK7x9QZATyKs2vpG8AKYGSC
   toIQpkaRk+0VfWR7odD4iFVlvcIoLS3dXkOh8R3qy48ZY80a1R/+UHXUKNWXXkrcZnwU1tq2VXAm
   7yl6fMWZuc1MSRmYkoDtgcHh4vGeQ4B/Au8As9xrJwEnRdQ5C3jDFQqnJ2knsBdnFCfZ+p87PGmG
   QuO1rKyqrT1noq11hUJVxKRbqWVlfTosHLxM1K2tqjffrFpVpXrhhapffZW4XrKEPuH2/YQGnzDh
   cK2oGF6UgkHVnM9eJvcpOKadjTg7klqBN/x05reYYDAyJRt24nSr45KSfgmvhyfdZG1mY8L58EPV
   yZNVx4xRXbYsdd1E76KkpF/bWDri0yiWlXWx4lcweMn5XIez9fRfqjoM5+TzC5mYqwyjkGloaGDi
   xKlMnDg1Kl/yvHk3uD6EaUB8DubRo0clzXfsZEyLZs6cORx66JEpc0AnG0v4t1Cohm23PZ+dd/6U
   qqq3efFFGDMm9fPV1s6gvHwmsABYQEnJmVxyyZkd9gVE5qqeMGGh+ReKiXSSA3jZ/bsc6OZ+fs2P
   FPJbMI3ByBCvq9lU9aJX0/UamVUtXM/JuBZtSoLeceYkx87fL+XqPNVY6uvrtXv3sQprFNYp7J+R
   ySqZpmKr/uKGAE1JjwK9cLaW3oWT//lZP535LSYYDD94MdukMqXETpplZX00FKpOOLn26jVIYaDC
   eE2UXtPpJ9LsVK8wTisqhicRRO1j2bpVdeDAqxSaPZmsOpKXwYRCceFXMHgJifEj4CvgTOBooDfw
   v9nSWAyjoaEhIlpn4qiefogM+eD3/uhIoncljTg6bty344LjxTMemImzz2IBcEVU6IpEbNq0A2PG
   bGD16vE4WXb/J+r3WJOVn220HX1PRhGSTnLgJOoZ4EfqZKtgGkPRkm9TRkf6j0ypOXz4GC0r6xO1
   a2natGltK/G6ujq3n1pXs4jXDOrq6hS2dTWL72pJybnau/fXOmDAfIWRCu3tJzNZtWsd9QqOlhIK
   jQ/q9RkFDgGakmbjbCd9GjgV6O+no44UEwzFSyGcMvVjSnF8C32i/A7duvVVkV7uxD5VoXeUwKmr
   q0u6zTMUqtbS0u0UKhTuVfhE4Wk97LBzIs4g1Cr0VeccQiqTVa1bp11ImYmoaxKYYND2yXk0MAfn
   XMJjfjrzW0wwFC+FIBj8EO8z0Da7v/PZu++ivLy/e2DuuwovRdxzlkJ4S2xdmwYAI5K+Ly9Obq+Y
   76Hz41cweNmuGuY/wL+B9Th5nw2jw8RupXTCUM/I97ACJey7CIVuoWfPWZSWlvHPf24D/AEY4NZq
   AG4F5gEnA5fjRLs/GViLcy40/n1NmjSJ0aNHdXiMYV9Fqq21RhGTTnLghMZuBN7EcTrHhawIumAa
   Q1HTGVem9fX1Wlq6bZTdv7S0n4r0jDD7RJuSYncyOaaoHRRedncc/UphuwjNI1YbObztc6rAddnw
   23RWTc6IhgB3JQ0CfqWqrwYjmoyuTlC7YoLa7QSwdOlStm7dBidSzBWIrGX27F9zzz2PsGzZdTgH
   4s4BrqOiYh133OHsDAqP6eWXl9PSciBwA05gYYAxlJR0o7W1lpISaG1N3v9ee41Omg7U8jIbHcaP
   NMl1wTQGI0OC2O0UGTdJJHKHUHtQunSH1Jzf/qywTJ2Dak8n8FE4ZyaiD871drWQ3OzcyvduMSM7
   UGj5GLJZTDAYmdJRU0jkVtRQqFrr6urciTrxdtPIw2bJTGPOmBZG3Pe/Gr0FtcLdaeS0GQqNj9ru
   GhZKiQ7ZBUFnNPEZ0ZhgMIwIOhocLnYrqsh26mxB7Z/Q/h8OShfbTnhiXbiwQYcN+7vCpqj7evbc
   USsqhrsnp2tTjtdW8UammGAwjAhiYxiJ9PGcKyH5VtSwplCvkecEoLdOmzYtrv/2SfxBFVmpu+++
   Unv02MMVAOO0pKSf1tXVJaifeNI3h7CRKX4FQ5AZ3Awjz2wGrgNAtZVly8a1hZ/w54zd5P6dhLO9
   djawDpjO2rWromo6kVmvwMloC6pLGDBgPkcffQQXXTSP1taraG2FOXNmsvfee5vD2Cgs0kkOnLDb
   v8SJkdQb+AVwqR8p5LdgGoORIYmT0xzetspOZT9PZEoqLe2nJSXbxPgEqtQ5eBYdDE9V9dvfnqXw
   Wdzq3vIfGLmEAENixIXYTnQtyGKCwciUVIIhFKpOO8HGOp/DIbZDofHaq9dgFenr+hyiQ3Hff/9i
   /dWvVPv2/UrLyk6I68OrYEgVJtscwoZXghQMz+How93ccjQWdtsoANKt+iMnf2cCr40IPxE5OU/V
   kpLttLR0ex0+fEzSCTeyv7q6ugQxjx7R8vK1etRRquvWJR5fulV/WPhE5mg2zcDwS5CCYRiwEGhy
   yz+AoX4681tMMBixeDGrtJ87qG7b+hm/aq9T+EbUyr+kpFfcltD43AxV7k6iSMGwUUeP/o2nscdq
   I9F9xDu+zcls+CEwwVAIxQSDEUv2bPUDNT6BTnymtsSmqRGueepzhbe1R49vek6Mk0iotfcR31co
   VO2pXTMzGZH4FQxpg+iJyG4i8piIvOF+31NELvDv7jaM/BKZq7i0tCXm1xuAK3AS7jgJb8I7haLp
   jpPM8HbgIioqJnP//fM97SSKziWdqI8ZwNmEAwvCWbzxxvKUQew6GvQuVa5powuSTnIATwL7Asvc
   7wK84UcK+S2YxmDEkK0dOk5ynEhTUmIzTnR/ixQ+UvjQl6knmbYT3ccIdyyHa6K8C17b9ILtdipe
   CNDHsNT9uyzi2qt+OvNbTDAUP36T5aS7x0uduro67dlzRy0t3V4rKnaIioMUOUmef/4VWlb2NxX5
   UEtKvu97Ik0XTylZMp9kE319fb1bf5xGhtTwKhjs4FzxEqRgeATYJUJj+DHwiJ/O/BYTDMWN1xWr
   nyT3mayE2+vHn0yePftlhTUKbyncrmVlVVEObT/PnOpZMnknyXZfeR2XCYbiJUjBMBx4DGjGyRDy
   DLYryciAdJOgl4kp0UQZDiyXbLtqslV0sl1BicZRXX2cHn20ann5xwqP5HTy9CIIE4059rCdl37M
   lFScBCkYurl/ewK9/XTS0WKCofOSrRhAieo4h8xShbeODIvt2OlDofFxp5rLyqqS7Dx6XMvK1usZ
   Z6geeOCRMb/VpkyW4/ddZWpOy9Zq33Y0FSdBCoYPcbZqHASIn046WkwwdF78agNehEdsmOpwu9F1
   612tYaCWlfXRUKha0zuY71R4T0Xe1unTb20Ld90elC91drbI5/I62fpdtdtq30hFkIJhW+AI4D7g
   A+B3wP5+OvNbTDB0XjoaAiLy98jTwOHkOJHpLsPtOpN/fBTUsF8gmWBobVWdNetVLSvboEOG3K8X
   X/x/MYfaHMHixTGc6YTdPubEbaY75W2rfSMRgQmGqMrQF7gN2OqnM7/FBEPnwE8IiExon9Tbt3CG
   J/lwu9EB5XxCAAAgAElEQVQB8OIFwPDhYxKaku6443E97DDVkSNVn3/e6S+ZUPNr+kq1q6ikpF/S
   +l52MZlQMBIRpMYgQA3wR2AVcDcw1U9nfosJhsInF5NXorAUsTuD2ifkeoX4lX14p1HY+bzzzmN0
   yJBLtLR0vfbv/2c96KCftgmYduf1NLetgTp8+MgMx5FeMDh1a6O0m5KSvmnbMjOSkY4gBcP7wP3A
   kUBPP510tJhgKHxyteUxVsjEfo/3LyTOzayqevvtT2hJSYM621APjJroHa3iVnUiqPaOaKO31tXV
   RcVhaq8brbl4nbSjhZmTJCgUGp/23do2UyMdgQgGnGiqF/lpOJvFBEPhk49JKtkW1shrIr3izE9j
   xlTryJG/U5EmhWUKP40bu3OPJtQ6Sku3TyKInN/D2kNs8L5MniM24mrq2EomGIzEBKkxvOSn4WwW
   EwyFTz7MGqlMLOEVdayg6N59jIosVmhSOMa9nmjHU3LBEOnXiA/hXesrZLYX53uQ/hujOPErGMS5
   NzkichVOxLC/AhvD11X1lZQ3ZhER0XTjNPJPQ0NDRGrKGZ4CymV6T2T9pqZPWLZsOk4wOgAnMN6i
   RfdE1auuHktj4yusWTOBt98+gi1bVgFjgQacVCPTcILVXQGAyK+ArajOwLGkLsYJmAdwFk7gPCe9
   Zyh0CytXrnSD4kFJSS2trfMSjikI/Lxzo+sgIqiqZHyfB8HQCMRVUtUDMu3MLyYYipNwRNDwpFpe
   PjNlPubY+mVlZwObaWm5Oup+IKreNtv8nt13f4Tu3ftRWvornnsuRPvEfRYVFfczZMgOfP75Rlat
   +ojW1p8BewBnADvixJB8Gvive98V7r3OpF9bOyNCWK1n2bLjyZVgMIxU+BUMpekqqGqNrxEZRhpm
   zZobEX4ampudkNTJBIMTrvoYnLxR0NJyHKHQUiorne+1tY5QmThxalS7X399BBs33s2bbx7HY48d
   4goNp83y8tu54472+95991e0T+phDeAA4ADKyn4F/JmWlj3ce2dSW7sgaoxTp05g5cqZEe07dWxl
   b3Qm0goGEdkBmAMMUNXJIjIS+I6q/snDvZOBq3Gc2Dep6mVJ6n0bJ4XoT1X13kwewOicNDQ0sHz5
   6xnd09T0CU4U+PCK/SxgtySr8V4Rnx9h6NB/UFp6XFsuhvZJOrmGAjB69KgIwXMXQNS9EK2dPP30
   TM4//zTuuedGPvjg3wwZMoKlS5cyZ861EXWmpdSMDCPvpHNCAPU4J59fc793B173cF834B1gqHvP
   q8DuSeo9DjxIkvMRmPM5rwRxiCrd3v1EY+jVa7DGBsUbPnxM1Ni2bFGdPv0t17n8osKCqC2k2Yho
   Gv8c8buSIttxnNG1tnvIyDn4dD6n1RiASlX9q4ic687Qm0Vki4f79gHeUdX3AUTkLuBHwFsx9U4D
   /g5820ObRo6JtesnW+1maipxVv9rgRHAjcBWRo8emfC+2DG0O4th1aqPXPMPPPnkpey88z5UVY3g
   5puf4s47r3bHk2hln/g5RozYhQ8+uJQhQwYyd66/Vf0HH/w7ypTV2gpwXcbtGEbeSCc5gEagH+35
   GMYBSzzc92PgxojvxwDXxtQZADyBc7r6FuDwJG0FI06NtGQrCF5s/faAdE4OgbKyPhmFloZxbkKd
   yJV4s+622426dWvmz+E3J0Si+xLFPXIC/tXatlIjpxCgxlALPADsLCLPAlXupJ9W5nioczVwrqqq
   iIgrIIyCpQG4jpdfXkdDQ0Pbajo6h7E3J3JLy2+JdPJ+61u3ZLg6fx/YCuwfce0hBg9+hJKSnye5
   ZwUw1f08LG5M6Z4hVnN57LEjGTZsR/r02Y7Nm8+hR49unHvuaey9995RDm7HF3I8JSU3M3r0SN+a
   iGHkCi+7kl4WkWpgN/fSP1V1s4e21wCDIr4PAlbH1NkLuMuRCVQCh4jIZlVdGNvY7Nmz2z7X1NRQ
   U1PjYQhGR6mtncHTT0+juXkF4b3+GzY4ZplsOlArK/ulHMOSJcfS0hK+chYAqtcD44FngbcT7hIK
   U109lsWLL6f9PMLpVFefk9EYY4VHa+sK3n33JuAaAL788iwuueQKFi68i/vuW8BRR53Chg1VhM89
   tLbuQWXlQhMKRmA0NjbS2NjY8YbSqRTAT3ET9AAXAvcCYz3cVwq8i+N8LiOJ8zmivpmSCpR2x2/m
   0T+TtZepkzc6suozrnP5Sd12233dTGzjU7aRDVNSfBvx0VthXJK8EOZ0NnIPAZqSLlTVu0VkP5xk
   PVfgeNL2SSNwtojIqTj2h27An1T1LRE5yf39es/Sy8g7GzduTPpbpltAM60PUFnZH/gf4Fj3yl+B
   GWzcOJ+NG6G5eWYGT+NvTO3aU/jK2ynbjK2fSqMxjIIineQAXnX//h9wtPt5mR8p5LdgGkNeyXRr
   abapr6/X3XY7QeF1hQ8V7lbom9FqPFYjKCnpm1bLSNZOKFTt5k+YGjEOx4kO39C6urqo+pYvwcgX
   BBhE7yGc1J6rgD5AD2C5n878FhMMicl00vE7SaULCx3UeFVVFy5cpKWlVytsUviLivRzo5ZWZ2ym
   iZ7UO7ZDKPwsPXvuqDBKnWB74zUyrLdh5JsgBcO2OFs5dnW/7whM9NOZ32KCIZ5c2PX93Jts8vfT
   //PPq/bo8Z7CEwonaDhsdkeS1GTb7m9+BKOQCUwwOG2zF05EsdPw4HjOdjHBEE+mE1JHcy97We2n
   mqwzGe+mTaq1tap9+36lIsdHmGqiE+340UCyPZFb6GujkAlSY7gIZwP4/wKXAMtxHNImGFIQtG05
   CMHQ0UkuVR9ex/vUU6q77qp6xBGq1dU/i7unpKRfVtODZmMiNz+CUagEKRj+BfSI+F4O/MtPZ35L
   ZxMMuVhFBmFK6uhqOtH9FRXDPaW6/PJL1dNPV91xR9V779WYnMvh2Ei12qvX4Kzkjk6U9MYmd6PY
   CFIwPAH0jfjeF3jcT2d+S2cTDEHanWOzk2XT+dzRccdO/s4undq0QeyeeEJ1551VjzlGtakpWTvR
   uZdLSvpG7f7pCGYOMoqVrAsG4Fq33I8T7exWt6wB7vPTmd9igsEh6AksW07m+JV+4uf//HPVX/xC
   dcAA1QceaL+eaBdUt279sm5Wiu8v+4LcMPKJX8GQ6oDbyzjxjl52hYPixDJqxFscpC5LUAebMo1J
   lCleD56lirg6adIk9tprNIsXT8FJf5mYxYth+nQ48EB4/XXo0ye2xgpgJuD0sXXrmXFttLbumtXn
   NwzDJZ3kwPEp7AGMIsLXkMtCJ9MYVIPMYZD/lW2yHATh5502bZp7VmBclClJVfWzz1R//nPVqqpN
   GgpdkvD91NfXu/dH9lGrsF2EealKYYRWVAxPalLz+t/ATElGsUIApqTuwOVAE/CKW5qA3wLd/XTm
   t3RGwRAEhTKBxQuGWjcZza3uBJ7YF/DQQ6oDB6oecsiH2qPH8JTPkegA2/Dhe7gCY1RUH87n6ANr
   fpzzfgW5Oa6NQiUIwXA1cBPQK+Jab5ysKtf46cxvMcHQjt9JKJunpOPDS0Su7uO1iZqaY/W441SH
   DlV99NGObZ1t92HEBq87PKqtXGlXhSKsDSMRfgVDKh/DD4BvqmprhNnpcxE5GfgnzoE3I8eE7fiZ
   4DULm9f6sb6IpqZRLFuWrPdBPPPMPN599xF22eV2tmw5ztOYU/k7HB+Gp2YCJ2i/j2HkhWQSgxRn
   FVL9FkTBNIYOEdQp6TDRq+awKelOhXcV3tPS0kOjVtR1dXVRGdzKyqoy1n6it7N6MyVlur3XC4Xi
   9zGMRBCAxvCWiExT1ajtNCJyLLAyIDll5InInM1NTeszujd2dd/SMocnnzyQsrIH2Wmn61i16hQi
   V9T33HMLra0tOAF7v6SlZTOnnPJrfv/7Kz2ttCdNmsT555/GlVdeCsAPf3gYa9euAlZFaRaRY6qu
   Po05c671rDV5xUJrG0VJMokBDAReBJYAV7plCfASMNCPFPJbMI2hQ6Szg8f+XlbWJ2pFH1k/le/h
   k09U99vvYxX5l8JDCfwPzoraSfoTHcYb+qTM+5zJ8yQiV4cOzb9gFBIEcfIZ59zCQcDpOAH0DvLT
   SUeLCQZvpJqg6urqtKJieNv2zkjSbT+NFAqJJuTWVtU77lDt31916NB7FW6LaCt6m2l5eX83VHXq
   7Gep8DPJJ35GJ3y3lwxwhtEZCUQwFEoxwZCeVKvodCtsrxNtonrf+94J+qMfqY4cqfrCC4nrwAgt
   KXHyKDg5EcYrVORUMNTV1WnsFteSkl4R3ys9ayyG0VkwwVCEZGKi8BPZNNx+KFStZWV90ppmnAk9
   MinNHdq9+2d6/vmqX32lbZN++5mGcKjs+igtxEmWs41GH1jzPjH7NyXVqrOt9XCFEb4Fk2F0FvwK
   Bi85n408kOkW00xpavokqv2ysrMJhW6ksrJ/wlAYDQ0NrFixEpgHfAMnmd8OTJ/+d+rqTowZ7wrg
   TGAAsAAnNMZZLF/+Jq2t093+3mTQoCr+858LAWGXXXZj7twLPTufM80Z7bAHTspygO94qG8YXRQ/
   0iTXhS6gMcRqBx3bMprelJRpasz48SxT2K9tvM6hs1EKI12NYpSKhE01tQlNR9kKLOgn7EVZWR8t
   Le2nZkoyihlMY+i8JNIORozYJaM2Uq2iE/0W/uyV5uYq4OCIK8uBLTGaxwqcg/Hz3TpnMHz45axa
   9TGtrbtl1J8XMtGq4t/BXQDMmjWXDz5YzZAh3jUWwyh6/EiTXBeKXGNIpB0MHz4mryG2wyvxgw+e
   qqee+oZut93X2q3bpQoLolbY0ZpH4kQ9zrV6jdyemo3nscNlhpEafGoMJXmWS0YSVq36iPPPP40J
   ExYyYcJCT/6FhoYGJk6cysSJU2loaIj7Hkl4BZ2o/Tlz5nDIIT9h8eJ/8+ijZ/P732/k8stf5KGH
   vk0odDMVFZcSCu3GwoV3UVnZL+WYPv/8i3CPOP6G66iouDSr/pJCItU7N4xOgx9pkutCkWsMTpjp
   2J08tRmtfjM5pJauHeilMEuhWeElhUoNhcZ7GHt0ZFXn81R1MrBlX/MptAB2hTYew8C2q3ZunK2g
   41xzTH3GZpH4rGcD48wsoVB1lKM2keP2u989xXUs/zvOJKSa2NkbPfY6dTKuVbmCQt0xjdNevQa3
   bVnNpnAolFPHZtoyCg0TDJ2cjq422/fph9tIdLK4XSuJ1Sh69NhRp09/S7t3/6/Cla4vIVqopAqF
   7WXXU3t4jOJcTZtgMAoNEwxFQCar39i68VnP6qNMONAnYgWvCQTHJ9q37wq9+eYl7mG39nu7deur
   odD4hHkQwhNfovFER1ztq8521tR5oDszZkoyCg0TDF2IZBNQ/Cq9VmF7VwiMj/ltXIxW8JwefPDU
   tvbDgmD48DERp6LjtZB0ZytSnYQuNsGgWlimLcMwwVDgZHPCiA/vUNvWbnSegm+4E/E41/bf/ltp
   6XdVZJnCGoW/RZmFkh+0i9ZCItN2ph5rfOgJv/kR0r1Hm5gNox0TDAVMtk0MjrM30kzUvmuoPf7R
   +LiTvTBVRfrrDjv8Sbfb7ms944zX9eCDk5l/wr6CWE1jqmsWGqdQm/ZZEgmGcITXTN9JpuHDzZRj
   dHVMMBQw2XZKJnLshkLVafvs3fuHuvPO/9XJk1U//NDbOEOh6hT5ndM/S7LJOluhsyPvMeevYUTj
   VzDYAbdOSKJDZZWV/aIOV7333nsxNcawadMCLrywNw8/DIMGee8r8iDc6NGjMhprqoN0hmEUKH6k
   Sa4LnVxjyLaJI1k+48hrIj3VCVz3gMIGhUd01KgpUW14TcQTxLP4acdMSYaRGZgpqbDJtlM0fTTW
   PyvcoPCpOucS2k9Sp4vEmm6c2XoWP+2Y89kwvONXMIhzb2EjItoZxplPJk6cyuLFU4Bp7pXPgCeB
   ZuArysrOZuHC25g0aVKCuo6pZ9GiewIdY0NDQ0R00xlmUjKMgBERVFUyvc98DAVKpsHYamtn0KPH
   bOANYBNQixMG+27gOr71rW/mdSIOh8hevHgKixdP4bDDpvkKMmdB6hzsPRiB4kfNyKQAk4GVwNvA
   zAS/H40T3P814BlgzwR1sq1hFTR+bOVPPqk6YMCX2r//U9qnz7c1+pRz9O6cfNjis7FjyHwIDvYe
   DK9QiD4GoBvwDjAU6A68CuweU+c7wHbaLkSeT9BOIC8tCLJh485kEv3yS9XTTlPdaSfV++5rH4MX
   J3IubfHZEAy2HdXB3oPhFb+CIegMbvsA76jq+wAichfwI+CtcAVVfS6i/gvAwIDHFBhB52mO5Ykn
   4MQTYb/9YMUKqKhwrnvJiTxp0qScmpZqa2fw9NPTaG52vpeXz6S2dkHO+jcMIwP8SBOvBfgxcGPE
   92OAa1PUPwu4IcH1LMvRYMjWSi4+t0JVVLjqzz9XPflk1YEDVR98MIAHCYiOailmQnGw92B4hQLV
   GDxvJRKRA4ATgPGJfp89e3bb55qaGmpqajo4tMIlcsXf1LSeN97YzLJl0wF44ok/ITKaqqrXmT9f
   +P73D+o0u306qqV40YT80lneIQT7HozOTWNjI42NjR1vyI808VqAcUB9xPdZJHZA74nji9glSTvZ
   F6UBEMRKLl4L+ULh9AQH22oVxmlJSb+0ge2MaGwFbhQrFKjGsBTYVUSGAmuBI4AjIyuIyGDgXuAY
   VX0n4PEESiYrucgVanX1WJYsecW9J9FqdUDE54XAauAampvhyisvpbn5GOB24DJaW+Gii84ESNNm
   /Dhqa2cAdJqVc7aYN+8G1y/knOtobnaudYVnN4yE+JEmmRTgEOCfOBrBLPfaScBJ7uebgPXAMre8
   mKCNQKRpOoLauRMfHru3u+KPXq2uX6960EGrVeR9DecwcKKk1qkTbnuc9uy5oybKk+AlW1oiX0Z7
   7oWus3K2XT5GsUIhblfNVsmHYAjSvJA4R8HhUZPSRRe9otts06SDBj2oF154mRtKu1pLS7fV2Oxq
   0DNBe+PSTnTJciV0tQnSTElGseJXMARtSuq05M+8sA0vvfQLHn10O1Rf5KOPmrjiiivbtr2OHVvD
   smXHt41r61aAK3A2dDmUlJxJa+sJAY+zeDBnrmFEY4IhD8Tu6YfTgelAI3AA3bs3ovpvnN29DTQ3
   D+Ooo07hjjt+nzDkNnwTRzjMpqJiHb/+dS1z5lxLc/MeQPIzA7HjKCs7G9hMS8uCpPd1pt07mZDr
   cx2GUdD4UTNyXSgyU1K4/bD/4txzr9Ltt39Wt932I73yyuciTDz1GpmOM1F4bcesVB9n+vHqH4mt
   l+o+M7kYRucC8zFkn6DDRrS2qt5+u+r226vOnKna3NzerzMBxzuVI8Nsh0LVOXUWm5PWMDoXfgWD
   mZJSEKR54eOP4eST4d134cEH4dvfju73vvsWcNRRp7BhQ+pxRZt2zDZuGEbHsXwMOUYV/vxnOPts
   OOkkuOAC2GabxHVjYy+Vl8/Ma2rMQhuPYRip8ZuPwQRDDlm9GmbMgLVrYfr0Z/nHP+YBmR1Ay/ck
   XGjjMQwjOSYYChhVuOkmOO88OP10GDt2ET/5yXG28jYMI1D8CgbzMQTM++/D9Onw6afw+OOwxx4w
   ceL1FoLBMIyCxVJ7BkRrK/zhD7D33nDQQfD8845QKGYs3aRhFAdmSgqAd9+Fn/8cvvoKbr4Zdt89
   +vdYJ25Z2dl861vfpLKyf6e125tj2jAKD7+mJNMYskhrK1xzDey7L/zgB/D00/FCAdq3o06YsJBQ
   6BbAybewePEUDjtsWl5X235X/dEhRBwBEXZSG4bRuTDBkCX+9S/43vfg73+HZ5+F2lro1i26TuSk
   C7Bo0T1UVvajpeVqIifUo446JS/mmPCqf/HiKQUhpAzDyBN+TsXlulDAiXq2bFH97W9V+/VTnT9f
   devWxPWShZNIHuE09yEnOnKy2cJlGEbhgZ18zj1vvgnHHw/bbgsvvgg775y8brupZQfgBpqbhzFr
   1qXMnXthTEC9s3CS7kzqVLuVLEKpYRQPJhh8sGULXH45XHUVXHqpc2itxJNR7gGcCKq7AuNZvvxm
   gLYJ9eWXl7NhwzQgPxNqbLTVZFFZk2ERSg2jOLBdSRny2muOlgBNfOMbF1Bevs7TTqI5c+ZwwQWX
   A/PdKzOBY5gwYRWLFt0DFMbOHjvZbBjFg518DpiWFpg7F373OzjmmNe57rqD+eorbxN4Q0ODGxDv
   QsKH2mABcB2hUDcqK/sDXTfnsmEYwWAnnwPklVccLWHgQFi2DE444WJXKKQ/udyuBQyL+03kn7zx
   RiktLdMBePrpadx334I2DcIwDCMfmGBIwddfOz6EG26AK66AY48FyVD2Rjudp7VdLyk5k2HDBvLu
   u7XkOzSGmY8Mw4jEzjEk4cUXYa+94PXXYflyOO64dqFQWzuD8vKZOOagBa6TdkaaFicRNh9VVFzK
   ww/fyc477xpXq6lpfU7DStjZBcMw4vCzxzXXhRyeY9i0SfXss52sanfc4WRZS0QmqTOT7e+P/a2s
   rCqnGdlULSubYRQz2DmGjvPss3DCCU6wu9deg/79k9f1ujUz1f7+2N+amr7JsmXTybdpyTCMro0J
   BmDTJieT2p13wrXXwo9/nN32UwmRyN/CoTJySUfPLhiGUXx0ecHw5JNw4olOzuUVK6CyMn9jycck
   bSeWDcOIpcueY/jyS5g1C+6918mb8KMfZbV539gOIcMwsoUdcMuAxx5z8iVUVzthLfr2zVrThmEY
   BYMdcPPA55/D2WfDww/D9dfDoYfme0SGYRiFR5c5x9DQ4Ow2am11ziaYUDAMw0hM0WsMn30Gv/41
   PP443HQTTJiQ7xEZhmEUNkWtMTz4IIwaBT16ODuOTCgYhmGkpyg1hg0b4IwznANrt90GBxyQ7xEZ
   hmF0HopOY7jvPkdL6NfPOb1sQsEwDCMzikZjWLcOTj3VCYt9992w3375HpFhGEbnpNNrDKrw1786
   O44GD3YioZpQMAzD8E+gGoOITAauBroBN6nqZQnqzAcOATYBP1PVZV7b//e/4Ze/hJUr4R//gH33
   zdbIDcMwui6BaQwi0g34HTAZGAkcKSK7x9Q5FNhFVXcFZgB/9NK2Ktx+O4weDSNGOBnWuopQaGxs
   zPcQCgZ7F+3Yu2jH3kXHCdKUtA/wjqq+r6qbgbuA2IhEU3Cy16CqLwB9RCRFsGtYswamTIHLL3dO
   MP/mN8521K6C/aNvx95FO/Yu2rF30XGCFAwDgI8ivq92r6WrMzBRY6pwyy0QCsHYsbB0qZNhzTAM
   w8guQfoYvEa9iw3wlPC+Qw6BTz6BRYtgzJiODcwwDMNITmDRVUVkHDBbVSe732cBrZEOaBG5DmhU
   1bvc7yuBalX9JKatwg8BaxiGUYAUWnTVpcCuIjIUWAscARwZU2chcCpwlytIPosVCuDvwQzDMAx/
   BCYYVHWLiJwKNOBsV/2Tqr4lIie5v1+vqg+LyKEi8g6wETg+qPEYhmEY3ugUiXoMwzCM3FFQJ59F
   ZLKIrBSRt0VkZpI6893fl4tIKNdjzBXp3oWIHO2+g9dE5BkR2TMf48wFXv5duPW+LSJbROTwXI4v
   V3j8/6NGRJaJyOsi0pjjIeYMD/9/VIpIvYi86r6Ln+VhmDlBRG4WkU9EZEWKOpnNm6paEAXH3PQO
   MBToDrwK7B5T51DgYffzvsDz+R53Ht/Fd4Dt3M+Tu/K7iKj3OPAgMDXf487Tv4k+wBvAQPd7Zb7H
   ncd3MRuYG34PwHqgNN9jD+h97A+EgBVJfs943iwkjSGQA3GdlLTvQlWfU9X/ul9fIMn5jyLAy78L
   gNOAvwPrcjm4HOLlPRwF3KOqqwFUtSnHY8wVXt7Fx0Bv93NvYL2qbsnhGHOGqj4FfJqiSsbzZiEJ
   hqweiOvkeHkXkZwIPBzoiPJH2nchIgNwJoZwSJVidJx5+TexK1AhIk+IyFIROTZno8stXt7FjcC3
   RGQtsBw4I0djK0QynjcLKex2Vg/EdXI8P5OIHACcAIwPbjh5xcu7uBo4V1VVRIT4fyPFgJf30B0Y
   CxwEfAN4TkSeV9W3Ax1Z7vHyLs4DXlXVGhEZDiwWkdGq+kXAYytUMpo3C0kwrAEGRXwfhCPZUtUZ
   6F4rNry8C1yH843AZFVNpUp2Zry8i71wzsKAY08+REQ2q+rC3AwxJ3h5Dx8BTaraDDSLyJPAaKDY
   BIOXd/FdYA6Aqr4rIquA3XDOV3U1Mp43C8mU1HYgTkTKcA7Exf6PvRA4DtpOVic8EFcEpH0XIjIY
   uBc4RlXfycMYc0Xad6GqO6vqMFUdhuNn+EWRCQXw9v/HP4D9RKSbiHwDx9H4Zo7HmQu8vIuVwMEA
   rj19N+C9nI6ycMh43iwYjUHtQFwbXt4FcBHQF/iju1LerKr75GvMQeHxXRQ9Hv//WCki9cBrQCtw
   o6oWnWDw+G/iN8AtIrIcZwF8jqpuyNugA0RE7gSqgUoR+Qi4GMes6HvetANuhmEYRhSFZEoyDMMw
   CgATDIZhGEYUJhgMwzCMKEwwGIZhGFGYYDAMwzCiMMFgGIZhRGGCwSgKRGSrG256mYi8IiJDROQZ
   97chInJkRN3RInKIjz4aRWSvLIw1K+0YRlAUzAE3w+ggm1Q1Ns58OH7UMJzIo3e630M4YTQeybAP
   JTuxubLVjmEEgmkMRtEiIl+6H/8P2N/VJs4B/hc4wv3+ExHZ1k128oKrbUxx7y8XkbtE5E0RuRco
   JyYYmZsw5u6I7zUi8oD7+Y8i8pKbKGZ2mjEiIj8WkVvcz1Ui8ncRedEt33WvV8doRj2z9LoMow3T
   GIxioVxElrmf31PVqbSvymcCZ6nqDwFE5BNgL1U93f3+G+AxVT1BRPoAL4jIo8DJwJeqOlJE9gBe
   IX6l/yhwvYiUu8HrjqBdMzlPVT8VkW7AoyKyh6rGZtnSJJ+vAa5S1WfcuFj1wEigFvilqj7nxkP6
   OkHuKOUAAAGySURBVNMXZRjpMMFgFAvNCUxJYWJDDseG5p4I/FBEznK/bwMMxsmMdQ2Aqq4Qkddi
   G3bj9tQDU0TkHpxsWeF2jhCR6Tj/n+0I7A4kTb8Yw8HA7m4cLIBeIrIt8AxwlYj8BbhXVYsxurCR
   Z0wwGF2RRPb9w2PzFriTspfcDncBpwIbgKWqulFEhuGs7vdW1f+6JqIeacZSHtk9sK+qtsTUv0xE
   HgS+DzwjIpNU9Z8exmgYnjEfg9EV+ALoleJ7A3B6+EtEsvQncZzWiMgoYM8k7S/BSZAznXYzUm+c
   SJafu2Gfk+2C+kRERohICXAY7YJiUcyYxrh/h6vqG6p6OfASTjhpw8gqJhiMYiGRFhC+thzYKiKv
   isgZwBPAyLDzGbgU6C4ir4nI6zjOaXBShfYUkTfdawmTvKhqK/AgMNn9i6ouB5bh5AX4C/B0knGf
   697zDLA24vrpwN4islxE3gBmuNfPEJEVbjjpFjLfWWUYabGw24ZhGEYUpjEYhmEYUZhgMAzDMKIw
   wWAYhmFEYYLBMAzDiMIEg2EYhhGFCQbDMAwjChMMhmEYRhQmGAzDMIwo/j/G71js/VPzEwAAAABJ
   RU5ErkJggg==
   "
   >
   </div>
   
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p>Plot yhat vs. Pearson residuals:</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[13]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">fig</span><span class="p">,</span> <span class="n">ax</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">subplots</span><span class="p">()</span>
   
   <span class="n">ax</span><span class="o">.</span><span class="n">scatter</span><span class="p">(</span><span class="n">yhat</span><span class="p">,</span> <span class="n">res</span><span class="o">.</span><span class="n">resid_pearson</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">hlines</span><span class="p">(</span><span class="mi">0</span><span class="p">,</span> <span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">set_xlim</span><span class="p">(</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">set_title</span><span class="p">(</span><span class="s">&#39;Residual Dependence Plot&#39;</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">set_ylabel</span><span class="p">(</span><span class="s">&#39;Pearson Residuals&#39;</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">set_xlabel</span><span class="p">(</span><span class="s">&#39;Fitted values&#39;</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt output_prompt">
       Out[13]:</div>
   
   
   <div class="output_text output_subarea output_pyout">
   <pre>
   &lt;matplotlib.text.Text at 0x7fd0d953b850&gt;
   </pre>
   </div>
   
   </div>
   
   <div class="output_area"><div class="prompt"></div>
   
   
   <div class="output_png output_subarea ">
   <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAYoAAAEZCAYAAACJjGL9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz
   AAALEgAACxIB0t1+/AAAIABJREFUeJztnXucHWV9/9/fJUQWSQhLEBARcEW5iLjiBRvbrNYkWDX+
   Qop3Gy8Vq61aXduI/JRU1grWeKn9/WyhaCLeSgU0WN0lXhaFVryFgCIqKiiK1BgU0fzKJd/fHzNn
   95yzM3Nm5sycmXP283695rU7Zy7Pd54z5/k8z/f7XMzdEUIIIeIYqtoAIYQQ9UZCIYQQIhEJhRBC
   iEQkFEIIIRKRUAghhEhEQiGEECIRCYXoKWb2QjObTjg+Y2YvLyCdcTP7abf3qQNmtsnMLq7ajjQM
   Ur6LOSQUIhYzu8XMfm9mvzWzX5jZxWa2tJt7uvtH3X1N0inhVipmttfM7g6fbZeZfd7MnlN2ujmp
   1WAnM3uJmd0f5t1vzGyHmT0jx322mNm5ZdgoikVCIZJw4JnuvgQ4GTgJ+N/VmlQojw6f7RHAFuCf
   zOyt1ZrUN1wT5t0y4CLgEjNbVrFNoiQkFCIV7n4HcCVwYuMzMzvVzP7TzO40s+vMbGXTsZeY2Q/N
   7C4z+5GZvaDp8680nbfKzG4ys1+b2fsBazrW4nIxs6PDlsBQuP9SM7sxTOOHZnZmzmfb7e4fAV4F
   nGVmI+H9DzSzi8zs52Z2m5md25T2S8zsGjN7f2j7d83sqU22drr2ajP7BzPbHebPaU3XHmNmV4XP
   dSWwvNneDvk+Y2ZvC+9/l5lNm9nBTcef3HTtT8xsQ/j5A8zsXWZ2a9h6/ICZ7ZeQbRbmnQMfAoaB
   h807yez40KY7zezbZvas8PMzgRcAfxu2TD6d6ssSlSChEJ0wADN7CHAacG24fwTwGeBt7n4Q8Ebg
   UjM72MweCLwPOM3dlwJPAq6bd2Oz5cClwJuBg4EfAiuaTunkcrkDeEaYxkuB95jZWN4HBbYBi4DH
   h/tbgHuAUWAMWA38edP5TwBuDm0/B7isqVad5tqbwmvfSVArb/Ax4OvhsXOBDYR5kZTvTdc/H3gJ
   8CBgcXgOZnYU8FmC72Y58BjmvpfzgIcTtBwfDhwBdGxdmdmi8Ll+C/yg7di+wBXAFHAI8Brgo2b2
   CHe/APgocL67L3H3Z3dKS1SHhEIkYcCnzOwu4CcEBflkeOxFwGfdfQrA3T8PfAN4BkGhthc4ycyG
   3f0Od78x4v5/Anzb3S9z9/vd/b3AL9rSj8XdP+vuPw7//zJBi+cPcz4r7n4vsAsYMbNDgacDr3f3
   Pe7+S+C9wPOaLvlvd39faPslwPeAZ6a89lZ3vyiskX8YONzMHmRmDwUeB7zF3e91968QFLYNkvId
   grz/kLvf7O7/D7iEQBAgqMFvd/d/C23e7e47zcyAVwBvcPdfu/vdwDva7G3nVDO7E7gdeC6wzt1/
   234O8EB3P8/d73P3LxGI3PPD40aH71jUg0VVGyBqjQPPdvcvmtkfERRYjwO+BhwFnNFwJYQsAr7o
   7r83s+cS1GQvMrNrgAl3/17b/R8M3Nb2WeoeM2b2dIKa/LEElZ79getTP938++1LUPPdTfB8+wK3
   B+UohGn8pOmSn7Xd4laCZ3poimtnBTHML4ADCFoBd7r7nrb7Hhn+H5vvUfcG9oT3JbzHj5jPIQR5
   980me43kiuRX3b2TKD+Y+d9nI4+gZkF6EY+EQqTC3b8cxhDOB55CUOhd7O6RcQF3vxK40sweALwd
   uBD4o7bTfg7MuhzCmu2RTcfvJijAGhzWdO4DCNxWLwI+7e73m9nldFdDfTZwH4EQ7gf8D3Cwu++N
   Of+Itv2jgE8TFI6dro3jduAgM9vf3X/fdN/7w/8T870DPyFwebWzi0BQTnD323PcN46fA0eamfnc
   NNVHEbjcQELRN8j1JLLwXuAJZvZE4CPAs8xstZntY2b7WdCH/ojQhfLsMFZxL/A75gq6Zj4LnGhm
   60Jf92tpEgMC//kfmdmRZnYgcFbTscXhtgvYG7YuVmd8nkb8ZcTMXgj8E3Ceu98ZFphXAu82syVm
   NmRmo2HLqsGDzOy1ZravmZ0BHEfgFvpFimsjcfdbCVxJfxfe98nAM5tOic339ueK4GPA08zsDDNb
   FMaTTg7F7ELgvWZ2SJgnR5hZ1vxs51rg9wQB633NbDx8lk+Ex+8gIgAu6oeEQqTG3XcBW4GN7n4b
   QQ38zcB/E9RWJ5hzWbyewDXzK4K4wasatwm3xv3OIAik7iIIol7dlN7ngX8jcCd9ncD11bj2twTC
   cgmBq+j5BLX5FpM7PNJOM2sEYV8G/LW7b2o6/mcEYnRjmMa/0ypk1xK4vX5JEHRe7+53prg2aqxI
   8/4LgCeG172VIM8Jnzsp36Pu1ZzfPyGIC00QfC87gEeH520kCMx/1cx+A2wn6DYcRaexLo307gGe
   RRCv+SWBEL/Y3b8fnncRcELYI+qyhPuJirGqFi4ysyMJgngPInixLnD3f7Sga+K/ETRRbwGe4+6/
   rsRIIWIws5cAL0/hpxei76myRXEvQa+QEwl6R/ylmR0PvImgZ8YjgC+E+0IIISqiMqFw91+4+3Xh
   /3cD3yUIDq5lrqm9Ffhf1VgoRCI9mWpEiDpQmeupxQizo4GrgEcBPwkHEjV6wexu7AshhOg9lQez
   zewAgm6Or2sfsBN2qateyYQQYgFT6TiKcIDTpQT9wj8VfnyHmR3m7r8ws8MJena0XyfxEEKIHLh7
   5rFGlbUoQrfSRcCN4dQNDbYRzG1D+PdT7dcCuLs2d84555zKbajLprxQXigvkre8VNmiWEEwqvZ6
   M9sRfnYWQZ/6SyxYvOYWoK5rBAghxIKgMqFw96uJb9E8rZe2CCGEiKfyYLbojvHx8apNqA3KizmU
   F3MoL7qnFt1js9I6x5gQQog0mBneT8FsIYQQ/YGEQgghRCISCiGEEIlIKIQQQiQioRBCCJGIhEII
   IUQiEgohhBCJSCiEEEIkIqEQQgiRiIRCCCFEIhIKIYQQiUgohBBCJCKhEEIIkYiEQgghRCISCiGE
   EIlIKIQQQiQioRBCCJGIhEIIIUQiEgpRGtPT06xevZ7Vq9czPT1dtTlCiJxozWxRCtPT06xbt4E9
   e84HYHh4I5dfvpU1a9ZUbJkQC5e8a2ZLKEQprF69nu3b1wIbwk+2smrVNq688tIqzRJiQZNXKCp1
   PZnZB83sDjO7oemzTWZ2m5ntCLfTqrRRCCEWOosqTv9DwPuBDzd95sC73f3d1ZgkimBi4kyuvnoD
   e/YE+8PDG5mY2FqtUUKIXFQqFO7+FTM7OuJQ5qaRqBdr1qzh8su3snnzBQBMTCg+IUS/UnmMIhSK
   K9z9pHD/HOClwG+AbwAT7v7rtmsUoxBCiIzkjVFU7XqK4gPA28L/zwU2Ay9vP2nTpk2z/4+PjzM+
   Pt4D04QQon+YmZlhZmam6/vUrkWR5phaFEIIkZ2+7PUUhZkd3rS7Drgh7lwhhBDlU6nrycw+DqwE
   lpvZT4FzgHEzewxB76cfA6+s0EQhhFjwVO56yoNcT0IIkZ2BcT0JIYSoFxIKIYQQiUgohBBCJCKh
   EEIIkYiEQgghRCISCiGEEIlIKIQQQiQioRBCCJGIhEIIIUQiEgohhBCJSCiEEEIkIqEQQgiRiIRC
   CCFEIhIKIURtmJ6eZvXq9axevZ7p6emqzREhmmZcCFELpqenWbduA3v2nA/A8PBGLr98K2vWrKnY
   ssEh7zTjEgohRC1YvXo927evBTaEn2xl1aptXHnlpVWaNVBoPYoFgJrlQogqqHQpVJGe9mb51Vdv
   ULNcDBQTE2dy9dUb2LMn2B8e3sjExNZqjRKAXE99g5rlYiEwPT3N5s0XAIFwqCJULHldT2pRCCFq
   w5o1ayQONURC0SeoWS6EqAq5nvoINcuFEN2g7rFCCCES6cvusWb2QTO7w8xuaPpsxMy2m9n3zexK
   M1tWpY1iYaCux0LEU/U4ig8Bp7V99iZgu7s/AvhCuC9EaTS6Hm/fvpbt29eybt0GiYUQTVTuejKz
   o4Er3P2kcP8mYKW732FmhwEz7n5c2zVyPYnCUNdjsVDoS9dTDIe6+x3h/3cAh1ZpTJHIvSGE6Edq
   3T3W3d3MIpsOmzZtmv1/fHyc8fHxHlmVD42sri/qeiwGlZmZGWZmZrq+T11dT+Pu/gszOxz40iC4
   nuTeqDfqeiwWAoM0MnsbQWl6fvj3U9WaIxYCGhHcHRLawaZSoTCzjwMrgeVm9lPgrcB5wCVm9nLg
   FuA51VlYHHJviEFFbtXBp3LXUx760fUEqnWJwURu1f6hNNeTmf01wXiHu4B/BR4LvMnd1W0nI3Jv
   CCH6kTTdY1/m7r8BVgMjwIsJ3ENCCMHExJkMD28EtgJbQ7fqmVWbJQokjVA0minPAC5292+XaI8Q
   s2jcSX+wZs0aLr88cDetWrVN8YkBpGOMwsy2AA8GHgacDOxD0GX1lNKti7epL2MUIj3tAdLh4Y0q
   gIToktJmjzWzIWAM+KG7/9rMDgaOcPfr85naPRKKwUcBUiGKp/BgtpmdAjRKYwceZmYQuKJUSgsh
   xAIhqdfTZpIF4SkF2yLELBp3IkR90DgKUVs07kSIYil1hTszOwk4Htiv8Zm7fzhrYkUhoRBCiOyU
   Ns24mW0C/hH4JwJ30zuBtVkTEkIMLurKPNikGUfxp8DTgNvd/aUEXWS1PKkQC5AoQdAKgYNPmkkB
   97j7/WZ2n5kdCPw3cGTJdgkhakbc5H+bN18QfhZ0Zd6zBzZvvkAxpQEijVB83cwOAi4EvgH8DvjP
   Uq0SA4GC0YNFnCCIwaejULj7q8N//9nMpoElVQ62E/2Bpp5eOKgr8+CTZmT2SiLGU7j7l8syqhPq
   9VR/NLJ68EiaVkWtx/6gzBXu/oY5odgPeALwTeCpWRNbiOgHJPqBNO9pY/K/ufPmWoiaQn/AcfdM
   G0Eg+7Ks1xW5BWbXn6mpKR8ePtRhi8MWHx4+1Kempqo2qycshGefmpryVatO91WrTu/rZ1sI35UI
   CMvO7GVu5guCuZ6+myexorZ+EYpVq04Pf3weblt81arTqzarZwxKQRrFIBWuC/09XUjkFYo0K9y9
   v2l3CHgMgetJiEQG2R2hLqFiIZEmRtEsCvcBH3f3q0uyZ6BQbxDRD+g9FZ3QpIAlU9dgdl3t6hfq
   vLBSnu9W78PCoPBJAc3shoTr3N0fnTWxougnoagjdS7k+ok6Fq76bkUSZQjF0eG/jQF3FxMEsl8I
   4O4bM1tZEBKK7lioYxyqKNh7neZC/W5FOgofR+Hut4Q3Xu3uj2k6dL2Z7QBKFQozuwW4C7gfuNfd
   n1BmeqJ76ljDblDFSPHqRqffAKwP/z+m5LTEgqBTtyhgJ/Dkpv0VwHV5ulhl2YAfAyMxx/L3D+sT
   yuxaWkbXzrp3F62iC2i3aeZ5ByYnJx2Wzn4PsNQnJyfzPsJAMMjdtLNCWeMogFOA64Fbw20n8Ng8
   iWUyLBCKg2OOFZx99aIXhW7RP56698XvN6GYnJz0oaGDMr8Ddf8eek3dKzC9pjSh8LnC+UDgwDyJ
   5DIMfgTsIJix9hVtxwrNvLoR9WMfGRlN9YJXVXuqQwGV9OxVFBh505yamvKhoYNz5Wcdvoc6ofxo
   Ja9QxMYozOzF7n6xmU3QNCmgmVmY2Lu7dXt1YIW7325mhwDbzewmd/9K4+CmTZtmTxwfH2d8fLxk
   c6pl9+5DWLcu2cdd5YytVffF7/TsjXmKzjrrXG699RccddRxpduUNDdSEps3X8DevcfmSrPq70HU
   i5mZGWZmZrq/UZyCAK8M/24CzmnaNgHn5FGlvFuY7kTTfnESW0Paa6JwqMNUx9pQ1bWnKn3BaZ69
   X9wQwbNMhN97YOvQ0EGpbU3zPSwUv32/fOe9gl7N9dSLDdifYN0LgAcC1wCrm44XmXe1ZGpqykdG
   Rh1ODUWic8FftVBUSZpnLzt/2gvfvIXxXOE24XCqDw0dnDkgXTc3XBaKFrGFIoppKE0ogHcCS4F9
   gS8Au4AX50kstVFBn77rwu3bwFltx4vPwRqS9Qdd9wKgTNI8e5lC0Z7+4sXLfPHiQ3J/F90Ubp3y
   Iks+9LqQXcjvcC8oUyh2hn/XARcRBLWvz5NYUdtCEQr37D/Ufq89TU1N+djYCh8ZGfWxsZWFFq5l
   FkLzC99TK2vddRKCTscb+Tg2tqIrsSvDdtEdeYUizaSAjXOeCXzS3X9jZp50gSiOtDOw1nmwW1qm
   p6dZu/Z53HPPIuBd7N4Na9e+mG3bLk71PJ3yKm9wedBICni3dgr4Z+Af0Ay5Ik3t/TzgJgI30GLg
   QcC1eVSpqI0F1KJIQ9nN9V61UoLaZHU18TTE5UXRrqdubez0PsQ9R2uNvve1e7meyoUyg9nACLBP
   +P8DgcPyJFbUJqFopZe+9zJ/uFUKRdqeQkl5UVQwuxsbo86dnJxMfV3ruzTlsLznhXa/u0/rTGlC
   EQrDW4ALw/1jgWfmSayoTULRSplC0Uuf8dTUlC9evKylcFq8+JCeD4xbvHiZj42t7FDbLi4vihCp
   tM+WtUNEXF4sRAZBwMoUiksIJgD8js8Jx848iRW1SShaKWvuplWrTg+76E70rHDsJpidl/m16KVh
   y+ZUX7x42awNZQhF2u8ub9p5rutVgVhkOmXbPCgusTKF4pvh3x1Nn0koCqCuP5T5A/6WhmLR2d2S
   5/5V/+haC9MVLS0aWO5jYysKs7s9v9IW5L0Uil5Q5DuQp9WU9Z2taz5mpUyh+E9guCEUwCjwtTyJ
   FbVVKRRFFci9KCzz2hr1oxgZGe0YwO3Xietan+Mhkc/efG6aPI06Lyq/xsZWpsqLXrmeekWR70DW
   cSGD8M7mpUyhWA1cBfwS+BjBDLJPyZNYUVtVQlHkj67IFy/KXRNl6+TkpI+NrQzPWxEbcK2qltsL
   90Hc/RvHliw5cp5tY2MrM6cT9Z5EPffY2IrU71Te/Kmjb70qocibbl0FNytl93paTjCO4pnh/w/J
   k1hRW1VCUdXLnURcAHhsbEXb/SfcbNnsObDcFy16YGQXzrL95nEiVnYX3zT3D/LzkJa8zGpHXL50
   K5BFFPh1EY2qXE/d/O7qknfdUIpQEKxFcQZwYrh/JHAB8JM8iRW1DYJQFPVDCWya36U0CEJv8SA4
   e7pHuVSiPms8T5k9caLuH+fuKuoHGQjnqWFeTCW2FDo9e6fjSYLQTX5164cvs9NDnsKzihjdoLQM
   8lK4UACTwHeBjwM/BDYTLCb018B+eRIrahsE11Pjfs0vd54CKk4oxsZWtrU05p+TJBRpbR8bW1FI
   98mowhVOLawwa14EKJiVdcKHhg7OVbilGcwWd06RcaOsfvii/exFVhR6ySC0DPJShlDc2BCEcMDd
   74Cj8yRS9DYIweyo+3YazBV1PM71FMQtVjYVDK2Dp2C5Dw09wM1GQhGZqDRA2n7PtFOrpyFahA52
   mMh8706F7ZyArvSxsRWFuJOmpqaaWoidC/msrq+8LHR/fz9ShlDsaNsvfZ3s1EYPSPfYZjr96JKO
   RwWzo6+Z8CVLHuojI6M+OnpCiy9+aOig1FNZl9UDZK5ATD+1el57YWWue3f6Hoqe7XfunAlPO0q6
   DNdX1rwo4zrRPWUIxW+AK5q2Xzf9vy1PYkVtEop0P66kdZe7+bFG+fuT3CBZWmC9aa0s79iCirO7
   /V7NI5fTdnVtkOY7aD1nyuHUjrGbMlxfWdNJot+Fop9dV2UIxXjCtjJPYkVtgygUeV1Pne8XvfhN
   N26D5pYILG8ZvZzlmZLSKDpA2uwSGh09IXHkd5rvohGfaW2VZVvnOrtQBOek6a7bq8IsT0VgbGxl
   mFfRgzjrTL+7zUrtHlu3bRCFwr373jbNpPGlF1UbjCu4el1z7DbQ3CBtyyDKtdfcBblT19q09qYV
   5rrT/rxDQwfNjufpF/q9NSShEC2keaGzCE/j3CxzP0XZEDXCu1uy2JZGQOe3DCYi7Z5/r/UOB3jQ
   m+xRvmjRganzNSlPsrj66ky/F7Lu/f8MEgrRQpFN5Pk+/vi5n4q4rjvblntSILxToRsUBBMe9Lra
   Ev6/NDIfW9NuPa/RBbe9tZXHJdTvhVODQXgOuZ76aJNQpKMoP3U3LYM8LZG0xPWSaowZiYotRLlx
   mtdrmBvRHj9QMaq32ZIlD513XiPw3Jz+IM3XVHZHhboGjRtxluapcPqFMud6eiRwIbAd+FK4fTFP
   YkVtC0Eo6vQjKaImWHRtsr1HV/O4izgRi7JhdPQxbb2YDgnHpcQHqEdGRn1ycjI8L5iO3OygeefB
   iI+OnjRbqETNI5VVcHvxPqQd21FmR4W6iqN7vW3rRJlCcT3wKuCJwOPC7ZQ8iRW1DbpQVPkiRv2Q
   i7CnaFdYVAHeaSR3XMsoKjjfyIP2OagaXWuDoHWzq2mpmx3YtL/Mh4b290WLmj8bibS5kR9ZVqIr
   i7TfUxlupOZ3b/5cZfVxU/WzC61MofhmnhuXuQ26UFT1Ihbd/779mjgRyuezj5qSZCR20GDDTdQ+
   riRNgZTs4pq77oADDm9xSURNzAitEzPO3W8idsxLN2TN37TvXtHvaFSPqDIWzCoCCUV0obwJ+Evg
   8HAqjxFgJE9iRW0SiuJorcWtLCzdorqqRjE/4LzF4UCH42Z9xu0C1Rx0Hho6OGI69ujxJq1pRrUG
   5vaXLDmypUCOvuY4HxkZjYjZFL9WeJ78TRuPKrrVG5XuXKuxXu4duZ6iC+VbwskAm7cf5UkstVFw
   GnAT8ANgY8TxwjOwTvTiRYyqYWcdMJZEGrHLK4jthXtQS1/vsMXNlrXEGOJaDY15mFatOt03bNjQ
   sTbf/p0sXnxIm1tpfo+upDEQ82vQ8XmfNz6RJ3+z9FRr2JV2Xqustja7AOtWENcphpiFgen1BOwD
   3AwcDewLXAcc33ZOoZlXR4p6EaNq13MjY49r+3Hmc39E2Tr/hz9/LEI3LafGcwR+/1YXRXtNf34c
   YiK1QDY/W3sMoSG2ixY9yOFRHtUtN6mHTPu94yZ9zFtp6EaI0/ZUK6pS08+19H6izBbFYuB1wKXA
   J4HXAPvmSSyVQfAkYKpp/03Am9rOKTwDB5H5teBlbTXcdj/whB9wwOGJ01u0s2HDhrCwfojD+pgC
   LnosQhGFw+joY+YVhkGhPbc/NrayQ+092u3Tyb75te9DHSY9aT6mpArA5OTkrFuq4f7qVkybe2Zl
   HdFdZqswzt5+rKX3E2UKxUXAVuCpwB8DW4B/zZNYKoPgT4ELm/ZfBLy/7ZzCM3AQmf8jjgoANwrN
   +IFlcUxOTnp7zx9YP682njRFdreFQ+BWap0+3WxJpCi1jnloFciollSnQjCq1dQcrM4yX1dcHGV0
   9ATPOyq73fWVdbW+NEKepuXYLb0SkIUgVHmFwoJr4zGz69390Z0+KwozWw+c5u6vCPdfBDzR3V/T
   dE6y0UIIISJxd8t6zVCKc+4zs4c3dsxsFLgva0IZ+BnBkqsNjgRuaz8pjyqWvU1NTbFq1emsWnU6
   U1NTpaYzPHwoQeNuC7AUmAC2MDx86Gza7ectXryMRYsOnt03W8bY2IpZu4PPPdy2sGrV6YnPNTIy
   Ou8aGJn37O12NNtYXB4sByZi7x31fCMjo5HfVdp7zz9vWWwextnQOD53bP45wWfB/6Ojj0mdR0np
   Ffn+Tk5OMjIyyqJFDwrfw3TplfEMdU+n6i03KW78x8BPgKvC7VbgqSU+yCKCpVePJoiP9EUwu5fB
   uOhul6dHukfaA6ZxPus4+6P6tzeCslGupw0bNsTmT5HN+qg8GB5eHrssaxZfelwX0UZ+tAe059xr
   6725y+7Q0EE5XE9R7sHTW+xIm79pYiydgumdmB+nSZ5rKytRXbbTTLOelX4eG5EFyuz1BOwHnAw8
   GnhAnoQyGQVPB75H0PvprIjjBWdf9/TyRcsiFM10GieRrvdS0KuoUZBEBWB7QXR8YE602qewziLk
   cd9l0j1WzY7tWOkw6o0xHe0kCWYjjtI6Nclc91RY7qOjJ0VcE7++Q1x6URWAPAPc4t6PoipLUTGo
   qHztloXS66o0oQCeAywN/38LcBnw2DyJFbUtdKGYX4tbGtZm4weMTU1FT3vRycZ4UZorPNO0FOJq
   4u014LStjuB5mgvUzlN6ZLE1qtBI+o6jWld5hbORViAA+3ujBQhLE1uBWdYZjy/guxeKIoPZcwJ8
   erhlX+c8LQpmdycUN4R/nwzMAM8EvpYnsaK2OgpFr2skzS91mgFjcz+4eNdI5+eaCAvklR5Mob0i
   1TO35017N93GPEdZ8691yvCViYLWTf4mta4a9y6rotD6jK1dbzuJeBJR1+YZCV32e79Qavq9okyh
   uC78ex7wwvD/HXkSK2qro1C4l1sjafcnN6eT5CppnDd/+uxTW5rwnVwio6MnzasxR41hiCqg5ts3
   3w8fzKyarRvofBFrXw8ifXfSNN/dVMK4hG6EIintuftOtYj88PCh4XeS3IqKIy6+lOf9LbsmvhBq
   +r2iTKH4D+ACgqk7loXxip15Eitqq6tQlEVSgRhMUbFyXoHRPsisffrs+D78afvLbwlHJBcjFHOD
   /071oNWyPrGgnXPPrJgNYE9OTuZajzltrXUqYVxC3ppvmoBzXJA7mNp8zn9vtiz1+gi9dOmI+lCm
   UOwPrAeODfcPB1bnSayobZCEIk1tqbWgjZoTZ74LKHp+o+i5c9LUhqPdHMfNE608rqfAbdbaYyjJ
   x5+lN0/6QrPxbFMeN7I6KoDe7I/PU/NNk/dTU1ORgxbnZrOd30IsIt1uUCugnuQVikVJXWfNbBHw
   LXc/rqk77e3A7UnXiXRMT0+zdu3zuOeeIHuvuup5bNv2CdasWZPpPsuXH8rll7+FzZsvYNeuXwEP
   5+abfxpx3sFceeWluWydmDiTq6/ewJ49jU82EgzY387IyLmccsrJTExsjbX9uOOO49Zbz+Woow7j
   He/4BACbN18AwK5dJ7Bjx8+A84ENs9dcddU2zj47+H96errp/DvYs2fu3D17gns10l6zZk3mPAxT
   Ce95PrsbDNTRAAAUCklEQVR3w9q1L+bEEx/B8uWHMjFxZsS5W9m9+11s3w5XX72Byy/fmjt/Yy0K
   n/uoow7j7rv/hnvuCT4fGno9e/e+LPd927/P4eGNTExsLcDiwOZ16zaE39Fc3uT7TkQt6KQkwKeB
   o/KoUFkbA9KiSNv1r5PrKdr90Xpe0vQNaWvoY2MrmmIJyX3lm69pD1xHxUA6TcyXtStn1h5U8WMY
   WhcWmrMjeVrwLD2s0o5tWbx4WYubLbhuuaf5juPSLqPW38segCIblOh6+gpwN/BF4Ipw25YnsaK2
   fheKxg808PG3FnhRA6oa/esbk/XFrYQWPbZgNJVbIm1/+6RYR/Q16dZZaF/adP4YhfQ9dPLEC+Ld
   O6fP5mUj/xvfRVphS0o/Ke7RySVW5PohRSKhqC9lCsV41JYnsaK2fhaK+X3fW0eyto86DQrQg8MC
   dyLzQLFuuojG3bPTOgGdYipxtsSJVZINzcHsPD2Q2nuTRX83Uy019+HhQ8MuydHfS5b003W5be3x
   tHjxIU1iVd4qcHlbHOrSWl9KE4o6bv0sFNGF+anzapPuUYPKgsV6kgq9+QVdsrjksTdbH/35hWxW
   W+IKnrjP09ocdX2jtTY2tjJsOU14MIV6s7utfT2Lg1oC70UJRRqXWJqFhfK4lrot7BXMridltiie
   BHw9dD/dC+wF7sqTWFHboAlF3EjWOFFJ1220+1XHGvfL48aJ863ntSVqHEmc+yetzZ0K9HZ32NzY
   jPkFd3vvp6LmWOrsEitnqVK5jwaTMoXim8CxwA6C1edeCpyXJ7Gitn4WiqTacXsNLM433+saWp7a
   YRBXiV7ZrZv00sY/snc7nl8Yxgl1VOC9fX6j9vSbBbw9xhMXc4rrEJBm4r1uC3oJxWBSqlCEf69v
   +uy6PIkVtfWzULhHz3uUdubWbibey9oTqBvXQdoabfukgp2uK9Kt1R7/SeNCatiZZcbUTnGpNG6x
   RqssSmjytJY6UXWcQa6rcihTKL4MPAC4GHgn8AY0MrtQktxRcbXNrGTtidNtIZEmCB41jUSwolt8
   KyFww7TGC/JMQje/8F7mzUu5dsqHRkGWJqAc38kgvgCPy7/mtJOeuQj3X1WFddUiNciUKRRHA8PA
   gcAm4N3Aw/MkVtS2EISiyKma49LIE2DtJr32Lq3BkqSt5+yzT6OG3zrnU9GB+rQ9xDoVlml6peX5
   fotwO8a5rupe8MrtVR6l9noimMbjkXkSKGMbNKHI45rISq+Fon18gFmjB9HcPaPmigrmfJrLh8YA
   tCQ3UJ5ab1qh6PSMaVyD7XkBI262xMfGVia2CMwO8LkpxoP8a3RQaI99dBOLqRv9Zm8/UWaLYi3B
   IkK3hPtjGnBXPFlcGXnv30vXU1A4zo04Nlsy78cfuJlaXU/tz95wt8S5YvLaGVV4N7cIii6AW6cL
   7zyrbSAUy5rsO8RhfUsvrKjp2ouIV1QdH5DrqTzKFIpvEcwau6Pps2/nSayobRCFokGZP5JeBrPn
   F04TkSOvG8HsYGqQA2ILtKh8iZr4MEtroFXIDpjtnVVU99qs5zbnedSzBWLW/Fm6Ue+9riQUQdVi
   NaiUKRTXhn+bheL6PIkVtQ2yULgPxo8kugWwIva5gvOTF1Zqz5e5a7JPlZ1UcKcVgDTB7k492+Lu
   FTWf1fyYTjqhiLInT76I/qdMofgg8ELghnA8xfuBf86TWFHboAvFINBeY2+e7C6KuQIq/bTZ3Sw/
   WoRQNJ6zvQBO6u6cLJTJwf/2brlpXU9ZkFAMNmUKxf7A3wPfCLe3A/vlSayobdCFom4tijz2tMcA
   Os1smsfl0U2h1qk10MmWuDyJ7sKbb2nSqDm1oloqRb4rRcWn6vT+ijkKF4qwS+zrgf8DvBLYN08C
   ZWyDIhRZaqNV2pjHnjyFeNYCptvab1Jh3zxbb5RIpBkgmWU51jp970n5Evf99GtX3IVGGUJxCfAR
   4C+ATwHvy5NAGdsgCEVcwVC3pn/eOEAvnqOMwjXNPeOeLW68RFRrKq6SUNeaePoWWPq4ieg9eYUi
   aYW74939JAAz+1eCiQFFQWzefEHkKm11Y9euOwgG578r/OSN7Nr1yI7XlbmCWoM1a9Zw+eVbZ/Mt
   aYW9tMR9L/nvexsnnvgIAFavXg/AypWP5e1vf3/kCnB1XQUuKV9aj22r0sxSaV5lcWLizNp+V2WQ
   JBT3Nf5x9/vMrAfmgJltAv4c+GX40VnuPtWTxGtALwrYbCwiEIkNTZ99qONVa9as4eyzX8O7330u
   AG94w2tK+WFVUbgmfUdf+MLz2bu3cebrgUO56647W5YG/cIXGsuYFiVGdeJM4EWze9W/v8Ww4Jd3
   jWtqAPcDv23a7mv6v7RpxoFzgDd0OKfY9lgFZOlaWSV5XUh18blnzcu0dscFlYNBhAd5YwAfbHGz
   9vEPDRdN/7hnsgT/i5hWvm7UzSWcF8qcwqOXWygUEx3OKSzjqqROghBHL4PZRZPX9m7Fxeygtmef
   77dPWso1y/P1agBlp3v0w7vcDXV4n4tg0ITiFmAncBGwLOKcIvNOdCBPIRD3w+plgVJWr6jO6bQL
   w/xR6d3OCpxFBOvSuutnBiUP8wpFUoyiNMxsO3BYxKGzgQ8Abwv3zwU2Ay9vP3HTpk2z/4+PjzM+
   Pl60mSIkTxwgyo+/cuVr+sbP251PegVDQ6+fjVUMD3+Es8+e4KqrgkBvI+h+9tn57csSdC8+QL/w
   KKPjRC+YmZlhZmam+xvlUZdebQRTnN8Q8XlhCiu6J22/+14337upBWYdnd2eTlHriBRh36C4TUT3
   MECup8Ob/n898LGIc4rLOdEVWQrjKgqsvK6urLb22kcfN/K9HwZxiuoYJKH4MHA9QYziU8ChEecU
   mXeiC7qpeS9efMjsRIF1KbjmRhinW3I06R5lPlcgFK1zabXPBVXXnnSiOgZGKFIZLaGoDXlr3t0U
   xHH3HRtbGU67saKQAHEeIetV7T0q34M5phqfTTmc6iMjoxIGMYuEQlRC3oIx60I+nWxoXYRoecfZ
   aqMowjXWK/daslBMefN07XI1iQZ5hWKo+3C4WMg0eoOsWrWNVau2peoZND09zc6dNxJMI7aWoDfO
   Dblt2Lz5Au655x/C+2wA3sU99xxXyylRimJi4kyGhzcCW4GtDA9v5A1veGn42Sag0csp6Lk1yHkh
   yqeS7rFiYXPWWe9g79730DwtyNDQBBMTH63OKIqZPqVXU7DEddd83OMexwte8Jfs3t19Ggt5biPR
   Rp5mSNUbcj3Vhqyup6mpqaZRyXNuk8ba2HltKML11LhX0SOYq+gRVcSaEuopNXigGIWogqw++blp
   y+OXPM1DEcHsMqiqwO1WnDT2YjDJKxRyPYkKOInAt34B8HNOPvmErt0azaPHGy6TzZsvqNxlUtWo
   6DpPWS76DwmF6IqsPvm5888H1jI8vJF3vKM4H/6Cnw66IOo33b2olDzNkKo35HqqFXlmWy3LZ183
   l0k/+/o1SG/wIKfryYJr+wsz8360W5TP6tXr2b690eUWIOi6e+WVl847t1e9etR7SNQFM8PdM69C
   J6EQA0W762l4eGOk6ynteUIMEhIKIULS1OCztDyEGBTyCoWC2WLgUI8fIYpFU3iIBUnUFBgTE2dW
   bVZXTE9Ps3r1elavXs/09HTV5ogBQq4nsWAZpCCzYi4iDYpRCLGAUcxFpCGvUMj1JIQQIhEFs4UY
   ADSSWpSJWhRiQTGoAd8864IUxaDmqZhDMQrRwiAFeNtRwLd4lKf9hYLZomsG/UevgG/xKE/7Cw24
   E11T1ZTYQoh6I6EQCwYFfItHebowkOtJzDLoricY7BhMVShP+4e+ilGY2RnAJuA44PHu/q2mY2cB
   LwPuB17r7ldGXC+hKAn96IUYXPpNKI4D9gL/Akw0hMLMTgA+BjweOAL4PPAId9/bdr2EQgghMtJX
   I7Pd/SZ3/37EoWcDH3f3e939FuBm4Ak9NU4IIUQLdRtw92Dgtqb92whaFkIIISqitF5PZrYdOCzi
   0Jvd/YoMt4r0MW3atGn2//HxccbHx7OYJ4QQA8/MzAwzMzNd36fSXk9m9iVaYxRvAnD388L9KeAc
   d7+27TrFKIQQIiN9FaNoo9nobcDzzGyxmR0DHAt8rRqzhBBCQEVCYWbrzOynwKnAf5jZ5wDc/Ubg
   EuBG4HPAq9V0EEKIatGAOyGEWCD0s+tJCCFEjZFQCCGESERCIYQQIhEJhRBCiEQkFEIIIRKRUAgh
   hEhEQiGEECIRCYUQQohEJBRCCCESkVAIIYRIREIhhBAiEQmFEEKIRCQUQgghEpFQCCGESERCIYQQ
   IhEJhRBCiEQkFEIIIRKRUAghhEhEQiGEECIRCYUQQohEJBRCCCESkVAIIYRIpBKhMLMzzOw7Zna/
   mT226fOjzWyPme0It/9bhX1CCCHmqKpFcQOwDvhyxLGb3X0s3F7dY7v6jpmZmapNqA3KizmUF3Mo
   L7qnEqFw95vc/ftVpD1o6Ecwh/JiDuXFHMqL7qljjOKY0O00Y2ZPrtoYIYRY6Cwq68Zmth04LOLQ
   m939ipjLfg4c6e53hrGLT5nZie7+27LsFEIIkYy5e3WJm30JmHD3b2U5bmbVGS2EEH2Mu1vWa0pr
   UWRg1mgzWw7c6e73m9nDgGOBH7VfkOdBhRBC5KOq7rHrzOynwKnAf5jZ58JDK4GdZrYD+Hfgle7+
   6ypsFEIIEVCp60kIIUT9qWOvp1nM7DQzu8nMfmBmG2PO+cfw+E4zG+u1jb2iU16Y2QvDPLjezK4x
   s0dXYWcvSPNehOc93szuM7PTe2lfL0n5GxkPexJ+28xmemxiz0jxG1luZlNmdl2YFy+pwMzSMbMP
   mtkdZnZDwjnZyk13r+UG7APcDBwN7AtcBxzfds6fAJ8N/38i8NWq7a4wL54EHBj+f9pCzoum874I
   fAZYX7XdFb4Xy4DvAA8J95dXbXeFebEJeEcjH4BfAYuqtr2EvPhDYAy4IeZ45nKzzi2KJxCM0r7F
   3e8FPgE8u+2ctcBWAHe/FlhmZof21sye0DEv3P2/3P034e61wEN6bGOvSPNeALwG+CTwy14a12PS
   5MULgEvd/TYAd9/VYxt7RZq8uB1YGv6/FPiVu9/XQxt7grt/Bbgz4ZTM5WadheII4KdN+7eFn3U6
   ZxALyDR50czLgc+WalF1dMwLMzuCoJD4QPjRoAbi0rwXxwIjZvYlM/uGmb24Z9b1ljR5cSFwopn9
   HNgJvK5HttWNzOVmHbrHxpH2x93eVXYQC4XUz2RmTwFeBqwoz5xKSZMX7wXe5O5uZsb8d2RQSJMX
   +wKPBf4Y2B/4LzP7qrv/oFTLek+avHgzcJ27j5vZKLDdzE72hTmgN1O5WWeh+BlwZNP+kQTKl3TO
   Q8LPBo00eUEYwL4QOM3dk5qe/UyavDgF+ESgESwHnm5m97r7tt6Y2DPS5MVPgV3uvgfYY2ZfBk4G
   Bk0o0uTFHwBvB3D3H5rZj4FHAt/oiYX1IXO5WWfX0zeAY8OpxxcDzwXaf+jbgD8DMLNTgV+7+x29
   NbMndMwLM3socBnwIne/uQIbe0XHvHD3h7n7Me5+DEGc4lUDKBKQ7jfyaeDJZraPme1PELy8scd2
   9oI0eXET8DSA0Cf/SCIG9C4AMpebtW1RuPt9ZvZXwDRBj4aL3P27ZvbK8Pi/uPtnzexPzOxm4HfA
   Sys0uTTS5AXwVuAg4ANhTfped39CVTaXRcq8WBCk/I3cZGZTwPXAXuBCdx84oUj5Xvw98CEz20lQ
   Sf5bd99dmdElYWYfJxi8vDwc2HwOgQsyd7mpAXdCCCESqbPrSQghRA2QUAghhEhEQiGEECIRCYUQ
   QohEJBRCCCESkVAIIYRIREIhBhIzuz+cWnuHmX3LzI4ys2vCY0eZ2fObzj3ZzJ6eI40ZMzulAFsL
   uY8QZVHbAXdCdMnv3b19nv3G/FfHEMyq+vFwf4xg2o/PkQ2nmLnFirqPEKWgFoVYMJjZ3eG/5wF/
   GLY2/hb4O+C54f4ZZvbAcPGXa8PWyNrw+mEz+4SZ3WhmlwHDtE2uFi6ec0nT/riZXRH+/wEz+3q4
   aM6mDjZiZn9qZh8K/z/EzD5pZl8Ltz8IP1/Z1nI6oKDsEmIWtSjEoDJswdrrAD9y9/XM1do3Am90
   92cBmNkdwCnu/tpw/++BL7j7y8xsGXCtmX0e+Avgbnc/wcxOAr7F/JbA54F/MbPhcCK+5zLXcnmz
   u99pZvsAnzezk9y9fRUyj/n/fcB73P2acF6vKeAEYAJ4tbv/VziX0/9kzSghOiGhEIPKngjXU4P2
   KZbbpyJfDTzLzN4Y7j8AeCjBymHvA3D3G8zs+vYbh3MOTQFrzexSgtXEGvd5rpm9guB3dzhwPBC7
   XGUbTwOOD+fxAlhiZg8ErgHeY2YfBS5z90GcPVlUjIRCiOj4wOntazaEhXSatS0+AfwVsBv4hrv/
   zsyOIaj9P87dfxO6lPbrYMtwc/LAE939nrbzzzezzwDPAK4xszXu/r0UNgqRGsUoxELkt8CShP1p
   4LWNnabF579MEATHzB4FPDrm/lcRLBb0CubcTksJZuq8K5ziOq6X1R1mdpyZDQHrmBOOK9tsekz4
   d9Tdv+Pu7wS+TjB1thCFIqEQg0pUK6Hx2U7gfjO7zsxeB3wJOKERzAbOBfY1s+vN7NsEwW4IllY9
   wMxuDD+LXPDG3fcCnwFOC//i7juBHQRrInwUuDrG7jeF11wD/Lzp89cCjzOznWb2HeDM8PPXmdkN
   4dTZ95C955YQHdE040IIIRJRi0IIIUQiEgohhBCJSCiEEEIkIqEQQgiRiIRCCCFEIhIKIYQQiUgo
   hBBCJCKhEEIIkcj/B7O++MHdETnRAAAAAElFTkSuQmCC
   "
   >
   </div>
   
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p>Histogram of standardized deviance residuals:</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[14]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="kn">from</span> <span class="nn">scipy</span> <span class="kn">import</span> <span class="n">stats</span>
   
   <span class="n">fig</span><span class="p">,</span> <span class="n">ax</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">subplots</span><span class="p">()</span>
   
   <span class="n">resid</span> <span class="o">=</span> <span class="n">res</span><span class="o">.</span><span class="n">resid_deviance</span><span class="o">.</span><span class="n">copy</span><span class="p">()</span>
   <span class="n">resid_std</span> <span class="o">=</span> <span class="n">stats</span><span class="o">.</span><span class="n">zscore</span><span class="p">(</span><span class="n">resid</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">hist</span><span class="p">(</span><span class="n">resid_std</span><span class="p">,</span> <span class="n">bins</span><span class="o">=</span><span class="mi">25</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">set_title</span><span class="p">(</span><span class="s">&#39;Histogram of standardized deviance residuals&#39;</span><span class="p">);</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   
   
   <div class="output_png output_subarea ">
   <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAW8AAAEKCAYAAADdBdT9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz
   AAALEgAACxIB0t1+/AAAHDFJREFUeJzt3Xv8HHV97/HXmwANECAgnBC5hQocihdCH4oXaF2k0Dyo
   cumpiLRClFJqe4DDqcjFWn7CURS5+HhQy6PniBBA4KAoEq2SCFkuakE8CbdAg5VwM/lhIEAiXoB8
   zh/z3TAse5n97e5vdzbv5+Oxj9/Mzsx3PjM7+/nNfmZ2v4oIzMysXDYadABmZtY5J28zsxJy8jYz
   KyEnbzOzEnLyNjMrISdvM7MS2mCTt6QHJP3xoOMYJElHSnpC0hpJ+ww6nhpJY5Ku6mF7V0g6Nw3/
   kaSHe9V2bh3rJP1+wXl7sn2SLpX0j922M4xavT8lVSQ90aP1VCUd34u2JttIJm9JyyUdVPfcXEl3
   1MYj4i0RcXubdmalN+VI7ifgAuDvImLLiLi36EL5ZNgnvf7yQdTajIg7ImKvHrc/kXi6byTi4xHx
   v3rR1rAp8v7s1aro/fE2KUY1KfX6BVEP23q1UWlKP9otuG4BuwBLBxVDP7T4R9uX19CaG+TxvSEY
   1eTdyGuSeTo7f18a3k/SPZKel7RS0gVpttp//udSaeGdyvxjWn5c0jxJW+XaPVbSY5JW5earrWdM
   0jckXSXpeeA4Se+Q9GNJqyX9QtIlkjbJtbdO0sclPSLpBUnnSHpTWuY5Sdfl56/bxoaxSvo9YA0w
   BbhX0iNNlr84Lfe8pPskvVnS3wDHAJ9M++Tbad4zJP0sxfigpCNy7cyVdKekL0p6VtLPJc3JTd9N
   0m1p2QXAdnVxfF3SirS9t0naOzftilQ++DdJa4GKpH0l/b/U3nXA1Nz86z9yS/pQ2oba47eSFqVp
   vyfpgvRarkzryLdzWnq9npT0sUb7r4Pte5ekH6VjYImk9+bi+0ndvKfm9nm+HLSNpO9Iejrt4/mS
   dswtV03Hzp0pjpslvSE3/YBcDI9LOq7IfqiLba6kH0q6SNIq4GxJmzZbXtJ2KebVkp6RdHuurfWf
   niVtlrb1WUkPAu+oW+9rSlad7Je6dnZPr9Nzkn6Zjp3hFREj9wAeBQ6qe24ucEfdPO9Lwz8G/jIN
   bw68Mw3vCqwDNsot9zHgEWAWsAVwA3BlmrY3WVJ8D7AJ8EXgd7n1jKXxw9L4VOAPgf3I/pHuSnYm
   fEpufeuAbwHTUvu/BW5N698KeBA4tsl+aBprru3fb7LsnwL3AFul8f8K7JCGLwfOqZv/L3LTjwLW
   AjNy+/53wPFkZ8B/CzyVW/bHZCWcTYA/Al6oi3Nuin8T4GJgcW7aFcBzwLvT+FbAY8ApZP+c/lta
   9zlpegV4osH2bpn2/Qlp/GLgRmB62vc3AZ9L0+YAK9PrsTlwTZt92XT7gB2BVcCcNP4nafwNqe0X
   gN1zbf0EOKr+dQC2BY4kO6amAdcD38otV03Hwu5pnkXAebnj/AXgQ2mfbQvs024/NNjOucBLwN+T
   Hc9T2+zH84BL0zqnAPs3eX9+HrgttbET8ADweLPjuMP9sgj4WBq+FjgzDW8KvGfQuaxlnht0AH3Z
   KFhOlkRX5x6/Am5vcnDcRpZYt6trZxavT963AH+bG9+TLDlMAf4J+Fpu2mZkyTafvKttYv8fwDfr
   Dsx358bvAU7LjV8AXNykrWaxbpRru1nCORD4D+Cd+e1P0y4Hzm2zHYt59Z/UXOCR3LTN07r/C1np
   5iVgs9z0rwFXNWl3elp2yzR+BXBFbvofk/vHkJ77IS2SN1mi+Q7w5TQusn8++YTwbuDnafir5BIY
   sEezfdli+2rJ+3Ry/6jSc98n/UMGrgI+nVvPC8DUdq8DMBt4Nje+CDgrN/5x4Htp+EzghgZttNwP
   DeafCzxWdHngM2SJ/U0N2sq/P/8TOCQ37YT8a1i/7yewX2rJex7wr8COrY7tYXmMatkkgMMjYpva
   A/g7mtc9jydLbA9JulvSn7VoeybZmV3N48DGwIw07cn1QUT8Gnimbvkn8yOS9kwf61YoK6V8luys
   K288N/zrBuPTJhBrSxGxCPhn4MvAuKR/lbRls/mVlYsWp4/Aq4G31G3HylzbL6bBacAbgdVpX9Ws
   j1nSFEmfV1aSeZ7sTQ2vlh6C1+7TNwJP1YX3GK19luzM/uQ0vj3ZP5if5rbne7l1zgTydzs83qLt
   ZttXOxZ3BT5YW09a1/7ADmn6NcCH0/AxZGeNv6lfiaTN02u0PO2n24CtJeWP+ZW54fxxszPw8wax
   t9sPjeT3S7vlvwj8DFgg6T8lnd6kzTdSfH+/RsH9UvNJstflbmV3u3y06HoGYVSTdyNNL1hFxM8i
   4piI2B74AvANSZvR+KLnL8jOyGt2AV4me2OsIPtYl60wa6M+Ede3eSnZx/XdI2Jr4FP07nVpFut4
   w7nrRMQlEfF2svLAnsBptUn5+STtCvxvso/L26Z/lg9Q7CLhCmAbSZvnnts1t45jgMPIymBbA7vV
   VpsPta69+prmrs1WLulosnLBX0TEK+npVWTJbe/cCcD0iKhd21hBti9r8sOdbt/jZJ8ytsk9toyI
   89P0HwDbK7uV82iyZJ5Xa+cfyF6j/dJ+ei/ZPiryGjwOvKnB8+32QyP516Ll8hGxNiI+ERFvInuN
   /6ekAxu02W5/v0j2T6JmJhPYLxExHhF/ExE7AicC/6KCt38OwoaUvJuS9FeStk+jz5O98OuAX6a/
   +QP7WuBUZbcRTgM+B1wXEevIasofkPRuSZuSlUnavXmmkZV4XpS0F9nH2bYhNxmu1yrW1iuQ3q7s
   Au0mZG+O3wC15DYO5A/qLcj22Spgo3TG8pYC20FEPEZWCvqMpE0kHQC8PzfLNLLS07OStkjb8JpQ
   68Z/BLws6eTU3p9Td4Ert437ApcAR0bE+k9Iaf/8H+BLteNC0o6SDkmzXA/MlfQHKSmf3cX2XU12
   zBySPmVMVXZRdce0/EvA18nKY9sAC+u2vbb908gS5fOStm0SU7Nj5RrgTyR9UNLGkt4gaZ8C+6Gl
   dstL+rN0kVBk5aBXyN5v9a4HzpQ0XdJOwEl105cAf5n23xyy0llNkf1CiueDqX3IrqPU8sBQ2pCS
   d9D89sE/BR6QtIbsAsvREfHb9PH+s8AP08e+/cjqnVeR3Ynyc7LEdhJARDyYhq8jO+tdAzxNlnya
   xfAJsrPLF8jOXq+rm6dRzPXTm21X01hbtF2zVYrnWbJrCKvIPuYCXAbsnfbJNyNiKXAh2YW5lWSJ
   +842MebHjyGrrT9Ldt1gXm7alWRlhqfIzuZ/TIvtT8nuz8nqr8+QXTy9ocm6Dyerod+pV+84+W6a
   djrZR/p/Tx+3F5KdwRER3we+RHbheBnZtYVW+7Lp9kXEkymOs8iOlcfJzhbz781rgIOAr9f9481v
   +5fIrrGsIvsH9r0GMTXcbxHxOHBoWu8zZNcr3tZuPzTQ6HVutfweaXxNivnLEXFbg3Y/Q3YMPEp2
   PeDKuvWcAnyA7NrWMWQX+GuK7Jeat6c41wDfBk6OiOVN5h04pUJ965my+zXvAZ6MiA9IGgP+muzM
   FLIrtN/vW5Qllc52V5OVRNrVXc3MCtu44HynkNVlaxesArgoIi7qS1QlJukDZGdiIvuoe58Tt5n1
   WtuySaoBHQp8hVdrZkUvhGyIDiP7iP8UWa386MGGY2ajqEjN+2Kyuwzqa20nSbpX0mWSpvcluhKK
   iBNyV9UPjoiG3140M+tGy+Qt6f3A0xGxmNeeaV9KdsvWbLLbeC7sW4RmZvY6LS9YSvoc8BGye4On
   kt2BcENEHJubZxYwPyLe2mD59ldDzczsdSKiZWm65Zl3RJwVETtHxG5ktdtbI+JYSTNzsx0J3N+i
   jaF/nH322QOPYRRidJyOc9gfZYmziKJ3m0BWNqm1en76xleQ3Xt5YgftmJlZlwon74iokv0yGRHx
   kT7FY2ZmBWxI37BsqlKpDDqEtsoQIzjOXnOcvVWWOIso9A3LCTcuRT/bNzMbRZKIbi5YmpnZcHLy
   NjMrISdvM7MScvI2MyshJ28zsxJy8jYzKyEnbzOzEnLyNjMrISdvM7MScvI2MyshJ28zsxJy8jYz
   KyEnbzOzEuqkMwaz15Fa/vBZQ/6lSbPuFTrzljRF0mJJ89P4tpIWSlomaYF7j9/QRQcPM+uFomWT
   U4ClvPruOwNYGBF7ArekcTMzmyRtk7eknYBDga+Q9WMJcBgwLw3PA47oS3RmZtZQkTPvi4HTgHW5
   52ZExHgaHgdm9DowMzNrruUFS0nvB56OiMWSKo3miYiQ1LSYOTY2tn64UqmMVB9yZma9UK1WqVar
   HS3Tsg9LSZ8DPgK8DEwFtgK+CbwDqETESkkzgUURsVeD5d2H5YjL7jbp5DWW7zYxa6PrPiwj4qyI
   2DkidgOOBm6NiI8ANwHHpdmOA27sRcBmZlZMp1/SqZ0yfR44WNIy4H1p3MzMJknLsknXjbtsMvJc
   NjHrva7LJmZmNpycvM3MSsjJ28yshJy8zcxKyMnbzKyEnLzNzErIydvMrIScvM3MSsjJ28yshJy8
   zcxKyMnbzKyEnLzNzErIydvMrIScvM3MSsjJ28yshJy8zcxKqG3yljRV0l2SlkhaKum89PyYpCcl
   LU6POf0P18zMoGBPOpI2j4gXJW0M3Al8AjgIWBMRF7VYzj3pjDj3pGPWez3rSSciXkyDmwJTgNW1
   dUw8PDMzm6hCyVvSRpKWAOPAooh4ME06SdK9ki6TNL1vUZqZ2WtsXGSmiFgHzJa0NXCzpApwKXBO
   muVc4ELg+Pplx8bG1g9XKhUqlUpXAZuZjZpqtUq1Wu1omY57j5f0aeDXEXFB7rlZwPyIeGvdvK55
   jzjXvM16ryc1b0nb1UoikjYDDgYWS9ohN9uRwP3dBGtmZsUVKZvMBOZJ2ogs2V8VEbdIulLSbLLT
   rkeBE/sYp5mZ5XRcNumocZdNRp7LJma917NbBc3MbLg4eZuZlZCTt5lZCTl5m5mVkJO3mVkJFfqG
   pW1YsjtIzGyYOXlbE0Vv53OiNxsEl03MzErIydvMrIScvM3MSsjJ28yshJy8zcxKyMnbzKyEnLzN
   zErIydvMrIT8JZ0NhL81aTZaWp55S5oq6S5JSyQtlXReen5bSQslLZO0wD3Hl0UUfJjZsGuZvCPi
   N8CBETEbeBtwoKQDgDOAhRGxJ3BLGjczs0nStuYdES+mwU2BKcBq4DBgXnp+HnBEX6IzM7OGivQe
   v5GkJcA4sCgiHgRmRMR4mmUcmNHHGM3MrE7bC5YRsQ6YLWlr4GZJB9ZND0lNC6VjY2PrhyuVCpVK
   ZcLBmpmNomq1SrVa7WiZjnqPl/Rp4NfAXwOViFgpaSbZGfleDeZ37/FDorNe3vs1bza/jwmz1rru
   PV7SdrU7SSRtBhwMLAZuAo5Lsx0H3Nh9uGZmVlS7sslMYJ6kjcgS/VURcYukxcD1ko4HlgNH9TdM
   MzPL66hs0nHjLpsMDZdNzMqj67KJmZkNJydvM7MScvI2MyshJ28zsxJy8jYzKyH/JKxNuk5+ntZ3
   ppg15uRtA9DJbYhm1ojLJmZmJeTkbWZWQk7eZmYl5Jp3SblPSrMNm5N3qXX2myJmNjpcNjEzKyEn
   bzOzEnLyNjMroSIdEO8saZGkByU9IOnk9PyYpCclLU6POf0P1zY0kjp6mG0o2nbGIGkHYIeIWCJp
   GvBT4Aiy3nPWRMRFLZZ1Zwx90lnnCtC/DhaGJY5sfh9vNgqKdMZQpPf4lcDKNLxW0kPAjrV1dB2l
   mZl1rKOat6RZwL7Av6enTpJ0r6TLah0Vm5lZ/xVO3qlk8g3glIhYC1wK7AbMBlYAF/YlQjMze51C
   X9KRtAlwA3B1RNwIEBFP56Z/BZjfaNmxsbH1w5VKhUqlMvFozcxGULVapVqtdrRMkQuWAuYBz0TE
   qbnnZ0bEijR8KvCOiDimbllfsOwTX7BsPL+PNxsFRS5YFkneBwC3A/fx6jvpLODDZCWTAB4FToyI
   8bplnbz7xMm78fw+3mwU9CR5dxmAk3efOHk3nt/Hm42CIsnb37A0MyshJ28zsxJy8jYzKyEnbzOz
   EnLyNjMrISdvM7MScjdoNlI6+VlY31ZoZebkbSOmk3vIzcrLZRMzsxJy8jYzKyEnbzOzEnLyNjMr
   ISdvM7MScvI2MyshJ28zsxJy8jYzK6G2yVvSzpIWSXpQ0gOSTk7PbytpoaRlkha493gzs8lTpBu0
   HYAdImJJ6kH+p8ARwEeBVRFxvqTTgW0i4oy6Zd2TTp+4J53u2/axacOqJz3pRMTKiFiShtcCDwE7
   AoeRdUxM+ntEd+GamVlRHdW8Jc0C9gXuAmbkOhweB2b0NDIzM2uqcPJOJZMbgFMiYk1+WqqN+DOo
   mdkkKfSrgpI2IUvcV0XEjenpcUk7RMRKSTOBpxstOzY2tn64UqlQqVS6CniUdfJzpmY2OqrVKtVq
   taNlilywFFlN+5mIODX3/PnpuS9IOgOY7guW3ensIuSwXCgcljg6b9vHpg2rIhcsiyTvA4Dbgft4
   9Z1xJnA3cD2wC7AcOCoinqtb1sm7A07ek9u2j00bVj1J3l0G4OTdASfvyW3bx6YNq57cKmhmZsPH
   ydvMrIScvM3MSsjJ28yshJy8zcxKyMnbzKyEnLzNzErIydvMrIScvM3MSsjJ28yshAr9qqDZKOrk
   Vxz9VXobNk7etgHr5DdWzIaLyyZmZiXk5G1mVkJO3mZmJeTkbWZWQm2Tt6SvShqXdH/uuTFJT0pa
   nB5z+hummZnlFTnzvhyoT84BXBQR+6bH93sfmpmZNdM2eUfEHcDqBpN8/5SZ2YB0U/M+SdK9ki6T
   NL1nEZmZWVsT/ZLOpcA5afhc4ELg+EYzjo2NrR+uVCpUKpUJrtLMbDRVq1Wq1WpHyxTqPV7SLGB+
   RLy1w2nuPb4D7j1+WNt2T/M2ufrWe7ykmbnRI4H7m81rZma917ZsIula4L3AdpKeAM4GKpJmk526
   PAqc2NcozczsNQqVTSbcuMsmHXHZZFjbdtnEJlffyiZmZjZYTt5mZiXk5G1mVkJO3mZmJeTkbWZW
   Qk7eZmYl5ORtZlZCTt5mZiXk5G1mVkJO3mZmJeTkbWZWQk7eZmYl5ORtZlZCTt5mZiXk5G1mVkIT
   7cNyg5X95rZtaDp93f3739Zvbc+8JX1V0rik+3PPbStpoaRlkhZseL3HR8FHJ/P6zT7c/DracClS
   NrkcmFP33BnAwojYE7gljZuZ2SRpm7wj4g5gdd3ThwHz0vA84Igex2VmZi1M9ILljIgYT8PjwIwe
   xWNmZgV0fcEyIkJS00Lf2NjY+uFKpUKlUul2lWZmI6VarVKtVjtaplDv8ZJmAfMj4q1p/GGgEhEr
   Jc0EFkXEXg2WG7ne48vZw3s/2x6WOPrZdudxjNpxb5Orn73H3wQcl4aPA26cYDtmZjYBbc+8JV0L
   vBfYjqy+/U/At4HrgV2A5cBREfFcg2V95j0UZ479bHtY4uhn2z7ztslV5My7UNmkiwCcvIci+fSz
   7WGJo59tO3nb5Opn2cTMzAbIydvMrIScvM3MSsjJ28yshJy8zcxKyMnbzKyEnLzNzErIydvMrISc
   vM3MSsjJ28yshJy8zcxKyMnbzKyEnLzNzErIydvMrIScvM3MSsjJ28yshLrqgFjScuAF4BXgpYjY
   rxdBmZlZa932Hh9kHRE/24tgzMysmF6UTVp21WNmZr3XbfIO4AeS7pF0Qi8CMjOz9rotm+wfESsk
   bQ8slPRwRNyRn2FsbGz9cKVSoVKpdLlKM7PRUq1WqVarHS3Ts97jJZ0NrI2IC3PPuff4oej9vJ9t
   D0sc/Wzbvcfb5Opr7/GSNpe0ZRreAjgEuH+i7ZmZWXHdlE1mAN/KzkTZGPhaRCzoSVRmZtZSz8om
   DRt32aSDeTudf1jaHpY4+tm2yyY2uYqUTbq9YDkS0qcHs57p5JhyoreJcPJer5OzMLN2fDxZf/m3
   TczMSsjJ28yshEpRNrn66mtYtOhHheffffddOPPMT/YxIjOzwSrF3SbHHHMC1167Fti/wNyPscce
   VZYt+0nh9vt3B8mw3C3Rz7aHJY5+tt3fODp5j3R6cd0XQ8tpxO42eR9Q5OdTfgJU+xuK2UD5Yqi5
   5m1mVkpO3mZmJeTkbWZWQiWqeRf3yCP3+FuTVhrDcKxOJAZfDB2skUzemU7vDjAblGE5VoclDivC
   ZRMzsxJy8jYzKyEnbzOzEuoqeUuaI+lhSY9IOr1XQZlZb0gq9Ohn252230m7nbTdabvDcCG5lW66
   QZsC/DMwB9gb+LCkP+hVYFavOugACqoOOoARU+1y+Sj46LbtRX1su5dxt4qzF3FPnm7OvPcDfhYR
   yyPiJeA64PDehGWvVx10AAVVBx3AiKkOOoCCqoMOoKDqoAPomW6S947AE7nxJ9NzZmbWZ93c5z1p
   nyumTIGpUy9h001vajvvK688x69+NQlBmZkN0IR/ElbSu4CxiJiTxs8E1kXEF3LzDH/hyMxsCLX7
   SdhukvfGwH8ABwG/AO4GPhwRD02oQTMzK2zCZZOIeFnSfwduBqYAlzlxm5lNjr72pGNmZv3R129Y
   SjpX0r2Slki6RdLO/VzfREn6oqSHUqzflLT1oGNqRNIHJT0o6RVJfzjoeOqV4Utbkr4qaVzS/YOO
   pRVJO0talF7vBySdPOiY6kmaKumu9P5eKum8QcfUiqQpkhZLmj/oWJqRtFzSfSnOu1vN2++vx58f
   EftExGzgRuDsPq9vohYAb46IfYBlwJkDjqeZ+4EjgdsHHUi9En1p63KyGIfdS8CpEfFm4F3A3w/b
   /oyI3wAHpvf324ADJR0w4LBaOQVYynB/AyeASkTsGxH7tZqxr8k7ItbkRqcBq/q5vomKiIURsS6N
   3gXsNMh4momIhyNi2aDjaKIUX9qKiDuA1YOOo52IWBkRS9LwWuAh4I2Djer1IuLFNLgp2bWvZwcY
   TlOSdgIOBb7C8P+ebaH4+v7DVJI+K+lx4Djg8/1eXw98DPi3QQdRQv7SVp9ImgXsS3ZiMVQkbSRp
   CTAOLIqIpYOOqYmLgdOAde1mHLAAfiDpHkkte1zvujMGSQuBHRpMOisi5kfEp4BPSTqDbAd+tNt1
   TkS7ONM8nwJ+FxHXTGpwOUXiHFLD/FG0tCRNA74BnJLOwIdK+sQ6O10nullSJSKqAw7rNSS9H3g6
   IhZLqgw6njb2j4gVkrYHFkp6OH1afJ2uk3dEHFxw1msY4BltuzglzSX7WHXQpATURAf7c9g8BeQv
   SO9MdvZtEyRpE+AG4OqIuHHQ8bQSEc9L+i7wdobvB0TeAxwm6VBgKrCVpCsj4tgBx/U6EbEi/f2l
   pG+RlSMbJu9+322yR270cGBxP9c3UZLmkH2kOjxdhCmDYavb3QPsIWmWpE2BDwHtf8/AGlL2e6SX
   AUsj4kuDjqcRSdtJmp6GNwMOZgjf4xFxVkTsHBG7AUcDtw5j4pa0uaQt0/AWwCFkNyk01O+a93mS
   7k81sQrwD31e30RdQnZBdWG6RedfBh1QI5KOlPQE2d0H35X0vUHHVBMRLwO1L20tBf7vMH5pS9K1
   wI+APSU9IWkgZbwC9gf+iuwOjsXpMWx3ycwEbk3v77uA+RFxy4BjKmJYS3wzgDty+/M7EbGg2cz+
   ko6ZWQm5GzQzsxJy8jYzKyEnbzOzEnLyNjMrISdvM7MScvI2MyshJ28zsxJy8jYzK6H/D6w8QxZ8
   wi+7AAAAAElFTkSuQmCC
   "
   >
   </div>
   
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p>QQ Plot of Deviance Residuals:</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[15]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="kn">from</span> <span class="nn">statsmodels</span> <span class="kn">import</span> <span class="n">graphics</span>
   <span class="n">graphics</span><span class="o">.</span><span class="n">gofplots</span><span class="o">.</span><span class="n">qqplot</span><span class="p">(</span><span class="n">resid</span><span class="p">,</span> <span class="n">line</span><span class="o">=</span><span class="s">&#39;r&#39;</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt output_prompt">
       Out[15]:</div>
   
   
   <div class="output_png output_subarea output_pyout">
   <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAYUAAAEPCAYAAACtCNj2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz
   AAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VPX1x/H3SUhICGFxRZE1IIoi7sWlEFGButZWRbGt
   +0rA/qy1KlpS1FZtbasE3Ou+VOtSNRWhrTFarStSEUHZFBBwAWQLgSTn98e9M5mEJCRkJpPl83qe
   PJm5cyf3XIxz8t3O19wdERERgJRkByAiIs2HkoKIiEQpKYiISJSSgoiIRCkpiIhIlJKCiIhEJS0p
   mFkPM3vVzD42s9lmNj48voOZzTCzT81supl1SVaMIiJtjSVrnYKZdQO6ufuHZtYReB/4IXAu8I27
   32pmvwK6uvvVSQlSRKSNSVpLwd1XuPuH4eP1wCdAd+Ak4KHwtIcIEoWIiDSBZjGmYGa9gQOAt4Fd
   3X1l+NJKYNckhSUi0uYkPSmEXUfPAJe7+7rY1zzo21IdDhGRJtIumRc3szSChPCIuz8fHl5pZt3c
   fYWZ7QZ8VcP7lChERLaDu1tdrydz9pEB9wNz3P3PMS+9AJwdPj4beL76ewHcvdV+TZw4Mekx6P50
   f23x/lrzvbnX72/pZLYUjgB+AvzPzGaGx64BbgaeMrPzgcXA6ckJT0Sk7UlaUnD3N6i9pXJMU8Yi
   IiKBpA80y9Zyc3OTHUJC6f5attZ8f6353uoraYvXGsPMvCXGLSKSTGaGN9eBZhERaX6UFEREJEpJ
   QUREopQUREQkSklBRESilBRERCRKSUFERKKUFEREJEpJQUREopQUREQkSklBRESilBRERCRKSUFE
   RKKUFEREJEpJQUREopQUREQkSklBRESilBRERCRKSUFERKKUFEREJKpdsgMQEZG6FRYWc8cd0ykt
   bUf79mWMHz+C448fmpBrJTUpmNlfgOOBr9x9UHgsH7gA+Do87Rp3n5acCEVEkquwsJjLL3+FBQtu
   ih5bsGACQEISQ7K7jx4ARlU75sAf3f2A8EsJQUTarDvumF4lIQAsWHATkyfPSMj1kpoU3P11YHUN
   L1lTxyIi0hyVltbcobNpU2pCrpfslkJtxpnZLDO738y6JDsYEZFkad++rMbjGRnlCblecxxovhOY
   FD6+AbgNOL/6Sfn5+dHHubm55ObmNkFoIiJNa/z4ESxYMKFKF1JOzrWMG1e9531rRUVFFBUVNeh6
   5u4NjTGuzKw38GJkoLk+r5mZJztuEZGmUlhYzOTJM9i0KZWMjHLGjTt2uwaZzQx3r7N7vtklBTPb
   zd2Xh4//DzjE3cdUe4+SgohIA9UnKSR7SuoTwDBgJzNbAkwEcs1sf4JZSIuAi5MYoohIm5L0lsL2
   UEtBRKTh6tNSaK6zj0REJAmUFEREJEpJQUREopQUREQkSklBRESilBRERCSqOZa5EBFptZpyb4Tt
   oaQgItJEmnpvhO2h7iMRkSbS1HsjbA8lBRGRJtLUeyNsDyUFEZEm0tR7I2wPJQURkSYyfvwIcnIm
   VDkW7I1wbJIi2poK4omINKF47Y2wPVrEfgrbQ0lBRKThVCVVREQaRElBRESilBRERCRKSUFERKKU
   FEREJEpJQUREolQQT0Qkjpp7FdRtUVIQEYmTllAFdVvUfSQiEictoQrqtiS1pWBmfwGOB75y90Hh
   sR2AvwK9gMXA6e6+JmlBiojUoKZuopZQBXVbkt199AAwGXg45tjVwAx3v9XMfhU+vzoZwYmI1KS2
   bqJOnVbXeH5zqoK6LUntPnL314Hq/4onAQ+Fjx8CftikQYmIbENt3USwudlXQd2WZLcUarKru68M
   H68Edk1mMCIi1dXWTdSp0x7ccMNwJk++PqYK6qgWM8gMzTMpRLm7m1mN5VDz8/Ojj3Nzc8nNzW2i
   qESkrak+frB27aoaz8vIKOf444cmLgmUl0Nq/ccnioqKKCoqatAlkl4628x6Ay/GDDTPBXLdfYWZ
   7Qa86u57VXuPSmeLSJOoafygW7fzgc6sWPHH6LGcnGu5/fYEtQo++wymToVnn4W5cyEzc7t+TH1K
   ZzfHlsILwNnALeH355Mbjoi0ZTWNH6xYcT8HHHABgwcnsJuovBxefhkKCuCDD+CCC+D117c7IdRX
   sqekPgEMA3YysyXAr4GbgafM7HzCKanJi1BE2rq6xg+mTcuP/wVXrYK//CVoGey4I4wbB88/DxkZ
   8b9WDZKaFNz9zFpeOqZJAxERqUX79mU1Ho/7NNOZM2HKFHjmGTjxRHjySTj00Pheox60ollEpA7j
   x49I3DTTzZvhiSfgyCPhpJOgb1+YNw8efjgpCQHqMdBsZv2Ape6+ycyOAgYBDydzlbEGmkWkKRUW
   FjN58oyY8YNjGzd+8OWXcPfdcM89MHAg5OUFrYN2ie28qc9Ac32SwizgIKA38A/g78A+7n5cnOJs
   MCUFEUmkhFQ6dYc33ggGjmfMgDPPhMsug332iU/Q9RCv2UcV7l5mZj8CJrv7ZDObGZ8QRUSSLzYJ
   rF27lOXLO1WZbtqoSqcbNsDjjwfJYNOmoFVwzz3QuXO8wo+r+iSFzWY2BvgZcGJ4LC1xIYmINJ2t
   1yFcB9xY5Zyg0un1DUsK8+cHM4gefjgYM/jDH+DooyGleQ/l1ie684AhwE3uvsjM+gCPJDYsEZGm
   sfU6hEZUOq2ogH/8A447Dg47DNLT4b33gimlxx7b7BMC1KOl4O4fm9nVQM/w+SKChWUiIi3e1usQ
   tmMK6urV8MADQcugc+dgbcEzzyR8oVkibDNtmdlJwExgWvj8ADN7IdGBiYg0ha3XIYwA6jkFddYs
   uOiiYCrpBx/Ao48GLYNzzmmRCQHqN/voA2A4QQ2iA8Jjs9193yaIr7aYNPtIRBolMri8bNnXLFxo
   lJTcFX2tW7fz2H33TLKzd956CuqWLUENooICWLQILr00KEGxa/Mv6Byv2Udb3H2NWZWfU9GoyERE
   kqD2RFBMZuZocnJ2o3v3bMaNO2frQeXly4NZQ3ffDQMGwM9/DiefnPC1BU2tPnfzsZmdBbQzs/7A
   eODNxIYlIhJfVWcZVZ9hNJSSkqF0734906bdUHnYHf7zn6D8xLRpcMYZMH067Ju0jpKEq89Q+Dhg
   H6AUeAJYC/w8kUGJiMRb1VlG25hhtHEj3HcfHHAAnHdeMJNo8WK4885WnRCgfrOPNgDXhl8iIi1S
   1VlGNc8w6lPxNVx5JTz4IBx+ONx6KxxzTIuYShovtSYFM3uxjve5u5+UgHhEROKi+irl+fNjd0uL
   zDC6CaOCEUznqg55HDHrGxhyEbz7LvTpk6TIk6uulsJtTRaFiEgcVR0/KAZeISjKECQCGEpn1nFR
   2v6MtcVsaZ/JhnN/RvvfTYQOHZIZetLVmhTcvagJ4xARiZuq4wfTCRJBYBAXMZb3GW0fseHwYXT/
   3Z0wZAhYnTM124y6uo+edvfTzGw2UH1RgLv7fokNTUSk4QoLi3nnnSUxR9rRji38kOfJo4B+zOcu
   LuG8IUfxbNEfkhZnc1VX99Hl4ffjgeopVCvHRKRZKSws5vrrH+aTT9LYtKkHALuyggt5lUu4l/n0
   o4A8nueHlJHGyE7XJzni5qnWIXV3/zJ8eJm7L479Ai5rkuhEROohMoYwc2Y3Nm2aymHszmPswyfs
   TXe6MIrjyOU1/sZplJEWv53TWqH6lLmYGSlvEXPsI3cflNDI6o5JZS5EJGrkyOsonj6BMzmDPJaQ
   zTqmMIoHyeA7ssjKepcBA3rXXLaiDWlUmQszu5SgRZBjZh/FvJQN/Cc+IYqINM6/7/8rx7/2PI9y
   N2/ThWuZzHRG4DEdIUceWW2lstSqrjGFx4GXgZuBX1E5rrDO3b9NdGAiIrWqqODtG//Iptsms9/a
   r/iAwQzh7yxkGcH001HRU4OuolG1/iipapvdRwBmlgrsSkwScfcvEhjXtuJR95FIGxJZiPbdF8sY
   uvBtLtzyORs8g8kcxhM8RQnvESSDyLqEGWRkfM7AgdlMmjS6TXYV1SQuVVLNbBwwEfgKiN1lIqFj
   Cma2mKDOUjlBpdZDE3k9EWk+IjOJPv10OSUlpQys6MRllHMG05hGP87hVN7kQeA3QAcg8qF/PZBK
   167zeeSRS5UMtkN9qqT+HBiQhC4jB3LdfdU2zxSRVqOwsJgLLniIb1ZUcDLtGcvHDGAzd9OPgVzC
   Cm4H8gl6tGNrGA0lkhwOPbSB+ylLVH2qPH1B8Bd7MmiJoUgbUlhYzJU/vY3zVyxgEc8wnm+4k0Po
   xedMYjgr6BqeGUkGDdglTeqlPi2FRcCrZlYIbA6Pubv/MXFhBdcA/mlm5cDd7n5vgq8nIsnizn0X
   XkXWg8/zZvlSnmYgxzOG/3EXwd4HaVRtFVQWtAtcrzGEOKlPUvgi/EoPv4ymWdF8hLsvN7OdgRlm
   NtfdX4+8mJ+fHz0xNzeX3NzcJghJROLp5Wdn8PGv72DEZ//hqM2lTOEGLmMZa4jd3ziSAEYCD1E9
   GaSlLWLQoM5MmnSBkkE1RUVFFBUVNeg99Zp9lGxmNhFY7+63hc81+0ikhYkdPN5l4xou9sWcyyre
   YzAFpDONo3B+QzB76KHwXd2InVGUmvoh6embSU3NJi0tiz59Oqpl0ADxmn20C3AVMBCi6dvdfXjj
   Q6z1mh2AVHdfZ2ZZBH8q/CZR1xORxCksLGb8+NtZuLCCYyjhMT7hCFbyEL04nNdYwIMEH0WR7qHI
   B/zDwIfAiWRlZTNgQFcmTfqlEkCC1af76DHgr8AJwMXAOcDXCYwJgjURz1lQyrYd8Ji7T0/wNUUk
   jiLJ4OuFmzmbZYzlC0pJpYCDGcNoNtIB6EdlQogdJwhmEmVmXszTT5+lRNCE6pMUdnT3+8xsvLu/
   BrxmZu8lMih3XwTsn8hriEji5OdP5bnfvsAVW+ZzJsuYQQ4XciJv0JPKSYVlMd9HECw+G0lkrUFK
   yiyuuqpt1ihKpvokhciMoxVmdgLwJUTnhYmIRP3jhX/z4kXXctrKBVxMCXdzMPtyMsujPc+x619j
   B5AjCWEGkEpm5idcddWx5OerIHNTq0+V1BOB14EewGSgE5Dv7i8kPrxaY9JAs0gzMuOxvzNr3LWM
   Xr2Iz8migGE8y95siU5UHEHl4HFENyoTwVeYLaNnz93Ya69ubbaKaaLVZ6C5Rcw+qk5JQST5CguL
   ue+iiZzy5cecyBqeYQ+mMJwP2YOqi8seojIBPAwsB0oASEnpSGZmx3AQWbOIEi0uScHMHqh2yAHc
   /bzGhbf9lBREkqOwsJgbJvyFfT9+m4vLvmRHyplCPx7gaFaTHZ4VGSOongw2AGXk5HTg9tvPVwJI
   grhMSQUKqVyslgmcQjCuICJtQGR9wfqP53Du5m94gSW8Txfy2Y9pfJ8K0sMzY1sHrwBnEySDO4Et
   5ORkcfvtKlLX3DW4+8jMUoD/uPthiQmpXjGopSCSYIWFxYwf92d6L/qKPJYwlOU8zAFM5SDmswuw
   GOhN7V1FwaBxevpsrrlmuAaNm4GEjCmY2V7AS+7erzHBNYaSgkjiFBYWc/m5v2HU198wloWUkUkB
   B/MYXdjAnuFZZcBnQH/UVdRyxGtF83oqu48cWEmwE5uItBL5+VO57baX2GP9UsaymndYyT/px8Uc
   wet8j2BtwVyqtgrmAJ8SfCyoq6i10OwjkTYqP38qt9zyNFs2reME0shjAfuylnvI4R5GsoxOBIkg
   0imwgsq/DyOtgtuB1UAWZun07dtRLYNmrNHdR2bWHhgD7EPQqvgIeMrd18Uz0IZSUhBpmEjJicWL
   V1NRUQpsZic6cj7fcinzWcYOFLAHz3AMm1kK7BG+MzYRRBjBtNJ2QDtycpQIWopGdR+Z2UDgBeBN
   4D2C34RhwDXhyuYfuPuf4hiviMRZYWEx559/EytXVhAUIujCwawljy84iZU8x978iKP5gIMIBo4h
   WENQfbFZJBGUAhlkZHRk4MCdtLagFaprTKEAuNTdZ8QeNLNjgH8TtBqUFESaofz8qfz2t4+yZUsW
   0J50MjidUvJ4g10oZyr9uIIjWMVewCcEYwWRZDAHeJzKcYLlgJGW1oFBg3orEbRytXYfmdk8dx9Q
   y2uLgH3dfUMig6uNuo9EajdmzK944olZQCZ7UMYlfM4FfMaHdKeAXvyDI6hgXnh2f2B3oIggMexC
   0GN8O0HZswzS0zdxzTU/0JTSVqCxs4/MzDLcfVO1gxnAlmQlBBHZWuWeBSuATuRSSh5fcBQLeIRe
   DCWXTzmIYBppBUFP8ItUzh7KJVinOhu4GcgiOzudK644SsmgjamrpXAdMATIc/fF4bE+wB3AO+5+
   Q1MFWUNsaimIUHXMoCMd+SkLGcvnOEYBQ3iUTDawA7CKyn0L2hG0ELoDfwt/UgegvQaNW7l4zD7K
   I9h1LSs8tAH4vbtPjluU20FJQSQYN5g06QX6+ybGsoyfsJR/swsF9OU1dgbygD8TJIFsYA3B9NFS
   gu3WO5KS0p4+fbKUCNqIuK1oNrNOAO6+Nk6xNYqSgrRV0amlC7/geJw8ljKY77iX/tzN0SxlOcGu
   ZS8CvagcHwjWEkC6WgNtmEpni7QiY8b8iulP/JfzWcGlLGE5WRSwN39jVzaTQtAltIJgN9tI11Am
   kInZBn796xM1PtDGKSmItHCRVccDNy0nj/X8kK94nu5MYSDvk00wSJxN5ZhBpHz17kAqUE6XLot4
   9NGL1DKQuJXOFpEmEjuLKJ0NnEoF/+YLdqeUqQykP0fyLRnAJqAjQTJYS7AwbQ2wKPxJK0lLy2TQ
   oF2YNEkJQeqvPpvsZAFXAD3d/UIz6w8McPeXmiLAWmJSS0FaldhZRN0p4xK+4AKW8hFdKaAfL9GX
   CjZTWXIiO+ZxZAA5C0ijW7d23HdfnhKBbCVeLYUHgPeBw8PnXxJ0ViYtKYi0FpXJoJxhbKCApQzn
   Kx6lJ7kcxzw6h2eWUNlNVEJlt5ERzCTqQUrKKq6//hiNG0ij1Ccp5Lj76WZ2BoC7bzCrM9GIyDYE
   q45fJYssfsI35PEFKRgF7MW5DGU9aQQf/iXhO4YRrDruSpAIlhEkg6CMRU5OOrff/ku1DqTRUupx
   TqmZZUaemFkOwUTnhDGzUWY218w+MzPt3SCtRn7+VFJSDuS9J97kT3zN57zJCNYyniHsw0jupBfr
   o3WIhhGMF6wF/kew6ngFwQK0XenatTMvvXQV7s8wf/5DSggSF/UZUxgBTAAGEuyvdwRwjru/mpCA
   zFKBecAxBH8OvQuc6e6fxJyjMQVpEWIHjlPYyHFAHkvYn/XcT3/uYgBLaEfVWUSR1kEvYD+C3tpS
   oDPQgczMLVx11Uh1E0mDxXPx2k4EJS8A/uvu38QhvtqudRgw0d1Hhc+vBnD3m2POUVKQZqsyESwC
   dmAHnPNYxmV8wVdkUkB/nqY/pWyh5imla6jcx6ALSgQSL43dT+Egqu6usTz83tPMerr7B3GIsSbd
   gSUxz5cC30vQtUTiojIRLCH4i34LB5DFWBbyI5byArszmqN4l52BbwkqkEaSQfUppVnArppFJElR
   10DzbWy95VKso+IcS0S9mgD5+fnRx7m5ueTm5iYoHJHa5edP5cYbH6C8vDOQRho7cSqrGMun9KCU
   O9mLPTmUb8igcuB4IDCL4Fe9ejLoomQgcVNUVERRUVGD3tPsVjSb2RAgP6b76Bqgwt1viTlH3UeS
   NJFVxps2fUfw1342u2NczEIuZBEf04UC+vISu1NObEs9Uq66F0Fr4g0qZxClk5OTrZpEklBxGVMI
   Zx5dBhxJ8KfN68Cd1fdZiBcza0cw0Hw0wZqId9BAszQDwTTSGcDOBBP3jO+zhjy+5BiW8zh9mEof
   PqETsJGgBlFtA8dBTaJ27UqYMOF4jRVIk4hXUniaoNPzUYIJ0mOAzu5+WrwCreGaPyCo+ZsK3O/u
   v6v2upKCNInKVsEqYDcglQ5kcBafk8dC0imjgP48TC/WkUblOMEqYCeCFoEGjqV5iFdSmOPuA7d1
   rCkpKUgiVe0eirQKoB/GZSzmZyykONy34N90BLpR2SLIpnIG0RKCzWuyUfkJaQ7iVebiAzM7zN3f
   Cn/oEIKyFyKtzrBh51JcvBzIADqQQgdGsYI85nMQa7if3hzIMXxBFsEYwVMECaCmlcY9tHeBtDj1
   aSnMBfYk+LPHgZ4Eff5lgLv7fokOsoaY1FKQuKi6t/F6YA+gI10p5VyWcRlzWU17JrMHf6UfpQwn
   GCyGyjGCx4jsZKYWgTRn8eo+6l3X65H9m5uSkoI0VmxV0qDBXAJkMZhU8pjLj/mcl9iNAvrxDjsA
   fQimkfZGq4ylpYrniuauQA9iupsSuHitPvEoKch2qZoM2gNZpFHOj1hIHsvoxTruJIf76MvXjKRq
   qyAyjTRIBCkpG7n++hOUCKTFiFdL4QbgHGAhUBE57u6JWry2TUoK0hCVA8frCT7QM4AsdqOCi/iM
   i/iUuWRTQD9e4AzK+Uf4TtUektYlXknhU2Bfd98cz+AaQ0lB6qNytfEOBDOIUoEOHMFa8pjHSJbx
   BD2ZQj/mkAb8gKAaqRKBtE7xmn30McHUipVxiUokgarWIOoK7AJ0pANljOFL8niTDMqYQg4X833W
   sp5gPGEdMBM4ACgmmE/xNfvsk87s2X9L1u2INLn6tBQOAf4OzKZyHwV395MSHFtdMamlIFFbF6Pr
   SvBB34G+bOEy5nE2C/gPO1JAP/5FP5zVbF2VdDPBgrOOpKRsZPTofXn88VtquapIyxOv7qNPgDsJ
   kkJkTMHd/bW4RLkdlBQEappBlApkYfRjJG+Tx1wO5Rv+Qm/uIofF7MbWJaq1t7G0HfFKCu+6+yFx
   jayRlBTatppmEAF0wTiXBVzGAr4DCujHk4xhE6+E74xdbaxkIG1PvJLCHwm6jV4gZhtOTUmVplR1
   kVkHgoVi6UAG+7GasczlNBZTyG4UMJy3ORh4ia1nEAWrjc3a07evVhtL2xKvpFBEDXscaEqqNIXK
   6aTtCMYKNgPtaceenMLr5DGHvqznLnK4l/34ihKC7qH+BPs1VVYkTU8v5ZprfqAZRNJmxW3xWnOj
   pNB6RVoEixevpqJiDcFK4nKCbqL+dGMeFzKfi5nPZ2RSwCD+TjZllKLuIZG6xXNF8wkE20VlRI65
   +6RGR7idlBRal/z8qdx220usX/8t0ImgRVBG0N1zAPAFh/MNeSxhFPN4kp5M4To+5gW2nkEUKUbX
   XsXoRKqJV/fR3QTt7+HAvcBpwNvufn68Am0oJYWWr+o00h4EPZRO8KHeH/iMTNI4k7Xk8TpZlDGF
   w3iIXfmODQTrD8YAt1PZKkhXIhCpQ7ySwkfuPsjM/ufu+5lZR2Caux8Zz2AbQkmhZQtWGr9IeXk2
   wXqCwcBnwN7AYvrQmUv5F+fyGW/RkwJymEEmHt3EBoKWwWaCQef1nHmm1hSIbEu8VjRH9hLcaGbd
   gW8JdhURqZfCwmKuv/5hPv10ORs2fEPw69OJoEWwFGiHkcEI5pHHa3yP73iQHhzKBBYxm6qJYFH4
   uANpaZkMGrQTkyaNVctAJE7qkxReDKuk/h74gKCNf29Co5JWo7CwmJ/85F7WrEknGCzenaBlsBho
   R2e+4xzeZCxvsZ5uTKYbpzMk/EvkSyAXKCQYX8ggNXUHBg/elUmTRisRiCRAg2YfmVl7IMPdv0tc
   SPWKQ91HLUBhYTGnnXYLJSUHhEciXURl7Mu7jOVrRjOXl9mJAvrwFnsBZ1E5ThAUpDPLomPHCq64
   4mhNJxVphEaNKZjZocASd18ePj8b+DHBn3j57r4qvuHWn5JC8xPpIpozZyGlpU7Q378LQWN0HwDa
   sZCT2UAeC+nPXO4mh3s5mhUMIGgNfEfQrZRFdraSgEi8NTYpzASOdvdVZjYU+CuQRzBHcC93PzXe
   AdeXkkLzUTmLKPLfIzJFdD3BmMF8dqUbFzKTi3mXBezJFPrwHOMoYwqRmUNm6VphLJJgjR1oTolp
   DYwG7nb3Z4BnzGxWvIKUlqdqq6ALkEZkGmnwHSCVISwlj/c4jiU8RR+O4zA+oj/BVNIZwCBgFmee
   uadmDok0E3UlhVQzS3P3LcAxwEX1fJ+0UlVbBV2BnQmSQOTXIZMMnDOYTR7v0JlSpjCGPNawhgxg
   OUGx3ZtjxgmOVReRSDNS14f7E8BrZvYNsBF4HcDM+hPMDUwIM8sHLgC+Dg9d4+7TEnU9qZ/8/Kn8
   7ndFbN4cqSsElb8+ZfRmDZfyX85lCW+zOxO4gul8gtMZOJ2gZbAH6emzueaa4UoEIs1UnbOPzOww
   gknl0919Q3hsT6BjoqqkmtlEYJ27/7GOczSmkECR7qHFi9dTWlpCWdlqNm/eAdiX2L8jjC0cwyLy
   +JzDeZ8H6cGd9GchuxH82owEHgY2AGXk5HTQmIFIEjV68Zq7v1XDsU8bG1g91Bm0JEbV7qG9gZ8B
   DxGUvNo7PKuMTmzibD5iLJ9QwgYKGMIZvEAJTwELw69FwCwyMrIYOHAnrSsQaSGaXZXUsKVwLsH8
   xPeAX7j7mmrnqKUQR5GWwezZ69iyBYLuoRuB68Iz2gFlDORrxrKSM3iZ6XSngCP5D+cDj6DWgEjz
   F68yF3FnZjOouVTGBIKtPyMVWG8AbgO2Kr6Xn58ffZybm0tubm68w2wT8vOncuut/6OkpBvBf5LY
   X4l2pFLBycwmjwUMYB53czD78CQrKCQYOP4DGRkd1BoQaYaKioooKipq0HuaXUshlpn1Bl5090HV
   jqulEAfBiuMplJT8FcgPj5YBsDPjuZDTuYQPWExvCsjhOfLYwhMErYIt5ORkqVUg0oLUp6WQ0lTB
   1JeZ7Rbz9BTgo2TF0loVFhYzcuR1jB79J0pKKscKYAuHsgcP8yTz6EVvOnIiwxnKiTzF/7GFIoIZ
   RFuYOHE48+c/pIQg0so0x/UGt5jZ/gSF9xYBFyc5nhatsLCYO+6YzrJlX7NixRpSUjbw7be7UlFx
   NvAFUEYGJYymnDzupSvlTGUMl7OW1ZQBq0lPf5v09HmkpWXRp09HJk0ar2Qg0ko16+6j2qj7qGZ1
   J4BXCKZV43C+AAAPM0lEQVSITiGoWHIdPfmOS1nFeTzPewyjgFymMQ/nK7KyshkwoKvGCURakWY7
   0CzxEZsEPv98CaWle7Bly08IEsBZBAngfoJZRDeF3/fiaP5JHn/j+yznIfbiCO5lPp8AGzFbzcRf
   a3N7kbZKSaGFiU0ECxcaJSVjCJLAzlROI40kgMh4QTuyWcvP+IA83mEzz1FAX87ifTbyPsFq41Sg
   nP3330EJQaQNa3YDzVK7wsJiLr/8FaZPv5GPP96ZkpK7gOkESSCS32O/l7E3cyjgJRbTm++zmIuY
   yGBO4F6uZiO/BYYSzPzNJyennBtu+FmT35eINB9qKTRjkVZBaWk71q5dyvz5q1i37tnw1epJoCz6
   PZUyTmQ2eSxiIH/mHk5lEIfzJacStCpGEbQOviEl5UR69OjGXnt1Y9y4URo/EGnjlBSama27h+4C
   igk+zPeIObOs2vcR7MTPuYCNXMqOLAk3vH+GW8OppF+RkvJ7dtyxPRUVC9htt250774L48adpUQg
   IlGafdSMRLqHFiyIjAncGL4SeRx7LJIoRnIID5BHOSfyHM/Qj3vadWJu5k7ssEMG69eXhwkgm3Hj
   jlUCEGnDNPuoBYjtIpo9+xO+/fav4StVy00ERhBUArmJ9hzK6UxjvJ3CLqml3J/ejxN7jqFjr12Y
   qA9/EdlOSgpJVLVlAJWlJqCyWyj28VB6sJJLOIIL+JCP0jrh1/ySnr/+Jb9JTeU3TRK1iLRmmn2U
   RHfcMT0mIUDVRBBpFQAcy3DG8Aw/YiaX0IFDGNPjJ2x67q8c8purITW16YIWkVZNLYUmFttdNGvW
   kmqvVnYPwVA6spEL0g4kL+VzKsx5auf9uLfXxXhWKv837kfqIhKRuFNSaEJbdxddV+2M4EP+e12O
   46oO33DsVx+x/tDvsdtNz8LQoUww7T0kIomlpNCEtu4uqmwZpFDOCbzELzN/zsGsJeO8y+DiZ8ne
   Y49afpqISPwpKSRQbFdR+/ZlfPnl+mpnDGVH1jAuczAXV3zOqsxsSs8/h4ybroP27ZMRsoi0cUoK
   CbJ1VxFkZo6OPj6I98ijgJP5O+/s0Ituf/8X3Q46KBmhiohEafZRgmzdVQTlJRdybtpw3mIIf+NU
   5jCQEb3PpuzuO0AJQUSaAbUUEqS0tPKfdg+WcAl3cQH3sSA9jWn7HkFx9gjSM9eRr1lEItKMKCnE
   QfWxg/HjR9A+fQu5vEoeBRzFqzzCTxlKMX2OfJRp025IdsgiIjVSUmik6mMHHVnH9z44jUdsFt+l
   3cMft0ziHB5kPdnk5FzLuHGjkhyxiEjtVBCvHmpqCUS6fEaOvI7p029kT+Yxlin8hEd5laN4++AO
   DJt4AZML/smmTalkZJSrIJ2IJJUK4sVBTbOIFiwIyk8cP+oIDlq2gF8wgsHM4l4uZDCzWEoPhmXl
   c+sJwzj+hGHJCl1EpMGUFLahpllEqxdcwfIrToexC/jpmi3cxC08zWlspnJtQUZGeVOHKiLSaG0q
   KdTVDVSb2FlEB/I+Y5nCKTzH29/1ghefZuFXJfz38lfYvKAyIWjsQERaqjaTFOrsBqojMXRM28QY
   HiOPAnbnS+7kUvbkUw7a/w5GHXIIx4fnTZ58fczYgba1FJEWyt2b/As4DfgYKAcOrPbaNcBnwFxg
   RC3v94YaMWKCg2/1NXLkdTW/YckS9+uu85IuO/jrmb39JJ73FMoc3HNyrvGXXnqtwTGIiCRT+NlZ
   5+dzsloKHwGnAHfHHjSzgcBoYCDQHfinme3p7hWNvWBsN1CsTZti9iJwh+JiKCiAf/0LzjqLjDff
   4LuFX1M6eQbf3zRTLQERadWSkhTcfS4E06OqORl4wt23AIvNbD5wKPDfxl6zffuyGo9nZJTD+vXw
   6KNBMigvh7w8uP9+6NQJgOP33ltJQETahOY2prA7VRPAUoIWQ6ONHz+CBQsmVBlTOLrHxdzZfhX0
   6gXDhsHtt8Pw4aB9C0SkjUpYUjCzGUC3Gl661t1fbMCPqnGVWn5+fvRxbm4uubm5df6QyF/6U+6Y
   wOAvF/Hj5W8zaO0q2u9zKUy+DXr2bEBIIiLNX1FREUVFRQ16T1JXNJvZq8Av3P2D8PnVAO5+c/h8
   GjDR3d+u9j7f7riLi+HKK4MuotNPh4yMxtyCiEiLUZ8Vzc0hKVzp7u+HzwcCjxOMI3QH/gn0q54B
   GpUU3NU9JCJtUn2SQlL2UzCzU8xsCTAEKDSzlwHcfQ7wFDAHeBm4bPs//Wu9eFx/nIhIa6KCeCIi
   bUSzbSmIiEjzpKQgIiJRSgoiIhKlpCAiIlFKCiIiEqWkICIiUUoKIiISpaQgIiJRSgoiIhKlpCAi
   IlFKCiIiEqWkICIiUUoKIiISpaQgIiJRSgoiIhKlpCAiIlFKCiIiEqWkICIiUUoKIiISpaQgIiJR
   SgoiIhKlpCAiIlFJSQpmdpqZfWxm5WZ2YMzx3mZWYmYzw6+pyYhPRKStSlZL4SPgFKC4htfmu/sB
   4ddlTRxXs1BUVJTsEBJK99eyteb7a833Vl9JSQruPtfdP03GtVuC1v6Lqftr2Vrz/bXme6uv5jim
   0CfsOioysyOTHYyISFvSLlE/2MxmAN1qeOlad3+xlrd9CfRw99XhWMPzZraPu69LVJwiIlLJ3D15
   Fzd7FfiFu3/QkNfNLHlBi4i0YO5udb2esJZCA0QDNLOdgNXuXm5mfYH+wMLqb9jWTYmIyPZJ1pTU
   U8xsCTAEKDSzl8OXhgGzzGwm8DRwsbuvSUaMIiJtUVK7j0REpHlpjrOP6sXMbjCzWWb2oZn9y8x6
   JDumeDKz35vZJ+E9PmtmnZMdUzzVtoCxJTOzUWY218w+M7NfJTueeDOzv5jZSjP7KNmxxJuZ9TCz
   V8PfydlmNj7ZMcWTmWWY2dvh5+UcM/tdree21JaCmWVHZiWZ2ThgsLtfkOSw4sbMjgX+5e4VZnYz
   gLtfneSw4sbM9gIqgLupY7JBS2FmqcA84BhgGfAucKa7f5LUwOLIzL4PrAcedvdByY4nnsysG9DN
   3T80s47A+8APW9l/vw7uvtHM2gFvAFe6+xvVz2uxLYVq01Q7At8kK5ZEcPcZ7l4RPn0b2COZ8cRb
   K1zAeCjBavzF7r4FeBI4OckxxZW7vw6sTnYcieDuK9z9w/DxeuATYPfkRhVf7r4xfJgOpAKrajqv
   xSYFADO7ycy+AM4Gbk52PAl0HvCPZAchdeoOLIl5vjQ8Ji2MmfUGDiD4Y6zVMLMUM/sQWAm86u5z
   ajqvOUxJrdW2FsC5+wRggpldDfwJOLdJA2yk+izwM7MJwGZ3f7xJg4uD7VzA2FK1zH5YqSLsOvob
   cHnYYmg1wp6H/cPxyVfMLNfdi6qf16yTgrsfW89TH6cF/iW9rfszs3OA44CjmySgOGvAf7/WYBkQ
   O9mhB0FrQVoIM0sDngEedffnkx1Porj7d2ZWCBwMFFV/vcV2H5lZ/5inJwMzkxVLIpjZKOCXwMnu
   vinZ8SRYa1iM+B7QPyz/ng6MBl5IckxST2ZmwP3AHHf/c7LjiTcz28nMuoSPM4FjqeUzsyXPPvob
   MAAoBxYAl7r7V8mNKn7M7DOCAaHIYNBbramUuJmdAtwB7AR8B8x09x8kN6rGMbMfAH8mGMS7391r
   nfbXEpnZEwQLTHcEvgJ+7e4PJDeq+AiLbxYD/6OyK/Aad5+WvKjix8wGAQ8RNARSgEfc/fc1nttS
   k4KIiMRfi+0+EhGR+FNSEBGRKCUFERGJUlIQEZEoJQUREYlSUhARkSglBUkqM9vRzGaGX8vNbGn4
   eLWZfdzEsZxsZnvHPP+NmTV4NXm4gK3G8tJmto+Z/Tsssf2pmV3XmJjriKGmexkePi4ys4MScV1p
   +ZQUJKnc/Vt3P8DdDwDuAv4YPt6foLR2XIUlrmtzCjAwJraJ7v6vOF47E/g78Ft33wsYDBxuZolY
   lFjTvfw78hTVapJaKClIc2Mx31PN7J5w05NXzCwDwMxyzOxlM3vPzIrNbEB4vHf4V/gsM/tnZOMl
   M3vQzO4ys/8Ct9T0fjM7HDgR+L2ZfWBmfcP3/Tj8GYeY2X/CTUreNrOO4fWKzez98OuwbdzbGOAN
   d/8ngLuXAHnAr8Jr5JvZL6L/EMF99wwfPxfGO9vMLow5Z72Z3RjG9ZaZ7bKte6nyj202wszeDON/
   ysyywuM3W7DhzCwzq3Hlq7ROSgrSnPUHCtx9X2ANEPlQuwcY5+4HE9SHmhoenww84O6DgccIymhE
   7A4c5u5X1vR+d3+ToFbRle5+oLsvJPyLOqxl9CQw3t33JyhQWEJQgvhYdz8IOKPa9WoykGDzlqjw
   Otlhdc7qf73HPj8vjPcQYLyZdQ2PdyAogbI/QZmGC+u6l9gfbmY7AROAo8N7eB+4wsx2INhgZp/w
   3/KGbdyXtCLNukqqtHmL3P1/4eP3gd7hX7KHA08HNcyAoEYUwBDgh+HjR4Fbw8cOPO3uHn74HlbL
   +2Hr4nxGUGNrubu/D9FNWAiTRYGZDSaowbVnPe6ptuJ/adt43+VmFrm3HgQJ8x2CsuqF4fH3CQqd
   betakdeGECSqN8N/i3TgTYJaVJvM7H7gpfBL2gglBWnOSmMelwMZBK3b1eG4Q01q+yCM7DqVAqyp
   4/019bXX1v/+fwTJ4qfhWMW2qtnOAYZWCdasL7DR3VebWRlVW++R7rJcgtbJEHffZGavRl4DtsSc
   X0HV/6frM24ww93HVD9oZoeG1zyVoIurRZZvl4ZT95G0JBZuw7rIzE6FoOSxme0Xvv4mQTcOwFkE
   3SlVuPvaOt6/DuhU/S0Eey/vZmYHh+/JDpNAJ2BFeN7PCKqj1uUx4MjIjKZw4Pl24Lbw9cXAgeFr
   BwJ9wuOdCBLhJgv2th6yjevUdi/V7+u/wBFmlhNeM8vM+oetsS7u/jJwBcGAuLQRSgrS3Hgtj2Of
   nwWcb8HWgrOBk8Lj44BzzWxWeM7ltfys2t7/JPDLcNC1b/SNwZ7Lo4HJ4XteAdoTjGWcHR4bQLCp
   fW2xE+6LcRLBboHzgK+Bz9z9T+EpzwA7mNlsYCxBMgKYBrQzsznA74C3arlO7LhBjfdSLZ5vgHOA
   J8J/szfD+8gGXgyPvU7QIpI2QqWzRZLEzE4Gfg0c7+4rtnW+SFNQUhARkSh1H4mISJSSgoiIRCkp
   iIhIlJKCiIhEKSmIiEiUkoKIiEQpKYiISNT/AyqzUd5HUH+2AAAAAElFTkSuQmCC
   "
   >
   </div>
   
   </div>
   
   <div class="output_area"><div class="prompt"></div>
   
   
   <div class="output_png output_subarea ">
   <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAYUAAAEPCAYAAACtCNj2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz
   AAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VPX1x/H3SUhICGFxRZE1IIoi7sWlEFGButZWRbGt
   +0rA/qy1KlpS1FZtbasE3Ou+VOtSNRWhrTFarStSEUHZFBBwAWQLgSTn98e9M5mEJCRkJpPl83qe
   PJm5cyf3XIxz8t3O19wdERERgJRkByAiIs2HkoKIiEQpKYiISJSSgoiIRCkpiIhIlJKCiIhEJS0p
   mFkPM3vVzD42s9lmNj48voOZzTCzT81supl1SVaMIiJtjSVrnYKZdQO6ufuHZtYReB/4IXAu8I27
   32pmvwK6uvvVSQlSRKSNSVpLwd1XuPuH4eP1wCdAd+Ak4KHwtIcIEoWIiDSBZjGmYGa9gQOAt4Fd
   3X1l+NJKYNckhSUi0uYkPSmEXUfPAJe7+7rY1zzo21IdDhGRJtIumRc3szSChPCIuz8fHl5pZt3c
   fYWZ7QZ8VcP7lChERLaDu1tdrydz9pEB9wNz3P3PMS+9AJwdPj4beL76ewHcvdV+TZw4Mekx6P50
   f23x/lrzvbnX72/pZLYUjgB+AvzPzGaGx64BbgaeMrPzgcXA6ckJT0Sk7UlaUnD3N6i9pXJMU8Yi
   IiKBpA80y9Zyc3OTHUJC6f5attZ8f6353uoraYvXGsPMvCXGLSKSTGaGN9eBZhERaX6UFEREJEpJ
   QUREopQUREQkSklBRESilBRERCRKSUFERKKUFEREJEpJQUREopQUREQkSklBRESilBRERCRKSUFE
   RKKUFEREJEpJQUREopQUREQkSklBRESilBRERCRKSUFERKKUFEREJKpdsgMQEZG6FRYWc8cd0ykt
   bUf79mWMHz+C448fmpBrJTUpmNlfgOOBr9x9UHgsH7gA+Do87Rp3n5acCEVEkquwsJjLL3+FBQtu
   ih5bsGACQEISQ7K7jx4ARlU75sAf3f2A8EsJQUTarDvumF4lIQAsWHATkyfPSMj1kpoU3P11YHUN
   L1lTxyIi0hyVltbcobNpU2pCrpfslkJtxpnZLDO738y6JDsYEZFkad++rMbjGRnlCblecxxovhOY
   FD6+AbgNOL/6Sfn5+dHHubm55ObmNkFoIiJNa/z4ESxYMKFKF1JOzrWMG1e9531rRUVFFBUVNeh6
   5u4NjTGuzKw38GJkoLk+r5mZJztuEZGmUlhYzOTJM9i0KZWMjHLGjTt2uwaZzQx3r7N7vtklBTPb
   zd2Xh4//DzjE3cdUe4+SgohIA9UnKSR7SuoTwDBgJzNbAkwEcs1sf4JZSIuAi5MYoohIm5L0lsL2
   UEtBRKTh6tNSaK6zj0REJAmUFEREJEpJQUREopQUREQkSklBRESilBRERCSqOZa5EBFptZpyb4Tt
   oaQgItJEmnpvhO2h7iMRkSbS1HsjbA8lBRGRJtLUeyNsDyUFEZEm0tR7I2wPJQURkSYyfvwIcnIm
   VDkW7I1wbJIi2poK4omINKF47Y2wPVrEfgrbQ0lBRKThVCVVREQaRElBRESilBRERCRKSUFERKKU
   FEREJEpJQUREolQQT0Qkjpp7FdRtUVIQEYmTllAFdVvUfSQiEictoQrqtiS1pWBmfwGOB75y90Hh
   sR2AvwK9gMXA6e6+JmlBiojUoKZuopZQBXVbkt199AAwGXg45tjVwAx3v9XMfhU+vzoZwYmI1KS2
   bqJOnVbXeH5zqoK6LUntPnL314Hq/4onAQ+Fjx8CftikQYmIbENt3USwudlXQd2WZLcUarKru68M
   H68Edk1mMCIi1dXWTdSp0x7ccMNwJk++PqYK6qgWM8gMzTMpRLm7m1mN5VDz8/Ojj3Nzc8nNzW2i
   qESkrak+frB27aoaz8vIKOf444cmLgmUl0Nq/ccnioqKKCoqatAlkl4628x6Ay/GDDTPBXLdfYWZ
   7Qa86u57VXuPSmeLSJOoafygW7fzgc6sWPHH6LGcnGu5/fYEtQo++wymToVnn4W5cyEzc7t+TH1K
   ZzfHlsILwNnALeH355Mbjoi0ZTWNH6xYcT8HHHABgwcnsJuovBxefhkKCuCDD+CCC+D117c7IdRX
   sqekPgEMA3YysyXAr4GbgafM7HzCKanJi1BE2rq6xg+mTcuP/wVXrYK//CVoGey4I4wbB88/DxkZ
   8b9WDZKaFNz9zFpeOqZJAxERqUX79mU1Ho/7NNOZM2HKFHjmGTjxRHjySTj00Pheox60ollEpA7j
   x49I3DTTzZvhiSfgyCPhpJOgb1+YNw8efjgpCQHqMdBsZv2Ape6+ycyOAgYBDydzlbEGmkWkKRUW
   FjN58oyY8YNjGzd+8OWXcPfdcM89MHAg5OUFrYN2ie28qc9Ac32SwizgIKA38A/g78A+7n5cnOJs
   MCUFEUmkhFQ6dYc33ggGjmfMgDPPhMsug332iU/Q9RCv2UcV7l5mZj8CJrv7ZDObGZ8QRUSSLzYJ
   rF27lOXLO1WZbtqoSqcbNsDjjwfJYNOmoFVwzz3QuXO8wo+r+iSFzWY2BvgZcGJ4LC1xIYmINJ2t
   1yFcB9xY5Zyg0un1DUsK8+cHM4gefjgYM/jDH+DooyGleQ/l1ie684AhwE3uvsjM+gCPJDYsEZGm
   sfU6hEZUOq2ogH/8A447Dg47DNLT4b33gimlxx7b7BMC1KOl4O4fm9nVQM/w+SKChWUiIi3e1usQ
   tmMK6urV8MADQcugc+dgbcEzzyR8oVkibDNtmdlJwExgWvj8ADN7IdGBiYg0ha3XIYwA6jkFddYs
   uOiiYCrpBx/Ao48GLYNzzmmRCQHqN/voA2A4QQ2iA8Jjs9193yaIr7aYNPtIRBolMri8bNnXLFxo
   lJTcFX2tW7fz2H33TLKzd956CuqWLUENooICWLQILr00KEGxa/Mv6Byv2Udb3H2NWZWfU9GoyERE
   kqD2RFBMZuZocnJ2o3v3bMaNO2frQeXly4NZQ3ffDQMGwM9/DiefnPC1BU2tPnfzsZmdBbQzs/7A
   eODNxIYlIhJfVWcZVZ9hNJSSkqF0734906bdUHnYHf7zn6D8xLRpcMYZMH067Ju0jpKEq89Q+Dhg
   H6AUeAJYC/w8kUGJiMRb1VlG25hhtHEj3HcfHHAAnHdeMJNo8WK4885WnRCgfrOPNgDXhl8iIi1S
   1VlGNc8w6lPxNVx5JTz4IBx+ONx6KxxzTIuYShovtSYFM3uxjve5u5+UgHhEROKi+irl+fNjd0uL
   zDC6CaOCEUznqg55HDHrGxhyEbz7LvTpk6TIk6uulsJtTRaFiEgcVR0/KAZeISjKECQCGEpn1nFR
   2v6MtcVsaZ/JhnN/RvvfTYQOHZIZetLVmhTcvagJ4xARiZuq4wfTCRJBYBAXMZb3GW0fseHwYXT/
   3Z0wZAhYnTM124y6uo+edvfTzGw2UH1RgLv7fokNTUSk4QoLi3nnnSUxR9rRji38kOfJo4B+zOcu
   LuG8IUfxbNEfkhZnc1VX99Hl4ffjgeopVCvHRKRZKSws5vrrH+aTT9LYtKkHALuyggt5lUu4l/n0
   o4A8nueHlJHGyE7XJzni5qnWIXV3/zJ8eJm7L479Ai5rkuhEROohMoYwc2Y3Nm2aymHszmPswyfs
   TXe6MIrjyOU1/sZplJEWv53TWqH6lLmYGSlvEXPsI3cflNDI6o5JZS5EJGrkyOsonj6BMzmDPJaQ
   zTqmMIoHyeA7ssjKepcBA3rXXLaiDWlUmQszu5SgRZBjZh/FvJQN/Cc+IYqINM6/7/8rx7/2PI9y
   N2/ThWuZzHRG4DEdIUceWW2lstSqrjGFx4GXgZuBX1E5rrDO3b9NdGAiIrWqqODtG//Iptsms9/a
   r/iAwQzh7yxkGcH001HRU4OuolG1/iipapvdRwBmlgrsSkwScfcvEhjXtuJR95FIGxJZiPbdF8sY
   uvBtLtzyORs8g8kcxhM8RQnvESSDyLqEGWRkfM7AgdlMmjS6TXYV1SQuVVLNbBwwEfgKiN1lIqFj
   Cma2mKDOUjlBpdZDE3k9EWk+IjOJPv10OSUlpQys6MRllHMG05hGP87hVN7kQeA3QAcg8qF/PZBK
   167zeeSRS5UMtkN9qqT+HBiQhC4jB3LdfdU2zxSRVqOwsJgLLniIb1ZUcDLtGcvHDGAzd9OPgVzC
   Cm4H8gl6tGNrGA0lkhwOPbSB+ylLVH2qPH1B8Bd7MmiJoUgbUlhYzJU/vY3zVyxgEc8wnm+4k0Po
   xedMYjgr6BqeGUkGDdglTeqlPi2FRcCrZlYIbA6Pubv/MXFhBdcA/mlm5cDd7n5vgq8nIsnizn0X
   XkXWg8/zZvlSnmYgxzOG/3EXwd4HaVRtFVQWtAtcrzGEOKlPUvgi/EoPv4ymWdF8hLsvN7OdgRlm
   NtfdX4+8mJ+fHz0xNzeX3NzcJghJROLp5Wdn8PGv72DEZ//hqM2lTOEGLmMZa4jd3ziSAEYCD1E9
   GaSlLWLQoM5MmnSBkkE1RUVFFBUVNeg99Zp9lGxmNhFY7+63hc81+0ikhYkdPN5l4xou9sWcyyre
   YzAFpDONo3B+QzB76KHwXd2InVGUmvoh6embSU3NJi0tiz59Oqpl0ADxmn20C3AVMBCi6dvdfXjj
   Q6z1mh2AVHdfZ2ZZBH8q/CZR1xORxCksLGb8+NtZuLCCYyjhMT7hCFbyEL04nNdYwIMEH0WR7qHI
   B/zDwIfAiWRlZTNgQFcmTfqlEkCC1af76DHgr8AJwMXAOcDXCYwJgjURz1lQyrYd8Ji7T0/wNUUk
   jiLJ4OuFmzmbZYzlC0pJpYCDGcNoNtIB6EdlQogdJwhmEmVmXszTT5+lRNCE6pMUdnT3+8xsvLu/
   BrxmZu8lMih3XwTsn8hriEji5OdP5bnfvsAVW+ZzJsuYQQ4XciJv0JPKSYVlMd9HECw+G0lkrUFK
   yiyuuqpt1ihKpvokhciMoxVmdgLwJUTnhYmIRP3jhX/z4kXXctrKBVxMCXdzMPtyMsujPc+x619j
   B5AjCWEGkEpm5idcddWx5OerIHNTq0+V1BOB14EewGSgE5Dv7i8kPrxaY9JAs0gzMuOxvzNr3LWM
   Xr2Iz8migGE8y95siU5UHEHl4HFENyoTwVeYLaNnz93Ya69ubbaKaaLVZ6C5Rcw+qk5JQST5CguL
   ue+iiZzy5cecyBqeYQ+mMJwP2YOqi8seojIBPAwsB0oASEnpSGZmx3AQWbOIEi0uScHMHqh2yAHc
   /bzGhbf9lBREkqOwsJgbJvyFfT9+m4vLvmRHyplCPx7gaFaTHZ4VGSOongw2AGXk5HTg9tvPVwJI
   grhMSQUKqVyslgmcQjCuICJtQGR9wfqP53Du5m94gSW8Txfy2Y9pfJ8K0sMzY1sHrwBnEySDO4Et
   5ORkcfvtKlLX3DW4+8jMUoD/uPthiQmpXjGopSCSYIWFxYwf92d6L/qKPJYwlOU8zAFM5SDmswuw
   GOhN7V1FwaBxevpsrrlmuAaNm4GEjCmY2V7AS+7erzHBNYaSgkjiFBYWc/m5v2HU198wloWUkUkB
   B/MYXdjAnuFZZcBnQH/UVdRyxGtF83oqu48cWEmwE5uItBL5+VO57baX2GP9UsaymndYyT/px8Uc
   wet8j2BtwVyqtgrmAJ8SfCyoq6i10OwjkTYqP38qt9zyNFs2reME0shjAfuylnvI4R5GsoxOBIkg
   0imwgsq/DyOtgtuB1UAWZun07dtRLYNmrNHdR2bWHhgD7EPQqvgIeMrd18Uz0IZSUhBpmEjJicWL
   V1NRUQpsZic6cj7fcinzWcYOFLAHz3AMm1kK7BG+MzYRRBjBtNJ2QDtycpQIWopGdR+Z2UDgBeBN
   4D2C34RhwDXhyuYfuPuf4hiviMRZYWEx559/EytXVhAUIujCwawljy84iZU8x978iKP5gIMIBo4h
   WENQfbFZJBGUAhlkZHRk4MCdtLagFaprTKEAuNTdZ8QeNLNjgH8TtBqUFESaofz8qfz2t4+yZUsW
   0J50MjidUvJ4g10oZyr9uIIjWMVewCcEYwWRZDAHeJzKcYLlgJGW1oFBg3orEbRytXYfmdk8dx9Q
   y2uLgH3dfUMig6uNuo9EajdmzK944olZQCZ7UMYlfM4FfMaHdKeAXvyDI6hgXnh2f2B3oIggMexC
   0GN8O0HZswzS0zdxzTU/0JTSVqCxs4/MzDLcfVO1gxnAlmQlBBHZWuWeBSuATuRSSh5fcBQLeIRe
   DCWXTzmIYBppBUFP8ItUzh7KJVinOhu4GcgiOzudK644SsmgjamrpXAdMATIc/fF4bE+wB3AO+5+
   Q1MFWUNsaimIUHXMoCMd+SkLGcvnOEYBQ3iUTDawA7CKyn0L2hG0ELoDfwt/UgegvQaNW7l4zD7K
   I9h1LSs8tAH4vbtPjluU20FJQSQYN5g06QX6+ybGsoyfsJR/swsF9OU1dgbygD8TJIFsYA3B9NFS
   gu3WO5KS0p4+fbKUCNqIuK1oNrNOAO6+Nk6xNYqSgrRV0amlC7/geJw8ljKY77iX/tzN0SxlOcGu
   ZS8CvagcHwjWEkC6WgNtmEpni7QiY8b8iulP/JfzWcGlLGE5WRSwN39jVzaTQtAltIJgN9tI11Am
   kInZBn796xM1PtDGKSmItHCRVccDNy0nj/X8kK94nu5MYSDvk00wSJxN5ZhBpHz17kAqUE6XLot4
   9NGL1DKQuJXOFpEmEjuLKJ0NnEoF/+YLdqeUqQykP0fyLRnAJqAjQTJYS7AwbQ2wKPxJK0lLy2TQ
   oF2YNEkJQeqvPpvsZAFXAD3d/UIz6w8McPeXmiLAWmJSS0FaldhZRN0p4xK+4AKW8hFdKaAfL9GX
   CjZTWXIiO+ZxZAA5C0ijW7d23HdfnhKBbCVeLYUHgPeBw8PnXxJ0ViYtKYi0FpXJoJxhbKCApQzn
   Kx6lJ7kcxzw6h2eWUNlNVEJlt5ERzCTqQUrKKq6//hiNG0ij1Ccp5Lj76WZ2BoC7bzCrM9GIyDYE
   q45fJYssfsI35PEFKRgF7MW5DGU9aQQf/iXhO4YRrDruSpAIlhEkg6CMRU5OOrff/ku1DqTRUupx
   TqmZZUaemFkOwUTnhDGzUWY218w+MzPt3SCtRn7+VFJSDuS9J97kT3zN57zJCNYyniHsw0jupBfr
   o3WIhhGMF6wF/kew6ngFwQK0XenatTMvvXQV7s8wf/5DSggSF/UZUxgBTAAGEuyvdwRwjru/mpCA
   zFKBecAxBH8OvQuc6e6fxJyjMQVpEWIHjlPYyHFAHkvYn/XcT3/uYgBLaEfVWUSR1kEvYD+C3tpS
   oDPQgczMLVx11Uh1E0mDxXPx2k4EJS8A/uvu38QhvtqudRgw0d1Hhc+vBnD3m2POUVKQZqsyESwC
   dmAHnPNYxmV8wVdkUkB/nqY/pWyh5imla6jcx6ALSgQSL43dT+Egqu6usTz83tPMerr7B3GIsSbd
   gSUxz5cC30vQtUTiojIRLCH4i34LB5DFWBbyI5byArszmqN4l52BbwkqkEaSQfUppVnArppFJElR
   10DzbWy95VKso+IcS0S9mgD5+fnRx7m5ueTm5iYoHJHa5edP5cYbH6C8vDOQRho7cSqrGMun9KCU
   O9mLPTmUb8igcuB4IDCL4Fe9ejLoomQgcVNUVERRUVGD3tPsVjSb2RAgP6b76Bqgwt1viTlH3UeS
   NJFVxps2fUfw1342u2NczEIuZBEf04UC+vISu1NObEs9Uq66F0Fr4g0qZxClk5OTrZpEklBxGVMI
   Zx5dBhxJ8KfN68Cd1fdZiBcza0cw0Hw0wZqId9BAszQDwTTSGcDOBBP3jO+zhjy+5BiW8zh9mEof
   PqETsJGgBlFtA8dBTaJ27UqYMOF4jRVIk4hXUniaoNPzUYIJ0mOAzu5+WrwCreGaPyCo+ZsK3O/u
   v6v2upKCNInKVsEqYDcglQ5kcBafk8dC0imjgP48TC/WkUblOMEqYCeCFoEGjqV5iFdSmOPuA7d1
   rCkpKUgiVe0eirQKoB/GZSzmZyykONy34N90BLpR2SLIpnIG0RKCzWuyUfkJaQ7iVebiAzM7zN3f
   Cn/oEIKyFyKtzrBh51JcvBzIADqQQgdGsYI85nMQa7if3hzIMXxBFsEYwVMECaCmlcY9tHeBtDj1
   aSnMBfYk+LPHgZ4Eff5lgLv7fokOsoaY1FKQuKi6t/F6YA+gI10p5VyWcRlzWU17JrMHf6UfpQwn
   GCyGyjGCx4jsZKYWgTRn8eo+6l3X65H9m5uSkoI0VmxV0qDBXAJkMZhU8pjLj/mcl9iNAvrxDjsA
   fQimkfZGq4ylpYrniuauQA9iupsSuHitPvEoKch2qZoM2gNZpFHOj1hIHsvoxTruJIf76MvXjKRq
   qyAyjTRIBCkpG7n++hOUCKTFiFdL4QbgHGAhUBE57u6JWry2TUoK0hCVA8frCT7QM4AsdqOCi/iM
   i/iUuWRTQD9e4AzK+Uf4TtUektYlXknhU2Bfd98cz+AaQ0lB6qNytfEOBDOIUoEOHMFa8pjHSJbx
   BD2ZQj/mkAb8gKAaqRKBtE7xmn30McHUipVxiUokgarWIOoK7AJ0pANljOFL8niTDMqYQg4X833W
   sp5gPGEdMBM4ACgmmE/xNfvsk87s2X9L1u2INLn6tBQOAf4OzKZyHwV395MSHFtdMamlIFFbF6Pr
   SvBB34G+bOEy5nE2C/gPO1JAP/5FP5zVbF2VdDPBgrOOpKRsZPTofXn88VtquapIyxOv7qNPgDsJ
   kkJkTMHd/bW4RLkdlBQEappBlApkYfRjJG+Tx1wO5Rv+Qm/uIofF7MbWJaq1t7G0HfFKCu+6+yFx
   jayRlBTatppmEAF0wTiXBVzGAr4DCujHk4xhE6+E74xdbaxkIG1PvJLCHwm6jV4gZhtOTUmVplR1
   kVkHgoVi6UAG+7GasczlNBZTyG4UMJy3ORh4ia1nEAWrjc3a07evVhtL2xKvpFBEDXscaEqqNIXK
   6aTtCMYKNgPtaceenMLr5DGHvqznLnK4l/34ihKC7qH+BPs1VVYkTU8v5ZprfqAZRNJmxW3xWnOj
   pNB6RVoEixevpqJiDcFK4nKCbqL+dGMeFzKfi5nPZ2RSwCD+TjZllKLuIZG6xXNF8wkE20VlRI65
   +6RGR7idlBRal/z8qdx220usX/8t0ImgRVBG0N1zAPAFh/MNeSxhFPN4kp5M4To+5gW2nkEUKUbX
   XsXoRKqJV/fR3QTt7+HAvcBpwNvufn68Am0oJYWWr+o00h4EPZRO8KHeH/iMTNI4k7Xk8TpZlDGF
   w3iIXfmODQTrD8YAt1PZKkhXIhCpQ7ySwkfuPsjM/ufu+5lZR2Caux8Zz2AbQkmhZQtWGr9IeXk2
   wXqCwcBnwN7AYvrQmUv5F+fyGW/RkwJymEEmHt3EBoKWwWaCQef1nHmm1hSIbEu8VjRH9hLcaGbd
   gW8JdhURqZfCwmKuv/5hPv10ORs2fEPw69OJoEWwFGiHkcEI5pHHa3yP73iQHhzKBBYxm6qJYFH4
   uANpaZkMGrQTkyaNVctAJE7qkxReDKuk/h74gKCNf29Co5JWo7CwmJ/85F7WrEknGCzenaBlsBho
   R2e+4xzeZCxvsZ5uTKYbpzMk/EvkSyAXKCQYX8ggNXUHBg/elUmTRisRiCRAg2YfmVl7IMPdv0tc
   SPWKQ91HLUBhYTGnnXYLJSUHhEciXURl7Mu7jOVrRjOXl9mJAvrwFnsBZ1E5ThAUpDPLomPHCq64
   4mhNJxVphEaNKZjZocASd18ePj8b+DHBn3j57r4qvuHWn5JC8xPpIpozZyGlpU7Q378LQWN0HwDa
   sZCT2UAeC+nPXO4mh3s5mhUMIGgNfEfQrZRFdraSgEi8NTYpzASOdvdVZjYU+CuQRzBHcC93PzXe
   AdeXkkLzUTmLKPLfIzJFdD3BmMF8dqUbFzKTi3mXBezJFPrwHOMoYwqRmUNm6VphLJJgjR1oTolp
   DYwG7nb3Z4BnzGxWvIKUlqdqq6ALkEZkGmnwHSCVISwlj/c4jiU8RR+O4zA+oj/BVNIZwCBgFmee
   uadmDok0E3UlhVQzS3P3LcAxwEX1fJ+0UlVbBV2BnQmSQOTXIZMMnDOYTR7v0JlSpjCGPNawhgxg
   OUGx3ZtjxgmOVReRSDNS14f7E8BrZvYNsBF4HcDM+hPMDUwIM8sHLgC+Dg9d4+7TEnU9qZ/8/Kn8
   7ndFbN4cqSsElb8+ZfRmDZfyX85lCW+zOxO4gul8gtMZOJ2gZbAH6emzueaa4UoEIs1UnbOPzOww
   gknl0919Q3hsT6BjoqqkmtlEYJ27/7GOczSmkECR7qHFi9dTWlpCWdlqNm/eAdiX2L8jjC0cwyLy
   +JzDeZ8H6cGd9GchuxH82owEHgY2AGXk5HTQmIFIEjV68Zq7v1XDsU8bG1g91Bm0JEbV7qG9gZ8B
   DxGUvNo7PKuMTmzibD5iLJ9QwgYKGMIZvEAJTwELw69FwCwyMrIYOHAnrSsQaSGaXZXUsKVwLsH8
   xPeAX7j7mmrnqKUQR5GWwezZ69iyBYLuoRuB68Iz2gFlDORrxrKSM3iZ6XSngCP5D+cDj6DWgEjz
   F68yF3FnZjOouVTGBIKtPyMVWG8AbgO2Kr6Xn58ffZybm0tubm68w2wT8vOncuut/6OkpBvBf5LY
   X4l2pFLBycwmjwUMYB53czD78CQrKCQYOP4DGRkd1BoQaYaKioooKipq0HuaXUshlpn1Bl5090HV
   jqulEAfBiuMplJT8FcgPj5YBsDPjuZDTuYQPWExvCsjhOfLYwhMErYIt5ORkqVUg0oLUp6WQ0lTB
   1JeZ7Rbz9BTgo2TF0loVFhYzcuR1jB79J0pKKscKYAuHsgcP8yTz6EVvOnIiwxnKiTzF/7GFIoIZ
   RFuYOHE48+c/pIQg0so0x/UGt5jZ/gSF9xYBFyc5nhatsLCYO+6YzrJlX7NixRpSUjbw7be7UlFx
   NvAFUEYGJYymnDzupSvlTGUMl7OW1ZQBq0lPf5v09HmkpWXRp09HJk0ar2Qg0ko16+6j2qj7qGZ1
   J4BXCKZV43C+AAAPM0lEQVSITiGoWHIdPfmOS1nFeTzPewyjgFymMQ/nK7KyshkwoKvGCURakWY7
   0CzxEZsEPv98CaWle7Bly08IEsBZBAngfoJZRDeF3/fiaP5JHn/j+yznIfbiCO5lPp8AGzFbzcRf
   a3N7kbZKSaGFiU0ECxcaJSVjCJLAzlROI40kgMh4QTuyWcvP+IA83mEzz1FAX87ifTbyPsFq41Sg
   nP3330EJQaQNa3YDzVK7wsJiLr/8FaZPv5GPP96ZkpK7gOkESSCS32O/l7E3cyjgJRbTm++zmIuY
   yGBO4F6uZiO/BYYSzPzNJyennBtu+FmT35eINB9qKTRjkVZBaWk71q5dyvz5q1i37tnw1epJoCz6
   PZUyTmQ2eSxiIH/mHk5lEIfzJacStCpGEbQOviEl5UR69OjGXnt1Y9y4URo/EGnjlBSama27h+4C
   igk+zPeIObOs2vcR7MTPuYCNXMqOLAk3vH+GW8OppF+RkvJ7dtyxPRUVC9htt250774L48adpUQg
   IlGafdSMRLqHFiyIjAncGL4SeRx7LJIoRnIID5BHOSfyHM/Qj3vadWJu5k7ssEMG69eXhwkgm3Hj
   jlUCEGnDNPuoBYjtIpo9+xO+/fav4StVy00ERhBUArmJ9hzK6UxjvJ3CLqml3J/ejxN7jqFjr12Y
   qA9/EdlOSgpJVLVlAJWlJqCyWyj28VB6sJJLOIIL+JCP0jrh1/ySnr/+Jb9JTeU3TRK1iLRmmn2U
   RHfcMT0mIUDVRBBpFQAcy3DG8Aw/YiaX0IFDGNPjJ2x67q8c8purITW16YIWkVZNLYUmFttdNGvW
   kmqvVnYPwVA6spEL0g4kL+VzKsx5auf9uLfXxXhWKv837kfqIhKRuFNSaEJbdxddV+2M4EP+e12O
   46oO33DsVx+x/tDvsdtNz8LQoUww7T0kIomlpNCEtu4uqmwZpFDOCbzELzN/zsGsJeO8y+DiZ8ne
   Y49afpqISPwpKSRQbFdR+/ZlfPnl+mpnDGVH1jAuczAXV3zOqsxsSs8/h4ybroP27ZMRsoi0cUoK
   CbJ1VxFkZo6OPj6I98ijgJP5O+/s0Ituf/8X3Q46KBmhiohEafZRgmzdVQTlJRdybtpw3mIIf+NU
   5jCQEb3PpuzuO0AJQUSaAbUUEqS0tPKfdg+WcAl3cQH3sSA9jWn7HkFx9gjSM9eRr1lEItKMKCnE
   QfWxg/HjR9A+fQu5vEoeBRzFqzzCTxlKMX2OfJRp025IdsgiIjVSUmik6mMHHVnH9z44jUdsFt+l
   3cMft0ziHB5kPdnk5FzLuHGjkhyxiEjtVBCvHmpqCUS6fEaOvI7p029kT+Yxlin8hEd5laN4++AO
   DJt4AZML/smmTalkZJSrIJ2IJJUK4sVBTbOIFiwIyk8cP+oIDlq2gF8wgsHM4l4uZDCzWEoPhmXl
   c+sJwzj+hGHJCl1EpMGUFLahpllEqxdcwfIrToexC/jpmi3cxC08zWlspnJtQUZGeVOHKiLSaG0q
   KdTVDVSb2FlEB/I+Y5nCKTzH29/1ghefZuFXJfz38lfYvKAyIWjsQERaqjaTFOrsBqojMXRM28QY
   HiOPAnbnS+7kUvbkUw7a/w5GHXIIx4fnTZ58fczYgba1FJEWyt2b/As4DfgYKAcOrPbaNcBnwFxg
   RC3v94YaMWKCg2/1NXLkdTW/YckS9+uu85IuO/jrmb39JJ73FMoc3HNyrvGXXnqtwTGIiCRT+NlZ
   5+dzsloKHwGnAHfHHjSzgcBoYCDQHfinme3p7hWNvWBsN1CsTZti9iJwh+JiKCiAf/0LzjqLjDff
   4LuFX1M6eQbf3zRTLQERadWSkhTcfS4E06OqORl4wt23AIvNbD5wKPDfxl6zffuyGo9nZJTD+vXw
   6KNBMigvh7w8uP9+6NQJgOP33ltJQETahOY2prA7VRPAUoIWQ6ONHz+CBQsmVBlTOLrHxdzZfhX0
   6gXDhsHtt8Pw4aB9C0SkjUpYUjCzGUC3Gl661t1fbMCPqnGVWn5+fvRxbm4uubm5df6QyF/6U+6Y
   wOAvF/Hj5W8zaO0q2u9zKUy+DXr2bEBIIiLNX1FREUVFRQ16T1JXNJvZq8Av3P2D8PnVAO5+c/h8
   GjDR3d+u9j7f7riLi+HKK4MuotNPh4yMxtyCiEiLUZ8Vzc0hKVzp7u+HzwcCjxOMI3QH/gn0q54B
   GpUU3NU9JCJtUn2SQlL2UzCzU8xsCTAEKDSzlwHcfQ7wFDAHeBm4bPs//Wu9eFx/nIhIa6KCeCIi
   bUSzbSmIiEjzpKQgIiJRSgoiIhKlpCAiIlFKCiIiEqWkICIiUUoKIiISpaQgIiJRSgoiIhKlpCAi
   IlFKCiIiEqWkICIiUUoKIiISpaQgIiJRSgoiIhKlpCAiIlFKCiIiEqWkICIiUUoKIiISpaQgIiJR
   SgoiIhKlpCAiIlFJSQpmdpqZfWxm5WZ2YMzx3mZWYmYzw6+pyYhPRKStSlZL4SPgFKC4htfmu/sB
   4ddlTRxXs1BUVJTsEBJK99eyteb7a833Vl9JSQruPtfdP03GtVuC1v6Lqftr2Vrz/bXme6uv5jim
   0CfsOioysyOTHYyISFvSLlE/2MxmAN1qeOlad3+xlrd9CfRw99XhWMPzZraPu69LVJwiIlLJ3D15
   Fzd7FfiFu3/QkNfNLHlBi4i0YO5udb2esJZCA0QDNLOdgNXuXm5mfYH+wMLqb9jWTYmIyPZJ1pTU
   U8xsCTAEKDSzl8OXhgGzzGwm8DRwsbuvSUaMIiJtUVK7j0REpHlpjrOP6sXMbjCzWWb2oZn9y8x6
   JDumeDKz35vZJ+E9PmtmnZMdUzzVtoCxJTOzUWY218w+M7NfJTueeDOzv5jZSjP7KNmxxJuZ9TCz
   V8PfydlmNj7ZMcWTmWWY2dvh5+UcM/tdree21JaCmWVHZiWZ2ThgsLtfkOSw4sbMjgX+5e4VZnYz
   gLtfneSw4sbM9gIqgLupY7JBS2FmqcA84BhgGfAucKa7f5LUwOLIzL4PrAcedvdByY4nnsysG9DN
   3T80s47A+8APW9l/vw7uvtHM2gFvAFe6+xvVz2uxLYVq01Q7At8kK5ZEcPcZ7l4RPn0b2COZ8cRb
   K1zAeCjBavzF7r4FeBI4OckxxZW7vw6sTnYcieDuK9z9w/DxeuATYPfkRhVf7r4xfJgOpAKrajqv
   xSYFADO7ycy+AM4Gbk52PAl0HvCPZAchdeoOLIl5vjQ8Ji2MmfUGDiD4Y6zVMLMUM/sQWAm86u5z
   ajqvOUxJrdW2FsC5+wRggpldDfwJOLdJA2yk+izwM7MJwGZ3f7xJg4uD7VzA2FK1zH5YqSLsOvob
   cHnYYmg1wp6H/cPxyVfMLNfdi6qf16yTgrsfW89TH6cF/iW9rfszs3OA44CjmySgOGvAf7/WYBkQ
   O9mhB0FrQVoIM0sDngEedffnkx1Porj7d2ZWCBwMFFV/vcV2H5lZ/5inJwMzkxVLIpjZKOCXwMnu
   vinZ8SRYa1iM+B7QPyz/ng6MBl5IckxST2ZmwP3AHHf/c7LjiTcz28nMuoSPM4FjqeUzsyXPPvob
   MAAoBxYAl7r7V8mNKn7M7DOCAaHIYNBbramUuJmdAtwB7AR8B8x09x8kN6rGMbMfAH8mGMS7391r
   nfbXEpnZEwQLTHcEvgJ+7e4PJDeq+AiLbxYD/6OyK/Aad5+WvKjix8wGAQ8RNARSgEfc/fc1nttS
   k4KIiMRfi+0+EhGR+FNSEBGRKCUFERGJUlIQEZEoJQUREYlSUhARkSglBUkqM9vRzGaGX8vNbGn4
   eLWZfdzEsZxsZnvHPP+NmTV4NXm4gK3G8tJmto+Z/Tsssf2pmV3XmJjriKGmexkePi4ys4MScV1p
   +ZQUJKnc/Vt3P8DdDwDuAv4YPt6foLR2XIUlrmtzCjAwJraJ7v6vOF47E/g78Ft33wsYDBxuZolY
   lFjTvfw78hTVapJaKClIc2Mx31PN7J5w05NXzCwDwMxyzOxlM3vPzIrNbEB4vHf4V/gsM/tnZOMl
   M3vQzO4ys/8Ct9T0fjM7HDgR+L2ZfWBmfcP3/Tj8GYeY2X/CTUreNrOO4fWKzez98OuwbdzbGOAN
   d/8ngLuXAHnAr8Jr5JvZL6L/EMF99wwfPxfGO9vMLow5Z72Z3RjG9ZaZ7bKte6nyj202wszeDON/
   ysyywuM3W7DhzCwzq3Hlq7ROSgrSnPUHCtx9X2ANEPlQuwcY5+4HE9SHmhoenww84O6DgccIymhE
   7A4c5u5X1vR+d3+ToFbRle5+oLsvJPyLOqxl9CQw3t33JyhQWEJQgvhYdz8IOKPa9WoykGDzlqjw
   Otlhdc7qf73HPj8vjPcQYLyZdQ2PdyAogbI/QZmGC+u6l9gfbmY7AROAo8N7eB+4wsx2INhgZp/w
   3/KGbdyXtCLNukqqtHmL3P1/4eP3gd7hX7KHA08HNcyAoEYUwBDgh+HjR4Fbw8cOPO3uHn74HlbL
   +2Hr4nxGUGNrubu/D9FNWAiTRYGZDSaowbVnPe6ptuJ/adt43+VmFrm3HgQJ8x2CsuqF4fH3CQqd
   betakdeGECSqN8N/i3TgTYJaVJvM7H7gpfBL2gglBWnOSmMelwMZBK3b1eG4Q01q+yCM7DqVAqyp
   4/019bXX1v/+fwTJ4qfhWMW2qtnOAYZWCdasL7DR3VebWRlVW++R7rJcgtbJEHffZGavRl4DtsSc
   X0HV/6frM24ww93HVD9oZoeG1zyVoIurRZZvl4ZT95G0JBZuw7rIzE6FoOSxme0Xvv4mQTcOwFkE
   3SlVuPvaOt6/DuhU/S0Eey/vZmYHh+/JDpNAJ2BFeN7PCKqj1uUx4MjIjKZw4Pl24Lbw9cXAgeFr
   BwJ9wuOdCBLhJgv2th6yjevUdi/V7+u/wBFmlhNeM8vM+oetsS7u/jJwBcGAuLQRSgrS3Hgtj2Of
   nwWcb8HWgrOBk8Lj44BzzWxWeM7ltfys2t7/JPDLcNC1b/SNwZ7Lo4HJ4XteAdoTjGWcHR4bQLCp
   fW2xE+6LcRLBboHzgK+Bz9z9T+EpzwA7mNlsYCxBMgKYBrQzsznA74C3arlO7LhBjfdSLZ5vgHOA
   J8J/szfD+8gGXgyPvU7QIpI2QqWzRZLEzE4Gfg0c7+4rtnW+SFNQUhARkSh1H4mISJSSgoiIRCkp
   iIhIlJKCiIhEKSmIiEiUkoKIiEQpKYiISNT/AyqzUd5HUH+2AAAAAElFTkSuQmCC
   "
   >
   </div>
   
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <h2 id="glm-gamma-for-proportional-count-response">GLM: Gamma for proportional count response</h2>
   <h3 id="load-data">Load data</h3>
   <p> In the example above, we printed the <code>NOTE</code> attribute to learn about the
    Star98 dataset. Statsmodels datasets ships with other useful information. For
    example: </p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[16]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="k">print</span><span class="p">(</span><span class="n">sm</span><span class="o">.</span><span class="n">datasets</span><span class="o">.</span><span class="n">scotland</span><span class="o">.</span><span class="n">DESCRLONG</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
   
   This data is based on the example in Gill and describes the proportion of
   voters who voted Yes to grant the Scottish Parliament taxation powers.
   The data are divided into 32 council districts.  This example&apos;s explanatory
   variables include the amount of council tax collected in pounds sterling as
   of April 1997 per two adults before adjustments, the female percentage of
   total claims for unemployment benefits as of January, 1998, the standardized
   mortality rate (UK is 100), the percentage of labor force participation,
   regional GDP, the percentage of children aged 5 to 15, and an interaction term
   between female unemployment and the council tax.
   
   The original source files and variable information are included in
   /scotland/src/
   
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p> Load the data and add a constant to the exogenous variables:</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[17]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">data2</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">datasets</span><span class="o">.</span><span class="n">scotland</span><span class="o">.</span><span class="n">load</span><span class="p">()</span>
   <span class="n">data2</span><span class="o">.</span><span class="n">exog</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">add_constant</span><span class="p">(</span><span class="n">data2</span><span class="o">.</span><span class="n">exog</span><span class="p">,</span> <span class="n">prepend</span><span class="o">=</span><span class="bp">False</span><span class="p">)</span>
   <span class="k">print</span><span class="p">(</span><span class="n">data2</span><span class="o">.</span><span class="n">exog</span><span class="p">[:</span><span class="mi">5</span><span class="p">,:])</span>
   <span class="k">print</span><span class="p">(</span><span class="n">data2</span><span class="o">.</span><span class="n">endog</span><span class="p">[:</span><span class="mi">5</span><span class="p">])</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
   [[  7.12000000e+02   2.10000000e+01   1.05000000e+02   8.24000000e+01
       1.35660000e+04   1.23000000e+01   1.49520000e+04   1.00000000e+00]
    [  6.43000000e+02   2.65000000e+01   9.70000000e+01   8.02000000e+01
       1.35660000e+04   1.53000000e+01   1.70395000e+04   1.00000000e+00]
    [  6.79000000e+02   2.83000000e+01   1.13000000e+02   8.63000000e+01
       9.61100000e+03   1.39000000e+01   1.92157000e+04   1.00000000e+00]
    [  8.01000000e+02   2.71000000e+01   1.09000000e+02   8.04000000e+01
       9.48300000e+03   1.36000000e+01   2.17071000e+04   1.00000000e+00]
    [  7.53000000e+02   2.20000000e+01   1.15000000e+02   6.47000000e+01
       9.26500000e+03   1.46000000e+01   1.65660000e+04   1.00000000e+00]]
   [ 60.3  52.3  53.4  57.   68.7]
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <h3 id="fit-and-summary">Fit and summary</h3>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[18]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">glm_gamma</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">GLM</span><span class="p">(</span><span class="n">data2</span><span class="o">.</span><span class="n">endog</span><span class="p">,</span> <span class="n">data2</span><span class="o">.</span><span class="n">exog</span><span class="p">,</span> <span class="n">family</span><span class="o">=</span><span class="n">sm</span><span class="o">.</span><span class="n">families</span><span class="o">.</span><span class="n">Gamma</span><span class="p">())</span>
   <span class="n">glm_results</span> <span class="o">=</span> <span class="n">glm_gamma</span><span class="o">.</span><span class="n">fit</span><span class="p">()</span>
   <span class="k">print</span><span class="p">(</span><span class="n">glm_results</span><span class="o">.</span><span class="n">summary</span><span class="p">())</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
                    Generalized Linear Model Regression Results                  
   ==============================================================================
   Dep. Variable:                      y   No. Observations:                   32
   Model:                            GLM   Df Residuals:                       24
   Model Family:                   Gamma   Df Model:                            7
   Link Function:          inverse_power   Scale:                 0.0035842831735
   Method:                          IRLS   Log-Likelihood:                -83.017
   Date:                Wed, 20 May 2015   Deviance:                     0.087389
   Time:                        21:52:42   Pearson chi2:                   0.0860
   No. Iterations:                     6                                         
   ==============================================================================
                    coef    std err          z      P&gt;|z|      [95.0% Conf. Int.]
   ------------------------------------------------------------------------------
   x1          4.962e-05   1.62e-05      3.060      0.002      1.78e-05  8.14e-05
   x2             0.0020      0.001      3.824      0.000         0.001     0.003
   x3         -7.181e-05   2.71e-05     -2.648      0.008        -0.000 -1.87e-05
   x4             0.0001   4.06e-05      2.757      0.006      3.23e-05     0.000
   x5         -1.468e-07   1.24e-07     -1.187      0.235     -3.89e-07  9.56e-08
   x6            -0.0005      0.000     -2.159      0.031        -0.001 -4.78e-05
   x7         -2.427e-06   7.46e-07     -3.253      0.001     -3.89e-06 -9.65e-07
   const         -0.0178      0.011     -1.548      0.122        -0.040     0.005
   ==============================================================================
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <h2 id="glm-gaussian-distribution-with-a-noncanonical-link">GLM: Gaussian distribution with a noncanonical link</h2>
   <h3 id="artificial-data">Artificial data</h3>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[19]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">nobs2</span> <span class="o">=</span> <span class="mi">100</span>
   <span class="n">x</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">arange</span><span class="p">(</span><span class="n">nobs2</span><span class="p">)</span>
   <span class="n">np</span><span class="o">.</span><span class="n">random</span><span class="o">.</span><span class="n">seed</span><span class="p">(</span><span class="mi">54321</span><span class="p">)</span>
   <span class="n">X</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">column_stack</span><span class="p">((</span><span class="n">x</span><span class="p">,</span><span class="n">x</span><span class="o">**</span><span class="mi">2</span><span class="p">))</span>
   <span class="n">X</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">add_constant</span><span class="p">(</span><span class="n">X</span><span class="p">,</span> <span class="n">prepend</span><span class="o">=</span><span class="bp">False</span><span class="p">)</span>
   <span class="n">lny</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">exp</span><span class="p">(</span><span class="o">-</span><span class="p">(</span><span class="o">.</span><span class="mo">03</span><span class="o">*</span><span class="n">x</span> <span class="o">+</span> <span class="o">.</span><span class="mo">0001</span><span class="o">*</span><span class="n">x</span><span class="o">**</span><span class="mi">2</span> <span class="o">-</span> <span class="mf">1.0</span><span class="p">))</span> <span class="o">+</span> <span class="o">.</span><span class="mo">001</span> <span class="o">*</span> <span class="n">np</span><span class="o">.</span><span class="n">random</span><span class="o">.</span><span class="n">rand</span><span class="p">(</span><span class="n">nobs2</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <h3 id="fit-and-summary">Fit and summary</h3>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[20]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">gauss_log</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">GLM</span><span class="p">(</span><span class="n">lny</span><span class="p">,</span> <span class="n">X</span><span class="p">,</span> <span class="n">family</span><span class="o">=</span><span class="n">sm</span><span class="o">.</span><span class="n">families</span><span class="o">.</span><span class="n">Gaussian</span><span class="p">(</span><span class="n">sm</span><span class="o">.</span><span class="n">families</span><span class="o">.</span><span class="n">links</span><span class="o">.</span><span class="n">log</span><span class="p">))</span>
   <span class="n">gauss_log_results</span> <span class="o">=</span> <span class="n">gauss_log</span><span class="o">.</span><span class="n">fit</span><span class="p">()</span>
   <span class="k">print</span><span class="p">(</span><span class="n">gauss_log_results</span><span class="o">.</span><span class="n">summary</span><span class="p">())</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
                    Generalized Linear Model Regression Results                  
   ==============================================================================
   Dep. Variable:                      y   No. Observations:                  100
   Model:                            GLM   Df Residuals:                       97
   Model Family:                Gaussian   Df Model:                            2
   Link Function:                    log   Scale:               1.05311425588e-07
   Method:                          IRLS   Log-Likelihood:                 662.92
   Date:                Wed, 20 May 2015   Deviance:                   1.0215e-05
   Time:                        21:52:42   Pearson chi2:                 1.02e-05
   No. Iterations:                     7                                         
   ==============================================================================
                    coef    std err          z      P&gt;|z|      [95.0% Conf. Int.]
   ------------------------------------------------------------------------------
   x1            -0.0300    5.6e-06  -5361.316      0.000        -0.030    -0.030
   x2         -9.939e-05   1.05e-07   -951.091      0.000     -9.96e-05 -9.92e-05
   const          1.0003   5.39e-05   1.86e+04      0.000         1.000     1.000
   ==============================================================================
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>

   <script src="https://c328740.ssl.cf1.rackcdn.com/mathjax/latest/MathJax.js?config=TeX-AMS_HTML"type="text/javascript"></script>
   <script type="text/javascript">
   init_mathjax = function() {
       if (window.MathJax) {
           // MathJax loaded
           MathJax.Hub.Config({
               tex2jax: {
               // I'm not sure about the \( and \[ below. It messes with the
               // prompt, and I think it's an issue with the template. -SS
                   inlineMath: [ ['$','$'], ["\\(","\\)"] ],
                   displayMath: [ ['$$','$$'], ["\\[","\\]"] ]
               },
               displayAlign: 'left', // Change this to 'center' to center equations.
               "HTML-CSS": {
                   styles: {'.MathJax_Display': {"margin": 0}}
               }
           });
           MathJax.Hub.Queue(["Typeset",MathJax.Hub]);
       }
   }
   init_mathjax();

   // since we have to load this in a ..raw:: directive we will add the css
   // after the fact
   function loadcssfile(filename){
       var fileref=document.createElement("link")
       fileref.setAttribute("rel", "stylesheet")
       fileref.setAttribute("type", "text/css")
       fileref.setAttribute("href", filename)

       document.getElementsByTagName("head")[0].appendChild(fileref)
   }
   // loadcssfile({{pathto("_static/nbviewer.pygments.css", 1) }})
   // loadcssfile({{pathto("_static/nbviewer.min.css", 1) }})
   loadcssfile("../../../_static/nbviewer.pygments.css")
   loadcssfile("../../../_static/ipython.min.css")
   </script>