1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
|
<!DOCTYPE html>
<html>
<head>
<meta charset="utf-8" />
<meta name="generator" content="pandoc" />
<meta http-equiv="X-UA-Compatible" content="IE=EDGE" />
<meta name="viewport" content="width=device-width, initial-scale=1" />
<meta name="author" content="Johannes Ranke" />
<title>Short introduction to mkin</title>
<script>// Pandoc 2.9 adds attributes on both header and div. We remove the former (to
// be compatible with the behavior of Pandoc < 2.8).
document.addEventListener('DOMContentLoaded', function(e) {
var hs = document.querySelectorAll("div.section[class*='level'] > :first-child");
var i, h, a;
for (i = 0; i < hs.length; i++) {
h = hs[i];
if (!/^h[1-6]$/i.test(h.tagName)) continue; // it should be a header h1-h6
a = h.attributes;
while (a.length > 0) h.removeAttribute(a[0].name);
}
});
</script>
<style type="text/css">detaiks.chunk-details > summary.chunk-summary {
text-align: right;
}
details.chunk-details[open] > summary.chunk-summary::after {
content: "Hide";
}
details.chunk-details[open] > summary.chunk-summary > span.chunk-summary-text {
display: none;
}
</style>
<style type="text/css">
code{white-space: pre-wrap;}
span.smallcaps{font-variant: small-caps;}
span.underline{text-decoration: underline;}
div.column{display: inline-block; vertical-align: top; width: 50%;}
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
ul.task-list{list-style: none;}
</style>
<style type="text/css">
code {
white-space: pre;
}
.sourceCode {
overflow: visible;
}
</style>
<style type="text/css" data-origin="pandoc">
pre > code.sourceCode { white-space: pre; position: relative; }
pre > code.sourceCode > span { display: inline-block; line-height: 1.25; }
pre > code.sourceCode > span:empty { height: 1.2em; }
.sourceCode { overflow: visible; }
code.sourceCode > span { color: inherit; text-decoration: inherit; }
div.sourceCode { margin: 1em 0; }
pre.sourceCode { margin: 0; }
@media screen {
div.sourceCode { overflow: auto; }
}
@media print {
pre > code.sourceCode { white-space: pre-wrap; }
pre > code.sourceCode > span { text-indent: -5em; padding-left: 5em; }
}
pre.numberSource code
{ counter-reset: source-line 0; }
pre.numberSource code > span
{ position: relative; left: -4em; counter-increment: source-line; }
pre.numberSource code > span > a:first-child::before
{ content: counter(source-line);
position: relative; left: -1em; text-align: right; vertical-align: baseline;
border: none; display: inline-block;
-webkit-touch-callout: none; -webkit-user-select: none;
-khtml-user-select: none; -moz-user-select: none;
-ms-user-select: none; user-select: none;
padding: 0 4px; width: 4em;
color: #aaaaaa;
}
pre.numberSource { margin-left: 3em; border-left: 1px solid #aaaaaa; padding-left: 4px; }
div.sourceCode
{ }
@media screen {
pre > code.sourceCode > span > a:first-child::before { text-decoration: underline; }
}
code span.al { color: #ff0000; font-weight: bold; } /* Alert */
code span.an { color: #60a0b0; font-weight: bold; font-style: italic; } /* Annotation */
code span.at { color: #7d9029; } /* Attribute */
code span.bn { color: #40a070; } /* BaseN */
code span.bu { color: #008000; } /* BuiltIn */
code span.cf { color: #007020; font-weight: bold; } /* ControlFlow */
code span.ch { color: #4070a0; } /* Char */
code span.cn { color: #880000; } /* Constant */
code span.co { color: #60a0b0; font-style: italic; } /* Comment */
code span.cv { color: #60a0b0; font-weight: bold; font-style: italic; } /* CommentVar */
code span.do { color: #ba2121; font-style: italic; } /* Documentation */
code span.dt { color: #902000; } /* DataType */
code span.dv { color: #40a070; } /* DecVal */
code span.er { color: #ff0000; font-weight: bold; } /* Error */
code span.ex { } /* Extension */
code span.fl { color: #40a070; } /* Float */
code span.fu { color: #06287e; } /* Function */
code span.im { color: #008000; font-weight: bold; } /* Import */
code span.in { color: #60a0b0; font-weight: bold; font-style: italic; } /* Information */
code span.kw { color: #007020; font-weight: bold; } /* Keyword */
code span.op { color: #666666; } /* Operator */
code span.ot { color: #007020; } /* Other */
code span.pp { color: #bc7a00; } /* Preprocessor */
code span.sc { color: #4070a0; } /* SpecialChar */
code span.ss { color: #bb6688; } /* SpecialString */
code span.st { color: #4070a0; } /* String */
code span.va { color: #19177c; } /* Variable */
code span.vs { color: #4070a0; } /* VerbatimString */
code span.wa { color: #60a0b0; font-weight: bold; font-style: italic; } /* Warning */
</style>
<script>
// apply pandoc div.sourceCode style to pre.sourceCode instead
(function() {
var sheets = document.styleSheets;
for (var i = 0; i < sheets.length; i++) {
if (sheets[i].ownerNode.dataset["origin"] !== "pandoc") continue;
try { var rules = sheets[i].cssRules; } catch (e) { continue; }
var j = 0;
while (j < rules.length) {
var rule = rules[j];
// check if there is a div.sourceCode rule
if (rule.type !== rule.STYLE_RULE || rule.selectorText !== "div.sourceCode") {
j++;
continue;
}
var style = rule.style.cssText;
// check if color or background-color is set
if (rule.style.color === '' && rule.style.backgroundColor === '') {
j++;
continue;
}
// replace div.sourceCode by a pre.sourceCode rule
sheets[i].deleteRule(j);
sheets[i].insertRule('pre.sourceCode{' + style + '}', j);
}
}
})();
</script>
<style type="text/css">
/* for pandoc --citeproc since 2.11 */
div.csl-bib-body { }
div.csl-entry {
clear: both;
}
.hanging div.csl-entry {
margin-left:2em;
text-indent:-2em;
}
div.csl-left-margin {
min-width:2em;
float:left;
}
div.csl-right-inline {
margin-left:2em;
padding-left:1em;
}
div.csl-indent {
margin-left: 2em;
}
</style>
<style type="text/css">body {
background-color: #fff;
margin: 1em auto;
max-width: 700px;
overflow: visible;
padding-left: 2em;
padding-right: 2em;
font-family: "Open Sans", "Helvetica Neue", Helvetica, Arial, sans-serif;
font-size: 14px;
line-height: 1.35;
}
#TOC {
clear: both;
margin: 0 0 10px 10px;
padding: 4px;
width: 400px;
border: 1px solid #CCCCCC;
border-radius: 5px;
background-color: #f6f6f6;
font-size: 13px;
line-height: 1.3;
}
#TOC .toctitle {
font-weight: bold;
font-size: 15px;
margin-left: 5px;
}
#TOC ul {
padding-left: 40px;
margin-left: -1.5em;
margin-top: 5px;
margin-bottom: 5px;
}
#TOC ul ul {
margin-left: -2em;
}
#TOC li {
line-height: 16px;
}
table {
margin: 1em auto;
border-width: 1px;
border-color: #DDDDDD;
border-style: outset;
border-collapse: collapse;
}
table th {
border-width: 2px;
padding: 5px;
border-style: inset;
}
table td {
border-width: 1px;
border-style: inset;
line-height: 18px;
padding: 5px 5px;
}
table, table th, table td {
border-left-style: none;
border-right-style: none;
}
table thead, table tr.even {
background-color: #f7f7f7;
}
p {
margin: 0.5em 0;
}
blockquote {
background-color: #f6f6f6;
padding: 0.25em 0.75em;
}
hr {
border-style: solid;
border: none;
border-top: 1px solid #777;
margin: 28px 0;
}
dl {
margin-left: 0;
}
dl dd {
margin-bottom: 13px;
margin-left: 13px;
}
dl dt {
font-weight: bold;
}
ul {
margin-top: 0;
}
ul li {
list-style: circle outside;
}
ul ul {
margin-bottom: 0;
}
pre, code {
background-color: #f7f7f7;
border-radius: 3px;
color: #333;
white-space: pre-wrap;
}
pre {
border-radius: 3px;
margin: 5px 0px 10px 0px;
padding: 10px;
}
pre:not([class]) {
background-color: #f7f7f7;
}
code {
font-family: Consolas, Monaco, 'Courier New', monospace;
font-size: 85%;
}
p > code, li > code {
padding: 2px 0px;
}
div.figure {
text-align: center;
}
img {
background-color: #FFFFFF;
padding: 2px;
border: 1px solid #DDDDDD;
border-radius: 3px;
border: 1px solid #CCCCCC;
margin: 0 5px;
}
h1 {
margin-top: 0;
font-size: 35px;
line-height: 40px;
}
h2 {
border-bottom: 4px solid #f7f7f7;
padding-top: 10px;
padding-bottom: 2px;
font-size: 145%;
}
h3 {
border-bottom: 2px solid #f7f7f7;
padding-top: 10px;
font-size: 120%;
}
h4 {
border-bottom: 1px solid #f7f7f7;
margin-left: 8px;
font-size: 105%;
}
h5, h6 {
border-bottom: 1px solid #ccc;
font-size: 105%;
}
a {
color: #0033dd;
text-decoration: none;
}
a:hover {
color: #6666ff; }
a:visited {
color: #800080; }
a:visited:hover {
color: #BB00BB; }
a[href^="http:"] {
text-decoration: underline; }
a[href^="https:"] {
text-decoration: underline; }
code > span.kw { color: #555; font-weight: bold; }
code > span.dt { color: #902000; }
code > span.dv { color: #40a070; }
code > span.bn { color: #d14; }
code > span.fl { color: #d14; }
code > span.ch { color: #d14; }
code > span.st { color: #d14; }
code > span.co { color: #888888; font-style: italic; }
code > span.ot { color: #007020; }
code > span.al { color: #ff0000; font-weight: bold; }
code > span.fu { color: #900; font-weight: bold; }
code > span.er { color: #a61717; background-color: #e3d2d2; }
</style>
</head>
<body>
<h1 class="title toc-ignore">Short introduction to mkin</h1>
<h4 class="author">Johannes Ranke</h4>
<h4 class="date">Last change 18 May 2023 (rebuilt 2023-05-19)</h4>
<div id="TOC">
<ul>
<li><a href="#abstract">Abstract</a></li>
<li><a href="#background">Background</a>
<ul>
<li><a href="#derived-software-tools">Derived software tools</a></li>
</ul></li>
<li><a href="#unique-features">Unique features</a></li>
<li><a href="#internal-parameter-transformations">Internal parameter
transformations</a>
<ul>
<li><a href="#confidence-intervals-based-on-transformed-parameters">Confidence
intervals based on transformed parameters</a></li>
<li><a href="#parameter-t-test-based-on-untransformed-parameters">Parameter
t-test based on untransformed parameters</a></li>
</ul></li>
<li><a href="#references">References</a></li>
</ul>
</div>
<p><a href="https://www.jrwb.de">Wissenschaftlicher Berater, Kronacher
Str. 12, 79639 Grenzach-Wyhlen, Germany</a><br /> Privatdozent at the
University of Freiburg</p>
<div id="abstract" class="section level1">
<h1>Abstract</h1>
<p>In the regulatory evaluation of chemical substances like plant
protection products (pesticides), biocides and other chemicals,
degradation data play an important role. For the evaluation of pesticide
degradation experiments, detailed guidance has been developed, based on
nonlinear optimisation. The <code>R</code> add-on package
<code>mkin</code> implements fitting some of the models recommended in
this guidance from within R and calculates some statistical measures for
data series within one or more compartments, for parent and
metabolites.</p>
<details class="chunk-details"><summary class="chunk-summary"><span class="chunk-summary-text">Code</span></summary>
<div class="sourceCode" id="cb1"><pre class="sourceCode r"><code class="sourceCode r"><span id="cb1-1"><a href="#cb1-1" aria-hidden="true" tabindex="-1"></a><span class="fu">library</span>(<span class="st">"mkin"</span>, <span class="at">quietly =</span> <span class="cn">TRUE</span>)</span>
<span id="cb1-2"><a href="#cb1-2" aria-hidden="true" tabindex="-1"></a><span class="co"># Define the kinetic model</span></span>
<span id="cb1-3"><a href="#cb1-3" aria-hidden="true" tabindex="-1"></a>m_SFO_SFO_SFO <span class="ot"><-</span> <span class="fu">mkinmod</span>(<span class="at">parent =</span> <span class="fu">mkinsub</span>(<span class="st">"SFO"</span>, <span class="st">"M1"</span>),</span>
<span id="cb1-4"><a href="#cb1-4" aria-hidden="true" tabindex="-1"></a> <span class="at">M1 =</span> <span class="fu">mkinsub</span>(<span class="st">"SFO"</span>, <span class="st">"M2"</span>),</span>
<span id="cb1-5"><a href="#cb1-5" aria-hidden="true" tabindex="-1"></a> <span class="at">M2 =</span> <span class="fu">mkinsub</span>(<span class="st">"SFO"</span>),</span>
<span id="cb1-6"><a href="#cb1-6" aria-hidden="true" tabindex="-1"></a> <span class="at">use_of_ff =</span> <span class="st">"max"</span>, <span class="at">quiet =</span> <span class="cn">TRUE</span>)</span>
<span id="cb1-7"><a href="#cb1-7" aria-hidden="true" tabindex="-1"></a></span>
<span id="cb1-8"><a href="#cb1-8" aria-hidden="true" tabindex="-1"></a></span>
<span id="cb1-9"><a href="#cb1-9" aria-hidden="true" tabindex="-1"></a><span class="co"># Produce model predictions using some arbitrary parameters</span></span>
<span id="cb1-10"><a href="#cb1-10" aria-hidden="true" tabindex="-1"></a>sampling_times <span class="ot">=</span> <span class="fu">c</span>(<span class="dv">0</span>, <span class="dv">1</span>, <span class="dv">3</span>, <span class="dv">7</span>, <span class="dv">14</span>, <span class="dv">28</span>, <span class="dv">60</span>, <span class="dv">90</span>, <span class="dv">120</span>)</span>
<span id="cb1-11"><a href="#cb1-11" aria-hidden="true" tabindex="-1"></a>d_SFO_SFO_SFO <span class="ot"><-</span> <span class="fu">mkinpredict</span>(m_SFO_SFO_SFO,</span>
<span id="cb1-12"><a href="#cb1-12" aria-hidden="true" tabindex="-1"></a> <span class="fu">c</span>(<span class="at">k_parent =</span> <span class="fl">0.03</span>,</span>
<span id="cb1-13"><a href="#cb1-13" aria-hidden="true" tabindex="-1"></a> <span class="at">f_parent_to_M1 =</span> <span class="fl">0.5</span>, <span class="at">k_M1 =</span> <span class="fu">log</span>(<span class="dv">2</span>)<span class="sc">/</span><span class="dv">100</span>,</span>
<span id="cb1-14"><a href="#cb1-14" aria-hidden="true" tabindex="-1"></a> <span class="at">f_M1_to_M2 =</span> <span class="fl">0.9</span>, <span class="at">k_M2 =</span> <span class="fu">log</span>(<span class="dv">2</span>)<span class="sc">/</span><span class="dv">50</span>),</span>
<span id="cb1-15"><a href="#cb1-15" aria-hidden="true" tabindex="-1"></a> <span class="fu">c</span>(<span class="at">parent =</span> <span class="dv">100</span>, <span class="at">M1 =</span> <span class="dv">0</span>, <span class="at">M2 =</span> <span class="dv">0</span>),</span>
<span id="cb1-16"><a href="#cb1-16" aria-hidden="true" tabindex="-1"></a> sampling_times)</span>
<span id="cb1-17"><a href="#cb1-17" aria-hidden="true" tabindex="-1"></a></span>
<span id="cb1-18"><a href="#cb1-18" aria-hidden="true" tabindex="-1"></a><span class="co"># Generate a dataset by adding normally distributed errors with</span></span>
<span id="cb1-19"><a href="#cb1-19" aria-hidden="true" tabindex="-1"></a><span class="co"># standard deviation 3, for two replicates at each sampling time</span></span>
<span id="cb1-20"><a href="#cb1-20" aria-hidden="true" tabindex="-1"></a>d_SFO_SFO_SFO_err <span class="ot"><-</span> <span class="fu">add_err</span>(d_SFO_SFO_SFO, <span class="at">reps =</span> <span class="dv">2</span>,</span>
<span id="cb1-21"><a href="#cb1-21" aria-hidden="true" tabindex="-1"></a> <span class="at">sdfunc =</span> <span class="cf">function</span>(x) <span class="dv">3</span>,</span>
<span id="cb1-22"><a href="#cb1-22" aria-hidden="true" tabindex="-1"></a> <span class="at">n =</span> <span class="dv">1</span>, <span class="at">seed =</span> <span class="dv">123456789</span> )</span>
<span id="cb1-23"><a href="#cb1-23" aria-hidden="true" tabindex="-1"></a></span>
<span id="cb1-24"><a href="#cb1-24" aria-hidden="true" tabindex="-1"></a><span class="co"># Fit the model to the dataset</span></span>
<span id="cb1-25"><a href="#cb1-25" aria-hidden="true" tabindex="-1"></a>f_SFO_SFO_SFO <span class="ot"><-</span> <span class="fu">mkinfit</span>(m_SFO_SFO_SFO, d_SFO_SFO_SFO_err[[<span class="dv">1</span>]], <span class="at">quiet =</span> <span class="cn">TRUE</span>)</span>
<span id="cb1-26"><a href="#cb1-26" aria-hidden="true" tabindex="-1"></a></span>
<span id="cb1-27"><a href="#cb1-27" aria-hidden="true" tabindex="-1"></a><span class="co"># Plot the results separately for parent and metabolites</span></span>
<span id="cb1-28"><a href="#cb1-28" aria-hidden="true" tabindex="-1"></a><span class="fu">plot_sep</span>(f_SFO_SFO_SFO, <span class="at">lpos =</span> <span class="fu">c</span>(<span class="st">"topright"</span>, <span class="st">"bottomright"</span>, <span class="st">"bottomright"</span>))</span></code></pre></div>
</details>
<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAwAAAAKgCAIAAACjiI2DAAAACXBIWXMAAA7DAAAOwwHHb6hkAAAgAElEQVR4nOzdeVxM+/8H8DNtWqZp1x7tC2nSTTshaxKVlChCSFFcu9zsrvVLt2vNnksk0eJeZMmeVEIlRJMRKa3am98f8/12+4WQmc4083r+8X2c+fTpc16Tuef7PmfO+XwoLBaLAAAAABAkQmQHAAAAAOhqKIAAAABA4KAAAgAAAIGDAggAAAAEDgogAAAAEDgogAAAAEDgoAACAAAAgYMCCAAAAAQOCiAAAAAQOCiAAAAAQOCgAAIAAACBgwIIAAAABA4KIAAAABA4KIAAAABA4KAAAgAAAIGDAggAAAAEDgogAAAAEDgogAAAAEDgoAACAAAAgYMCCAAAAAQOCiAAAAAQOCiAAAAAQOCgAAIAAACBgwIIAAAABA4KIAAAABA4KIAAAABA4KAAAgAAAIGDAggAAAAEDgogAAAAEDgogPjZy5cvx40b5+fnt3fvXrKzAAAA8BAKi8UiOwNwy+LFi729vc3NzQcOHHjjxg2y4wAAAPAKXAHiZ6tXr6bT6UVFRYqKitzby/Xr18ePH9+Vv/iTPn78+OrVq67fLwAA8A4UQPxMQkLizJkzQUFBf/75J9lZfsr+/fuNjIxUVFSmTJlSW1vb2r527VrR/xEREbGwsCAIYtu2bVpaWhoaGqGhoSwW6/Xr1zY2NkZGRrt27WL/1sqVK+vr68l5JwAAwBtQAPGzxMTEu3fvxsXFqaiofE//xsbGL2533JPb0tPTN2zYcO3atWfPnpWVlW3evLn1R2FhYY3/s2DBgvnz52dlZW3bti0jIyMvL+/mzZtxcXEbN24MDw9/8uTJkSNHKisrnz17RhCEoaFhl+UHAAAehAKIfyxZssTCwqKuro4giIiIiICAgLi4uLy8PE9PzylTprTtGRERoaOjo6+vv3DhwpaWllu3bvn4+Hh7e69Zs6btNkEQ69ev19XV1dHRWbx4cbuenwdoN+zIkSPPnTvH/pGlpeWNGzfadfjO91VYWDh37lwVFRUajebu7v78+fPP+zx48CA3N9fX11dISEhUVFRCQkJMTExcXFxISEhGRubp06cMBqOurk5UVHT9+vUrVqz4oT8sAADwIRbwET8/vwMHDtTW1urq6hYVFX2xz7Vr18zMzN6/f19dXe3l5bVu3bqbN2/KyMjk5uayWKy22xcvXjQ0NKyoqKirq7O3tz969Gjbn7YdcNy4cZ8Pe+jQoSlTprBYrBcvXujo6Fy9erVdB/Yvto6TkJBg1UZAQEC75KWlpZaWlnFxcZ+/KQcHhydPnrC3f/31VwkJCRkZGfbgFRUVs2bNcnJy+ueff9LS0pYtW9bpPy8AAPANEbILMOCk2bNnL1y4sLKy0sPDQ11d/Yt9rly5UlFR4enpSRBESUlJVVWVo6MjnU5v/VaodTslJcXDw4NGoxEE4e3tnZKSoqOj07Znx8NGR0cvXbq0sbHx1KlTfn5+KSkp7TrY29u3HcHZ2dnZ2flrby06Ojo8PDw8PHzcuHHtfpScnNyzZ08TExOCIG7dupWQkHD9+nVpaekZM2ZER0f7+Pjs2bOH3XPChAn79++/cOFCZGSkqanpunXrevTo8T1/WAAA4DMogPiKtbV1aWlpREREenr61/pQqdRZs2YtXbqUIIi6urqGhobs7GxpaenWDq3bLBaLQqGwt4WFhZubm9v+9JvD0mg0a2vrq1evxsTEnD179vTp0+06ZGRktB0hMTFx7dq1rS/NzMzY0xe1tLRMnTq1oqIiNTX1izczHT9+3N3dnb2dkJDg5uZmaWlJEERAQMC5c+d8fHzYP/rnn38GDBjQ3NwcHh4eFxd34MCByMjIBQsWdPD3BAAAfoV7gPhN7969bW1t5eTkvtZh+PDhR44cKSkpqa+vd3d3v3Xr1td6DhkyJDY2trq6uqGhITo6esiQIR3s94vDenl5bdy4UU5OTltb+5v7dXZ2vttG6+SN8fHxTCbz3LlzrdUPi8XKy8tramoiCKK+vv7ixYsjR45k/8jU1DQpKYl9hSk2NrZfv36tvxIZGRkcHFxTUyMnJ6epqWlqalpaWtrxHxMAAPgVCiC+cuvWrdLS0pSUlA4e86bT6SEhITY2Nrq6urq6uqNGjfpaz5EjR06aNMnc3NzY2HjAgAG+vr4d7PqLw7q4uNy/f9/Pz++H9ttOampqamoqjUaTlpaWlpaeOnVqQ0ODkZHRmzdvCIL4559/zMzMWgs+b2/vMWPGWFhYGBkZKSsrL1q0iN0eExMzduxYcXFxLS0te3v7Pn36REREzJs37zszAAAAn8FM0HzFzs5u06ZNO3funDhx4oQJE8iOAwAAwKNwBYh/nDt3jkajOTg4TJ8+vfW2XwAAAPgcCiA+0dzcvHz58nXr1hEEMXLkSCUlJTqdXlJSQnYuAAAAXoSvwAAAAEDg4AoQAAAACBwUQAAAACBwUAABAACAwEEBBAAAAAIHBRAAAAAIHC6uBdbS0pKVlVVcXNzU1KSurk6n04WEUG8BQFeoqKig0WgUCiU9PT0nJ8fCwsLY2JjsUADAQ7j1GPzVq1fnzJmjqqrKXpOcyWQyGIyoqKiBAwd+5wi3bt26cOECN7IBQOdoaWkFBgaSneLb9u3bt3379uzs7C1bthw4cMDBweHatWsrV66cOXPmd46A4w8Ar+H48YdbBZCpqem5c+d0dXVbWxgMhoeHx717975zhMWLFz958uT7CyYA4KqKiopjx44xGAyyg3xbr169srOzaTSapqZmVlaWvLx8SUmJra1tfn7+F/uPGjWq3aGpqqpKXV19zpw5XZIXAL6BG8cfbn0F1tzcrKqq2rZFWVn5R4stR0fH1sUsAYBcRUVFx44dIzvFd1FSUmpoaCAIQlZWlv3Nu7CwcAf94+Liamtr27Y4OTmpqqouWbKEqzkB4Dtx4/jDrQIoKCjIwsLC1dVVQ0ODQqEwmczz588HBQVxanwGg7F27dqioiIrK6tevXppa2sPGjSIU4MDQLe2fv16W1vbsWPHGhgYODo6Dh8+PDExMTg4+Gv9xcXFxcXF27aIiHDx/kgA4AXc+o88MDDQxcUlOTmZyWQSBKGlpZWUlKSpqcmRwRkMhqGhoaioqJiY2MWLF7W0tGRlZXV0dM6cOYP7rAFgxIgR9+/fj4+Pl5WVNTIyUlJSio2NNTIyIjsXAPAQLp7lqKurW1patj4Fxr4bmiN+/fVXRUXF9PR0IyOjEydOTJo0KScnZ8iQIcnJyc7OzpzaCwB0X7Kysn5+fmSnAADexa0C6OefAuvAy5cvzc3NCwsLtbW1vby8fHx8Xr16ZWdnl5ubiwIIAAAAvolbBdC8efMSExO//ymwGzduvHv3rm1LRkZGnz59vtjZ2tr66NGj8+fPf/z4cUhICIVC6d2795UrV7Zs2cLBtwDQgSdPniQkJJCdggMGDRpkbW1NdgoAgK7GK0+BXblyJScnp23Lw4cPP3369MXOv//++7lz50aMGNHc3Lxr1y5DQ0NjY+MRI0Y4OTlxJDzAN505c+bSpUv29vZkB/kpjx49ys/PRwEEAAKIV54CW716dbsWKysrOTm5L3aWlJR8+fLl+fPnb9y4cfLkyXXr1hkbG5uYmHDyDQB0iMViOTk5hYeHkx3kp0RFRd25c4fsFAAAJOiWT4ERBCEqKuru7u7u7l5cXPz48WN3d3dOjQwAAAB8j4tPgWlqagYEBHBvfLZt27aZmZn5+Pjo6elxe18AAADAH7r3rDk3btywt7dvbGw0MTH5888/yY4DAAAA3QO3rgDt2bPn+vXrn7f/9ddfnNpFfX29t7f30aNHHRwcTExMVq9ebW1t3b9/f06NDwAAAPyKWwXQlClT7t27V1lZyb21o3NycpSVlYcOHUoQxJ49e9zd3S9fvowCCADge+Tk5PznP/95//790KFDZ8+ejdU/QNBw6xMvJSUVGhp68uRJdoHCDcrKykwms66uTlxc3MnJSVpa+u7du1zaFwBAV6qpqTlx4sT79+8HDx5sa2vL8fFfvHgxZMiQX3/9ddSoUZGRkU+fPsVdBCBouFjy9+vXr1+/ftwbX1VV1dnZediwYRMmTEhPT5eWlk5NTc3LyzM0NOTeTgEAuK2iomLAgAHm5uZ6enpTpkwJDg4OCQnh7C6io6OnTZu2cOFCgiCcnJzU1NR27NjRo0cPzu4FgJd175ug9+/fHxgY+OLFC3Nz87S0tFWrVs2aNauD6RYBulJdXV1hYWFLSwvZQQiCICorK/fv3092CvguBw8etLe3Z09yduPGjfDw8ObmZs7uoqqqSkFBgb0tJSUlIiJSV1fH2V0A8LjuXQAJCQl5e3vv3LkzJCSESqXOnTu3rq5u27ZtqIGAdFu3blVRUXFwcNDW1r5x4wZX99XY2PjNPlVVVUePHuVqDOCUwsJCU1NT9ra6urqYmFhZWRlnd+Hs7Lxnz57Hjx/X1dWtWbOmX79+MjIynN0FAI/r3gVQO1u3bn3y5MmSJUt69+59+/ZtsuOA4Lpz587evXtzc3Nfv3596NAhb2/v76lR2tq5c+fEiRMNDAz09PQ2bdpEEASLxZoxY4a+vn6vXr3Yy979/fffs2fPHjZsWEJCQnh4uJ6eHp1OP3z4MEEQV65ccXNzs7OzMzExWb58OUEQixcvfvz48fr16zn/boHTbGxsTp8+XVNTQxBEYmKijIyMkpISZ3fh6Oi4dOlSZ2dneXn59PT048ePc3Z8AN7HP7f9p6amRkVFPX/+fOfOnbdu3fLy8iooKBAWFiY7F/C5nTt3Hjt2rF1jcXFxc3PzmDFj2C8/fvxIp9MlJCTa9hEVFd21a5elpeXXRr5//352djZBEObm5sOGDaNSqc3Nzc+ePausrDQ2Nl60aBFBEBcuXMjMzExLS7t58+ajR4/q6+utra3pdDpBEJcuXSosLKRSqfr6+qGhoZs3by4sLFyxYgVn3z5ww4QJE27cuNGrVy8VFZXq6upTp05xYy/Tp0+fPn06N0YG6Bb4qgCaMGGCsrJyWFiYmZkZi8V6/vw5bogGbvP29v58SdSkpKSUlJStW7cSBPHp0ycXF5eIiIjPv2Lo27dvByOPHj2aSqUSBDF27NjU1NSQkJDly5cfPHgwIyPj48eP7D4ODg5KSkqpqanFxcUTJ04kCKK2tjYjI0NLS2vw4MHs1fT09fUrKyvFxcU59I7hv1gs1vnz5+/fv9+rV68pU6a0K3B/BoVC+eOPP1atWlVSUmJgYCAqKsqpkQGgFf8UQCoqKklJSYcOHSoqKpo5c+bSpUvxRAN0gZ49e/bs2bNdo7Gx8enTp7dv306n00+fPu3r6ztkyJAfHZlCobA3hIWFm5ubU1JSQkND58+fP3v27ISEBPaPpKWlCYKQkJBYuHDhtGnTCIKorq7u0aPHjRs3Wm9xhX379nFjWZ6goKAHDx64urpeunRp7969t27d4myV+cWPFgBwCv8UQC4uLnPnzk1LSzM0NLx586aSktLq1asPHTpEdi4QRJKSknfu3Dl06FBhYeHy5ctdXV07MUhycjL7LpD4+PijR49ev3594sSJ/v7+d+7cef/+fduHywYPHrxixQofH5/a2lozM7PU1NQvDsjxJ4l4VlJSUklJSevLNWvWsE+H/Pz8OLWLd+/enT59+tWrV5KSkgRBjB49+vz5856enpwaH0BwvHz58vjx4/X19e7u7l05mzH/FEDJyclOTk5OTk4MBsPX13fu3LnXr19PTk4eNWoU2dFAEElJSQUFBf3MCHQ63cnJqbS01M/Pz8rKSkNDY+rUqfHx8f369Rs8ePDKlSsHDRrE7ung4DBs2LB+/fo1NjauWLFCU1Pz2bNn7UZTVFSsrq5etWrVmjVrfiZVt/DkyZOVK1d6eXnJy8sTBFFXV5eZmUl8vQDau3dvQUFB25aioqKOL+cwmUxNTU129UMQhJGR0Zs3bziTHkCQZGZmjhgxYtq0aVJSUi4uLn/++Wfnzhg7gX8KoMLCQjqdPn/+fPbLDRs2LF68ePbs2Y8ePcLjndAd2drasuepY1NXV7906VK7PiNGjGBvhIWFhYWFtbYPHTq0dRL21t969OgRF+PykkWLFtna2i5btszX13fo0KGXL1/esWNHB/2lpaXZ90u1EhISav0K8ouMjY2ZTObNmzft7e2Li4vj4+NPnjzJmfQAgmTr1q1r165lf0ltb2+/ZMkSFEA/zMbGZv78+QsXLpSVlb1582Z5efnkyZPv378fEhKCL8IABI2dnV1CQkJwcHBCQkJ9fX3HnSdNmtSu5ezZs+w70L9GXFz82LFjXl5e4uLiZWVlK1eu7OCBPgD4mrdv3+rr67O3DQwMmExml+2afwqgoUOHjh8/XldXV01Nrays7NixYyIiIlu3bqXT6bGxse7u7mQHBPgBc+bMwXyeP4lGox05cuTkyZO5ubncGN/JyenVq1evX79WVVVt/S4MAH7IwIED9+3bZ2trKyoq+scff7R+s98FuFgAtbS0ZGVlFRcXNzU1qaur0+l0ISHuzru4du3a0NBQJpOpr6/PvudRUlLyyJEjbm5u9vb2ysrKXN07AAeJiYmRHYFPeHl5eXl5cWlwERERXV1dLg0OIAiWLFni5+enoqIiKipqYmISExPTZbvmVgF09erVOXPmqKqqqqurEwTBZDIZDEZUVNTAgQO5tEc2eXl59m2PrWxsbPz8/ObMmePm5nbp0iVFRcW5c+fq6OhwNQYAAAB8k7i4+KlTpz5+/NjU1MTxGc87xq0CaN68eYmJiW3PjRgMhoeHx71797i0xw6sXr1aS0vrzp07JiYmr169sra2vn//fu/evbs+CQAAALTT7imErsGt76Sam5tVVVXbtigrK5N1T4OIiEh5efmHDx8GDhxoa2tbW1u7a9cuUpIAAAAAL+DWFaCgoCALCwtXV1cNDQ0KhcJkMs+fP/+T06J0WnV1dWNj44YNG2JjY2/fvp2WloalUuHnffz48eXLl2Sn+CkfPnwgOwIAADm4VQAFBga6uLgkJyezH2nT0tJKSkrS1NT8Wv/o6OiioqK2LUwmk1Nr6wgJCbHvv+7Zs+fChQufPn1qYWHBkZFBYPXu3fvYsWOtS1J0X1gOEwAEExefApOXl/fx8ZGSkmIymampqWVlZR0UQLW1ta3rO7K1tLRw6iszaWnpfv36HTx4sKio6OLFi3JycjNmzODIyCCwpk6dOnXqVLJTAABAJ3GrADp27FhwcLCIiEh4ePiOHTtsbGwWL168aNGir30L9nlFcvXqVfZCjxxx+vTpKVOmvH37VlxcvKGhwcrKilMjAwAAQLfDrQIoPDw8NzdXXFxcU1Pz4sWLdnZ2Hz58GDBgAFm3Aenq6t6+fbu6ulpCQiI8PHzSpEmXL1/m9rxEAAAAwJu4VQG0tLTIy8vTaDQpKSlFRUWCIKhUasdr63QBKpUqLCz822+/NTU1bd68mdwwAAAAQBZuXQHy9PS0srISExMbNmzY1KlTfXx8kpOTR44cyaXd/RAREZETJ04MGDDAxsamK2fdBgAAAB7BrQLo999/d3Z2FhYWtrOzu3LlSnJy8rhx4/z9/bm0ux+loaERHR09adKktLQ0NTU1suMAAABAl+LiU2Ctq14MHTp06NCh3NtR5wwePHj27NkTJky4du2aqKgo2XEAAACg6wj0XcArVqyQkZFZsWIF2UEAAACgSwl0ASQkJHT8+PHY2NhTp06RnQUAAAC6Dhe/AusW5OXl4+Lihg4dqq+v379/f7LjAAAAQFcQ6CtAbP369du7d6+HhwfWRQIAABAQKIAIgiDc3Nw8PDy8vb2bmppaG5OSkuzt7Q0NDWfPnl1WVkZiPAD4Uenp6fHx8RUVFa0tfLBwGwBwEAqg/9q0aZOEhERoaCj7ZXp6ekBAwMqVKy9cuCAkJOTr60tuPAD4fhs3bnRzc4uJiaHT6RkZGezGmTNnkpsKAHiKoN8D1EpISCg6OtrW1nb37t1z5sw5e/bsnDlz2DM3RkREKCgoVFVVcXBtMgDgnsjIyKysLAUFhezsbHd39/T09I7/4/306VN9fX3blqamJk4txgwAvAkF0L+kpaXj4+Pt7OwMDAyEhYVbvw5raWkhCAILhwF0F9LS0lJSUgRBmJqazpkzJyQkJCoqqoP+np6et2/fbttSVVVF+tI9AMBV+D/1/0dHR+fUqVOTJk2ysLDYu3fvX3/9lZaWNm3atBEjRrCPpwDA+7y9ve3s7Pbt20cQREhISFlZmZeX16dPn77WPyEhoez/69+/v4qKShdGBoCuhitA7Q0cOHDXrl2hoaF79uzZt2/f+/fvhw4diskSAbqRVatW2dnZlZaWEgRBoVDOnDkTExMjJydHdi4A4CEogL5g4sSJeXl569evv3r1qqSkJNlxAPjfnDlzvti+e/fuzg3YdvkdYWFhb29vb2/vzg0FAHzpCwVQZmZmZmamj49PcXGxpqZm12fiBWFhYQUFBV5eXnFxccLCwmTHAeBzY8aMITsCQLfR0tKC21J/Xvu/4OHDhydOnLhq1ar6+noHB4cDBw6QEot0FArlwIEDQkJCs2fPJjsLAP9z/syoUaNu3bpFdi4A3nL69GltbW1xcXFra+usrCyy43Rv7QugHTt23L1718zMjEqlZmZmbtq0qdNDt7S0ZGRkJCcnX7hw4eHDh+xnqboRYWHh6OjorKys8PBwsrMACITNmzerqalJSEj07t2bSqXi+A7Q1tOnT0NDQ2NjY2tra4OCgsaPH9928l74Ue0LoJqaGnFxcfa2pKSkqKho58a9evWqiYnJggULoqOjT5069euvvxoaGt64ceOnwnatysrKQ4cODRo0KCoqKjIykuw4APzvjz/+yMnJ8ff3T0xMvHLlirq6OtmJvqyqqmr37t3h4eHd65gG3d2VK1fc3d379+8vLCw8efJkaWnpnJwcskN1Y+0LIB8fn7FjxxYUFBw8eNDJycnHx6dz486bNy8xMfHq1avHjx8/fvx4SkpKSkrKokWLfjpwF/n48WP//v3T0tKkpaXFxMTCwsIOHTpEdigAPtejRw8ajWZubn7v3j0bGxvevAJUXl7ev3//27dvUyiUGTNmbNmyhexEIChkZGTev3/P3m5qaiotLZWVlSU3UrfW/ibo1atXJyYm3rx5882bNxs2bLC3t+/cuM3Nzaqqqm1blJWVu9HMqgcPHhw6dOjevXsJgpg5c6aBgcHKlStlZGTc3NzIjgbAtwYOHOjv779s2bIJEya8e/dORkaG7ERfcOjQIUdHx/379xMEERAQYGRkFBoaKiKCJ2qB68aOHRseHr5kyRILC4uTJ0/+8ssvAvugEke0/4+2tLTU2tra2tq69aWCgkInxg0KCrKwsHB1ddXQ0KBQKEwm8/z580FBQT+bt6u8efPG2NiYva2qqiohIXH48OHJkyfLysoOGTKE3GwA/OrAgQPPnj0zMDDYvn17eno6b152LSoqMjExIQiCxWKxDw5lZWU9e/bs+LdaWlpev34tJyeHU3boNFlZ2du3b+/YsSM2NtbBwSEwMJDsRN1b+wKIvfoVQRAVFRX5+flmZmaZmZmdGDcwMNDFxSU5OZnJZBIEoaWllZSU1EGtunHjxlevXrVtefXqFYknVXZ2dlu2bPH396fRaHFxcfLy8sOGDYuLixs/fnxMTMygQYPICgbAx8rKyhQVFUtLS+l0Op1OJzvOl9na2i5atGjz5s1lZWUmJiYSEhLfrH6ysrI8PT0bGxs/fvw4ZcqUnTt3Yp0N6BwVFZXff/+d7BR8on2FkZaW1rqdnZ39M6ug02i0mTNnUiiU9PT0nJyc6urqDjpbWFi0u9R06dIlCQmJTu/9J7m7u9++fbtXr16KioosFuvUqVMEQdja2p44ccLT0/PChQsDBgwgKxsAv+LUCRhXycvLf/z4saWlpXfv3oWFhWpqat/8FV9f3/DwcG9v76qqqlGjRp08eRKzMgKQj9WhPn36dNzha/bu3WtoaNjQ0LB+/XptbW1fX18tLa19+/Z9/wgDBgxwdnbu3N45paysLD8/n70udKvExERlZeX79++TlQqAFAwGQ0NDo8t29+jRIzqd3mW7a6eD48/ixYs3b95cWlqam5vb2NiooKBQWlrawVAVFRUyMjKtLyMjIwMDAzkcF4DfceP40/4KUGhoaOt2bm6uoaFh5+qq9evXZ2dni4qK7t69OysrS15evqSkxNbWdubMmZ0v1rqcnJzc5+sHjR49+vDhwy4uLvHx8VZWVqQEA+B7pqamjY2NZKf4AvaTOPLy8vLy8tXV1Q0NDR0vmEOlUoWEhIqKijQ0NAiCyMnJYW8AALnaF0BtH/tydnYePHhw58ZVUlJqaGggCEJWVpY9Yzc/LSgxcuTIgwcPurq6ogYC4CBOnYBx1eTJk62trWVlZfX09Pbu3evr69s6d9oXCQkJrVq1avDgwb6+vq9evUpJSWl7pwEAkOXfAujevXsEQbQ7NXnw4EHn/g9+/fr1tra2Y8eONTAwcHR0HD58eGJiYnBw8E/G5R2jR48+dOjQ2LFjz5w54+DgQHYcAH7AqRMwrtLS0kpNTd25c2d2draXl9f06dO/+SshISH9+/e/dOlSv379tm3bhgfBAHgBhfW/uXnYN+UxmczU1NR+/fpRqdR79+6xH3rq3NDl5eXx8fEMBqO2tlZJSWnkyJFGRkbf/+tWVlZKSkoJCQmd23vXuH79+sSJE48ePTp8+HCyswBwV1FRkY2NDYPB4Mbg7BOwz5F1hbVbHH8ABAc3jj//XgH666+/CIIYNmxYTk4O+8pzQUFBQEBAp4eWlZX18/P7+Yi8bNCgQbGxsW5ubvv37x87dizZcQC6q//85z8ER0/AAAA61v4eoNevXxsYGLC3e/fuzZ7FBzpgZ2eXlJTk4uJSXl7+M7MGAAgyjp+AAQB0rH0BZGNj4+Hh4ePjw2Kx/vrrL0tLS7QSkdEAACAASURBVFJidS8WFhYpKSkjR4589+5dN1rvDIDX8OsJWF1d3datWy9duqSoqLhw4UJbW1uyEwHAZwXQvn379uzZc+LECSEhoUGDBuEM7DsZGRmlpqaOGDGitLR048aNmOYVoBM4eALW0tKSlZVVXFzc1NSkrq5Op9PZj6OSIjAw8OPHj+vWrXv58qWHh8fFixf79etHVhgAYPu3AAoODl67dm1YWBhBEOx1TJ89e/brr79GRESQlq5b0dTUvHXrlqurq6en57Fjxzp+MhYAPsepE7CrV6/OmTNHVVVVXV2dIAgmk8lgMKKiogYOHMjRvN+loaHhzJkz7969k5CQcHBwKCoqOnnyJAogANL9WwA5OTmJi4u3TkUPnSAnJ/f3339PmjTJxcUlNjaWRqORnQige+DsCdi8efMSExN1dXVbWxgMhoeHx9eeNVuxYkV+fn7blpycnBcvXnh6erbGY892kZiYeOTIkbaxv9l+8ODB2tpaX19fCoUSHBzMXjy1E+OgHe2C3L5nz56KigqCo/4tgFxdXQmC0NXVvXXr1vTp01euXHnt2rVNmzZxdn98T0JC4syZMyEhIfb29gkJCVpaWmQnAugGOHsC1tzczK6iWikrK7dO+fG5UaNGtVt7NS0tTVpaesKECeyX+vr67I2+ffu2Nn5nu5eX1+vXr6uqqsaNG1dWVrZr167jx49ramr+6DhoR7sgtzs7O3/tBKbz2i2NMXz48JMnT967d2/QoEGpqakjR47k7NIb348X1gL7GXv37tXQ0MCSYcA3umAtsJycnAMHDrBYrBUrVtjZ2aWmpnZunMjISCMjoyVLlkRERPzxxx/Lly/v27fvnj17vn8Ezh5/Pnz4MG3aNA0NDXNz85iYGE4NCyA4uHH8aX9XYHFx8cSJE8+ePTt58mR7e/uXL19yuOASGAEBAbt37x4zZgx7Jfmvyc3NTU1Nraqq6rJgADxr/vz5VCr1/v37N2/e3LRp0/r16zs3TmBg4D///KOjo/Phw4eSkhItLa2kpKRZs2ZxNu33U1BQOHjwIIPBePjwYdsTXAAgUfunwAwNDRcuXBgTE5OWlrZjxw7cxfIzxowZc/ny5XHjxmVkZGzYsKHdQyjNzc3e3t5paWmampr5+fnR0dFDhgwhKyoAL2CfgC1duvTnT8A0NTXxECsAdKD9FaDIyEh1dfVz586pqKiUlpayZyeDTjM1Nb179+6dO3fc3d3bXeY5duxYRUVFfn7+jRs3YmJipk2bRlZIAB7BPgGLjo4eM2YMTsAAgKvaF0BKSkpDhgzJzs5ubGycNWuWnp4eKbH4iZKS0qVLl1RUVKysrPLy8lrbMzIyxowZIyIiQhCEg4NDbW1tSUkJeTEByIcTMADoMu0LoMOHD0+cOHHVqlX19fUODg4HDhwgJRafERMT271796+//jpw4MBz586xG3V1dVvvac/Pz29paVFUVCQvIwD5cAIGAF2mfQG0Y8eOu3fvmpmZUanUzMxMPAbPQf7+/gkJCSEhIYsXL25qapoxY8bTp0+HDRs2Z84cR0fH7du3Y/5oEHA4AQOALtO+AKqpqWmdwlhSUlJUVLTTQ7e0tGRkZCQnJ1+4cOHhw4ctLS2dj8kvLC0t09PTHz9+PGTIkPLy8nv37s2ZM8fMzOzKlStYSBUAJ2AA0GXaPwXm4+MzduzYt2/fHjx48PDhwz4+Pp0bl6emoucpCgoKCQkJGzdu/OWXX6Kiotzc3MhOBMArOHgCBgDQsfYF0OrVqxMTE2/evPnmzZsNGzbY29t3btwfnYpeoAgJCa1YscLJyWnSpEmjRo1avHgxJowGIDh3AgYA8E3tCyCCIJydnZ2dndnbVVVV0tLSnRj3R6eiF0BWVlbjxo3btWvXvn37evXqFRcX17dvX7JDAZCJUydgAADf9G8B9PDhw2XLlpWVlY0ePTo4OHjZsmWFhYVZWVnFxcWdGDcoKMjCwsLV1VVDQ4NCoTCZzPPnzwcFBXEuebd37ty569evv3v3LiEhISgoaMSIEUVFRbgPGgQcR07AAAC+6d8CaNq0aePHj3dycjpy5Ei/fv0CAgJ8fHx69+7duXEDAwNdXFySk5OZTCZBEOyp6DU1Nb/Wf+7cuc+ePWvbkpub27lddxfXr1+fMmWKvLy8r6+vlZVV3759R4wYceTIkXZXzgAEAWdPwAAAvunfAujjx4/h4eEEQRgbG8fHx7O3f0bbqej37dvXQfVDEMTcuXPZpVKroKAgOTm5n8zAy5SVlZ8/f87elpOTk5SUtLa2Njc337Fjh7e3N7nZALoYZ0/AAFpVVVVRqVRcXIfP/VsACQsLszcUFBSkpKR+ctykpKS28xqvWbOmR48eBEH4+fl9sb+JiYmJiUnbFhkZGfYsyfzK39/f0tKypaVFT0/v4MGDv/76a1hYmKurq5+f35kzZyIjI1VUVMjOCNBFOH4CBpCRkeHn5/fq1SthYeH169cHBgaSnQh4S/t5gDjlyZMnAQEBKSkpmZmZmZmZdXV17A0u7a476tmzZ3p6uqamJoPB+P3338PCwgiCsLCwSE9PNzExMTMzO3r0KNkZAboIZ0/A4HucPXvWwsJCU1Nz6tSpHz58IDsOh7W0tHh4eKxcubKysjIzM3Pbtm137twhOxTwln8LoNevX0v/T9vtzo27aNGilJSUgoIC9qKGqqqqO3bs2LFjB4di8wlFRcWlS5du37599OjRrY09evRYu3bt33//vXPnzuHDh7948YLEhADAl+7evRsaGsouC2RlZflvItaCggIhISFPT0+CIHr16uXt7Z2SkkJ2KOAt/xZAtbW1H/6n7Xanh7azs0tISDh69GhoaGh9fT0n0goQOp1+7969kSNH2tjYbNq0qbGxkexEAFzE2RMwgiDS09Pj4+MrKipaWxISEjiR9KsKCwvnz5/v7u6+fft23j/inTt3LjAw0NHRUUNDY/v27bdv366qqiI7FCcpKCiUlpZ++vSJ/bKwsFBJSYncSMBr/i2AenzFz4xOo9GOHDliZWXVdkZE+E4iIiILFixIS0u7ffu2mZnZ1atXyU4EwC2cPQHbuHGjm5tbTEwMnU7PyMhgN86cOZNzedt7//69nZ2dvLz8pEmTrl+/zvsXVMTExOrq6tjbjY2NLBaLz+65lJWVnTBhgqOjo7+///jx42/evMm+GgTQqis+8V5eXl5eXl2wI77Uq1ev8+fPnz9/3t/f38bGZvPmzRoaGmSHAuCwnzzXaicyMjIrK0tBQSE7O9vd3T09PZ3b8wmdPXt2xIgRv/32G0EQLi4uGhoaJSUlvHzJwcvLa+jQob169dLT0/vPf/7j6uoqISFBdigOo9PpMTExFRUVNTU1BgYGVCqV7ETAW7h1EzRw1tixY588eaKrq0un0zds2NB66gYAn5OWlmbfSW1qajpnzpyQkBBu77GiokJRUZG9LSYmJiMjU1lZye2d/gwTE5Nz585duHBh2bJlpqame/fuJTsRh9XV1S1fvjw7OzsvL4/BYAgJCcXGxpIdCngLX13z5G+SkpJr16719/dfuHChiYnJpk2bJkyYgMktAD7n7e1tZ2c3a9asgICAkJAQNzc3Ly+v1ttBPvf06dN285BVVFTIysp+/x5HjBjh4uLi6elpamp64MABMTExHR2dzr+BLmFlZcXHNcHr169VVFTY18spFIq9vT3fT64LPwoFUDejra199uzZa9euLVy4cOfOnVu3brWxsSE7FABvWbVqlZ2dXWlpKUEQFArlzJkzMTExHUys+ueff+bl5bVtefv27Q8VQHQ6fevWrRMmTGAwGPb29mfPnsXJCbm0tbVLSkpycnKMjY0bGxv/+eefxYsXkx0KeAuFZxcotbKyUlJS4vaDG91XS0vLsWPHwsLCLC0tN2zYYGhoSHYi4HNFRUU2NjYMBoPsIF3he44/LBbr9OnTd+/e7d279/Tp0zF9Ea+JiYkJDg62sLDIy8uztbU9evQoqtLuixvHH9wD1F0JCQn5+fnl5eVZWVk5ODjMmjWrqKiI7FAAvMjR0ZEbwwYEBOzcuVNTU/P+/fs2Nja1tbXc2At0mqen56NHj4KCguLj448dO4bqB9pBAdS9SUhILF68OC8vT0FBgU6nL1iw4P3792SHAuAtzc3NHB+zuLj4/Pnzly9fDg0NPX78eK9eveLj4zm+F/hJysrKo0eP7tu3L9lBgBehAOIHcnJyGzZsePLkSXNzs4mJyaJFi1AGAbTixqw8b9++VVdXb310XF9fHwvXA3QvKID4h7Ky8s6dOx89etTQ0GBiYrJgwYJ2D7YACCZuTIFoYmLy7t079uoKRUVFcXFxDg4OHN8LAHAPCiB+o6amxi6DKBQKexIUrCYGwHE9evQ4ceLEtGnTNDU1TU1NFy5caGFhQXYoAPgBKID4k5qa2rZt2/Ly8hQVFW1sbLy8vB4+fEh2KAC+MmjQoIKCgjt37rx79y4oKIjsOADwY1AA8TNFRcW1a9e+fPnSyspq3LhxQ4YMSUhIaGlpITsXAJ8QEhLS0NAQExMjOwgA/DAUQPyPSqWGhoa+ePFi5syZq1evNjY2joiI4LOVnwEAAH4ICiBBISoq6u3tnZaWFhUVlZqa2rt37/nz52NueACA7uXp06dJSUl4xuXncbEAamlpycjISE5OvnDhwsOHD/HNC4+wt7ePiYnJysqi0WhDhgxxcnI6c+ZMY2Mj2bkAAOAbZs2aNXr06MjISHNz80OHDpEdp3vj1lpgV69enTNnjqqqqrq6OkEQTCaTwWBERUUNHDiQS3uEH6KhobF27dqwsLCzZ89GRkYGBwdPnTp1+vTpenp6ZEcDAIAvuHTp0v379/Py8nr06FFQUPDLL794eHhIS0uTnau74lYBNG/evMTERF1d3dYWBoPh4eFx7949Lu0ROkFMTMzLy8vLy+vZs2f79+93cHAwMDDw9/d3d3enUqlkpwMAgH89evRo6NChPXr0IAhCW1tbW1v72bNnmH+h07j1FVhzc7OqqmrbFmVlZZ5deBUMDAy2bNlSWFi4YMGCs2fPamlp+fn5XblyBV9cApDi7t277u7u9vb2v/3226dPn8iOAzzB0NDw9u3b7MPyu3fvCgoK2l5lgB/FrStAQUFBFhYWrq6uGhoaFAqFyWSeP38eU2XwOFFRUVdXV1dX15KSkhMnTixZsuTt27fe3t4+Pj7m5uZkpwMQFE+fPh0/fvymTZv09PR27do1Y8aMEydOkB0KyDd69OiDBw/279/f3Nz88uXLK1eulJWVJTtUN0bh3lUZBoORnJzMvlNdVVV19OjRmpqaX+vs7Ox8586dti1VVVVmZmYPHjzgUjz4Hjk5OSdOnDhx4oSoqOjEiRM9PT379OlDdiggR1FRkY2NDYPBIDtIV7CyslJSUkpISCBl7+Hh4SwWa/Xq1QRBNDQ0KCkpvX37VlJSkpQwwGuuX79eUFAwYMAAExMTsrN0HW4cf7h1BYggiPfv3ysrK0+cOFFGRobdkpCQMGbMmC92Pn36dH19fdsWJycnZWVl7sWD72FsbLx27dq1a9empaWdPHly1KhRVCrVw8PD3d3dzMyM7HQA/KmhoaH1zlYREREhIaGmpiZyIwHvGDRo0KBBg8hOwQ+4dQ/Qxo0b3dzcYmJi6HR6RkYGu7GDJQklJSXl/j8REREKhcKlePCjLC0tt23b9vr164MHD9bW1rq5uenq6i5atOjmzZvNzc1kpwPgK+PGjdu9e3dKSgqDwViwYIG1tTWNRiM7FAC/4VYBFBkZ+fDhw+jo6PPnz0+cOBHzDvMHCoVibW29ZcuWFy9enD17VkpKat68eaqqqv7+/rGxsZWVlWQHBPiCffv2kR3hxwwYMCAiImLJkiV2dnaVlZXHjh0jOxEAH+LWV2DS0tJSUlIEQbAXJA8JCYmKiuLSvoAUZmZmZmZm4eHhhYWFFy5c2L9/v7+/v6Wl5ahRo0aOHIlbhYBESUlJJSUlrS/XrFnDfnLYz8+PvFA/hv04AtkpAPgZt64AeXt729nZsU+8QkJCysrKvLy88DAnX9LS0po7d+7FixeZTOb8+fNfvHjh4uKipaU1Y8aMmJiY0tLS4uJizDQNXenJkycBAQEpKSmZmZmZmZl1dXXsDbJzAQAP4dYVoFWrVtnZ2ZWWlhIEQaFQzpw5ExMTIycnx6XdAS+QkpJycXFxcXEhCOLZs2f//PPPH3/84eXlRaFQhISEJk2aFBkZifkVoQssWrTI1tZ22bJlvr6+Q4cOvXz58o4dO8gOBQC8hYuPwf8kch9DhZ/X0tIiKyurpqY2evTolJSUR48eiYuLm5ubDxo0aODAgXZ2dpjBvXvpdo/BV1ZWBgcHy8vLJyYmPnv2rIOefn5+T548aduSm5trbGyclpbG5YwA8F262WPwIOCysrJqamrS09Nb7wbT1tYODQ29fv365s2bHzx4oK+vb29vb2dnZ2dnx14zDoCDaDTakSNHTp48mZub23HP1atXsy9Xt5o6daqSkhI30wEAyVAAAXe1vcQoLCw8ePDgwYMHEwTR0NDw4MGDW7duRUdHBwUFSUhI2NraWltbW1lZmZubs29ZBfh57NXuOu7Tu3fv3r17t22RlJQUEuLWLZIAwAtQAAG3mJmZUalUCwuLoKCgO3fu5OXlbdu2rfWnYmJitra2tra2ixYtIggiPz//zp079+7dO3r0aG5ubp8+fQYMGGBpaWlhYWFsbCwsLEze+4Buz9HR8dq1a2SnAADeggIIuEVISOjy5cvjxo1bsmQJQRBr164dPnz41zrr6+vr6+v7+voSBFFbW5uRkZGWlnbp0qVNmzYVFRWZmpqa/0/fvn1xfQh+CObqBIDPoQACLrK0tHzz5s379+/l5eVFRL73w8b+OszW1pb9sqqqKiMjIyMjIzU1NSIiIj8/X09Pr9//9O3bV0NDg2vvAPgBu7DmtitXrvz2228MBsPOzm7Lli24rQ2Ax6EAAq7r2bPnz/y6tLT0wIEDBw4cyH7Z0NDw+PHjR48eZWdnb9u27fHjx3V1dX379u3Tp0+fPn1MTEyMjY3V1NQ4ERz4RAeL8HBKXl7e5MmTDxw40Ldv36ioKDc3t3v37nF7pwDwM1AAQTcjJibWv3///v37t7aUlpZmZ2fn5OQ8fvw4Li7u6dOndXV1hoaGxsbGhoaGBgYG7O/XJCQkOBIgIyNj9erVr169srKyWrNmDZbsBYIgEhMTfXx8nJ2dCYJYs2bN0aNHX79+3atXL7JzAcBXoQCCbk9BQcHR0dHR0bG15ePHj7m5ubm5uXl5ecePH3/27NnLly979uypp6enp6enq6vL/l8dHZ0fnYvozZs3zs7O69ev79+///Hjx11dXe/cucMHq/ZWVVUlJSU1NDQMHz4cJV0niIiI1NfXs7dbWloaGxtFRUXJjQQAHUMBBHxITk7OxsbGxsamtaW5ubmwsPD58+fPnz9/8eLFnTt3nj9//vLlSykpKR0dHW1tbfaD0L169WL/79cuFyUlJTk7O0+bNo0gCDMzM319/efPn+vr63fN++KSwsJCe3t7CwsLSUnJRYsWxcfHW1lZkR2qmxk3bpyVlZWpqWm/fv2ioqL69u2L72EBeBwKIBAIwsLC2tra2traw4YNa9teXFxcUFBQUFDw6tWr9PT02NjY169fFxYW0mg0TU1NTU1NLS0tTU1NDQ0NDQ0NTU3NhoaGts/kCwkJtbS0dPm74bB169YFBgYuXbqUIIjY2NglS5bgofEfpaWldfHixQ0bNhw6dMje3n7Lli1kJwKAb0ABBAJNRUVFRUWl7bUituLiYsb/FBYWPnjwoKioiMFgvH37trm5+eLFi9ra2hUVFVVVVSkpKXl5ecrKympqaj179uyOj+jn5+f7+Piwt21tbYODg8nN002ZmZmdOnWK7BQA8L1QAAF8AbswsrS0bNfOYrEuXry4bt26/Px8LS2tSZMmZWVlXbx48d27d+wH/qlUqrKycs+ePVVUVHr27KmkpNSzZ09lZWXFNkh5Rx2g0+nnz58fNGgQQRDx8fHm5uZkJwIA4DoUQAA/gEKhVFdXM5nM9+/fa2hoTJkyxczMrG2H0tLSd+/elZSUMJnMkpKSkpKSBw8elJSUfPjwgf2ysrJSUVFR4Uvk5eXl5eXl5OTk5eVlZWWpVGrXvKmwsLBhw4bR6XQqlcpkMi9evNg1+wUAIBEKIIAf8PTp09DQ0PPnz5uZmf3111/jx49/9uxZ2zke2aVMByM0NTV9+PChtLT0w4cPZWVl7I2SkpK8vLyy//n48WN5eXlDQ4OsrKycnJxsG3JycjQaTUZGRkZGhkaj0Wg0aWlpOTk5aWlpaWlpcXHxTrwpeXn5+/fvP3jwoL6+fsCAAZ0bBACge0EBBPADrly54u7uzp6FaPLkyVu2bMnJyTE1Nf3+EURERNjfr32zZ0NDA7sS+tybN28qKioqKysrKyurqqrKy8vZGy0tLVQqVU5OjkqlSklJSUlJycnJSUlJSUpKSktL02g0CQkJKSkpGo0mLi5OpVKlpaV79OhBo9EkJSUNDAxkZWX54JF+AIDvwcUCqKWlJSsrq7i4uKmpSV1dnU6nY3Vl6O5kZGTev3/P3m5qaiotLZWVleXSvsTExJSVlX9oVp6Ghobq6ury8vLq6urq6uqampry8vKampqamprq6uqKiooPHz7U1NRUVFTU1dXV1NRUVlbW19dXVVV9+vSpvr4+LCwsNDSUS28HAICncKsAunr16pw5c1RVVdkL4jCZTAaDERUV1bqgAUB3NHbs2PDw8CVLllhYWJw8efKXX37R1NQkO9S/xMTE2DcSkR0EAIDXcasAmjdvXmJioq6ubmsLg8Hw8PD42vo4dXV1tbW1bVuampq4lA2g02RlZW/fvr1jx47Y2FgHB4fAwECyEwEAQGdwqwBqbm5WVVVt26KsrMxisb7Wf9y4cffv32/bUlVVhZsxgQepqKj8/vvvZKeAb6ioqKDRaBQKJT09PScnx8LCwtjYmOxQAMBDuFUABQUFWVhYuLq6amhoUCgUJpN5/vz5oKCgr/X//MnbxYsXKykpcSkeAPCxffv2bd++PTs7e8uWLQcOHHBwcFixYsXKlSu7YFl4AOguuFUABQYGuri4JCcnM5lMgiC0tLSSkpJ46m4JAOBX69evz87OFhUV3b17d1ZWlry8fElJia2tLQogAGjFxafANDU1AwICuDc+AMAXKSkpNTQ0EAQhKyvLfvi07QpuAAAE5gECAP6zfv16W1vbsWPHGhgYODo6Dh8+PDExsYM1zi5evFhYWNi2paSkRFpamvtJAYA0vFsACQsL79q1q4PFBaurq1++fCkqKtqVqTrQ0NAgJiZGdor/4rUwoqKiPDLDXmNjo7CwMI9MSdXU1EShUHjk4kRzczOVSu3du/fXOjQ0NPBI1G8aMWLE/fv34+PjZWVljYyMlJSUYmNjjYyMvtY/MzOzoKCgbYuQkNDDhw9/+eWXr/0Kjj8d4LUwOP58EY4/lA6ezCJXeXn5ixcvOujw8OHDrVu3hoWFdVmkDpSXl4eFhUVERJAd5L+mT5++e/duHjkGLVu2bO7cuRoaGmQHIQiC2LFjx6BBg9hTOZMuOjpaQUFh5MiRZAchCIK4dOlSZWXlypUrO+gjJyeno6PTZZFIhOPPz8Dx52tw/Pkaco4/rG4rNTXV3t6e7BT/xWQyVVVVyU7xL0lJyZqaGrJT/Fffvn2zs7PJTvFfrq6u586dIzvFf4WGhm7fvp3sFP8VERERFBREdgquGDRoEMfHxPGnAzj+fA2OP19DyvGHJy7EAQBwT3NzM9kRAIDnoAACAD7n6+tLdgQA4DkogACAz2H6HwD4HAogAAAAEDjC4eHhZGfoJCqVqqCg0LdvX7KDEARBiIuLy8nJWVhYkB3kv+Tl5W1sbHjkyU8qlWpjY8Mjj4RISEj88ssvNBqN7CAEQRCSkpImJibKyspkByEIgqBSqb169ergMVRoC8efDuD48zU4/nwNKccf3n0MHgAAAIBL8BUYAAAACBwUQAAAACBwUAABAACAwEEBBAAAAAIHBRAAAAAIHBRAAAAAIHBQAAEAAIDA6a4F0J9//mliYtKvX7+rV6+SlWHt2rUmJiaampobN27knVRbtmxJSkoiPc/ff/9taWnZq1evffv2kR5m6dKl2traurq6x48fJzFMTU3Nnj17Wl9+nqErU7ULw5sfZp7FC38c3vwnw/Hnczj+fDMMaR/mLl59niMKCwuNjIxqamry8/MNDQ2bm5u7PsPly5f79+9fW1tbUlKipaWVnp7OC6lu3rwpJiZ2+PBhFql/pY8fP+rr65eUlFRUVBgaGlZUVJAY5vr163Q6va6u7u3btzIyMtXV1aSEefTokb+/v5eXF/vl5xm6MlW7MLz5YeZZvPDH4c1/Mhx/PofjzzfDkPhh7pZXgJKTk11dXSUlJfX09FRUVLKysro+Q2lp6axZs8TFxRUVFe3t7RkMBumpKisrly1b5uXlxX5JYp74+PixY8cqKirSaLTs7GxpaWkSwwgLC4uJiYmKioqLi4uIiBAk/WX++uuv8vLy1pefZ+jKVO3C8OCHmZfxwh+HB//JcPz5Ihx/vhmGxA9ztyyA3r59q6Ghwd7W0NAoLi7u+gyenp4BAQEEQWRnZ9++fdvR0ZH0VMHBwb/99pu8vDz7JYl5GAzGq1evTE1NtbS0tm7dSqFQSAxjZ2dnYGCgpqamqakZHh4uJSVFSpgNGzYEBga2vvw8Q1emaheGBz/MvIwX/jg8+E+G488X4fjzzTAkfpi7ZQHU0tLSdpm9pqYmUmKwWKwdO3a4u7ufPXtWRkaG3FSnTp2Sl5cfOnRoawuJeWpra/Pz82/cuJGZmbl///709HQSw6SkpOTn51+7di0pKWnHjh1MJpMXPj+fZyA3FU99tmXlVgAAIABJREFUmHkcj/xxeOqfDMefr8Hx53uQ9WEW4dK4XKWmplZYWMjeZjKZampqXZ+hubl5woQJsrKyDx48YC/tS26qw4cPv3jx4vLly2/fvj19+nR9fT2JeZSUlEaMGCEnJ0cQhIODQ25uLolhkpKSJk+ebGRkZGRkZG1tfePGDV74/HyegcRUvPZh5nG88MfhtX8yHH++BsefbyLzw8yle4u46vXr13379q2vr2cwGHp6eqTc7nfixImJEyfyWioWixUSEsK+CZHEPDk5OX369CktLX3//r2WltbTp09JDLNv377Ro0e33mGXkZFBVpjLly+33vf3eYYuTtU2DM9+mHkTL/xxePafDMefdnD8+WYYEj/M3fIKkJaW1ty5c+3t7QmCOHDggJAQCV/kpaamJicnq6qqsl/u379/zJgxpKdqi8S/kpGRUUBAgJWVFYvFWr58ubGxMUEQZIXx9/fPyMgwMTEhCGLRokV0Op3EMK0+/9ch8d+L9z/MPAXHn++B4w8bjj/fROKHmcJisbg0NAAAAABvwokdAAAACBwUQAAAACBwUAABAACAwEEBBAAAAAIHBRAAAAAIHBRAAAAAIHBQAAEAAIDAQQEEAAAAAgcFEAAAAAgcFEAAAAAgcFAAAQAAgMBBAQQAAAACBwUQAAAACBwUQAAAACBwUAABAACAwEEBBAAAAAIHBRAAAAAIHBRAAAAAIHBQAPGzly9fjhs3zs/Pb+/evWRnAQAA4CEUFotFdgbglsWLF3t7e5ubmw8cOPDGjRtkxwEAAOAVImQHAC5avXq1uLh4UVGRoqIi2VkAAAB4CL4C42cSEhJnzpwJCgr6888/ubeX69evjx8/vit/8Xs0NTV9+PDhiz/6+PHjq1evuLRfAADoFlAA8bPExMS7d+/GxcWpqKiQneWnZGRkWFtbq6mpzZgxo92PDA0NRUREREVFRUVFz54929q+YsWK0NBQgiBev35tY2NjZGS0a9cu9o9WrlxZX1/fZeEBAIAHoQDiZ3FxcXl5eZ6enlOmTPme/o2NjV/c7rgnt9XV1bm5ue3fv7+goKCwsPDixYutP2KxWGVlZU1NTY2NjY2NjW5ubuz269evHzx4kL29cePG8PDwJ0+eHDlypLKy8tmzZwRBGBoadll+AADgQSiA+MeSJUssLCzq6uoIgoiIiAgICDhw4EBCQsLp06ePHTvWtmdERISOjo6+vv7ChQtbWlpu3brl4+Pj7e29Zs2attsEQaxfv15XV1dHR2fx4sXten4eoN2wI0eOPHfuHPtHlpaWN27caNfhO99XQkLCgAEDTE1Ne/ToceHCBUdHx9YfvXnzRlVVtV3/8vLyRYsWrVq1iv1SRkbm6dOnDAajrq5OVFR0/fr1K1as+M5dAwAAv0IBxD9+//13U1PT6Ojourq6nTt3/vbbb1/sdv369aioqHv37mVmZjKZzI0bNxIEkZiYGB4evnbt2rbbf//997FjxzIyMnJycu7cuRMdHd2uZ8fDenl5sb+TevnyZVlZWUtLy+f7bSsxMdG6jVmzZrHbCwoK6uvrLSwsVFRUpk2b1tzc3PorBQUF5eXlTk5OKioq/v7+tbW1BEEEBgaGhYX17NmT3WfFihV5eXkzZ878z3/+8+TJE3V1dTU1tZ/+YwMAQPeGp8D4yuzZsxcuXFhZWenh4aGurv7FPleuXKmoqPD09CQIoqSkpKqqytHRkU6nt34r1LqdkpLi4eFBo9EIgvD29k5JSdHR0Wnbs+Nho6Ojly5d2tjYeOrUKT8/v5SUlHYd7O3t247g7Ozs7Oz8+cg1NTX379+/efNmz549fX19N2/evHr1avaPZGRkQkJCgoODKysrfXx8Nm/erKenR6VSXVxcTp06xe5Do9H27NnD3p4wYcL+/fsvXLgQGRlpamq6bt26Hj16/OgfGQAA+AAKIL5ibW1dWloaERGRnp7+tT5UKnXWrFlLly4lCKKurq6hoSE7O1taWrq1Q+s2i8WiUCjsbWFhYfall7Y9Ox6WRqNZW1tfvXo1Jibm7Nmzp0+fbtchIyOj7QiJiYltLyyZmZmx529UUVEZPny4jo4OQRDu7u5nzpxp7WNiYtKnTx9hYWEFBQUfH5+4uLjMzMx79+5pa2vX1NR8+vTJ2dk5MTGR3fmff/4ZMGBAc3NzeHh4XFzcgQMHIiMjFyxY8L1/XAAA4CP4Cozf9O7d29bWVk5O7msdhg8ffuTIkZKSkvr6end391u3bn2t55AhQ2JjY6urqxsaGqKjo4cMGdLBfr84rJeX18aNG+Xk5LS1tb+5X2dn57tttM5ePWbMmMuXL798+bK6uvrkyZO2trYsFisvL6+pqSkyMtLZ2bmurq6ysvL06dP29vZxcXFMJrOgoCAiImL8+PGt1Q+LxYqMjAwODq6pqZGTk9PU1DQ1NS0tLf3OvyoAAPAZFEB85datW6WlpSkpKR085k2n00NCQmxsbHR1dXV1dUeNGvW1niNHjpw0aZK5ubmxsfGAAQN8fX072PUXh3Vxcbl//76fn98P7bcdDQ2N7du3jxgxwtjYWFNTc/78+Q0NDUZGRm/evJk7dy6dTu/Tp0/fvn319PTmzZv3tUFiYmLGjh0rLi6upaVlb2/fp0+fiIiIDvoDAAB/w1IYfMXOzm7Tpk07d+6cOHHihAkTyI4DAADAo3AFiH+cO3eORqM5ODhMnz699bZfAAAA+BwKID7R3Ny8fPnydevWEQQxcuRIJSUlOp1eUlJCdi4AAABehK/AAAAAQODgChAAAAAIHBRAAAAAIHBQAAEAAIDAQQEEAAAAAgcFEAAAAAgc3l0L7NOnT8XFxWSnAIB/UanUnj17kp2iK+D4A8BrOH784ZXH4AcPHnzt2rV2jeLi4mpqamTEAYD2mpqahIWFX758SXaQrrB8+fK9e/fKysqSHQQACII7xx9euQJ09erVdi1WVlZKSkoJCQmk5AHoAKu+oeZmOnWoDdlBulRRUZGNjaC85aampqVLly5atIjsIABAENw5/uAeIIAfVnHu0ofd0XXZz8gOAgAAnYQCCODHNH34WHUxVX6qe2lUDKu5hew4AADQGSiAAH7MxyNnpZ0H0cYMFlGQq750k+w4AADQGSiAgOuqLt1qqf5EdgrOqM99WZ//SsZlKEEQ8tPcy08nt1TVkB0KAAB+GAog4K66nBel+05+PMkXN7OzWGWHzsj5uVF6iBEEIaqhImlNLz+TTHYsAAD4YbzyFBjwJxbr4+FYhYCJ5acSpYfbi2l170kN6vNf1b8oLNkWVdKmUUi8h5zveIqwMGmxAADgx6EAAi6qunKbEBGRdrIjmlvKDp1R+W0e2Yl+Sg8D7d5n/iA7BQAAcAC+AgNuaamtK49Jkvf3ICgU6eH2LRXVn+4/IjsUAAAAQaAAAu4pP50s2b9vD10tgiAIISH5aR5lR86yGhvJzgUA/8fencdFVfV/AD+zDzMM+zYM4IILbohbiCi5lbuWUEGWlrmkoWFlpdmqaTs9j9nj405lmvvCUqbinqgIhAviggoMyDIbzL7c3x/0Q55RUJGZOwyf96s/Zg6Hez7AdP3ee889FwBwCQxshKLUx86aFaqagyfvNjIY+oIifq8u9MUCAAAgBAUQ2AqDEbxuOd0hAFoxi1ZnuFHM79GZ7iAAzgmXwAAAHJHit7Q7n/9oqpLTHQTAOdnwDJDFYsnLyysvLzeZTBKJJCIigslEvQUA8GDG0jvqY2dFIwfJf97tu2A63XEAnJCtCqDMzMw5c+aIxWKJREIIkUqlxcXF69evj4mJsdGI4FAog9GsUJkVKrOy1lKrNteoLTVqi0Zr0Wgprd6iN1jUWspgoIwmSq+njGZCiEWnI408WovB4zLYbEII04VHWCwGh83gcpkuPAaHzXThMwR8Jo/HFPAZLnym0IUpFLBcBUxXIdNVwBQJmXyeXX9ygJYg27TTffIo0dODS5OW6S5d5XfHhTCAFmarAmj+/PlpaWmhoaH1LcXFxXFxcVlZWTYaEWhhVtSY7lQay6tMldXmSrmpSm6SKcwyBWU0sTzcWB5uLDdXpkjIFAlZrkK2twdT6MLg85l8LlPowuBwGFwOg8dlcNiEEAaf19hygpTeQJlMhBCLVk/MZspgpIxGi0ZHGU0Wnd6i0VI6g0WrtShrTNIKs1pjUWssNZq6wotYLEyRK8tNyPJwY7mLmG6ubC8PpruI7enO8nRjebkzBS52/ZUBPIgm+4Kpolo0JobBYnlOmSTbsDPwq3cJzqADtChbFUBms1ksFjds8ff3pyjKRsOBPVCUqVJmuCU1FpcZS8qM0gpj6R0Gm80O8GEH+LD9vLmd2wmi+rC9PVhe7kyhoAVHZvC4dU+faMZmKaPRrKy1qGrNyhqzqsairDXJFJaiYpNcaZarzNUKymJhe3uwPN3ZPp4sb0+2lzvL14vt68X29mQKURuBvVEmszxlt9ersXXHA8LovjV/HK85/JdoZDTd0QCciq0KoMTExH79+k2aNCkoKIjBYEil0n379iUmJtpoOLAJijKWVeqv3TJcu2W4WWIoKmEKXTghgdwQMT88TDR2KEfs5/glAoPDYft4Eh/PxjpQeoOpSm6WKU0yublKYbgtNedcMlXKTJUywmCwfb3Yfl5sP2+2nw/bz5vt783x96mrxgBsQXMmz1hWUb1uW32LRaM17fkTBRBAy7JVATR37twJEyZkZGRIpVJCSEhISHp6enBwcGP9i4qKZDJZwxa1Wu3p2eg/WmAjlNGkv3pTf/m6ruC6vrCIKRTwOrfjhrbzGNCL2yGY6dqS53UcBIPH5Uj8ORL/e79kUWtMFTJTpcxUUWW6U6n7u8B4p8p0p5op5LP9fTkBvmyxDyfAlx3gyxH7tuxJL2izBJERQT98bNXIwFQ2gJZmw7vAAgICpk6dyufzlUrl2bNnjU0uAfzpp59euHChYcvNmzcFAvyLYhcWi/5GsTb3su5CoeH6bU5wAC8sVPTUYJ/El1nuIrrD0YkpFHA7CLgdgqzazTKlsbzSdKfKWFapOfO3sbzSVFbJYLPZYl+OxJ8j9uME+rED/Thiv7rpTUCvNWvWzJo1i+4UD4vBYrL9fehOAeD8bLV3/uWXX+bMmePm5vbWW2+tWbOmU6dO58+fX7Zs2WuvvXbf/ps2bbJqiYyM9PX1tVE8IIRYatSanEvac/na/CtsLw+X3mHuE0fwuoUyXfh0R3N0LC93lpc76d6pYaNZUWOU3jGVVRrL7tQeO2OUVpgqqlle7hxJAFfiz5b4cyUBnKAApzyL5mjS09MrKyvr33722Wc8Ho8QMm3aNPpCAYBjsVUB9OGHH+bl5bVr1y4yMvKdd96ZOXNmVVVVVFRUYwUQ2I2pSq45k6c5nWcoKuaHdxX07en1aizL053uXK0ey0PE8hA1rIoos8VUUW0sKTeW3tFfKao9dMpYeofBYXOCxBxJADc4gBMUwAkJZLm50hjbKV28eHHJkiXx8fFeXl6EEJ1Ol5ubSxovgMaNG/fXX381bKmtrY2Kilq4cKEd0gIALWxVALHZ7I4dOxJCXnnllZEjRxJCPDw8bDQWPAyzQqU+laM+cc5UVunSv6fbpBEu4V0ZHA7duZwZg8XkiH05Yl8yoFd9o1mmNJSUG0vKDbek6pPZhuIyBoPJaRfICQrghgRyQsTcYDGmEz2mhQsXDho0aNGiRVOnTh0xYsTBgweTk5Ob6L99+3a9Xt+wZeTIka6uKEwBnJmtCqCYmJhJkya9++67dXd+Xbp0acWKFYMGDbLRcNAYymjUnPm79kiWvvCmoH8vj+fG8MPDGCwsKEIblpe7i5e7S3jX+hazQmUsLjPcLjPcKK49ctpYXM4U8DnBYm47CSdEzA2RcIIC7DCXSLnnT2NZhc+cKbYeyD6io6NTU1PnzZuXmppqVdzcSyAQWM04ZLPZDAbDlgEBmklXcL3yu42S5A8c/yZcB2erveratWt3795tsfyzsO+NGzf69+//+uuv22g4uJfhtrT2z5O1J87xOoa4Do30WziTwcX5npZBmc26C1ddeoe1yNbqVozk97pbEpkqqo3FZYbbUm3uZdXeQ8bySrafN7ddILedhBMSyG0nYft6tcjQd0esVqj2HiIctr7gBi+sY8tunC5ubm4pKSlbt24tKCigOwtAC6Eo2YYdLFeBYnu61yuxdKdp3WxVADGZzNjYu3+b8ePH22ggsEKZzZrTuaqMo6YKmWhEVODX77MbXwIHmqcm45gsZVfgl+9yOza6ssPjYPt5s/28Xfr1rHtLmc3GknLjbanhlrTmj+OGW6WU3sBtF8RtH8hpJ+G2D+IGix+zupX/vEc0JoYj9qtevy3wq/eIE538iI+Pj4+PpzsFQMuoOfwXg8Px/zCxNGmZaMQgTrD4wd8DjcA9us7DUqupOXBc9fsxTqC/+4QRLv174VKXLZhrapW7/nCfMFy2cUfAZ0l2qBUYLBa3nYTbTiIccjeD8Wap4Wap/vL1moxjRukdtp83t30Qt72E2z6I2z6I5fEI6xfoC4v0Bdc9nhtDGYxMPq/2SJbrsIE2+UmcjkWnL1/8rU/iS9yOIXRnAedn0eoUv6X5vTebKRK6P/u0bOMO/4/m0R2qFUMB5AxMVXJV6uHaI1mCJ3r7fzCX205CdyJnpvg1VRgzwPPlZ8re/1p9KkcY3df+GVgiV1avrvVXzSiz2VhcbrhVYrhZqtz7p6GohMFmczsEcTsEcTsEc9sHcQJ8Gi3UKEq2YYfnS5Oq12w1Vcp8k16p+HqdILI3HpH2MJS7DlgMhuoNO8RLFzjTaTNwTMrtGYI+PXihIYQQ0Zgnaw6d0mZfqD9VDI8KBVDrZqqoVu78Q30mTzQ8SvLdBywvx7ub3WKpTN7o8fxY5zhVaywu05zJk3y/hDAYXtPjKpM3Cvr3pP3JGAwWi9tewm0vIU/+02KqkhuKSgw3S9Qnzsl/3mOpUXPbSbgdg7jtg7gdgjnBYgb7n+fO1h7JMpbe0WZfNErvMEVCxbZ0Boet3P2n55SJtP08rYTpTlXtnyfF37xf+fVa9Yls4ZD+dCcCZ2a6U6VKPyIcMkD+y966FpaHmyxlt6RPdzwot3lQALVWpiq5Ynu69ky+aPSQoJUfO+zyejWH/9LmXzHXagI+doZTtbJNOz1eGMcUCQkhvK4deV07KPce9Hh+LN25rLF9PNk+noL/v/3eotYYikoMRSW6C4Wq1ExjeSVH4s/tEMztEETMZtHwgTWZWS59ezCFLurMLOGwSOeoVm1N9tNutwnD2d4eXtPjKr9dL3ginPZSGJwYg8vxiP+f2bQuvcMYbDZOPTYbCqDWx6ysUe78o/b4WbdRMZJVHzvypQqLVqfYlu7/wdzq1Vs0WXmCyN50J3os+sIibV6BUVqh3HuwroXS67XZF91jR9U9uNthMYUCfs8u/J5d6t5SBqPhttRQVGy4UWwoKtHfLGFyOYRBhNH9mK4Cc5XcNWYAvYEdn+5CoaGoxDfpFUIIr0sHXrdOjlkKg9Ngebq7P/MU3SmcCgqg1oQyGFWpmar9h4VPDpD860PHX0FYsS1d0K8nr3N7r1diq1b/6tK3e6teepHXuX3Qqk+sGhkctoNXP/dicDm8Tu14ndoRQswyRelbX/i++bJZUcPgcjwmjyp9c6mu4Do/LJTumA6MomQbd3KCAlRpR+oaWCKhcs+frsOjcNMlQGuBAqjVUJ88L/95N69ze/GKd9gBreBZicbySvWRM4HJiwkh/F5duO0CVamZ7s8+TXeux8BgON9TKlW/H6P0+up12+tbLGqtan8mCqCmUJTgiXDKaLLUauoaGDyu29ihtGYCgEeDAqgVMNyWyjbssKg1PkmvtKJ/luQbd/K6dtBdvl73lte5vXLXAdehkXjumEPxiB0tGmG9RHvdJCdoFJPp8cI4ukMAwGNBAeTQKL1BsS299kiWx/NjRU9Ft6ap/hTFFAkpvUFz6nx9m0tEN7OiBgWQQ2HwuM53WgsA4IFQADkubc6l6rW/8buFBiZ/4PjTfawxGD6JL9MdAgAA4P5QADkii1oj27BTV3Dd+/UXGz41EwAAAFoECiCHo82+UPXfrcLI3pLvFmNZEQAAAFtwlAJo//79ZWVlDVsqKiqEwjYwE5OiLBotUygghFB6gyxllzb3sm/SNH73znQnA4AHqzlwQhgzgMnn0R0EAB6NoxRAV69evXLlSsMWjUZjNBrpymM3NX8cV+4+IPn3R8bS8srvU3id2wV+u4jpwqc7FwA8mDbvcvXa30xVMs8XH/TkEIq6s2K125gnXfp0t0s0AHgARymA3nrrLauW3Nxcd3cnv13IUqtRbM/gSPwrvl5ruFHs9drztDxZEwCagTJb5Cm7vV9PUPyyz3VYFEfs20Rn9clsw81S2YbtgclL6h/EBgA0aj23VTsjxW9pggHhDB5Xl3/F750ZqH4AWpGa348x3UWiEYPcJgyrfz7lfVEGo3zzPt8Fr3IkAaq0TLslBIAmoACijbGkvPbEOd3fBRyxn9vkUaqMo3QnAoCHZanVKHf/4fVqLCHEbeJI4+1Sbd7lxjordx/ghXXkdwv1eiVWteegWa60Y1IAuD8UQLSp/H4jMZjcnx/rNT3OY/LT+uu3dZeu0h0KwElkZ2fv3btXqbxbaqSmprbg9uVb9gsH9eOGBBJCGGyW50uTZBt3UmbLvT1N1Yqa34/XTRJiB/i4Dh8o39qSSQCgeRxlDlDbQlGV320w3CwVRPbWX7mhv3KDEMJ0Fch/2iP+YiHd4QBavRUrVqxevXrw4MFJSUm7du3q06cPIWTmzJlWt5o2G2U2q4+dtWh1qvQjd1sZDMP1W7wuHaw6y3/azZb4aXP/OT/EdBcpUw+7jYrhdgxukTAA0DwogOzNotNX/TvFLFN6Tpvc8NZZbscQlocbjcEAnMaqVavy8vK8vb3z8/NjY2Ozs7NFIlET/bVarU6na9hiMpkoimqsP4PFCvn5m4cMw5H4MwV8w43b9S2i4VGUyfSQ3w4ANoICyK5MVfKKFat5XTr4vv0ag4U7QQBsQiQS1a0i1qtXrzlz5iQlJa1fv76J/pMnT87KymrYUlNT01JhPJ4f21KbgmYwSu9wxH6EwaA7CDgczAGyH8PN0vIPvnUdFuk9Ox7VD4DtJCQkREdHr1mzhhCSlJQkk8ni4+M1Gk1j/TMyMmT/q2/fvmKx2I6RwSbMMoV04Zc1f56kOwg4IpwBshNd/pXK7zd5TY8TRvdr2G4sr9ScynGf/DRdwQCcz0cffRQdHV1dXU0IYTAYO3bs2LZtm6enJ925wN7kv+wV9O2h2JYmjO7HFLrQHQccCwoge1CfOCfbtNN34Qx+WKjVl2Trt+vyr/C7d+KFdaQlG4BTGjFiRP1rFouVkJCQkJBAYx6wP31hke7iVcm/P5Kl7FJsT/d6JZbuROBYcAnM5moOnJD9tMf/w8R7qx9t9gVTRbX36y9Wr99GGp9xCQAAj4aiZBt2eL40icHjer44UX38nLG4Ze4BBKeBM0C2pdz5R23mafGyBWw/b6svUWazLGW316uxLhHdag+dqs087To8ipaQALSbM2fOfdv/85//2DkJOIfaI1mEyRAO7k8IYboK3J99WrZxh/9H8+jOBQ7kPgVQbm5ubm7ulClTysvLg4OxUkXzybfs157ND1i24L73t9ekHWGLfeuejOg1/bk7y/8jGBjBFOAqNbRF48ePpzsCOBGKkm9JZbrw7iz94Z8Gk0V36aru0jV+9070RgPHYV0Abdq0acWKFVqtNi4ubsiQIUuWLJkxY0bzNm2xWPLy8srLy00mk0QiiYiIYDLb0BU3xdZUbfZF/0/ns0Su937VrKxR7PzDN+kV050qQghTwOd1DFbu+sPzpWfsnhSAfuPGjbNqsVgsS5Ysubcd4MEYDL93XrPoDA3b3MloLD4JDVkXQMnJyadPn546daqrq2tubm7//v2bVwBlZmbOmTNHLBZLJBJCiFQqLS4uXr9+fUxMTAukdniylF26i9cCPpnPdBXct4P+6k2mq6B63baGjRa94b6dAdqIr7766vvvv5fL5f7+/hUVFcOGDaM7EbRW9y7JDWDFugBSq9V8Pr/utUAg4HA4zdvu/Pnz09LSQkPvTvstLi6Oi4uzWm3MwWlzLyt3Hwj4ZP4jLaIl37xPd/FawMfzmrjrUtC/l6B/r5bICOA8fvjhh8uXLy9evHju3LkqlWrjxo10JwIAe9BmXzBVyUWjhthzUOtrUlOmTJk4cWJRUdGGDRtGjhw5ZcqU5m3XbDZbLSPm7+/fxNLyDogym2UbdpgqqmsO//Xw36XYlq7NvhDw4RtYcwLgUfF4PDc3tz59+mRlZUVFReXl5dGdCABsjtIbqtf+Jv91v6lSZs9xrc8Affrpp2lpaSdOnCgtLV2+fPngwYObt93ExMR+/fpNmjQpKCiIwWBIpdJ9+/YlJiY+dmD7qZuk7Bk//s7nPwqj+jzM9GTl3oPqU+cDPn2TKRLaISGAk4mJiZk+ffqiRYuee+65O3fuuLu7050IAGxOuftPXrdOnEA/+c97fN+abrdxrQug6urqgQMHDhw4sP6tt7f1/dsPY+7cuRMmTMjIyJBKpYSQkJCQ9PT0Ju4pW7JkSWFhYcOWa9euseh7XoRZWaPc82fAsgWcQH9Bv57KHb97Tn226W+pOXSq5sAJ8dIFLPemHrsIAI1Zt25dYWFhly5dvvvuu+zsbFwCA3B6pip5zR/HxV+/x3JzLU1aprt0ld+9s32Gti6ARo8eXfdCqVRevXq1d+/eubm5zdt0RUWFv7//Cy+8UH8Yl5qa2tjNrqNHj+7du3fDlpycHIHg/jOI7UCxJVU4NJIT6E8I8XhxgnTB564jozmBfo3115z9W7E1LeDT+SwvHLMCNJNMJvPx8amuro6IiIiIiKA7DgDYnDzySz8SAAAgAElEQVRll2jck2wfT0KI55RJsg07A796l9jlnnHrAujs2bP1r/Pz86dOndq87a5YsWL16tWDBw9OSkratWtXnz59CCEzZ84sK7v/Wpz3Xmv75ptv6qdj25nhZmlN5l+eCRPrH6HHaRco/2m33/uz79tfd6GwevUW/w/fqCuYAKB5WvAADAAcn+7SNW1egWjcUMON24QQdoAPZTDUHP5LNDLaDqM3tRJ0r169jEZj87a7atWqvLw8b2/v/Pz82NjY7Oxskaj1XBiiKNHwKNOdyvoGToDvfRczJIQYi8sqv9/o+/Z0bvsge+UDcE4tdQAGAK2C4VYpO8BXtmFnfQuDzzdJK+wzunUBtGDBgvrXBQUFXbt2bd52RSKRUCgkhPTq1WvOnDlJSUnr169vdko743YI8p79UM9NNMuVd1as9nz5WbtdswRoIx7nAAwAWgW3MU+6jXmSrtGtC6CGl6LGjRvX7IXIEhISoqOjZ8+ePWvWrKSkpMmTJ8fHx2s0muYndTwWnf7O8v+Inh7s+uQTdGcBcAYtdQBGG4oyVcnZvl505wCAB7tbANUtURgU9D/Xcc6dOxcZGdmM7X700UfR0dHV1dWEEAaDsWPHjm3btnl6ej5eWkdCUVXfb+J1auf+zFN0RwFwEi11AEaX2szTsg07JP/+kOXlQXcWAHiAuwXQ999/TwiRSqXHjx8PDw93dXXNysp69tlnt23b1vi3N2XEiBH1r1ksVkJCQkLCQ11XahXkm/da1Brfd5r5oDQAaKhlD8BoYdHq5FtS+b3D5Jv3+czD7CUAR3e3ANqyZQsh5Kmnnrp8+XLdmeeioqJZs2bRFs2B1R49o/4rV/zFOww2bSsVATiTFj8Asz/l9gyXvt29pz9X+uZSXcF1fljog78HAAghhBhulpiVtS69w+w5qPWt9rdu3erSpUvd6/bt29ctYwgN6QuL5D/v9l/0+n0f8w4AzbBly5YtW7ZwudzLly/n5uaeOHGisLBQLpfTnethmcqrajOzPOPHM3hcj4QJsg07yWM/+cdSo9ac/btF4gE4NIul6odfqv6dYlFr7TmsdQEUFRUVFxe3a9eunTt3PvfccwMGDLBnGsdnVtRUfrfBZ84UTlAA3VkAnE0LHoBZLJacnJyMjIz9+/efP3/eYrG0UMb7k23a6fbMSJanOyHENWYAk8upPXrmMbcp37yv8tv1RnvdEgxAl5qDp5h8nmBAuGJHhj3Htb4LbM2aNatXr/7111+ZTOaTTz6JS2ANUWZL5XcbXEcOcunXk+4sAE6o7gBsypQpFEVt2bKl2QdgmZmZc+bMEYvFEomEECKVSouLi9evXx8TE9Oief+hv3JDcy7fotHq8i7XtZhVtfJf97k++QRhMJq3TcPNEs25fPdnnpal7PJf9HrLhbVGmc0M+h46BGDR6hTbM/wXzWb5eEnfXCYaHsUJFj/421rC3QJo3rx5S5cu/fDDDwkhdQ9yLywsfOedd1auXGmfKI5PnrKTKeB7xI6mOwiAc2qpA7D58+enpaWFht6dhVNcXBwXF1c31fpeRUVFMtn/PIZapVKxWKzs7Oy6t2FhYXULmykUiuvXr9d3q2vnhATyF0y7UVxc395l1BCRuxthMO7bv7HtNGyvXr3FpV+PPmOHqJf8S5tzSd8hsHnbabo9WKE17j8S+M0iZW2NLbaPdrQ/sF21/3CnXl24HUMIIWTUoCPLv/ea8fy9/fPz800mE2lRdwugkSNH8vn8+qXowYr6ZLY255L4y3ebfUgHAI1p2QMws9lct5F6/v7+VOOTcj799NMLFy40bCkqKiouLp49+5+n33z22Wdjx44lhKSlpSUnJ9d3q2tnuvAP3Sy0bo+MaKz/A9sttRqzQsXJCfisnf/Qac/KNu3M7N/++3/961G384B2iprv32VkWK+ajKNpqrLmbwftaG9uO2U0Gksrlv9nZd2T0o8alV/v+oV1NI0p4Fv1//LLL1UqFWlRDKudQkFBwcmTJ1977bUlS5YcOXLkiy++uPcpXfYRGRnp6+ubmppKy+j1LGotU+hiLK8s/+A7/yVvcDvgeRfQRpWUlERFRRU3OM/Rgvbu3Ttq1KhDhw5ZtY8bN64ZW/vxxx9Xrlw5adKkoKAgBoMhlUr37duXmJhYX9A8kC32P4Ybt9mB/kw+r+lulMFYmrTMJ/Gl+vXl73z+o0vvbm7jW3hVJMXWVFNFtfvzY8sXfRuYvLixp/0A2I5ia6pix+9WjcKoPr5vv2bVaIv9j/UcoDfffHP69Olnzpw5ceLEF1988fnnn2dk2HVSkkOpOXBCufP3wG8XVX673iN+PKofABuZNGkSISQ0NNTqAKx5W5s7d+6ECRMyMjLqplGHhISkp6cHBwe3ZOJHZFbWlH+y0nX4QK9XYpvuqTmXb6qUVXy5tr6FMpmMZRUtWwCZquQ1f5wQf/0e28dTOPQJxW9pD/n8H4AW5BE/3iN+PF2jWxdA5eXlL7zwwvvvv//SSy8NHjz4xo0btMRyBBa1RrEtjSP2u/PZDxxJgOgpezycFqAta8EDsODgYIe6h0Px635BZG/1sbOipwZzJP5N9BQO6nvvaigMDqdl88hTdonGDWX7eBJCPJ4fW/rmUtfrt3mhIS07CoAjs74NvmvXrm+//fbmzZvHjx+fnJzs5tZ2T4oqfksXRPYRxgww3Cxxn/w03XEAnF/dAdiuXbuc7ADMUFSiyb7g9Wqs++SnZesfvLQjUyiw+o/BbckCSFdwXX0616Kqlf+yV/7LXuXOP9heHvKf97TgEACOz7oAWrVqlUQi2bNnT0BAQHV1dd3y0G2QsaRcfeKc69PRit/ShCMGKbe33euAAHbjrAdgso07PBLGMwUuojFDTXKVNvvCg7/HllhuIs8pE1le7kxXQd1/goERrkOw6hu0LdaXwHx9fYcPH56bmxseHj579mx6r5rTSJayyz12lHzTTrfxw93GDi1NWqa7dLV+TiIA2MKqVat+/vlnJzsAU5/Mtmh1omEDCSEMFtPr1TjZmq2B4WEMjvXu1244gX54ijOA9f+BmzZtWrFihVarjYuLGzJkyJIlS2bMsMfzPmNjY3Nzcxu2lJaW9uxJz3qD2rwCbe5li7LGpFBxAnxlG3cwXQXylN3iL9+lJQ9AG+GUB2CKHb+bZcqSeZ/Vt5gqqtV/5bjG4IwLAJ2sC6Dk5OTTp09PnTrV1dU1Nze3f//+9imA1q5dq1AoGrZMnjzZ37+pqYK2ww0Wu08eVZN+xC12FEsoIIRwO4aw3EW0hAFoO+g6ALOpgI/nU3q9VSPLy52WMABQz7oAUqvVfD6/7rVAIOC09K0HjfHy8vLy8mrYwuPxGDQtOcgUCbVn//aeFS8c0p+WAABtE10HYDbF8hARgsMnAIdjPQl6ypQpEydOLCoq2rBhw8iRI6dMmUJLLHoptqayA3xQ/QDYGV0HYADQBlmfAfr000/T0tJOnDhRWlq6fPlyupaBppH+yo3ao2cCv11EdxCANqfuAKysrGzDhg2bNm1qmwdgAGAf97kNYdy4cfXLz9fU1IhEbejkLWUwVq36xXt2PGb8ANgfDsAAwG7uFkDnz59ftGiRTCYbO3bsvHnzFi1adPv27by8vPLychrz2ZnitzRuxxDBgHC6gwC0UW35AAwA7OnuHKBXX301KioqOTlZKpWGh4dLJJJFixadPn2axnB2ZigqqT16xuvVBzypBwBa3Pnz50eNGjVgwICPP/64qqpq5syZo0aN6twZK28BgK3cPQMkl8s/+eQTQki3bt327t1b97rtoMyWqh9/8XplMi5+Adjfq6+++uyzz44cOTIlJSU8PHzWrFlTpkxp37493bkAwGndLYBYLFbdC29vb6FQSFMe2ih3H2B7ewoH484vABq08QMwALA/69vg2yZjeaUqLdNrehzdQQDaqDZ+AAYA9nf3DNCtW7fq5xuq1er61zU1NTTksi/Z2t88Ykex/bzpDgIAAAD2cLcA0mq1NOawD/kve4SD+nE7/s8DhtTHz5kVNaIxT9KVCgDa8gEYANDibgHE4/FozGEHuguFqoxjukvXxJ+/Tf7/IRsWjVb+yx7ft19j/P8ZeACwvxY/AMvOzi4pKRk6dKi7+z9P3UpNTR0/fnzLjgIArZcN5wBZLJacnJyMjIz9+/efP3/eYrHYbqyHSSPbuNN33lRCkdpjZ+ub5b/ud+kfzuvSgcZoAMBrRPO2tmLFismTJ2/bti0iIiInJ6eucebMmS2XFwBavfusBN0iMjMz58yZIxaLJRIJIUQqlRYXF69fvz4mJsZGIzat5s+TTJFQMDCC5eVe+e16QWRvJp9nuC3VnM6RJC+hJRIA2MiqVavy8vK8vb3z8/NjY2Ozs7OxoCIAWLFVATR//vy0tLTQ0ND6luLi4ri4uKysLBuN2ASLWqPYnuH/4RuEEF6XDvwenVV7/vSIHy9bt90jfjxThFtOAJyKSCSqu5WsV69ec+bMSUpKWr9+Pd2hAMCx2KoAMpvNYrG4YYu/vz9FUTYarmmKbenc9hJiNhtu3CaECKL6Vn6/kSkSWnQ60fAoWiIBgO0kJCRER0fPnj171qxZSUlJkydPjo+P12g0jfU/fvy41TN/ZDIZThoBODdbFUCJiYn9+vWbNGlSUFAQg8GQSqX79u1LTEy00XBNM8tVZpW6avXW+ha22F+583e/918nTKyEBOBsPvroo+jo6OrqakIIg8HYsWPHtm3bPD09G+t/4MCBK1euNGyRyWR+fn42DwoA9GHY7qxMcXFxRkaGVColhIjF4rFjxwYHBzfWWSaTKRSKhi2TJ0+WSCRpaWm2yCbfst9crfBJfNkWGwdwSiUlJVFRUcXFxXQHsYfIyEhfX9/U1FS6gwAAIbbZ/9jqDBAhJDg4eNasWQ/ZeebMmbm5uQ1bSktL2WybxDNVK2oPnBR//Z4tNg4Ajmbo0KFHjhyhOwVAizFVyti+XnSnaPVsWAA9kp07d1q11B2B2WIs+U+7ReOeZPs0ej4cAJyJ2WymO4KjsljMChXLy4PuHPAI1Cezq374RfKvJXh6wWOyVQG0evXqo0eP3tu+ZcsWG434MCijyVBUrL9yw2fuFBpjAIA9TZ06le4IDkqZmqnae1Cy8iOmwIXuLPBQKL1B/vMel77d5T/v8X37NbrjtG62KoBefvnlrKwslUo1d+5cGw3xqCizWbpwBWUyeSRMYPC4dMcBADvBEoj3ZVbUqHb/ye0Uotjxu9fUZ+mOAw9FufcgLyzU540ppW8u0+YVuPQOoztRK2arAkgoFC5YsGDr1q0jRoyw0RCPqibjGKU3mmRKQf+edGcBAKCZ/Nd9riMGuk0YIV3wuWjkIE6gP92J4AFM1YqajGPir95lcDieLz8j27Qz8JtFDBbuZW4mG/7iwsPDly9fbrvtPxJLrUa5+wAhlEt4mHLnH3THAQCgk6GoRJtzyX3yKJa7yO2Zp+Qpu+lOBA8m/2m3aExM3fRnYVQflptr7eFTdIdqxRxlErStyTfv4waLCZvlk/iSdMHnrjjcAYA2TLZxh2f8+LqpP25jh0oPntLmXHLp053uXNAow80S9anzvK4d9AXX61rMyhrFllTX4VF4mHfztIkCyFhcps7KZTCZ/ovn1B/u+C16ne5cAAA00JzN112+znRz1eZe+qeJxZRv3osCyJFxxH7+H1ovJsxgs1D9NFubKIBkG3cwWCymUKA+dV596jxlMmvOX9T+fcUlvCvd0QAA7I0X1sH3rekNWwSEsNxd6coDD4PB4+LfrJbVJgogl8hw3ZUi1ycHMF0FdS2eUyayfbD0BQC0RSyRqzCqD90pAGjWJgogc4XcdehAz5eeoTsIAAAAOATnv33OrFDVHDrlPvlpuoMAAACAo3D+M0DKXQdchw5ke+OCF9Bmx44dX3zxBd0p7u/VV19944036E4BAGBvTl4AmWVK9fFzgckf0B0E2rQLFy7069fv4Z8NbDd79+7NycmhOwUAAA2cvABSbM9wHTmI5SGiOwi0dWKxuF+/fnSnsJabmyuVSulOAQBAA2eeA2SqkqtP57hPdJRncQAAAICDcOYCSLn7gOipaKZISHcQAAAAcCy2KoDMZnP965ycnM2bN1+5csUWA1Fmc8WX/zWWV1oHkCvVp867jRtmi0EBAACgVbNVARQQEFD3Ijk5efz48b///vuoUaN++umnFh+oJuOo7vJ1ecquf95TlFF6h9Sd/hkexXLH7B+Atm7NmjV0RwAAh2PzSdDffvvtuXPnxGJxZWVldHT01KlT79vt+PHj5eXlDVtkMplI9IDyxayqVe46EPBZUuU367S5l10iutUcOiXbuDNg2Vu1R89KvsfNXwBtUXp6emXl3bPCn332GY/HI4RMmzaNvlAA4FhsXgD5+PiIxeK6F6zGn9l24MABq2tkMpnMz8+v6Y0rtuwXPvkENyTQ8+VnZZt2BixNUmxNE/TtUfX9RtdhkSxP9xb5EQCgdbl48eKSJUvi4+O9vLwIITqdLjc3l6AAAoAGbHUJjMPh9OjRY8yYMSqVKiUlhaKoGTNmDB48uLH+S5cu3fa/OnXq5Onp2cQQhpslmrP57pOeKv/4XxyxL9vHs/LbDYIBvbxeiTWWVfC7dbbBjwVgK6Y7VZXJG4nF0ozvzcrKYjAYqamp9S1Lly5lMBh6vb7urVqtXr16dcsEbQ0WLlx4+PDhoqKi8ePHJycni8Xi5OTk5ORkunMBgAOx1RkgqVSq1Wpv3bp18+bNgIAAiqI6duz41ltvteAQsg07+d1CZSk7DTdLKr/fxOseWpt+zHvG8zWHTrl076zYkS54ohdhMFpwRADbkaXs0uUX1hz6S/RUdDO+3cfHZ8eOHePHj697m5qaWn8FOT8///vvv9doNK+//nqLxXV40dHRqamp8+bNS01NrS8EG7N58+aSkpKGLVKp1MXFxZYBAYBmNrwN3sXFJSwsbPTo0REREUwm84MPPmjBHQplNrN9PCijSfNXLq9ze1NFlfrYWU67QMpoqjlw3GvmC0w+rzbzdEsNB2BTuvxCwy2p/0eJit/SLGptM7bQt2/fnJwco9FICLlx44avr69Q+M8CEFu2bFEoFC0Zt5Vwc3NLSUmJjIwMDQ1tuqdWq5X/Lw8Pjw4dOtgnJ9iU/Je9pio53SnAEbXWlaAZLJbP/GlVP/zsNnG454sTlWmZ8o07mcHiyu82EEJk67eZVbXyLamuwwbiJBA4DkpvMCtU1q0WS/Xare7PjGS5ufJ7dpb9tMtj8ijrPkwm28eziQ8zk8kcNmzY4cOHR40atWvXrtjY2Ozs7LovLV++/NChQ+vWrWvJn6T1iI+Pj4+Pb7rPjBkzrFosFouvr6/NQoGdaLLyVKmHTZXVvgum050FHE5rLYAIIYYbt7V5BZJ/f0gIEY0YpDlxniUJMEorPJ4fyxb7EUKYPA6qH3AoqrTMmkN/WTVaNFpKb1DuOajcc5CyUJZqmTavgHHPHQO+86bywjo2sfG4uLiUlJRRo0bt3bt3z549ixcvbuH0rdbQoUOPHDlCdwqwN8pklv+yx/ft12Qbd+ouXeV3x8RQ+B+ttgCiqOp12/k9O2tzL9c1CPr3VOz8nde5vWjUEHqjATTGffIo9/89u2Op1ZTM+9TzxYmckMC6ltrDf1k0Wv/Fcx5144MGDZozZ87t27d5PJ63t3fLJHYKDddlhbZDte8QJyRQMCCcMphkG3YGfvUuYTrzww/gUbXWAoiyUByJP6XTa06d//8misFmC4dF0ZoL4NGYquS8DkHa8xe05y/UNzK4XGKxPOrOmslkxsTEvPnmm5MnT27pmK1bY8uPgRMzK2pU+w8HLH+bECKM7lvzx7GazNOiEYPozgUOpLUWQAwW0+eNlxq2aHMuGaUVoicH0BUJoBm47SX+H81rqa3FxcWNGDFi1apVLbVB5zBz5ky6I4C9yTfv5ffsTGm1hhu3CSHCmAGKLfuFUX2YAtzcB/9orQXQvVT7D7lNGolJP9AGRUZGZmRkEEKGDRtm+f+VhMrKyuo7jBgxYsSIEfSEA7A/ijLdqbbo9FWrt9a3sbw8jbfLmp5IB22KkxRAhttSQ8kdv+h+dAcBAAC6MRgBn71JdwhwdE4yI0y195Db2CcZ7EYftQEAAABQzxkKILNMqTmXj9ltAAAA8JCcoQBS/XHMNWYAUySkOwgAAAC0Dq2+AKKMptpDf4nGxNAdBAAAAFqNVl8AqU+c43YM5gT60x0EAAAAWo1WXwCpMo65jXmS7hQAAADQmtjwNniLxZKXl1deXm4ymSQSSd0z4Vt2CN2la5RW6xLRrWU3C9DiTp48+eWXX9KdwtrZs2c9PDzoTgEAQANbFUCZmZlz5swRi8USiYQQIpVKi4uL169fHxPTkpN1ajKOisYOxeKH4ODGjBmj0+nkcjndQax16tRp+PDhdKcAAKCBrQqg+fPnp6WlhYaG1rcUFxfHxcVlZWXdt/9333135cqVhi03b97k8XhNDGGWK7X5V7znTmmRwAC2ExkZGRkZSXcKAAC4y1YFkNlsFovFDVv8/f0pimqsf+fOnV1dXRu2XLhwoUePHk0MwXDh+709g+nCf8yoAAAA0NbYqgBKTEzs16/fpEmTgoKCGAyGVCrdt29fYmJiY/0nTJhg1XLt2jVfX98mhmDyefxeXVomLgAAALQltiqA5s6dO2HChIyMDKlUSggJCQlJT08PDg620XAAAAAAD8+Gd4EFBwfPmjXLdtsHAAAAaB6Hfhr8tWvXDh482NhXTSbTqVOnrGYO0Ugul3t6etKd4h+OFsbDw4PhGDfrqVQqoVDIYjnEc3M1Gg2LxWp6sr/dGAwGiUTSxGnaqqoqe+ahHfY/zeZoYbD/uS/sfxhNTEym19atW9evX99Eh8rKykuXLolEIrtFagJFUUql0nGWVHGoHZBSqXR1dXWQ/+dra2t5PB6Hw6E7CCGEaDQaJpPJ5zvERH69Xs/n8/v06dNEn44dO/73v/+1WyQaYf/zOLD/aQz2P42hZ/9DtVrHjx8fPHgw3Sn+IZVKxWIx3SnuEggEarWa7hT/6NmzZ35+Pt0p/jFp0qQ9e/bQneIfCxYs+O677+hO8Y+VK1cmJibSnaLVwP6nCdj/NAb7n8bQsv9p9Y/CAAAAAHhUKIAAAACgzUEBBAAAAG0OCiAAAABoc1pxAcTn8x1k+johhMvlCgQCulPcJRKJHOSuB+JgfylHC+Pi4kJ3in+4uLg4zm/G8TnUBwn7nyY41F/K0cK08f2P494G/zD0er2DrGFAEKZxCNMYo9HIZDId5J8Ki8ViNpsd5AbdVsGhPksI0xiEaQz2P627AAIAAABohlZ8CQwAAACgeVAAAQAAQJuDAggAAADaHBRAAAAA0OagAAIAAIA2BwUQAAAAtDmttQD68ccfu3fvHh4enpmZSVeGpUuXdu/ePTg4eMWKFY6T6uuvv05PT6c9zx9//DFgwIB27dqtWbOG9jDvv/9+hw4dQkNDf/nlFxrDqNXq1atX17+9N4M9U1mFccwPs8NyhF+OY/7JsP+5F/Y/DwxD24fZzk+fbxG3b98OCwtTq9VXr17t2rWr2Wy2f4aDBw/27dtXq9VWVlaGhIRkZ2c7QqoTJ05wudxNmzZRtP6W5HJ5586dKysrlUpl165dlUoljWGOHj0aERGh0+nKysrc3d1ra2tpCfP3339Pnz49Pj6+7u29GeyZyiqMY36YHZYj/HIc80+G/c+9sP95YBgaP8yt8gxQRkbGpEmTBAJBp06dAgIC8vLy7J+hurp69uzZfD7fx8dn8ODBxcXFtKdSqVSLFi2Kj4+ve0tjnr17906cONHHx8fNzS0/P18kEtEYhsVicblcDofD5/PZbDah6TezZcsWhUJR//beDPZMZRXGAT/MjswRfjkO+CfD/ue+sP95YBgaP8ytsgAqKysLCgqqex0UFFReXm7/DM8///ysWbMIIfn5+adOnRo6dCjtqebNm/fxxx97eXnVvaUxT3Fx8c2bN3v16hUSEvLNN98wGAwaw0RHR3fp0iUwMDA4OPiTTz4RCoW0hFm+fPncuXPr396bwZ6prMI44IfZkTnCL8cB/2TY/9wX9j8PDEPjh7lVFkAWi4XBYNS/NZlMtMSgKCo5OTk2NnbXrl3u7u70pvrtt9+8vLxGjBhR30JjHq1We/Xq1WPHjuXm5q5duzY7O5vGMIcPH7569eqRI0fS09OTk5OlUqkjfH7uzUBvKof6MDs4B/nlONSfDPufxmD/8zDo+jCzbbRdmwoMDLx9+3bda6lUGhgYaP8MZrP5ueee8/DwOHfunJubG+2pNm3adP369YMHD5aVlW3fvl2v19OYx9fXd9SoUZ6enoSQIUOGFBQU0BgmPT39pZdeCgsLCwsLGzhw4LFjxxzh83NvBhpTOdqH2cE5wi/H0f5k2P80BvufB6Lzw2yjuUU2devWrZ49e+r1+uLi4k6dOtEy3e/XX3994YUXHC0VRVFJSUl1kxBpzHP58uUePXpUV1dXVFSEhIRcunSJxjBr1qwZO3Zs/Qy7nJwcusIcPHiwft7fvRnsnKphGIf9MDsmR/jlOOyfDPsfK9j/PDAMjR/mVnkGKCQk5I033hg8eDAhZN26dUwmDRfyjh8/npGRIRaL696uXbt2/PjxtKdqiMbfUlhY2KxZsyIjIymKWrx4cbdu3QghdIWZPn16Tk5O9+7dCSELFy6MiIigMUy9e/86NP69HP/D7FCw/3kY2P/Uwf7ngWj8MDMoirLRpgEAAAAcEw7sAAAAoM1BAQQAAABtDgogAAAAaHNQAAEAAECbgwIIAAAA2hwUQAAAANDmoAACAACANgcFEAAAALQ5KIAAAACgzUEBBAAAAG0OCiAAAABoc1AAAQAAQJuDAggAAADaHBRAAAAA0OagAAIAAIA2BwUQAAAAtEyf67oAACAASURBVDkogAAAAKDNQQEEAAAAbQ4KIAAAAGhzUAA5sxs3bjzzzDPTpk3773//S3cWAAAAB8KgKIruDGAr7777bkJCQp8+fWJiYo4dO0Z3HAAAAEeBM0DO7NNPP42IiCgpKfHx8bHdKEePHn322Wft+Y0PQ6lUNvYluVx+8+ZNG40LAACtAgogZ+bi4rJjx47ExMQff/yR7iyP5b333gsODg4MDHznnXfqWvbs2RMWFhYaGvrvf//73v4qlSomJqbu9a1bt6KiosLCwup7LlmyRK/X2yc5AAA4JhRAziwtLe306dO7d+8OCAh4mP5Go/G+r5vuaWsHDhw4fPhwYWHhpUuXUlNT//jjD5lMNm/evAMHDuTk5Pz4448XL15s2H/58uWRkZEqlaru7YoVKz755JOLFy+mpKSoVKrCwkJCSNeuXe2WHwAAHBAKIOfx3nvv9evXT6fTEUJWrlw5a9as3bt3X7ly5fnnn3/55Zcb9ly5cmXHjh07d+789ttvWyyWkydPTpkyJSEh4bPPPmv4mhDy+eefh4aGduzY8d1337XqeW8Aq82OHj16z549dV8aMGDAsWPHrDo85M/l6uq6dOlSFxcXDw+Pnj17ajSatLS0oUOHhoSEuLm5JSQkbN++vWH/UaNGLV++vP6tu7v7pUuXiouLdTodh8P5/PPPP/jgg0f5vQIAgDOiwIlMmzZt3bp1Wq02NDS0pKTkvn2OHDnSu3fvioqK2tra+Pj4ZcuWnThxwt3dvaCggKKohq9///33rl27KpVKnU43ePDgn376qeFXG27wmWeeuXezGzdufPnllymKun79eseOHTMzM6061H1j/XZSU1MjG5g1a1bDUXbu3Dl27Njx48frdLq6IqaufcOGDTNmzLD6GUtKStq3b1/3WqlUzp49e+TIkQcOHDh79uyiRYse5zcMAADOgU13AQYt6fXXX3/77bdVKlVcXJxEIrlvn0OHDimVyueff54QUllZWVNTM3To0IiIiPqrQvWvDx8+HBcX5+bmRghJSEg4fPhwx44dG/ZserObN29+//33jUbjb7/9Nm3atMOHD1t1GDx4cMMtjBs3bty4cY39aFFRURaLZdmyZZmZmVZnj5q+Hufm5rZ69eq6188999zatWv379+/atWqXr16LVu2jMfjNfG9AADgrFAAOZWBAwdWV1evXLkyOzu7sT6urq6zZ89+//33CSE6nc5gMOTn54tEovoO9a8pimIwGHWvWSyW2Wxu+NUHbtbNzW3gwIGZmZnbtm3btWvX9u3brTrk5OQ03EJaWtrSpUvr3/bu3btu+aK6OUxRUVFxcXG3b9/eu3dvZGTkkSNH6rqVlJQEBwc/zC/nwIEDTzzxhNls/uSTT3bv3r1u3bpVq1a99dZbD/O9AADgZDAHyNm0b99+0KBBnp6ejXV4+umnU1JSKisr9Xp9bGzsyZMnG+s5fPjwnTt31tbWGgyGzZs3Dx8+vIlx77vZ+Pj4FStWeHp6dujQ4YHjjhs37nQD9Ys3lpaWfvDBByqVSqVSpaend+vWbdy4cQcPHqysrNTpdL/99tuzzz5LUdSVK1dMJlNj8SiKWrVq1bx589RqtaenZ3BwcK9evaqrq5v4iQAAwImhAHIqJ0+erK6uPnz4cBO3eUdERCQlJUVFRYWGhoaGho4ZM6axnqNHj37xxRf79OnTrVu3J554YurUqU0Mfd/NTpgw4cyZM9OmTXukca3MnTs3PDw8PDy8R48eYWFhc+fO9fX1TU5Ofuqpp/r37z9jxoy+ffsaDIawsLDS0tLGNrJt27aJEyfy+fyQkJDBgwf36NFj5cqV8+fPf8gMAADgZLAStFOJjo7+4osv/vWvf73wwgvPPfcc3XEAAAAcFM4AOY89e/a4ubkNGTLktddeq5/2CwAAAPdCAeQkzGbz4sWLly1bRggZPXq0r69vREREZWUl3bkAAAAcES6BAQAAQJuDM0AAAADQ5qAAAgAAgDbHVgVQ3aJ5dXJycjZv3nzlyhUbjQUAAADwSGxVANU/fjw5OXn8+PG///77qFGjfvrpJxsNBwAAAPDwbDUJ2tfXt+4WpKCgoLNnz4rF4srKyujo6MLCQlsMBwAAAPDwbD4HyMfHRywW171gsVi2Hg4AAADggWz1MFQOh9OjR4+QkBCVSpWSkjJ16tQZM2ZYPf27aQcPHty+fbuN4gFAM7Rr127x4sV0p7AH7H8AHE2L739suA6QVqu9devWzZs3AwICwsPDV6xY8dZbb7m4uNy3c2xsbG5ubsOWsrKyoKCgd955x0bxAOCRyOXyH374obi4mO4g9vDuu+/euHHj6aefpjsIABBim/2PoyyEKJfL5XJ5w5bJkydLJJK0tDS6IgG0KZdr/5YZqqK9hjfWoaSkJCoqqu0UQL6+vgsXLqQ7CAAQYpv9j60ugT0qT09PT0/Phi08Ho/BYNCVB6CtkeqKy/XSaNJoAQQA4ExsVQCtXr366NGj97Zv2bLFRiMCAAAAPCRbFUAvv/xyVlaWSqWaO3eujYYAgMdXUJt/R19GCLmmLlCaFEerDxBCRGy3vu4D6Y4GAGBDtiqAhELhggULtm7dOmLECBsNAQCPr1wvvaW9TgipNlZqzZqb2muEEB+uH925AABsy4ZzgMLDw8PDw223fQB4fEO9R9W9OFSVVq6XTpHMpDcPAIB94GGoAAAA0OY4yl1gAEAvN7aH1qyhOwUAgJ2gAAIAQggZ4BFNdwQAAPvBJTAAAABoc1AAAQAAQJuDAggAwBFJdcVrbyfTnQLAaaEAAgBwRLXmmmpDJd0pAJwWCiAAAABoc3AXGACAA7lQk5Ot/IsQojTJ7+jLUkp+JITwmPxY8UscBpfudADOAwUQAIADcWN7tHfpRAi5Yygr05bUvWYymCzCojsagFNBAQQA4EBCXDqEuHQghBSqL91QX3nS+2m6EwE4J8wBAgAAgDYHBRAAOLk1a9bQHaE5PDneXV170p0CwGnhEhgAOJv09PTKyrs3kH/22Wc8Ho8QMm3aNPpCPTJfrv+zAS/SnQLAaaEAAgBnc/HixSVLlsTHx3t5eRFCdDpdbm4uabwAmj59+t9//92w5fr165GRkQsXLrRDWgCgBQogAHA2CxcuHDRo0KJFi6ZOnTpixIiDBw8mJze1pPKSJUvkcnnDltZ1rggAmgEFEAA4oejo6NTU1Hnz5qWmpur1+qY7d+zY0apFKBSy2dg9AjgzTIIGAOfk5uaWkpISGRkZGhpKdxYAcDgogADAmcXHx2dkZNCdAgAcDgogAAAAaHNQAAEAAECbgwIIAAAA6GSmzAbLA25WaHEogAAAAIBOJ2SHdpVvtvOgKIAAAACATibKaKbMdh4UBRAAAAC0OVjpC6Ctq9CX3dAU3tBc7erao597FN1xAKCt2CbdVFCbTwhRmZQmynhdXUAI8ecFzm73th1GRwEE0OZozZobmsIbmsLrmsIbmkIek9dR0KWjoEs7F+sFkQEAbGek7/gnPIcQQs4qTlQZKsf4PUsIcWG62Gd0FEAAbUKFvuyapuCauuCauqDScKedS2gnYdcYr6deDXrDneNJdzoAaIu8OD5eHB9CyHV1gcFiaO9i10XbUQABOCcLZbmtK7pae6lQfemauoDJYHYWduskDIvxeirEpQOTwaI7IAAAnVAAATgPE2Uq0lwtVF8srL10TVPgxfHpLOze1z0yPnC6N9eX7nQAAPfX1bWnmB9s50FRAAG0bmbKfENTWFCbf0V98Ya6MIAv6SrsMcxn9EzBAle2iO50AAAPFsRvZ/9BUQABtD4WynJLe/1ybf7l2r+vq68E8CVhwp5P+0zo3K67C0tAdzoAgFbAhgWQxWLJy8srLy83mUwSiSQiIoLJxLJDAM1Xpi+5VJN3ufbvK7UXPThe3V17j/AZO6fdQgFLSHc0AIBWxlYFUGZm5pw5c8RisUQiIYRIpdLi4uL169fHxMTYaEQAp1RrqrlUm3exJvdiTS6Dweju2ru/e/TUoDlubA+6owEAtGK2KoDmz5+flpYWGnr3lrbi4uK4uLisrCwbjQjgNCyU5brmyoWanAs1OeX60q7Cnj1Evcf4PRvAk9AdDQDASdiqADKbzWKxuGGLv78/RVE2Gg7ACSiN8vya8/k15y/V5Hlz/XqJ+j4vfqWTMIyFW9adgokysRmYdgngKGz1f2NiYmK/fv0mTZoUFBTEYDCkUum+ffsSExNtNBxAK0UR6oam8G9V9t+q7CrDnR6iiHBRvymSmbjC9Ziys7NLSkqGDh3q7u5e15Kamjp+/HgaI/1486txfrGhwq40ZgCAerYqgObOnTthwoSMjAypVEoICQkJSU9PDw5u9C5/vV6v0WgatphMJhtlA6CdzqK9UJOTpzr3t+qcO9sz3K1/guS1ToIwJgM3CrSAFStWrF69evDgwUlJSbt27erTpw8hZObMmWVlZTSmMlB6vUVHYwAAaMiG52ODg4NnzZr1kJ0nTZp05syZhi01NTU2CAXQAjRm9Rbp+teC5z/qN8qMVbnKMzmqMzc0haGCrhFuAyb5x/tw/WwRsi1btWpVXl6et7d3fn5+bGxsdna2SNTUkkjffvttYWFhw5abN29yuVwbxwQAOjnKBenff//dqiUyMtLXF2vXgiOq0lecVZx8+AKoRHcrR5mVo8yqMlb2dus/1HtUYvv3eUy+TUO2ZSKRSCgUEkJ69eo1Z86cpKSk9evXN9G/a9euVhXSgQMHOBzO4ycxUcYNxSvNlJkQUqK7tf/OtqOyA4SQCLcBUZ5DH3/7ANBstiqAVq9effTo0Xvbt2zZYqMRAexGT+nM1AMu0VKEuqYuOK88naPMshCqn/vAeMn0ToJuuMhlBwkJCdHR0bNnz541a1ZSUtLkyZPj4+OtLrI3dO/coPXr1wsELbCkJJvBifIcqjNrCSFl+pIuwh5BLu0IIR0EnR5/4wDwOGxVAL388stZWVkqlWru3Lk2GgLAzsr0JTnKM4SQKkOFhbKkV+wihHAYnCHeI/lMl7o+FspcUHvhnPKvHFWWiOXez33gG+3fD3ZpT2PsNuijjz6Kjo6urq4mhDAYjB07dmzbts3Tk56H3vcS9a17cVR2oKtrj+6i3rTEAAArtiqAhELhggULtm7dOmLECBsNAWBnRZprOcrThBCdRUsIqXtNEUa4qD+Hx71c+/dZxckc1Rlfrn8/96hFnVb4cQNoTtyGNdzzsFishISEhIQEGvMAgKOx4Ryg8PDw8PBw220fwM46Cjor3AcSQqoMFWW6kj7uAynKUm2s3F+x7W/VOX9+YH/3QRP9X8Bz15thzpw5923/z3/+Y+ckttPXLdKfF0h3CgD4x30KoNzc3Nzc3ClTppSXlzdx4zpAWxPAk4z1m0wIuVJ78ZjsQKWh/LzytB9PPMA9erJ4ihfHh+6ArRi9K/TYx3CfsXRHAIC7rAugTZs2rVixQqvVxsXFDRkyZMmSJTNmzKAlGYCjoQh1Q12YpTh2Wn6cIpQ/L/Cjzt/ifE+LGDdunFWLxWJZsmTJve0AAC3C+oaU5OTk06dP9+7d29XVNTc394svvqAlFoBDKdXd3ln2y3uXZ28s+UHEdn+74yfRnsNH+z6D6qdlffXVV4GBgS4uLu3bt3d1dc3Ly6M7EQA4LeszQGq1ms//Z3kSgUDQIithALRScmN1luL4X/IjalPtE55D5rVfXH8/13TBPFqjOacffvjh8uXLixcvnjt3rkql2rhxI92JAMBpWRdAU6ZMmThxYllZ2YYNGzZt2jRlyhRaYgHQSGfRZiv++kt+5La2qJ9H1BTJzM7C7gzCoDuX8+PxeG5ubn369MnKypo+fXpSUhLdiQDAaVkXQJ9++mlaWtqJEydKS0uXL18+ePBgWmIB2J+FslyqzTslP/K36lxX157DfMb0duvPZuAkqP3ExMRMnz590aJFzz333J07d+qfY+ocUit2POEe7ccT0x0EAAi5twCqrq4eOHDgwIED6996e3vbPRWAXUl1xSflh/+SH/Xi+AzyHPpi4AxXdlOPjgIbWbduXWFhYZcuXb777rvs7GwnuwRWUJvf0aUzCiAAB2FdAI0ePbruhVKpvHr1au/evXNzc+2eCsAeNGZ1luL4CdkhhVEW5Tl0YehnYl4Q3aHaNJlM5uPjU11dHRERERERQXccAHBm1gXQ2bNn61/n5+dPnTrVvnkAbI4i1OXa/OPVf+bXnO8p6vNMQEIP1wg8oqvKcKfWVNOe1mdUOd8BGEWov+RHjRYDIURurP5blV1puEMI6SzsFsjHKmsAdGpqJehevXoZjUa7RQGwNZmx6oTs0AnZISFLNMRrxEtBs4UsV7pDOYo81blyvZTeAqg1HoBtlW54JiCh/mFwVsyU+abmqpEyEkJ0Fm25vlRP6QghIrY7CiD7uKYu6CDozGKw6A4CDse6AFqwYEH964KCgq5du9o3D0DLM1PmXNWZY9V/FmmuRnrGJLZfFOLSge5Q8ACt5QDsvPL0SJ/xfO79CyA2g/2iZGbd629ufPy070Q8DNXOUkp+fL3dOxJ+CN1BWpLWrHFhCehO0epZF0ANb/saN27csGHD7JsHoCVV6MuOyf48KcsU8yUxXk8ldnifw+DSHcqxGCmDwWIghBgseiNlVJtrCSFsBpvH5Ns/DA7AAB5IYZR99X/t3XlYVNfdB/Az+7DvyLCpgIBLKG6vC2g0mpBGq41i1OprqlGMRqJtlooaNTFq2tiSJ75aqzXVp31jYxLqwhIR3CBGDIhIcENl02Fghp1h1jv3/WN80aKgRWbOnZnv5687x8udLzOX4+/OnHvO7Q3bonfTDmLzHhRABQUFhJDg4H8bBFpYWDhmzBhrhwJ4NkbWWNxScLbhxF1t9Xjvyb+L2BqARSi78fe7fyluKSCEGFg9y5Ki5vOEkCBp6NqIbdYPYysXYHqTvqD5nIk1EUK0Jk1BU575tsFhbsN7mBzcU+jlLMRXrvCsWg0tjQYV7RT24EEB9NlnnxFC5HJ5Xl5eTEyMq6trQUHBq6++evjwYXrxAP4zSn3d2Ybs7xtPyaTBk3wSRniMFfJ6GugGS0KSSUgyISRXlaHQyRf8//c1VmZbF2BaU0dFRzlLWEKI0WS8q6uUGpwIIYHSkB4KoKWhmNfRSlqNze9dTTKwevPDD268bd54PXjl8z4v0cvVN3SslmEZ2inswYP/Gw4dOkQIefHFF69du2b+5LmioiIpKYlaNICnZmJNV9qKTquyKjW3xntN/l3ExwGSINqh4D9gWxdg7kLPRcErzNtlbZdnB/y3r9ifEKIzac83nR7vxdEPrhyHu9DzLzH3z5wPbrxtH2OASloLjyoOEUK0Jq2JZT66+Q4hRMAXruj/nrfIl3Y6m9T14riqqioyMtK8PWDAALlcbvVIAP+BVmPLucaTZxtOeIq8J/v8HKN8eo3PE1Bc7sM+LsAa9MrM+jQUQGAJEr4kyKk/IaTV2FKnk5u3CSHCHu/mhh50feHGjRuXmJi4YMEClmUPHTo0evRoKrEAnuiW+vqphszS1kujPMcnD1iHG7ue0VjPiXqTjm6GPrwAM5lMJSUlCoXCaDQGBQXFxsby+RaZ6slb5EtlwDg8JbtZxc9L5GOeqVXEE/MIMW8LeUKJAKdfL3UtgPbu3btnz54vv/ySz+c///zzNncFBnbPPP40V5WpN2kn+/x8YdByZ4EL7VD2wEngTP3G2r66ADt9+vSKFStkMllQUBAhRC6X19TU7N+/f+LEiX2alxBC1kZs05m0lZrbhBClrs7A6s3bIp7IDr52sQNLQpJldvGFeD9J4Cv+swgh5epreY0nzdvwLB4UQMnJyVu2bPnggw8IITKZjBBy8+bNd999d+fOndTSATxEpa8/3ZCV35gb7hw1R7ZoiNvP7ObaDsz66gLs7bffzsjICA8P72ypqalJTEw0D7V+VHZ2dmVl5cMtcrm8tbV179695ofTpk0z11K3bt06depU527m9sutPx768cCdC9UMy7QxrdeEb0ZOCvMIcH97wDplVcOj+3d3HLSj/WnaXQXu+nvM3gt7OZLHOu3ffvutVqslfepBATR16lSpVNo5FT0Ad1xvL81RZdxsL4vzfmHDoE/9xP1oJ4I+1rcXYAzDmA/SqV+/fizLdrd/UVFRlwKovb3dYDAUFRWZH8bExJg74oqKis7GzvYxnhNaifabu990MO21bfWDPcMSRYvHDhpLCLlYUfjo/t0dB+1of5p2T5GXX1Mwd/JYp720tLTPZ0bldekUrl+//v3337/xxhsbNmw4c+bMJ5988vDMHNY0ZswYPz+/9PR0Ks8OXGBg9Reazp1UppsIM9V32jivSRhsQdHdu3fHjRtXU1NjiYMfPXo0ISEhNze3S/u0adN6cbTdu3fv3Llz5syZwcHBPB5PLpcfO3Zs1apVy5cvf8oj9K7/kWtrdlf94eMom/nUXMN0lKuvxbiPpB0E4Aks0f90HRW4evVqV1fXixcv5ufnf/LJJ1u3bu3DJwN4Ss2GxjTF/753NelSy4V5gYu3RH0+yedlVD92bObMmVKpNDw8XKFQTJs27Ycffti+fbuHh0fvjrZy5crs7OywsDCVSqVUKkNDQzMzM5+++uk1Pk/AJ7a05lSNtiKzPo12CgA6ug6CVigUc+fOXbt27cKFC+Pj4+/cuUMlFtiT399e/0bIavNEKU9U2XErW3W8tLVonNfz6yK2+0tkT/4ZsBerV69esmTJwxdgWVlZvTtUSEiI9e/hCJAEvh/+kZWfFAB6p2sBFBUV9c477xw+fPjHH39MTU11d3enEgvsSbuxTWd6wuA1lrDFLQXZymMNBuVU3+n/HbSc+h1JYH12cAHmKrSBPlNv0rcYmwghzYYmA6tX6usIIRK+xF3oSTsagPV0LYB27dr197///ciRIwEBAQ0NDebZyQAsR2fS5jfmnlSluwncX/KbMdJjLJ9nS18iQB/CBZh1nGn4LleVQQgxsHq1Ub3j9kZCiJAn2hi5w/6+aN5XnTpH9rqnyJt2EOCcrgWQn5/fCy+8cPny5ZiYmOXLl4eEhFCJBbZOzbT/7719DGskhDQaVF/J/2b+ROd575eGuP3MvE+LoSlHlXG2MTvKZejSkNURLtE0EwMH4ALMOl7ym/GS3wxCyE11WVrtl2sj7HmsZ7WmQs20owCCR3UtgA4cOLB9+3aNRpOYmDhhwoQNGzYsXbqUSjKwaVK+00iPseblsm+rbwx1izWvVuMvCSCE3NNWf6c8crnl4liv5z/Abe3w/3ABBgBW07UASk1NvXDhwqJFi1xdXS9fvjxq1CgUQNALAp5gpMc48/axuq+GuQ03z4p7rf3K3+/+pVpbMcV32ieD97gIXKnGBG7BBZiV8Qifz7PD2UR1Jm1mfZr5E+hWY3O28pib0J0QMsozboBT+JN+GhxF1wJIrVZLpfe/A3Z2dhaJRFaPBHbIxJouNud/V/8vPatP8JuZ7JUi5OHUgq5wAWZl4c6RS0PX0E5hEU58JxNhCSE8wncSODsLXAnBzPHwb7oWQAsWLJgxY0Ztbe0XX3xx4MCBBQsWUIkFdkNv0vmI/D6v2Ooj9p8ZMD/GfSTWr4Du4ALMyvg8gfm7aTsj4Utf9n/VvP1D05kJ3lOxLhs8qmsB9OGHH2ZkZOTn59+7d2/btm20poEGO9BubDvVkHlKlTXIZfCbA94Nd46inQi4DhdgAGA1XQsgQsi0adM6p59va2tzc3OzbiSweY0G1Qnl0fONp0d6jFsbsS1AEkg7EdgGXIBBnwtzjnQV4H8xeIwHBdClS5dSUlIaGxtfeeWV5OTklJSU6urqkpIShUJBMR/YFrm2Jkv5r5KWH+N9pm6J+tx862mDXukt9sU3X/A0cAEGfWtxyCraEYCjHqwFtnjx4nHjxqWmpsrlcvMSrCkpKRcuXKAYDmzI7Y4bOyu3f3p7Yz9J4CeD97z20Mxjn1duVWjv0Y0HHHfp0qWEhITRo0dv2rRJpVItW7YsISFh0KBBtHMBgN168AlQU1PT5s2bCSGDBw8+evSoedvBaU2ayo5b0a7P0Q7CaWVtlzPqv1Xp6172++Xy0HfEfHGXHUysyURMVLKBrVi8ePGrr746derUgwcPxsTEJCUlLViwYMCAAbRzAYDdelAACQT31x/w8fFxcXGhlIdbqjUVRxSH1kagAHoMlrCXWi5k1H9rMOl/7j9rrOcELGEBvYYLMACwsscMgu4TDMN0VlTFxcVXr14dNWpUVBTuA7IHDMsUNJ/LqP9Wynea7j9nuMd/PTq+p9XYvP3WOhPLEEIaDQ2pdz4U8ISEkBn95sZ5v0AhNHAbLsAAwMoeFEBVVVWd4w3VanXndltbWy+OGxAQoFQqCSGpqak7dux44YUX1q9f/9FHHy1atOiZM1scwzJak4YQojVpGJZRM+2EED7hY31yI2vIa8zNqv+Xn7jfgqCkIa4x3e3pLvR8N+xD8zdff7rz4YKgZf0kgYQQT5GX9eICAAB040EBpNFoLPEEf/zjHwsLC2UymVKpjIuL664AqqqqUqlUD7d0dHSYTHQGjuSo0tPrviaEMIQxMPq1194khEgFTh9H7bS/pZKfks6kPdNwIlt5rL9T+Jv93wlzjnzij/iI/cwbQp7QS+SDBb+gB317AUYIKSoqunv37qRJkzw8PMwt6enp06dPf/aoAGAfHhRAEonEEk/g6+srk8nMG52fcj9q48aNZWVlD7dUVFR0zglrZQl+MxP8ZhJCbqqvptX+Y23ENioxOELNtJ9SZeaqMga7xqwJ+yBEOoB2IrBDfXsBtn379j179sTHx69ZsyYtLW348OGEkGXLltXW1vbhswCATbPUGCCRSDR06NDQ0NDW1taDBw8uWrRo6dKlPUxrdvDgwS4tY8aM8fPzs1A8eBptxtZs5bGzDSeGe4xJidjer7fzA6ziRgAAFRZJREFUGb4a8Cs/cUDfZgM707cXYLt27SopKfHx8SktLZ09e3ZRURHmEwKALixVAMnlco1GU1VVVVlZGRAQwLJsWFjYb3/7Wws9nYXwCI9nj0slP1GrsSVbefRsQ/bP3EdtGPQHf4nsWY42wmNsXwUDeBpubm7mkdTPPffcihUr1qxZs3//ftqhAIBbLFUAEUKcnJyio6Ojo6PND9evX2+557KQMOfI14Pfop3CqlT6+qz6tIvN+fHeUz6O2umBMctgg+bPnx8XF7d8+fKkpKQ1a9bMmjVr3rx5HR0d3e1fX1/f3t7+cItWq6U1BhEArMOCBZAdEPAEjrOOVZ1OnlH/7eWWi5N8ErZF73YTutNOBNBLGzdujIuLa2hoIITweLxvvvnm8OHDXl7dVvPJycmFhYUPt9y7d08s7jqlJwDYExRAQO5pqzPqv7naVjLFd9ong/c4CzALC9i8KVOmdG4LBIL58+fPnz+/u52/+uqrLi0Ygwhg91AAObRqzZ3jdV/fUl9/yW/GouAVUr4T7UQAfW/SpElnzpyhnQIAuAUFkIO603HzeN3hak3Fy/6/XBa6Rsy3yCQIAFzAMAztCADAOSiAHM5Nddnxuq/rdPJX/Ge9NeB3Qp6IdiIAy7KJCegBwMpQADmQsrbLx+u+bjE2TfOfPc5rkgBrl4JjWLZsGe0IAMA5KIDsH0vYK61Fx+sOa02a6f6J/+U5gc/j0w4FAABAEwoge8YS9lLLheN1XxPCTu83Z6THuEeXbQcAAHBAKIDsE0vYktbCI4pDhLDT/BNHeY5H6QMAANAJBZC9MbLGH5rOZNaneQi95sgWDXWLpZ0IAACAc1AA2Q+9SZ/XePI75RGZJHhxyKpIlyG0EwEAAHAUCiB7oDNpT6uyTqiOhTlHruz//kDnQbQTAQAAcBoKINumZtpzVRmnVJmDXWPeDfswSBpKOxEAAIANQAFkq1qMTdnKY3kNOcM9xqyL+MRfIqOdCLp15cqVrKws2ikeLy4uLj4+nnYKAHBo9XqF2thm5a8vUADZHpW+Lqv+Xxeb88d5Pb85KtVb5Es7ETxBWlpaTk4OB+uMK1eulJeXczAYADiU0tYihU6OAgi6dU9bnVn/bWnrpUk+CVujd7kLPWgngqc1derUzZs3007R1f79+3/44QfaKQAAKEABZBtuqa9n1qdVdJS/6Dd9YdByJ4Ez7UQAAADPRMN0tDNthJB2pk3DdCj1dYQQZ4GLi8DVCs+OAojTWMKWthZl1qc1Gxpf9v/ligHvinhi2qEAAAD6wNe1B8vaLhNCNEwHQ5hy9VVCiL9E9k7YZis8u2MVQEUtP0S6DHUTutMO8mQMyxQ0531X/y8+T/Bz/1dHe8RhAS8AALAni4JXmDdyVRkKnXxBkFXXLXasAuh0w3dOfOchbj+jHaQnWpPmbEN2jiq9nyRwbuBiTOUMAADQ5xyrAOK4JkNDjio9ryFnqFts8oCUUKcw2okA7MHevXuTkpJopwCAbvF5AusvWOkQBZBCd09r0hJCNCZNnU7uLHQlhPQTy7gzlLhaU5GtPFbS+uN4r8kbI//oK/annQjAhmVmZiqVys6HH330kUQiIYS8/vrr9EIBQLfGe00ymPRWflL7L4CMrGF/9ecMaySEKPTyk8rj0kYnQsh478lTfafTzWYe45ytPFaruzfVd/qvgpY6C1zoRgKK2o1trkI32insQVlZ2YYNG+bNm+ft7U0I0Wq1ly9fJiiAALhKwpdK+FIrP6n9F0BCnmj9oN+bt3fc2fSK3ywujAHSmbTnm86cVB7n8/hCnvAPg/cKeALaoYAmlb4utWLL1qj/6cXPFhQUjB079vjx49On36/pt2zZsnHjRq1WK5FItmzZcujQoba2tpUrV6akpPRpao567733xo8fn5KSsmjRoilTpuTk5KSmptIOBQDcYv8FENeo9PWnVJn5TaeiXIYuDlnFElNa7ZeofoBhGRNr6vWP+/r6fvPNN50FUHp6upubGyEkNzf3yJEjly5dam9vHzlyZEJCwogRI/omMbfFxcWlp6cnJyenp6frdLqed3755ZcvXrz4cEtbW9vIkSMtGRAAKHOsAkhABHxKpQZL2GttV041ZN5UX4v3nrJp0B99xH6EkJvqMip5wM6MGDGiuLjYYDCIRKI7d+74+flVV1cTQhoaGpYvXy6VSqVSaXx8fE1NjYMUQIQQd3f3gwcP/vOf/7x+/XrPex45ckSj0Tzcsnnz5uDgYEumAwDKHKsAeiP0bTerLx+hYTrON50+1ZAlIMIpvq8sC/2NhC+92Jx/uPYAIaSdaZFr7/656lNCiIAn+HXwKjEfUx3arRPKowVN5x5uUTPtrYZmQghLWANrWHFlLiGE5ZF+4gAB79/+PBcGLw9zjuzuyHw+f/LkyadOnUpISEhLS5s9e3ZRUREh5LXXXjPvUFpaev78+d27d/f5L8Vx8+bNmzdvXs/7mAvEh1vMg6YBwI45VgHkLvS05tNVam6fbThR2Hx+mNvw14NXRLoM7fynAc4R5lv+anU1zfqWUR7jCSFSgZOIL7JmQrCyeO8pUa7DHm4xmPRKQz0hpMXQlK08Oifw14QQEU/oK+738E2hfMIPlob2fPDExMSDBw8mJCQcPXr0yJEj69atM7ezLPvZZ5/9+c9/TktL8/BwxPXjJk2adObMGdopAIBbHKsAOqXK/Jn7aPN3T5ajNWkKmvLONma3G1sn+ry0Nfp/Hi28/MUB/uIAQshNddnVttLRnnEWjQQc4SJwdXHqusbNIDKYEFKnk59ryB7v+XyvDz5+/PgVK1ZUV1dLJBIfHx9zI8Mwc+bM8fT0LCwsdHe3gTnQLYFhGNoRAIBzHKsAutRaECAJslwBdFt941zjyUstF6Jch70a8KthbsOtP7MTOCw+nz9x4sTVq1fPmjWrs/Hw4cNisfiLL76gGIy6RYsW0Y4AAJzjWAWQhTQbGs83nfm+8RRL2AneUx/7kU93/MQBIzzGWDQe2AQvke90/8RnPEhiYuKUKVN27drV2ZKXl5eVlSWTycwP9+3b13mnmONYtsyqCwwBgE1wiAKotO2SltEQQlqNzTfay9RMOyFkoHOEr7jfsxxWZ9IWtxScbzpT0VE+ynP84pBVES7R/+lBvEQ+L/nNeJYYYB/EfHGc9wu9+9kxY8ZkZWURQiZPnmwy3b+Xvra2lhCye/duBxz4DADwRPZQAKmZdhdB13EVnYys4YemMwzLEEJajS031WUK/T1ze+8KIIZlytqKLzTnXWktjHAZHO89JXlgioiHW7cAAABshj0UQFvK3/td+MdeIp/H/quQJ0oK/a15+1lmgmZY5lr7lR+bvy9uvSiTBI3xnDg/8A03oYOOKgUAALBp9lAAmVjG/AGPJehNup/aii+1XChpLZRJg0d7xP0yYH53xRYAAADYBHsogJ7e0y+31qBXlrZdKmn98ab6aphz5HD3MYmyRZ4ib0snBAAAACuw1QKoQa/cVfV7wrKEkGZD4+cVW4U8ISEkwf+XWkbzvM9Lj/2plf3f72HVLZ1Je6P9p7L2krK2y+3G1mFuw8d7TU4K/a2TwNlCvwUAAABQYasFkLfY99fBb5mIiRDyecXWWbKFniJvPuGxLNlXk9pdAfRo9aNhOm51XL/RXnZTXVajqQxzjhzsFrM0ZHV/53BM4QMAAGCvLFgAmUymkpIShUJhNBqDgoJiY2P5fH5fHZxHeKFOA83bQp4wWNrfV+xPCJFra3r+Qb1Jf1dbWaW5U9lx63bHjUa9aqDzoEjXIbMCFoQ5R2EdLrCQpqamO3fu0E7RlUqloh0BAIAOSxVAp0+fXrFihUwmCwoKIoTI5fKampr9+/dPnDjREk+nYToqNbcJIUpdnYHVm7f5hOcp8lHp6+p0coVOXqu9e1db1WhQyiQh/Z3CBjhHTPGdFiwNpbU+PDiOiIiITZs2paen0w7yGG+++SbtCAAAFFiqAHr77bczMjLCw8M7W2pqahITEwsKCh67f35+vnnetk6NjY09L11kZA0lrYUm1hTlOixXlfFT6yUTMRlZQwfT8YdbGxjWyBCjE9/FTxLQTyzzl8hGe8bNlM4LkAT1MAwIwBIWLly4cOFC2ikAAOABSxVADMN0zr5v1q9fP5Zlu9s/KyurvLz84RaNRuPq2u30hoQQtbH9x+bvWcJK+U48Hu+/vCZIBFKjyXi+6fTS0NWeQm8fsZ+YL3nGXwQAAADsj6UKoFWrVo0cOXLmzJnBwcE8Hk8ulx87dmzVqlXd7b9169YuLe+//76fX0+rlnqIvN7s/26XRrm2pri1YLBrTK+TAwAAgN2zVAG0cuXKX/ziF1lZWXK5nBASGhqamZkZEhJioafrJOQJzffDAwAAAHTHgrVCSEhIUlKS5Y7/WP4S2e/Cu36YBAAAAPAwTn9Y0tjY2POdwzdv3nRycrJanp6p1WoXFxfaKe5DmO50dHQ4OTnxeJyY5Emn0wkEAqGQE3+GDMN4eHh4eXl1t4NCobBmHurQ//QawnQH/U93qPQ/nPjNH2vgwIE7duw4fPhwdzt0dHTU1dVx5M1jWZZhGI6EIYQYDAahUMiRPzOj0SgQCLgThs/n9+GUVM+CYRgej8eRMCaTSSwWd7l3oYuYGEcZXYf+51mg/+kO+p/u0Ol/WJuVl5cXHx9PO8V9crlcJpPRTvGAs7OzWq2mneK+YcOGlZaW0k5x38yZM48cOUI7xX2/+c1v/vSnP9FOcd/OnTtXrVpFO4XNQP/TA/Q/3UH/0x0q/Q8nSj8AAAAAa0IBBAAAAA4HBRAAAAA4HBRAAAAA4HBsuAASCoUCAVdW9eJUGEKISCTiyNh+Qgh37rQkCNM9rp3DHMepl4tTYQj6n+4hTHeonMM8tvv1uTjOZDIpFIrAwEDaQe67d+9eUFAQ7RT3cS1MYGAgR25Draur8/Hx4ciffVNTk1gs5sgkJR0dHVqt1tvbm3YQ24D+pwdcC4P+57HQ/9hwAQQAAADQO1z5lBIAAADAalAAAQAAgMNBAQQAAAAOBwUQAAAAOBwUQAAAAOBwUAABAACAw7HVAmj37t1DhgyJiYk5ffo0rQxbtmwZMmRISEjI9u3buZPq008/zczMpJ7nxIkTo0eP7t+//969e6mHWbt27cCBA8PDw//xj39QDKNWq/fs2dP58NEM1kzVJQw3T2bO4sKLw823DP3Po9D/PDEMtZPZyqvP94nq6uro6Gi1Wl1eXh4VFcUwjPUz5OTkjBgxQqPRKJXK0NDQoqIiLqTKz88Xi8UHDhxgqb5KTU1NgwYNUiqVLS0tUVFRLS0tFMOcPXs2NjZWq9XW1tZ6eHi0t7dTCXPlypUlS5bMmzfP/PDRDNZM1SUMN09mzuLCi8PNtwz9z6PQ/zwxDMWT2SY/AcrKypo5c6azs3NERERAQEBJSYn1MzQ0NCxfvlwqlfr6+sbHx9fU1FBP1drampKSMm/ePPNDinmOHj06Y8YMX19fd3f30tJSNzc3imEEAoFYLBaJRFKp1DwBK5Uwhw4dam5u7nz4aAZrpuoShoMnM5dx4cXh4FuG/uex0P88MQzFk9kmC6Da2trg4GDzdnBwsEKhsH6G1157LSkpiRBSWlp6/vz5SZMmUU+VnJy8adOmzqnEKeapqamprKx87rnnQkNDd+zYwePxKIaJi4uLjIwMDAwMCQnZvHmzi4sLlTDbtm1buXJl58NHM1gzVZcwHDyZuYwLLw4H3zL0P4+F/ueJYSiezDZZAJlMpocXdjEajVRisCybmpo6e/bstLQ0Dw8Puqm++uorb2/vKVOmdLZQzKPRaMrLy8+dO3f58uV9+/YVFRVRDHPq1Kny8vIzZ85kZmampqbK5XIunD+PZqCbilMnM8dx5MXh1FuG/qc76H+eBq2TmRNLsv2nAgMDq6urzdtyuZzKeoQMw8yZM8fT07OwsNDd3Z16qgMHDty+fTsnJ6e2tvbrr7/W6XQU8/j5+SUkJHh5eRFCJkyYcP36dYphMjMzFy5cGB0dHR0dPXbs2HPnznHh/Hk0A8VUXDuZOY4LLw7X3jL0P91B//NENE9mC40tsqiqqqphw4bpdLqampqIiAgqw/2+/PLLuXPnci0Vy7Jr1qwxD0KkmOfatWtDhw5taGior68PDQ29evUqxTB79+595ZVXOkfYFRcX0wqTk5PTOe7v0QxWTvVwGM6ezNzEhReHs28Z+p8u0P88MQzFk9kmPwEKDQ1966234uPjCSF//etf+XwKX+Tl5eVlZWXJZDLzw3379k2fPp16qodRfJWio6OTkpLGjBnDsuy6desGDx5MCKEVZsmSJcXFxUOGDCGEvPfee7GxsRTDdHr03aH4fnH/ZOYU9D9PA/2PGfqfJ6J4MvNYlrXQoQEAAAC4CRd2AAAA4HBQAAEAAIDDQQEEAAAADgcFEAAAADgcFEAAAADgcFAAAQAAgMNBAQQAAAAOBwUQAAAAOBwUQAAAAOBwUAABAACAw0EBBAAAAA4HBRAAAAA4HBRAAAAA4HBQAAEAAIDDQQEEAAAADgcFEAAAADgcFEAAAADgcIS0A4DNmzVrVnFxsVarbW5uDggIIITMnz+/ra1t586dtKMBgJ1D/wO9xmNZlnYGsAdnz579+OOPT548SQjRarVardbT05N2KABwCOh/oBfwFRj0vdzc3A8++CA3N3fq1KnPPfdcSEjIp59+OmPGjMGDB2/ZssW8z+bNmyMiImJjYw8cOEA1LADYFfQ/8JTwFRhY0I0bN65evXrt2rWxY8f+9NNP/fv3Dw8PX7du3YkTJ/Lz869cuaLT6caOHRsbGxsbG0s7LADYFfQ/0DMUQGBBkydPdnNzi4qKioyMHDJkCCHEy8vLaDTm5eUpFIq5c+cSQjQaTXFxMTogAOhb6H+gZyiAwIIEAoF5QyQSPdzu5OT0zjvvLF68mBDS3t4ukUgohAMAu4b+B3qGMUBAweTJk//2t7/p9fqWlpZhw4YpFAraiQDAUaD/ATN8AgQUTJgw4cUXX4yJiTEYDOvXrw8JCaGdCAAcBfofMMNt8AAAAOBw8BUYAAAAOBwUQAAAAOBwUAABAACAw0EBBAAAAA4HBRAAAAA4HBRAAAAA4HBQAAEAAIDDQQEEAAAADgcFEAAAADgcFEAAAADgcFAAAQAAgMNBAQQAAAAOBwUQAAAAOBwUQAAAAOBwUAABAACAw/k/uBfMiaBLr+IAAAAASUVORK5CYII=" /><!-- --></p>
</div>
<div id="background" class="section level1">
<h1>Background</h1>
<p>The <code>mkin</code> package <span class="citation">(J. Ranke
2021)</span> implements the approach to degradation kinetics recommended
in the kinetics report provided by the FOrum for Co-ordination of
pesticide fate models and their USe <span class="citation">(FOCUS Work
Group on Degradation Kinetics 2006, 2014)</span>. It covers data series
describing the decline of one compound, data series with transformation
products (commonly termed metabolites) and data series for more than one
compartment. It is possible to include back reactions. Therefore,
equilibrium reactions and equilibrium partitioning can be specified,
although this often leads to an overparameterisation of the model.</p>
<p>When the first <code>mkin</code> code was published in 2010, the most
commonly used tools for fitting more complex kinetic degradation models
to experimental data were KinGUI <span class="citation">(Schäfer et al.
2007)</span>, a MATLAB based tool with a graphical user interface that
was specifically tailored to the task and included some output as
proposed by the FOCUS Kinetics Workgroup, and ModelMaker, a general
purpose compartment based tool providing infrastructure for fitting
dynamic simulation models based on differential equations to data.</p>
<p>The ‘mkin’ code was first uploaded to the BerliOS development
platform. When this was taken down, the version control history was
imported into the R-Forge site (see <em>e.g.</em> <a href="https://cgit.jrwb.de/mkin/commit/?id=30cbb4092f6d2d3beff5800603374a0d009ad770">the
initial commit on 11 May 2010</a>), where the code is still being
updated.</p>
<p>At that time, the R package <code>FME</code> (Flexible Modelling
Environment) <span class="citation">(Soetaert and Petzoldt 2010)</span>
was already available, and provided a good basis for developing a
package specifically tailored to the task. The remaining challenge was
to make it as easy as possible for the users (including the author of
this vignette) to specify the system of differential equations and to
include the output requested by the FOCUS guidance, such as the <span class="math inline">\(\chi^2\)</span> error level as defined in this
guidance.</p>
<p>Also, <code>mkin</code> introduced using analytical solutions for
parent only kinetics for improved optimization speed. Later, Eigenvalue
based solutions were introduced to <code>mkin</code> for the case of
linear differential equations (<em>i.e.</em> where the FOMC or DFOP
models were not used for the parent compound), greatly improving the
optimization speed for these cases. This, has become somehow obsolete,
as the use of compiled code described below gives even faster execution
times.</p>
<p>The possibility to specify back-reactions and a biphasic model
(SFORB) for metabolites were present in <code>mkin</code> from the very
beginning.</p>
<div id="derived-software-tools" class="section level2">
<h2>Derived software tools</h2>
<p>Soon after the publication of <code>mkin</code>, two derived tools
were published, namely KinGUII (developed at Bayer Crop Science) and
CAKE (commissioned to Tessella by Syngenta), which added a graphical
user interface (GUI), and added fitting by iteratively reweighted least
squares (IRLS) and characterisation of likely parameter distributions by
Markov Chain Monte Carlo (MCMC) sampling.</p>
<p>CAKE focuses on a smooth use experience, sacrificing some flexibility
in the model definition, originally allowing only two primary
metabolites in parallel. The current version 3.4 of CAKE released in May
2020 uses a scheme for up to six metabolites in a flexible arrangement
and supports biphasic modelling of metabolites, but does not support
back-reactions (non-instantaneous equilibria).</p>
<p>KinGUI offers an even more flexible widget for specifying complex
kinetic models. Back-reactions (non-instantaneous equilibria) were
supported early on, but until 2014, only simple first-order models could
be specified for transformation products. Starting with KinGUII version
2.1, biphasic modelling of metabolites was also available in
KinGUII.</p>
<p>A further graphical user interface (GUI) that has recently been
brought to a decent degree of maturity is the browser based GUI named
<code>gmkin</code>. Please see its <a href="https://pkgdown.jrwb.de/gmkin/">documentation page</a> and <a href="https://pkgdown.jrwb.de/gmkin/articles/gmkin_manual.html">manual</a>
for further information.</p>
<p>A comparison of scope, usability and numerical results obtained with
these tools has been recently been published by <span class="citation">Johannes Ranke, Wöltjen, and Meinecke
(2018)</span>.</p>
</div>
</div>
<div id="unique-features" class="section level1">
<h1>Unique features</h1>
<p>Currently, the main unique features available in <code>mkin</code>
are</p>
<ul>
<li>the <a href="https://pkgdown.jrwb.de/mkin/articles/web_only/compiled_models.html">speed
increase</a> by using compiled code when a compiler is present,</li>
<li>parallel model fitting on multicore machines using the <a href="https://pkgdown.jrwb.de/mkin/reference/mmkin.html"><code>mmkin</code>
function</a>,</li>
<li>the estimation of parameter confidence intervals based on
transformed parameters (see below) and</li>
<li>the possibility to use the <a href="https://pkgdown.jrwb.de/mkin/reference/sigma_twocomp.html">two-component
error model</a></li>
</ul>
<p>The iteratively reweighted least squares fitting of different
variances for each variable as introduced by <span class="citation">Gao
et al. (2011)</span> has been available in mkin since <a href="https://pkgdown.jrwb.de/mkin/news/index.html#mkin-0-9-22-2013-10-26">version
0.9-22</a>. With <a href="https://pkgdown.jrwb.de/mkin/news/index.html#mkin-0-9-49-5-2019-07-04">release
0.9.49.5</a>, the IRLS algorithm has been complemented by direct or
step-wise maximisation of the likelihood function, which makes it
possible not only to fit the variance by variable error model but also a
<a href="https://pkgdown.jrwb.de/mkin/reference/sigma_twocomp.html">two-component
error model</a> inspired by error models developed in analytical
chemistry <span class="citation">(Johannes Ranke and Meinecke
2019)</span>.</p>
</div>
<div id="internal-parameter-transformations" class="section level1">
<h1>Internal parameter transformations</h1>
<p>For rate constants, the log transformation is used, as proposed by
Bates and Watts <span class="citation">(1988, 77, 149)</span>.
Approximate intervals are constructed for the transformed rate constants
<span class="citation">(compare Bates and Watts 1988, 135)</span>,
<em>i.e.</em> for their logarithms. Confidence intervals for the rate
constants are then obtained using the appropriate backtransformation
using the exponential function.</p>
<p>In the first version of <code>mkin</code> allowing for specifying
models using formation fractions, a home-made reparameterisation was
used in order to ensure that the sum of formation fractions would not
exceed unity.</p>
<p>This method is still used in the current version of KinGUII (v2.1
from April 2014), with a modification that allows for fixing the pathway
to sink to zero. CAKE uses penalties in the objective function in order
to enforce this constraint.</p>
<p>In 2012, an alternative reparameterisation of the formation fractions
was proposed together with René Lehmann <span class="citation">(J. Ranke
and Lehmann 2012)</span>, based on isometric logratio transformation
(ILR). The aim was to improve the validity of the linear approximation
of the objective function during the parameter estimation procedure as
well as in the subsequent calculation of parameter confidence intervals.
In the current version of mkin, a logit transformation is used for
parameters that are bound between 0 and 1, such as the g parameter of
the DFOP model.</p>
<div id="confidence-intervals-based-on-transformed-parameters" class="section level2">
<h2>Confidence intervals based on transformed parameters</h2>
<p>In the first attempt at providing improved parameter confidence
intervals introduced to <code>mkin</code> in 2013, confidence intervals
obtained from FME on the transformed parameters were simply all
backtransformed one by one to yield asymmetric confidence intervals for
the backtransformed parameters.</p>
<p>However, while there is a 1:1 relation between the rate constants in
the model and the transformed parameters fitted in the model, the
parameters obtained by the isometric logratio transformation are
calculated from the set of formation fractions that quantify the paths
to each of the compounds formed from a specific parent compound, and no
such 1:1 relation exists.</p>
<p>Therefore, parameter confidence intervals for formation fractions
obtained with this method only appear valid for the case of a single
transformation product, where currently the logit transformation is used
for the formation fraction.</p>
<p>The confidence intervals obtained by backtransformation for the cases
where a 1:1 relation between transformed and original parameter exist
are considered by the author of this vignette to be more accurate than
those obtained using a re-estimation of the Hessian matrix after
backtransformation, as implemented in the FME package.</p>
</div>
<div id="parameter-t-test-based-on-untransformed-parameters" class="section level2">
<h2>Parameter t-test based on untransformed parameters</h2>
<p>The standard output of many nonlinear regression software packages
includes the results from a test for significant difference from zero
for all parameters. Such a test is also recommended to check the
validity of rate constants in the FOCUS guidance <span class="citation">(FOCUS Work Group on Degradation Kinetics 2014,
96ff)</span>.</p>
<p>It has been argued that the precondition for this test, <em>i.e.</em>
normal distribution of the estimator for the parameters, is not
fulfilled in the case of nonlinear regression <span class="citation">(J.
Ranke and Lehmann 2015)</span>. However, this test is commonly used by
industry, consultants and national authorities in order to decide on the
reliability of parameter estimates, based on the FOCUS guidance
mentioned above. Therefore, the results of this one-sided t-test are
included in the summary output from <code>mkin</code>.</p>
<p>As it is not reasonable to test for significant difference of the
transformed parameters (<em>e.g.</em> <span class="math inline">\(log(k)\)</span>) from zero, the t-test is
calculated based on the model definition before parameter
transformation, <em>i.e.</em> in a similar way as in packages that do
not apply such an internal parameter transformation. A note is included
in the <code>mkin</code> output, pointing to the fact that the t-test is
based on the unjustified assumption of normal distribution of the
parameter estimators.</p>
</div>
</div>
<div id="references" class="section level1">
<h1>References</h1>
<!-- vim: set foldmethod=syntax: -->
<div id="refs" class="references csl-bib-body hanging-indent">
<div id="ref-bates1988" class="csl-entry">
Bates, D., and D. Watts. 1988. <em>Nonlinear Regression and Its
Applications</em>. Wiley-Interscience.
</div>
<div id="ref-FOCUS2006" class="csl-entry">
FOCUS Work Group on Degradation Kinetics. 2006. <em>Guidance Document on
Estimating Persistence and Degradation Kinetics from Environmental Fate
Studies on Pesticides in EU Registration. Report of the FOCUS Work Group
on Degradation Kinetics</em>. <a href="http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics">http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics</a>.
</div>
<div id="ref-FOCUSkinetics2014" class="csl-entry">
———. 2014. <em>Generic Guidance for Estimating Persistence and
Degradation Kinetics from Environmental Fate Studies on Pesticides in EU
Registration</em>. 1.1 ed. <a href="http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics">http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics</a>.
</div>
<div id="ref-gao11" class="csl-entry">
Gao, Z., J. W. Green, J. Vanderborght, and W. Schmitt. 2011.
<span>“Improving Uncertainty Analysis in Kinetic Evaluations Using
Iteratively Reweighted Least Squares.”</span> Journal. <em>Environmental
Science and Technology</em> 45: 4429–37.
</div>
<div id="ref-pkg:mkin" class="csl-entry">
Ranke, J. 2021. <em>‘<span class="nocase">mkin</span>‘:
<span>K</span>inetic Evaluation of Chemical Degradation Data</em>. <a href="https://CRAN.R-project.org/package=mkin">https://CRAN.R-project.org/package=mkin</a>.
</div>
<div id="ref-ranke2012" class="csl-entry">
Ranke, J., and R. Lehmann. 2012. <span>“Parameter Reliability in Kinetic
Evaluation of Environmental Metabolism Data - Assessment and the
Influence of Model Specification.”</span> In <em>SETAC World 20-24
May</em>. Berlin. <a href="https://jrwb.de/posters/Poster_SETAC_2012_Kinetic_parameter_uncertainty_model_parameterization_Lehmann_Ranke.pdf">https://jrwb.de/posters/Poster_SETAC_2012_Kinetic_parameter_uncertainty_model_parameterization_Lehmann_Ranke.pdf</a>.
</div>
<div id="ref-ranke2015" class="csl-entry">
———. 2015. <span>“To t-Test or Not to t-Test, That Is the
Question.”</span> In <em>XV Symposium on Pesticide Chemistry 2-4
September 2015</em>. Piacenza. <a href="https://jrwb.de/posters/piacenza_2015.pdf">https://jrwb.de/posters/piacenza_2015.pdf</a>.
</div>
<div id="ref-ranke2019" class="csl-entry">
Ranke, Johannes, and Stefan Meinecke. 2019. <span>“Error Models for the
Kinetic Evaluation of Chemical Degradation Data.”</span>
<em>Environments</em> 6 (12). <a href="https://doi.org/10.3390/environments6120124">https://doi.org/10.3390/environments6120124</a>.
</div>
<div id="ref-ranke2018" class="csl-entry">
Ranke, Johannes, Janina Wöltjen, and Stefan Meinecke. 2018.
<span>“Comparison of Software Tools for Kinetic Evaluation of Chemical
Degradation Data.”</span> <em>Environmental Sciences Europe</em> 30 (1):
17. <a href="https://doi.org/10.1186/s12302-018-0145-1">https://doi.org/10.1186/s12302-018-0145-1</a>.
</div>
<div id="ref-schaefer2007" class="csl-entry">
Schäfer, D., B. Mikolasch, P. Rainbird, and B. Harvey. 2007.
<span>“<span>KinGUI</span>: A New Kinetic Software Tool for Evaluations
According to <span>FOCUS</span> Degradation Kinetics.”</span> In
<em>Proceedings of the XIII Symposium Pesticide Chemistry</em>, edited
by Del Re A. A. M., Capri E., Fragoulis G., and Trevisan M., 916–23.
Piacenza.
</div>
<div id="ref-soetaert2010" class="csl-entry">
Soetaert, Karline, and Thomas Petzoldt. 2010. <span>“Inverse Modelling,
Sensitivity and Monte Carlo Analysis in <span>R</span> Using Package
<span>FME</span>.”</span> <em>Journal of Statistical Software</em> 33
(3): 1–28. <a href="https://doi.org/10.18637/jss.v033.i03">https://doi.org/10.18637/jss.v033.i03</a>.
</div>
</div>
</div>
<!-- code folding -->
<!-- dynamically load mathjax for compatibility with self-contained -->
<script>
(function () {
var script = document.createElement("script");
script.type = "text/javascript";
script.src = "https://mathjax.rstudio.com/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML";
document.getElementsByTagName("head")[0].appendChild(script);
})();
</script>
</body>
</html>
|
