Rsnake/Blog.html at master · mokjpn/Rsnake · GitHub

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
<!DOCTYPE html>
<!-- saved from url=(0014)about:internet -->
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8"/>

<title>Rでスネークゲーム</title>

<style type="text/css">
body, td {
   font-family: sans-serif;
   background-color: white;
   font-size: 12px;
   margin: 8px;
}

tt, code, pre {
   font-family: 'DejaVu Sans Mono', 'Droid Sans Mono', 'Lucida Console', Consolas, Monaco, monospace;
}

h1 {
   font-size:2.2em;
}

h2 {
   font-size:1.8em;
}

h3 {
   font-size:1.4em;
}

h4 {
   font-size:1.0em;
}

h5 {
   font-size:0.9em;
}

h6 {
   font-size:0.8em;
}

a:visited {
   color: rgb(50%, 0%, 50%);
}

pre {
   margin-top: 0;
   max-width: 95%;
   border: 1px solid #ccc;
   white-space: pre-wrap;
}

pre code {
   display: block; padding: 0.5em;
}

code.r, code.cpp {
   background-color: #F8F8F8;
}

table, td, th {
  border: none;
}

blockquote {
   color:#666666;
   margin:0;
   padding-left: 1em;
   border-left: 0.5em #EEE solid;
}

hr {
   height: 0px;
   border-bottom: none;
   border-top-width: thin;
   border-top-style: dotted;
   border-top-color: #999999;
}

@media print {
   * {
      background: transparent !important;
      color: black !important;
      filter:none !important;
      -ms-filter: none !important;
   }

   body {
      font-size:12pt;
      max-width:100%;
   }

   a, a:visited {
      text-decoration: underline;
   }

   hr {
      visibility: hidden;
      page-break-before: always;
   }

   pre, blockquote {
      padding-right: 1em;
      page-break-inside: avoid;
   }

   tr, img {
      page-break-inside: avoid;
   }

   img {
      max-width: 100% !important;
   }

   @page :left {
      margin: 15mm 20mm 15mm 10mm;
   }

   @page :right {
      margin: 15mm 10mm 15mm 20mm;
   }

   p, h2, h3 {
      orphans: 3; widows: 3;
   }

   h2, h3 {
      page-break-after: avoid;
   }
}

</style>

<!-- Styles for R syntax highlighter -->
<style type="text/css">
   pre .operator,
   pre .paren {
     color: rgb(104, 118, 135)
   }

   pre .literal {
     color: rgb(88, 72, 246)
   }

   pre .number {
     color: rgb(0, 0, 205);
   }

   pre .comment {
     color: rgb(76, 136, 107);
   }

   pre .keyword {
     color: rgb(0, 0, 255);
   }

   pre .identifier {
     color: rgb(0, 0, 0);
   }

   pre .string {
     color: rgb(3, 106, 7);
   }
</style>

<!-- R syntax highlighter -->
<script type="text/javascript">
var hljs=new function(){function m(p){return p.replace(/&/gm,"&amp;").replace(/</gm,"&lt;")}function f(r,q,p){return RegExp(q,"m"+(r.cI?"i":"")+(p?"g":""))}function b(r){for(var p=0;p<r.childNodes.length;p++){var q=r.childNodes[p];if(q.nodeName=="CODE"){return q}if(!(q.nodeType==3&&q.nodeValue.match(/\s+/))){break}}}function h(t,s){var p="";for(var r=0;r<t.childNodes.length;r++){if(t.childNodes[r].nodeType==3){var q=t.childNodes[r].nodeValue;if(s){q=q.replace(/\n/g,"")}p+=q}else{if(t.childNodes[r].nodeName=="BR"){p+="\n"}else{p+=h(t.childNodes[r])}}}if(/MSIE [678]/.test(navigator.userAgent)){p=p.replace(/\r/g,"\n")}return p}function a(s){var r=s.className.split(/\s+/);r=r.concat(s.parentNode.className.split(/\s+/));for(var q=0;q<r.length;q++){var p=r[q].replace(/^language-/,"");if(e[p]){return p}}}function c(q){var p=[];(function(s,t){for(var r=0;r<s.childNodes.length;r++){if(s.childNodes[r].nodeType==3){t+=s.childNodes[r].nodeValue.length}else{if(s.childNodes[r].nodeName=="BR"){t+=1}else{if(s.childNodes[r].nodeType==1){p.push({event:"start",offset:t,node:s.childNodes[r]});t=arguments.callee(s.childNodes[r],t);p.push({event:"stop",offset:t,node:s.childNodes[r]})}}}}return t})(q,0);return p}function k(y,w,x){var q=0;var z="";var s=[];function u(){if(y.length&&w.length){if(y[0].offset!=w[0].offset){return(y[0].offset<w[0].offset)?y:w}else{return w[0].event=="start"?y:w}}else{return y.length?y:w}}function t(D){var A="<"+D.nodeName.toLowerCase();for(var B=0;B<D.attributes.length;B++){var C=D.attributes[B];A+=" "+C.nodeName.toLowerCase();if(C.value!==undefined&&C.value!==false&&C.value!==null){A+='="'+m(C.value)+'"'}}return A+">"}while(y.length||w.length){var v=u().splice(0,1)[0];z+=m(x.substr(q,v.offset-q));q=v.offset;if(v.event=="start"){z+=t(v.node);s.push(v.node)}else{if(v.event=="stop"){var p,r=s.length;do{r--;p=s[r];z+=("</"+p.nodeName.toLowerCase()+">")}while(p!=v.node);s.splice(r,1);while(r<s.length){z+=t(s[r]);r++}}}}return z+m(x.substr(q))}function j(){function q(x,y,v){if(x.compiled){return}var u;var s=[];if(x.k){x.lR=f(y,x.l||hljs.IR,true);for(var w in x.k){if(!x.k.hasOwnProperty(w)){continue}if(x.k[w] instanceof Object){u=x.k[w]}else{u=x.k;w="keyword"}for(var r in u){if(!u.hasOwnProperty(r)){continue}x.k[r]=[w,u[r]];s.push(r)}}}if(!v){if(x.bWK){x.b="\\b("+s.join("|")+")\\s"}x.bR=f(y,x.b?x.b:"\\B|\\b");if(!x.e&&!x.eW){x.e="\\B|\\b"}if(x.e){x.eR=f(y,x.e)}}if(x.i){x.iR=f(y,x.i)}if(x.r===undefined){x.r=1}if(!x.c){x.c=[]}x.compiled=true;for(var t=0;t<x.c.length;t++){if(x.c[t]=="self"){x.c[t]=x}q(x.c[t],y,false)}if(x.starts){q(x.starts,y,false)}}for(var p in e){if(!e.hasOwnProperty(p)){continue}q(e[p].dM,e[p],true)}}function d(B,C){if(!j.called){j();j.called=true}function q(r,M){for(var L=0;L<M.c.length;L++){if((M.c[L].bR.exec(r)||[null])[0]==r){return M.c[L]}}}function v(L,r){if(D[L].e&&D[L].eR.test(r)){return 1}if(D[L].eW){var M=v(L-1,r);return M?M+1:0}return 0}function w(r,L){return L.i&&L.iR.test(r)}function K(N,O){var M=[];for(var L=0;L<N.c.length;L++){M.push(N.c[L].b)}var r=D.length-1;do{if(D[r].e){M.push(D[r].e)}r--}while(D[r+1].eW);if(N.i){M.push(N.i)}return f(O,M.join("|"),true)}function p(M,L){var N=D[D.length-1];if(!N.t){N.t=K(N,E)}N.t.lastIndex=L;var r=N.t.exec(M);return r?[M.substr(L,r.index-L),r[0],false]:[M.substr(L),"",true]}function z(N,r){var L=E.cI?r[0].toLowerCase():r[0];var M=N.k[L];if(M&&M instanceof Array){return M}return false}function F(L,P){L=m(L);if(!P.k){return L}var r="";var O=0;P.lR.lastIndex=0;var M=P.lR.exec(L);while(M){r+=L.substr(O,M.index-O);var N=z(P,M);if(N){x+=N[1];r+='<span class="'+N[0]+'">'+M[0]+"</span>"}else{r+=M[0]}O=P.lR.lastIndex;M=P.lR.exec(L)}return r+L.substr(O,L.length-O)}function J(L,M){if(M.sL&&e[M.sL]){var r=d(M.sL,L);x+=r.keyword_count;return r.value}else{return F(L,M)}}function I(M,r){var L=M.cN?'<span class="'+M.cN+'">':"";if(M.rB){y+=L;M.buffer=""}else{if(M.eB){y+=m(r)+L;M.buffer=""}else{y+=L;M.buffer=r}}D.push(M);A+=M.r}function G(N,M,Q){var R=D[D.length-1];if(Q){y+=J(R.buffer+N,R);return false}var P=q(M,R);if(P){y+=J(R.buffer+N,R);I(P,M);return P.rB}var L=v(D.length-1,M);if(L){var O=R.cN?"</span>":"";if(R.rE){y+=J(R.buffer+N,R)+O}else{if(R.eE){y+=J(R.buffer+N,R)+O+m(M)}else{y+=J(R.buffer+N+M,R)+O}}while(L>1){O=D[D.length-2].cN?"</span>":"";y+=O;L--;D.length--}var r=D[D.length-1];D.length--;D[D.length-1].buffer="";if(r.starts){I(r.starts,"")}return R.rE}if(w(M,R)){throw"Illegal"}}var E=e[B];var D=[E.dM];var A=0;var x=0;var y="";try{var s,u=0;E.dM.buffer="";do{s=p(C,u);var t=G(s[0],s[1],s[2]);u+=s[0].length;if(!t){u+=s[1].length}}while(!s[2]);if(D.length>1){throw"Illegal"}return{r:A,keyword_count:x,value:y}}catch(H){if(H=="Illegal"){return{r:0,keyword_count:0,value:m(C)}}else{throw H}}}function g(t){var p={keyword_count:0,r:0,value:m(t)};var r=p;for(var q in e){if(!e.hasOwnProperty(q)){continue}var s=d(q,t);s.language=q;if(s.keyword_count+s.r>r.keyword_count+r.r){r=s}if(s.keyword_count+s.r>p.keyword_count+p.r){r=p;p=s}}if(r.language){p.second_best=r}return p}function i(r,q,p){if(q){r=r.replace(/^((<[^>]+>|\t)+)/gm,function(t,w,v,u){return w.replace(/\t/g,q)})}if(p){r=r.replace(/\n/g,"<br>")}return r}function n(t,w,r){var x=h(t,r);var v=a(t);var y,s;if(v){y=d(v,x)}else{return}var q=c(t);if(q.length){s=document.createElement("pre");s.innerHTML=y.value;y.value=k(q,c(s),x)}y.value=i(y.value,w,r);var u=t.className;if(!u.match("(\\s|^)(language-)?"+v+"(\\s|$)")){u=u?(u+" "+v):v}if(/MSIE [678]/.test(navigator.userAgent)&&t.tagName=="CODE"&&t.parentNode.tagName=="PRE"){s=t.parentNode;var p=document.createElement("div");p.innerHTML="<pre><code>"+y.value+"</code></pre>";t=p.firstChild.firstChild;p.firstChild.cN=s.cN;s.parentNode.replaceChild(p.firstChild,s)}else{t.innerHTML=y.value}t.className=u;t.result={language:v,kw:y.keyword_count,re:y.r};if(y.second_best){t.second_best={language:y.second_best.language,kw:y.second_best.keyword_count,re:y.second_best.r}}}function o(){if(o.called){return}o.called=true;var r=document.getElementsByTagName("pre");for(var p=0;p<r.length;p++){var q=b(r[p]);if(q){n(q,hljs.tabReplace)}}}function l(){if(window.addEventListener){window.addEventListener("DOMContentLoaded",o,false);window.addEventListener("load",o,false)}else{if(window.attachEvent){window.attachEvent("onload",o)}else{window.onload=o}}}var e={};this.LANGUAGES=e;this.highlight=d;this.highlightAuto=g;this.fixMarkup=i;this.highlightBlock=n;this.initHighlighting=o;this.initHighlightingOnLoad=l;this.IR="[a-zA-Z][a-zA-Z0-9_]*";this.UIR="[a-zA-Z_][a-zA-Z0-9_]*";this.NR="\\b\\d+(\\.\\d+)?";this.CNR="\\b(0[xX][a-fA-F0-9]+|(\\d+(\\.\\d*)?|\\.\\d+)([eE][-+]?\\d+)?)";this.BNR="\\b(0b[01]+)";this.RSR="!|!=|!==|%|%=|&|&&|&=|\\*|\\*=|\\+|\\+=|,|\\.|-|-=|/|/=|:|;|<|<<|<<=|<=|=|==|===|>|>=|>>|>>=|>>>|>>>=|\\?|\\[|\\{|\\(|\\^|\\^=|\\||\\|=|\\|\\||~";this.ER="(?![\\s\\S])";this.BE={b:"\\\\.",r:0};this.ASM={cN:"string",b:"'",e:"'",i:"\\n",c:[this.BE],r:0};this.QSM={cN:"string",b:'"',e:'"',i:"\\n",c:[this.BE],r:0};this.CLCM={cN:"comment",b:"//",e:"$"};this.CBLCLM={cN:"comment",b:"/\\*",e:"\\*/"};this.HCM={cN:"comment",b:"#",e:"$"};this.NM={cN:"number",b:this.NR,r:0};this.CNM={cN:"number",b:this.CNR,r:0};this.BNM={cN:"number",b:this.BNR,r:0};this.inherit=function(r,s){var p={};for(var q in r){p[q]=r[q]}if(s){for(var q in s){p[q]=s[q]}}return p}}();hljs.LANGUAGES.cpp=function(){var a={keyword:{"false":1,"int":1,"float":1,"while":1,"private":1,"char":1,"catch":1,"export":1,virtual:1,operator:2,sizeof:2,dynamic_cast:2,typedef:2,const_cast:2,"const":1,struct:1,"for":1,static_cast:2,union:1,namespace:1,unsigned:1,"long":1,"throw":1,"volatile":2,"static":1,"protected":1,bool:1,template:1,mutable:1,"if":1,"public":1,friend:2,"do":1,"return":1,"goto":1,auto:1,"void":2,"enum":1,"else":1,"break":1,"new":1,extern:1,using:1,"true":1,"class":1,asm:1,"case":1,typeid:1,"short":1,reinterpret_cast:2,"default":1,"double":1,register:1,explicit:1,signed:1,typename:1,"try":1,"this":1,"switch":1,"continue":1,wchar_t:1,inline:1,"delete":1,alignof:1,char16_t:1,char32_t:1,constexpr:1,decltype:1,noexcept:1,nullptr:1,static_assert:1,thread_local:1,restrict:1,_Bool:1,complex:1},built_in:{std:1,string:1,cin:1,cout:1,cerr:1,clog:1,stringstream:1,istringstream:1,ostringstream:1,auto_ptr:1,deque:1,list:1,queue:1,stack:1,vector:1,map:1,set:1,bitset:1,multiset:1,multimap:1,unordered_set:1,unordered_map:1,unordered_multiset:1,unordered_multimap:1,array:1,shared_ptr:1}};return{dM:{k:a,i:"</",c:[hljs.CLCM,hljs.CBLCLM,hljs.QSM,{cN:"string",b:"'\\\\?.",e:"'",i:"."},{cN:"number",b:"\\b(\\d+(\\.\\d*)?|\\.\\d+)(u|U|l|L|ul|UL|f|F)"},hljs.CNM,{cN:"preprocessor",b:"#",e:"$"},{cN:"stl_container",b:"\\b(deque|list|queue|stack|vector|map|set|bitset|multiset|multimap|unordered_map|unordered_set|unordered_multiset|unordered_multimap|array)\\s*<",e:">",k:a,r:10,c:["self"]}]}}}();hljs.LANGUAGES.r={dM:{c:[hljs.HCM,{cN:"number",b:"\\b0[xX][0-9a-fA-F]+[Li]?\\b",e:hljs.IMMEDIATE_RE,r:0},{cN:"number",b:"\\b\\d+(?:[eE][+\\-]?\\d*)?L\\b",e:hljs.IMMEDIATE_RE,r:0},{cN:"number",b:"\\b\\d+\\.(?!\\d)(?:i\\b)?",e:hljs.IMMEDIATE_RE,r:1},{cN:"number",b:"\\b\\d+(?:\\.\\d*)?(?:[eE][+\\-]?\\d*)?i?\\b",e:hljs.IMMEDIATE_RE,r:0},{cN:"number",b:"\\.\\d+(?:[eE][+\\-]?\\d*)?i?\\b",e:hljs.IMMEDIATE_RE,r:1},{cN:"keyword",b:"(?:tryCatch|library|setGeneric|setGroupGeneric)\\b",e:hljs.IMMEDIATE_RE,r:10},{cN:"keyword",b:"\\.\\.\\.",e:hljs.IMMEDIATE_RE,r:10},{cN:"keyword",b:"\\.\\.\\d+(?![\\w.])",e:hljs.IMMEDIATE_RE,r:10},{cN:"keyword",b:"\\b(?:function)",e:hljs.IMMEDIATE_RE,r:2},{cN:"keyword",b:"(?:if|in|break|next|repeat|else|for|return|switch|while|try|stop|warning|require|attach|detach|source|setMethod|setClass)\\b",e:hljs.IMMEDIATE_RE,r:1},{cN:"literal",b:"(?:NA|NA_integer_|NA_real_|NA_character_|NA_complex_)\\b",e:hljs.IMMEDIATE_RE,r:10},{cN:"literal",b:"(?:NULL|TRUE|FALSE|T|F|Inf|NaN)\\b",e:hljs.IMMEDIATE_RE,r:1},{cN:"identifier",b:"[a-zA-Z.][a-zA-Z0-9._]*\\b",e:hljs.IMMEDIATE_RE,r:0},{cN:"operator",b:"<\\-(?!\\s*\\d)",e:hljs.IMMEDIATE_RE,r:2},{cN:"operator",b:"\\->|<\\-",e:hljs.IMMEDIATE_RE,r:1},{cN:"operator",b:"%%|~",e:hljs.IMMEDIATE_RE},{cN:"operator",b:">=|<=|==|!=|\\|\\||&&|=|\\+|\\-|\\*|/|\\^|>|<|!|&|\\||\\$|:",e:hljs.IMMEDIATE_RE,r:0},{cN:"operator",b:"%",e:"%",i:"\\n",r:1},{cN:"identifier",b:"`",e:"`",r:0},{cN:"string",b:'"',e:'"',c:[hljs.BE],r:0},{cN:"string",b:"'",e:"'",c:[hljs.BE],r:0},{cN:"paren",b:"[[({\\])}]",e:hljs.IMMEDIATE_RE,r:0}]}};
hljs.initHighlightingOnLoad();
</script>


</head>

<body>
<h1>Rでスネークゲーム</h1>

<p>R Advent Calendar 2012 の16日目の記事です。
昨年は <a href="http://mokjpn.blogspot.jp/2011/12/r-advent-calendar-rrthe-stanford-3d.html">Rで年賀状</a>というのをやりました。今年も何か、干支にちなんだものを&hellip;来年は巳年、ヘビですね。ヘビといえばスネーク。スネークといえば、最近は段ボールに隠れたり光学迷彩を使っていろいろなところに潜入する人が有名ですが、私にとっては<a href="http://ja.wikipedia.org/wiki/%E3%83%98%E3%83%93%E3%82%B2%E3%83%BC%E3%83%A0">スネークゲーム</a>かなぁ、ということで、Rでスネークゲームを作ってみました。</p>

<p>まずは遊んでみてください。</p>

<pre><code class="r">source(&quot;http://www.okada.jp.org/~mokada/Rsnake.R&quot;)
</code></pre>

<pre><code class="r">rsnake()
</code></pre>

<p>ただし、MacOS Xの場合は、Quartzデバイスではキー入力をリアルタイムにとれないので、X11デバイスをアクティブにしておく必要があります。WindowsでもRstudioから実行している場合は同様です。</p>

<pre><code class="r">x11()
</code></pre>

<p>で、X11デバイスを起動してから rsnake() を呼んでください。</p>

<p>起動すると、8x8のフィールドに四角形がいくつか連なっています□□□■。これがヘビです。ヘビに見えないかもしれませんが、断じてヘビです。グレーの■が頭で、白い□が体です。</p>

<p>ヘビは移動させることができます。h,j,k,lキーでそれぞれ左,下, 上, 右に動きます。vi と同じキーバインドですが、私はEmacs使いなので、これはviではなくNetHackのキーバインドだと思っています。
NetHackに慣れていない方のために、カーソルキーでも操作できるようにはしてあります。</p>

<p>さて、画面上に小さな赤い点が見えると思います。これはヘビのエサです。エサを食べると、ヘビは成長して体の長さが１増えます。そしてスコアも１増えます。こうして、できるだけたくさんエサを食べるのがゲームの目的です。</p>

<p>ヘビは成長していくので、フィールドの中で動くのがだんだん窮屈になってきます。上下左右どの方向にも動けなくなったらゲームオーバーです。</p>

<p>rsnake()関数は、戻り値としてスコアを返します。detailed=TRUE をつけると、総キー入力回数(Inputs)と、ヘビの動いた回数(Movements)もつけてベクタにして戻り値とします。</p>

<p>また、xおよびyをパラメータとして与えると、フィールドを広げることができます。rsnake(x=20, y=20)などとすると20x20の広いフィールドになります。</p>

<p>途中でギブアップするときは、q を押してください。</p>

<p>では、お楽しみください！。。。。。というだけではちょっとつまらないので、「オートモード」をつけてみました。rsnake()はautoパラメータに関数を与えることで、ヘビを自動操縦することができます。autoに与える関数は、たとえば名前をhebi() とすると、</p>

<pre><code class="r">hebi(snake, food, x, y)
</code></pre>

<p>として呼ばれ、&quot;h&quot;, &ldquo;j&rdquo;, &ldquo;k&rdquo;, &ldquo;l&rdquo; のうちどれから一文字を返すことが期待されます。
snakeパラメータは、ヘビの現在座標、foodはエサの座標、xとyはフィールドの大きさ(rsnake関数呼び出しの際に与えられた値と同じ、デフォルトでは8x8)です。</p>

<p>例として、単に一様分布に従う乱数を用いて進む方向を決めるヘビを作ってみましょう。下記のような関数を定義します。</p>

<pre><code class="r">random &lt;- function(s, food, x, y) {
    return(c(&quot;h&quot;, &quot;j&quot;, &quot;k&quot;, &quot;l&quot;)[floor(runif(1) * 4) + 1])
}
</code></pre>

<p>autoパラメータに与えられる関数の値は、キー入力のイベントハンドラ関数に直接渡されますから、実際には進めない方向の値を返しても問題ありません。</p>

<p>では、この関数を使ってヘビを走らせた場合のスコアを10回とってみましょう。</p>

<pre><code class="r">unlist(sapply(1:10, function(x) {
    rsnake(auto = random)
}))
</code></pre>

<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAYAAACmKP9/AAAACXBIWXMAAAsSAAALEgHS3X78AAAWVUlEQVR4nO3df4yV9b3g8c/BAccfgIAwgIjVWkVq9A/SbeuP9hoxoTX+giKNNl26118ga9vN7vamvdbUTbRNo93WtLvrbvauv1NMRasN/gHVBKzb7b25d7VrqFZJEToDHRAQmXEKfPePc2aYGYZhgJFz5sPrlZzAc57nnOczBzjvec7zjEYAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACDOSEi/l1EvB4RH0bEBxHxjxHxtxFRGcb9NEdEqd2GUyUivhIR/zuqs38QEb+JiK/G8M4PACPKD6Ia3b9GNZL/HPtD/N1h3M9HFfj/3Ot5/xwRW3ot/2iY9wUAI0IlIrZHNYaf7nX/dbX7tsXwHQV/FIGfV3u+7RFxRVRnrUTE9VH9NKJExEXDuD8AGBGaYn90z+t1/6iIWFy7jYq+cb4hIt6O6hH/byNidu0xLRHxDxHRGhF7I2JHRKyMiE/W1vcP/KiI+C8RsTYi/k3tvvMi4ldR/Zh9e0Q8EhGTBpn/hdrzLRtg3e0R8V8j4m8Oc7a/q+3/zdpj/0NEbKh9vf8cEZ/qtY/DnRcAjplXohq23RHxfER8JyKuimr0uvUO4K6oRnlrbfmV2jara8sbI+KpqMa/RMS6AZ4jIuK+2u+fjeo3GqdFxKaI6IqIR6Ma4FJ7nhMOMvt7tW1mHOJrPJzZ2qIa9+7TFh/Wvt7uj/5frz3mSOYFgGPmjIhYHhEdsT9yJapB+5vaNr0D2H3fRbE/ghER/ysiHo+IWbXlmb0ec2K/5/jXtV/XRMRJte3/be2++3rNtqJ23xcGmLv3pw9Nve7/h6iGu/u25DBnuyCq8e5enlt7zFm15X1HOC8A1MWYiPhXEfGNiPh/UQ3Vm7V1vQM4qnbfSb3u637830bEkxHxL9H3G4bmfs/RVfv19l77f7TX+v63uweYtxLVTx1KREzvdf+qfo/9+8Ocrfuag+7lkwd4DY5kXjjujDr0JsBHZFRUzzn/XVTPTf+fqF6VfkVt/Sei79FxxP4j2P4Xyz0aEf8jqp8IPBIR1w6y3021X++NiPG133cfybdGxB/63boGeI5Smzci4sZe98+NaqR7R/ZwZuv/de0bcKvDnxcAjqm2qEbty73uuzT2n4+OGPgK+P73dR+VT60tXx4HP0qeFhHPRd8fZes+J//ve+1jTlSviD/nILPfGPuvC7gm9l9F/8WI2Bn7j+APZ7ZuvdcP9PUeybwAcMz8x9gfrv8b1f9ITHcQv13bZiiBf7vXc/wiqleqd68/aYDtZ0f1U4M9Ub2a/ZyoXo3eERH/MyIeq83xXlTPmQ+kEtVz7t3Pu7l26774rzvwhztbxKEDfyTzAsAxU4nq+el/jGqsuqJ6Dn5Z7D+FNpTAXxoRv4/qRXf/FNWfUW+vrV9wkOd4uLa8ujbHZ6N64V1nRLwf1aP8Cw4x/6ionsv/l6he8d4e1R+/GxfVUwF/f4SzHSrwcYTzAgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAxdpd4DcFD/PSL21XsIgAa0MSL+U72HaHQC37j+FBHz6z0EQAN6IiJm1XuIRtdU7wE4qA0R8U/1HgKgAf2l3gOMBAI/QpRSSr1ngONFpVLx6SYj3qh6DwAADD+BB4CEBB4AEhL4Eeydd96JefPmxfjx42PChAlx3XXXxbvvvlvvsQBoAAI/gn3lK1+Jq666KrZu3Rqtra1x0UUXxeLFi+s9Vlx66aVRqVR6bnfccUe9RwI47rhStHGtiYjLuxcGuop+zJgx0d7eHuPGjYuIiM7Ozvja174WTz311LGbsp9SSpx++unx2muvxfjx4yMioqmpKZqbm+s2ExwuV9E3vD7vjwzMEfwI9tnPfjbmz58fzz33XOzatSuam5t74r5ly5a4+uqrY/z48XHWWWfF448/HhER7e3tceONN8aECRNi4sSJsWjRomhvb+95zkqlEg8//HBMmDAhNm3aFNdcc01MmjQpJk+eHIsWLYrW1tY+2/a3efPm6Orqii9+8YsxceLEWLBgQWzfvv0jfiUA6E/gR7Bnn302rrjiirjvvvti8uTJcd1118U777wTERF33nlnfPKTn4xt27bFCy+8EEuWLImdO3fGkiVLYty4cfHnP/85Nm3aFOPGjYs777yzz/P+8Y9/jDfeeCMWL14cixcvjvb29ti0aVNceOGFcdNNN/Vs99JLLx0wU1tbW3zmM5+Jn//859HW1hZTp06NpUuXfrQvBACMIGt6L5QB/PWvf+35fVtbW1m6dGmZM2dOKaWUU045pbS3t/es37p1a9mzZ0855ZRTyubNm/s8buzYsT3LEVG2b99eSill9OjRJSL63EaPHj3QKAfV2tpaTj755MN6DNRbHf69c3jWHHoTHMGPYGPHjo0tW7ZERERLS0vcf//98fvf/z4iIvbt2xejRu3/421ra4vOzs4opfT5aH3UqFGxd+/ePs/bfe583LhxsW3btiilRCkldu3aFW+99dagM/3ud7+L3/72tz3LY8aMidGjRx/dFwrAYRP4EWzevHnxrW99K3bs2BG7d++OH//4x/H5z38+IiKuvPLK+OEPfxh79+6NN954Iy699NLo7OyML3zhC/Gd73wnOjs7o6OjI7797W/H1VdfPeDzz58/P+65557o6uqK999/P26++eb43ve+17P+5ZdfPuAxbW1tcdNNN8WmTZuiq6sr7r333rjhhhs+kq8fAEaiQ35E397eXhYtWlROPfXUctJJJ5V58+aVDRs2lFJK2bhxY5k7d24ZO3ZsOfPMM8sTTzxRSinlL3/5S1mwYEEZP358Oe2008rChQv7fJTfe1c7duwoX/7yl8uECRPKxIkTy1e/+tWyY8eOAbfttm/fvvL973+/TJs2rYwdO7YsXLiw5yN/GCnq80+ew+Aj+iHwoyCN65A/Jgd8NPyYXMPzY3JD4CN6AEhI4AEgIYEHgIQEHgASaqr3AAxNo1/04yLAo+PPt7Ecb19vo2v0fx+NyhE8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEJN9R6AoSmllHrPwEfHny8w3BzBA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJNRU7wEAqK9KpVKp9wwMP0fwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPx5Hf/OY3cfnll0dzc3Occsopcdlll8WKFSv6bLNnz55oaWmJCy64IEopfdZVKpW46aabDnjem2++OXr/H0crlUo0NTUdcBvIli1b4s4774wzzzwzxowZEzNnzoxly5bFli1bIiLi/vvvj2uvvfaAx82dOzd+9rOfDWl/A62/6KKLYu3atUN85QCGz5reCwWO0muvvVamTJlSnnjiibJ79+7y3nvvlccee6ycdtpp5fHHH+/ZbtWqVWX27Nnl9NNPL6+//nqf54iIMm3atNLZ2dlz34cfflimT59eev81Hepf2c2bN5ePfexj5Rvf+EZZv3592bNnT3nnnXfKXXfdVc4+++yyZcuW8vbbb5cTTzyx7Nixo+dx27ZtK2PGjCmtra1D2l//9V1dXeXBBx8s55577pDmzK4O729Ha82hN4HGJfAMq4ULF5af/vSnB9z/3HPPlXvuuadn+bbbbis/+MEPyh133FG++93v9tk2IsrChQvL888/33PfypUry8KFC48o8EuXLi133XXXgOuWLVtWli1bVkop5VOf+lSfb0Iee+yx8rnPfW7I+xto/fvvv19Gjx49pDmzq8P729ESeEY0gWdYnX766WXz5s2DbtPV1VWmTJlSNm7cWF566aUya9assm/fvp71EVGeeuqpsnjx4p77brvttvLkk08eUeDPOOOM8tZbbw247g9/+EOZOXNmKaWUBx54oFx//fU96xYsWFAeeuihIe+v//oPPvigPPjgg+XKK68c0pzZ1eH97WgJPCOawDOsTjjhhPLhhx/2LEdEn1sppbz44otl7ty5pZRS9uzZU6ZNm9bnY/qIKDt37izTp08vXV1dZe/evWXmzJll+/btBwR+oFt/TU1Npaura8B5Ozs7e46wN2zYUE466aSyc+fO0tHRUcaOHVs2bdo05P0NtK5SqZSXXnrpCF7JfD7yd7PhJ/BD4CI7OE5MnTo11q9f37NcSolSSuzevbvnvuXLl8eqVat6LkprbW2Np59+us/zjB07NubMmRMvv/xyvPrqq3HhhRfG+PHjD9hf9/P3vvXX0tISGzZsGHDejRs3xtSpUyMi4swzz4w5c+bEr371q1i9enVcfPHFMX369MPaX+/7Ozo64oEHHoilS5ce4lWDkUvg4Thx1VVXxZNPPnnA/c8//3xERHR1dcWKFSviT3/6U08IX3zxxVi+fPkBsZw/f3784he/iBUrVsQNN9xwxDNde+218dBDDw247ic/+Ulcf/31PcuLFi2Kp59+Op599tn40pe+dMT7jIhobm6OJUuWxJtvvnlUzwNwJHxEz7B66623yqRJk8rDDz9cPvjgg7Jz587yyCOPlPPPP7+ccMIJ5YUXXigXX3xxn8d0dHSUk08+uedj+u6/ilu3bi1nnHFG+cQnPtFzXj/i8M/Bt7W1lbPPPrt885vfLOvXry979+4t69evL1//+tfLOeecU7Zs2dKzbWtrazn11FPLlClTyrvvvtvneQ61v4Ot90+r6ti/vR01H9EPgSN4OE6ce+658eKLL8ajjz4aEydOjHPOOSd+/etfx6pVq2L27NmxfPnyuOaaa/o8prm5Oa688soDPqafOHFizJ49O6ZNmxZTpkwZcH9D+Tn4lpaWePXVV6OzszMuu+yyOPHEE+OSSy6Jjo6OeOWVV2Ly5Mk9206dOjU+/elPx8c//vGYMWPGEe2vv0mTJh1yGxipKofehDpZExGXdy+M0O+ygRGg0vu/UjQy9Hl/ZGCO4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASaqr3AADUVyml1HuGwYzA/199Q3AEDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAk1FTvAQDq7dZbb41Row5+vDNjxoy4++67j+FEcPQq9R6Ag1oTEZd3L5RSSh1ngdTOqlTimUHW33z++bFu3bpjNg99VSqV/q3q8/7IwBzBA8e9mRExZ5D1kydPPlajwLBxDh4AEhJ4AEhI4AEgIYEHgIQEHgASchU9cNxbGxG3D7Lej8gxEgk8cNxbt25d7Nq166Dr725pOYbTwPAQeOC4d/7559d7BBh2zsEDQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQk31HgA40KxZs2LUqPp9/z1t2rRYvXr1Qdebb3AjfT5yEHhoQB0dHfHMM8/Ubf+XXHLJoOvNN7iRPh85CDw0oJkzZ8acOXPqtv/m5uZB15tvcCN9PnJwDh4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASaqr3AMCB1q5dG7fffnvd9r9r165B15tvcCN9PnIQeGhAK1eujPb29rrt/5Zbbhl0vfkGN9LnI4dKvQfgoNZExOXdC6WUUsdZAOqmUqn0b1Wf90cG5hw8ACQk8ACQkMADQEICDwAJCTwAJOTH5KABzZo1K0aNqt/339OmTYvVq1cfdH2jz3frrbfWdb4ZM2bE3XfffdD1jf76kYPAQwPq6OiIZ555pm77v+SSSwZd3+jz/fKXv4x77733GE1zoB/96EeDBr7RXz9yEHhoQDNnzow5c+bUbf/Nzc2Drm/0+VpaWuK88847RtMcaPLkyYOub/TXjxycgweAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEmuo9AHCgtWvXxu233163/e/atWvQ9Y0+3+uvvx4PPvjgMZrmQOvWrRt0faO/fuQg8NCAVq5cGe3t7XXb/y233DLo+kafb926dXWNWEtLy6DrG/31I4dKvQfgoNZExOXdC6WUUsdZAOqmUqn0b1Wf90cG5hw8ACQk8ACQkHPwI8QAH1EBwEE5ggeAhAQeABLyEX3juiwi/lu9hwBoQLPqPcBI4Lxu4zo/Ik6t9xAADWhzRGys9xAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADHp/8PBX1J7QjoHzcAAAAASUVORK5CYII=" alt="plot of chunk unnamed-chunk-7"/> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAYAAACmKP9/AAAACXBIWXMAAAsSAAALEgHS3X78AAAW4ElEQVR4nO3df4zV9b3n8dcZfnRQAUFxwCpW2ipag3+QXS0/2mvExNaoFX9g0LTaq0H8kW43u9ub9rpWm2g3TevaJt1du9m7emubYqtoNZpWFAPW2+29uXdLa1jt1YhSfjgqUHRwBL/7xzkzzAzDDAzCmfnweCQn8D2fc77f94w4zznf7xlIAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACAgYxK8u+TrEnyXpJ3kvxjkr9OUvsQj9OapGrcPky1JFcn+YfUZ38nyW+SfDEf7vwAMKL8l9Sj+37qkfzn7A7xf/4Qj3OwAv9fe+z3z0k299i++0M+FgCMCLUkW1KP4Vk97r+4cd9b+fBeBR+MwJ/f2N+WJOekPmstyRdSPxtRJZn1IR4PAEaE0dkd3VN63N+S5JrGrSW943xJkn9N/RX/b5Oc3nhOW5K/S7Ihya4kW5M8keRTjfW+gW9J8t+SrE7y5cZ9pyR5PPXT7FuS3JfkmAHmf6yxv5v7WVuS5L8n+av9nO1vGsd/sfHc/5hkXePj/eck/6bHMfZ3XgA4ZJ5LPWzvJvllkm8kOS/16HXpGcDtqUf5zcb2c43HrGhsv57kp6nHv0qytp99JMmdjd8vT/0bjaOTrE/SmeT+1ANcNfYzai+zv914zAmDfIz7M9vG1OPeddnivcbH23Xqf03jOUOZFwAOmY8mWZakI7sjV6UetL9qPKZnALvum5XdEUyS/53kx0lmNran93jOR/rs40uNX1clGdd4/C2N++7sMdvDjfs+18/cPc8+jO5x/9+lHu6u29L9nO201OPdtb2g8ZyTGtsfDHFeAGiKsUn+bZJ/l+SPqYfqxcZazwC2NO4b1+O+ruf/dZKfJPmX9P6GobXPPjobvy7pcfz7e6z3vd3az7y11M86VEmO73H/U32e+7f7OVvXew66to/o53MwlHnhsNMy+EOAg6Ql9WvOf5P6ten/k/q70s9prH8yvV8dJ7tfwfZ9s9z9Sf5n6mcE7kty0QDHXd/49Y4kExu/73olvyHJ/+tz6+xnH1Vj3iS5osf9C1KPdM/I7s9sfT+uD/p91P7PCwCH1MbUo3Zlj/vmZvf16KT/d8D3va/rVfnUxvb87P1V8rQkj6T3j7J1XZP/Dz2OMTv1d8TP2MvsV2T3+wIuzO530X8+ybbsfgW/P7N16bne38c7lHkB4JD5T9kdrv+b+l8S0xXErzcesy+B/9ce+/hF6u9U71of18/jT0/9rMHO1N/NPiP1d6N3JPlfSf6+McfbqV8z708t9WvuXfvd1Lh1vfmvK/D7O1syeOCHMi8AHDK11K9P/2PqsepM/Rr8zdl9CW1fAj83yR9Sf9PdP6X+M+rtjfVL97KPexvbKxpzfDr1N97tSPKX1F/lnzbI/C2pX8v/l9Tf8d6e+o/fTUj9UsDfDnG2wQKfIc4LAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAsO9qzR6AvfpRkg+aPQTAMPR6km81e4jhTuCHr1eTLGz2EADD0ANJZjZ7iOFudLMHYK/WJfmnZg8BMAy90ewBRgKBHyGqqqqaPQMcLmq1mrObjHgtzR4AAPjwCTwAFEjgAaBAAj+Cvfzyyzn//PMzceLETJo0KRdffHFee+21Zo8FwDAg8CPY1VdfnfPOOy9vvvlmNmzYkFmzZuWaa65p9liZO3duarVa9+2GG25o9kgAhx3vFB2+ViWZ37XR37vox44dm/b29kyYMCFJsmPHjlx77bX56U9/euim7KOqqhx77LH5/e9/n4kTJyZJRo8endbW1qbNBPvLu+iHvV5fH+mfV/Aj2Kc//eksXLgwjzzySLZv357W1tbuuG/evDkXXHBBJk6cmJNOOik//vGPkyTt7e254oorMmnSpEyePDmLFi1Ke3t79z5rtVruvffeTJo0KevXr8+FF16YY445JlOmTMmiRYuyYcOGXo/ta9OmTens7MznP//5TJ48OZdeemm2bNlykD8TAPQl8CPY8uXLc8455+TOO+/MlClTcvHFF+fll19Oktx000351Kc+lbfeeiuPPfZYli5dmm3btmXp0qWZMGFC/vznP2f9+vWZMGFCbrrppl77/dOf/pQXXngh11xzTa655pq0t7dn/fr1OeOMM7J48eLuxz3zzDN7zLRx48acffbZ+dnPfpaNGzdm6tSpufHGGw/uJwIARpBVPTeqfrz//vvdv9+4cWN14403VrNnz66qqqqOPPLIqr29vXv9zTffrHbu3FkdeeSR1aZNm3o9b/z48d3bSaotW7ZUVVVVY8aMqZL0uo0ZM6a/UfZqw4YN1RFHHLFfz4Fma8L/7+yfVYM/BK/gR7Dx48dn8+bNSZK2trbcdddd+cMf/pAk+eCDD9LSsvs/78aNG7Njx45UVdXr1HpLS0t27drVa79d184nTJiQt956K1VVpaqqbN++PS+99NKAM/3ud7/Lb3/72+7tsWPHZsyYMQf2gQKw3wR+BDv//PPzta99LVu3bs27776be+65J5/97GeTJOeee26+853vZNeuXXnhhRcyd+7c7NixI5/73OfyjW98Izt27EhHR0e+/vWv54ILLuh3/wsXLsxtt92Wzs7O/OUvf8lVV12V22+/vXt95cqVezxn48aNWbx4cdavX5/Ozs7ccccdueSSSw7Kxw8AI9Ggp+jb29urRYsWVUcddVQ1bty46vzzz6/WrVtXVVVVvf7669WCBQuq8ePHVyeeeGL1wAMPVFVVVW+88UZ16aWXVhMnTqyOPvro6vLLL+91Kr/nobZu3VpdeeWV1aRJk6rJkydXX/ziF6utW7f2+9guH3zwQfXtb3+7mjZtWjV+/Pjq8ssv7z7lDyNFc/6XZz84Rb8P/CjI8DXoj8kBB4cfkxv2/JjcPnCKHgAKJPAAUCCBB4ACCTwAFGh0swcAhj9vOoORxyt4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQKObPQAw/FVVVTV7hoHUarVas2eA4cYreAAokMADQIGcoh+hZs6cmZaW5n1/Nm3atKxYsWKv68N9PoDSCfwI1dHRkYceeqhpx58zZ86A68N9PoDSCfwINX369MyePbtpx29tbR1wfbjPB1A61+ABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACjS62QMwNKtXr86SJUuadvzt27cPuD7c5wMoncCPUE888UTa29ubdvzrrrtuwPXhPh9A6WrNHoC9WpVkftdGVVVVE2eBYa1Wq/ladnjp9fWR/rkGDwAFEngAKJBr8MCI5xJW2VyCGRqv4AGgQAIPAAVyin6Euv7669PS0rzvz0444YTceuute12fOXNmU+ebNm1aVqxYsdf14T4fwIES+BHq0UcfzR133NG04999990DBr6joyMPPfTQIZyotzlz5gy4PtznAzhQAj9CtbW15ZRTTmna8adMmTLg+vTp0zN79uxDNM2eWltbB1wf7vMBHCjX4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKNLrZAzA0a9asyfe+972mHX/t2rUDrq9evTpLliw5RNPsafv27QOuD/f5AA6UwI9Qa9eubWok2traBlx/4okn0t7efoim2dN111034Ppwnw/gQNWaPQB7tSrJ/K6NqqqqJs4C0DS1Wq1vq3p9faR/rsEDQIEEHgAKJPAAUCCBB4ACCTwAFMiPyY1QM2fOTEtL874/mzZtWlasWLHXdfMNbLD5AA6UwI9QHR0deeihh5p2/Dlz5gy4br6BDTYfwIES+BFq+vTpmT17dtOO39raOuC6+QY22HwAB8o1eAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwMNh5De/+U3mz5+f1tbWHHnkkZk3b14efvjhXo/ZuXNn2tractppp6Wqql5rtVotixcv3mO/V111VWq1Wq/HjR49eo9bfzZv3pybbropJ554YsaOHZvp06fn5ptvzubNm5Mkd911Vy666KI9nrdgwYL88Ic/3Kfj9bc+a9asrF69eh8/czDyCDwcJtasWZNLLrkkS5cuzdtvv53169fnhhtuyJe//OU88MAD3Y979tlnc+yxx6a9vT1//OMf99jPypUr895773Vvd3Z2ZuXKlXs8bufOnXvc+tq8eXPOOuusjB07NqtWrUpHR0eeffbZjBo1KmeffXbeeOONLFq0KL/61a+ybdu27ue9/fbbWbVqVRYuXLjPx+t5f0dHR6699tpce+21+/U5hJFE4OEw8a1vfSu33XZbFi9enHHjxuXoo4/O1Vdfnfvuuy8vvfRS9+OWLVuWL33pS7nsssvy4IMP7rGfefPm5de//nX39tNPP525c+cOaabbb789F110Ue6+++587GMfy6hRo3LyySfnnnvuyQUXXJA77rgjM2bMyKxZs/LLX/6y+3mPP/54zj777EydOnVIxx0zZkyuv/76vPrqq0N6PowEAg+HiWeeeSaXXXbZHvdfdNFF+eY3v5kkef/997N8+fJcddVVWbRoUZYtW7bHafqFCxfmF7/4Rff2ww8/nEsuuWRIMz3yyCO55ZZb+l275ZZb8uijjyZJrrzyyvz85z/vXlu+fHkuv/zyIR0zSd5999386Ec/ymc+85kh7wNgqFb13Kj6mDdvXt+7DqkJEyYMuG6+gQ0238EwatSo6r333uveTtLrVlVV9eSTT1YLFiyoqqqqdu7cWU2bNq1as2ZNr+ds27atOv7446vOzs5q165d1fTp06stW7ZUPf+Y9t13z2P0NHr06Kqzs7PfeXfs2FGNGTOmqqqqWrduXTVu3Lhq27ZtVUdHRzV+/Phq/fr1+3y8/tZqtVr1zDPPDOEzyaE22NdH+ucVPBwmpk6dmldeeaV7u6qqVFWVd999t/u+ZcuW5amnnup+U9qGDRv2OE0/fvz4zJ49OytXrszzzz+fM844IxMnTtzjeF3773nrq62tLevWret33tdff737FPyJJ56Y2bNn5/HHH8+KFSty5pln5vjjj9+v4/W8v6OjI9/97ndz4403DvJZg5FL4OEwcd555+UnP/nJHvd3Xdvu7OzMww8/nFdffbU7hE8++eSAp+kP5PR8Ur888IMf/KDfte9///v5whe+0L29aNGiPPjgg1m+fHm/lxr2R2tra5YuXZoXX3zxgPYDMBRO0R+AkT7fwfDSSy9VxxxzTHXvvfdW77zzTrVt27bqvvvuq0499dRq1KhR1WOPPVadeeaZvZ7T0dFRHXHEEd2n6bv+KL755pvVRz/60eqTn/xktWnTpl5rfX8/kI0bN1Ynn3xy9dWvfrV65ZVXql27dlWvvPJK9ZWvfKWaMWNGtXnz5u7HbtiwoTrqqKOq4447rnrttdd67Wew4+1tfV/npLkG+/pI/7yCh8PEJz7xiTz55JO5//77M3ny5MyYMSNPP/10nnrqqZx++ulZtmxZLrzwwl7PaW1tzbnnnrvHafrJkyfn9NNPz7Rp03Lcccf1e7x9+Tn4tra2PP/889mxY0fmzZuXj3zkI5kzZ046Ojry3HPPZcqUKd2PnTp1as4666x8/OMfzwknnDCk4/V1zDHHDPoYGKlqgz+EJlmVZH7XRt/vYufPn59Vq5r3TezEiROzdevWva6bb2CDzQfsVuv5tyjV9fr6SP+8ggeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFGjwv6yZYWn16tVZsmRJ046/ffv2AdfNN7DB5gM4UAI/Qj3xxBNpb29v2vGvu+66AdfNN7DB5gM4UP6xmeFrwH9sBuBw4R+bGRrX4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAo5s9AIfG9ddfn5aWvX8/d8IJJ+TWW289hBMBcDDVmj0Ae7UqyfyujaqqqgPZ2Um1Wh4aYP2qU0/N2rVrD+QQAAdFrVbr26peXx/pn1fwh4npSWYPsD5lypRDNQoAh4Br8ABQIIEHgAIJPAAUSOABoEACDwAF8i76w8TqJEsGWPcjcgBlEfjDxNq1a7N9+/a9rt/a1nYIpwHgYBP4w8Spp57a7BEAOIRcgweAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgfxzsRwWarVardkzDKSqqqrZMwBl8QoeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAo0OhmDwCHQlVVVbNnADiUvIIHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKNLrZA7BvarVardkzADByeAUPAAUSeAAokFP0w9e8JP+j2UMADEMzmz3ASOC67vB1apKjmj0EwDC0KcnrzR4CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADg8PT/Ae9ZtJqaf8ukAAAAAElFTkSuQmCC" alt="plot of chunk unnamed-chunk-7"/> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAYAAACmKP9/AAAACXBIWXMAAAsSAAALEgHS3X78AAAWKUlEQVR4nO3df4xW9b3g8c8zgAwiv4UBRFQqitRgWrLbFjG9RkzwGtGiFaOmJXu1BnX7Y7Pd3rTXNdpEs2nUVZPuLrvde7WKKaYiLUaTgtqIdXdzb+7d0jVUqkSE8qMDCkUGp+B3/3ieGWaGYWZAnGfmw+uVPIHznPOc85kB5s055xmIAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAejIkIv5dRGyIiI8i4sOI+MeI+JuIqJzE4zRGRKk9TqZKRNwaEf8rqrN/GBG/iYivxcmdHwAGlf8U1ej+JaqR/Oc4EuL/eBKP82kF/j932O8fI2JXh+VHTvKxAGBQqETEB1GN4Rc6PH9t7bk9cfLOgj+NwC+s7e+DiLg8qrNWIuK6qF6NKBEx5yQeDwAGhaFxJLoXdHi+ISKW1h4N0TnOX4mIt6N6xv+/I2J27TVNEfH3EbE9Ig5HxN6IeDEiPltb3zXwDRHxXyJifUT8m9pzF0TEC1G9zP5BRDwRERN6mH9NbX93d7Pujoj4rxHxV8c529/Wjv9W7bXfjYgttY/3nyPiX3U4xvHOCwD95vWohu1ARPwyIn4QEVdGNXptOgZwf1SjvLu2/Hptm3W15a0R8UxU418iYmM3+4iIeKD28+ej+heNsRGxLSJaI+LJqAa41PYz5Bizv1/bZlovH+PxzLYjqnFvu23xUe3jbbv0v6H2mhOZFwD6zVkRsTIiWuJI5EpUg/ZXtW06BrDtuTlxJIIREf8QEU9FxKza8vQOrxneZR9fr/34WkSMqG3/b2vPPdBhtlW1567qZu6OVx+Gdnj+76Ma7rbHsuOc7aKoxrtteUHtNefUlj8+wXkBoC5Oi4h/HRHfjoj/F9VQvVVb1zGADbXnRnR4ru31fxMRKyLiX6LzXxgau+yjtfbjHR2O/2SH9V0f93QzbyWqVx1KREzt8PzaLq/9u+Ocre09B23Lp3fzOTiReeGU09D7JsCnpCGq95z/Nqr3pv9PVN+Vfnlt/czofHYcceQMtuub5Z6MiP8R1SsCT0TEoh6Ou6324/0RMab287Yz+e0R8fsuj9Zu9lFq80ZE3Njh+QVRjXTHyB7PbF0/ro+73er45wWAfrUjqlG7qcNzl8aR+9ER3b8DvutzbWflk2vLl8Wxz5KnRMTq6PytbG335P99h2PMjeo74mccY/Yb48j7Aq6JI++i/+uI2BdHzuCPZ7Y2Hdd39/GeyLwA0G/+QxwJ1/+N6j8S0xbE79e26Uvg3+6wj59H9Z3qbetHdLP97KheNTgU1Xezz4jqu9FbIuJ/RsRPa3O8H9V75t2pRPWee9t+d9YebW/+awv88c4W0XvgT2ReAOg3lajen/7HqMaqNar34O+OI7fQ+hL4SyPid1F9090/RfV71Jtr668/xj6W15bX1eb4UlTfeHcwIv4c1bP8i3qZvyGq9/L/JarveG+O6rffjY7qrYC/O8HZegt8nOC8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAH1XqfcAHNN/j4iP6z0EwAC0NSJ+WO8hBjqBH7jejYjF9R4CYAB6OiJm1XuIgW5ovQfgmLZExD/VewiAAehP9R5gMBD4QaKUUuo9A5wqKpWKq5sMeg31HgAAOPkEHgASEngASEjgB7F33nknFi5cGGPGjIlx48bFtddeG++99169xwJgABD4QezWW2+NK6+8Mnbv3h3bt2+POXPmxNKlS+s9VqxevTpmzpwZ48aNi5kzZ8aaNWvqPRIADBivdVwo3Rg2bFjZu3dv+3JLS0u56aabutu037S2tpaxY8eW119/vZRSyquvvlomTpxY15ngeNXrDz199lrvm+AMfhD70pe+FIsXL47Vq1fH/v37o7GxMZ555pmIiNi1a1dcffXVMWbMmDjnnHPiqaeeioiI5ubmuPHGG2PcuHExfvz4WLJkSTQ3N7fvs1KpxPLly2PcuHGxbdu2uOaaa2LChAkxceLEWLJkSWzfvr3Ttl19/PHHsWLFipg3b14cOnQoDh06FCNHjvyUPxMAdCXwg9jzzz8fl19+eTzwwAMxceLEuPbaa+Odd96JiIi77rorPvvZz8aePXtizZo1sWzZsti3b18sW7YsRo8eHX/84x9j27ZtMXr06Ljrrrs67fcPf/hDvPnmm7F06dJYunRpNDc3x7Zt2+Liiy+Om2++uX27V1555aiZhg8fHldddVXs378/hg0bFgsWLIjHHnvs0/1EAMAg0usl+r/85S/tP9+xY0e58847y9y5c0sppYwcObI0Nze3r9+9e3c5dOhQGTlyZNm5c2en140aNarjpcnywQcflFKqtwAiotNj2LBhfb7M2draWn7yk5+Uiy66qM+vgYGgDn/eOT4u0feBM/hBbNSoUbFr166IiGhqaooHH3wwfve730VE9VJ5Q8ORX94dO3bEwYMHo5TS6dJ6Q0NDHD58uNN+x4wZExERo0ePjj179kQpJUopsX///ti0aVOPM23atCm++93vRkTEsGHDYsmSJb2+BoCTT+AHsYULF8b3vve92Lt3bxw4cCAeffTR+PKXvxwREVdccUX86Ec/isOHD8ebb74Zl156aRw8eDCuuuqq+MEPfhAHDx6MlpaW+P73vx9XX311t/tfvHhx3HvvvdHa2hp//vOf45Zbbon77ruvff2rr7561GsmTZoUy5cvj1//+tdRSomf/exn8bnPfe5T+fgBYDDq9RJ9c3NzWbJkSTnjjDPKiBEjysKFC8uWLVtKKaVs3bq1LFiwoIwaNaqcffbZ5emnny6llPKnP/2pXH/99WXMmDFl7Nix5atf/WqnS/kdD7V3795y0003lXHjxpXx48eXr33ta53etX+Mscrq1avLzJkzy+mnn14+//nPlw0bNnyy66XQz+rzR57j4BJ9H/gPFQau1yLisrYFX3Sg//jPZga8Tl8f6Z5L9ACQkMADQEICDwAJCTwAJDS03gMAfFID/U1x3iT7yQz0X9+Byhk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEJD6z0AwCdVSin1ngEGGmfwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQv67WOCUd/vtt0dDw7HPd6ZNmxb33HNPP04En1yl3gNwTK9FxGVtC/6/a/j0nFOpxHM9rL/lwgtj48aN/TYPnVUqla6t6vT1ke45gwdOedMjYm4P6ydOnNhfo8BJ4x48ACQk8ACQkMADQEICDwAJCTwAJORd9MApb31E3NHDet8ix2Ak8MApb+PGjbF///5jrr+nqakfp4GTQ+CBU96FF15Y7xHgpHMPHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4OIX85je/icsuuywaGxtj5MiRMX/+/Fi1alWnbQ4dOhRNTU1x0UUXRSml07pKpRI333zzUfu95ZZbolKpdNpu6NChRz26s2vXrrjrrrvi7LPPjtNOOy2mT58ed999d+zatSsiIh588MFYtGjRUa9bsGBB/PjHP+7T8bpbP2fOnFi/fn0fP3MAJ89rHRcKfEK//e1vy6RJk8rTTz9dDhw4UN5///3y05/+tIwdO7Y89dRT7dutXbu2zJ49u5x55pllw4YNnfYREWXKlCnl4MGD7c999NFHZerUqaXjb9O+/pbduXNnOffcc8u3v/3tsnnz5nLo0KHyzjvvlG9+85vlvPPOK7t27Spvv/12GT58eNm7d2/76/bs2VNOO+20sn379j4dr+v61tbW8vDDD5fzzz+/T3NSX719faR7zuDhFPHDH/4w7r333rj55ptjxIgRMXbs2Lj11lvjiSeeiE2bNrVvt3Llyvj6178eN9xwQzz77LNH7Wf+/Pnxq1/9qn355ZdfjksvvfSEZrrvvvti0aJF8cgjj8S5554bQ4YMifPOOy8effTRuPrqq+P++++PGTNmxJw5c+KXv/xl++teeOGF+OIXvxiTJ08+oeMOGzYsbr/99nj33XdP6PUAn4QzeE6qM888s+zcubPHbVpbW8ukSZPK1q1byyuvvFJmzZpVPv744/b1EVGeeeaZsnTp0vbnvvGNb5QVK1ac0Bn8WWedVTZt2tTtut///vdl+vTppZRSHnrooXLddde1r7v++uvL448/3ufjdV3/4YcflocffrhcccUVfZqT+urt6yMMNgLPSTVkyJDy0UcftS9HRKdHKaW89NJLZcGCBaWUUg4dOlSmTJnS6TJ9RJR9+/aVqVOnltbW1nL48OEyffr08sEHHxwV+O4eXQ0dOrS0trZ2O+/BgwfLsGHDSimlbNmypYwYMaLs27evtLS0lFGjRpVt27b1+XjdratUKuWVV145gc8k/a23r490zyV6OEVMnjw5Nm/e3L5cSolSShw4cKD9uZUrV8batWvb35S2ffv2oy7Tjxo1KubOnRuvvvpqvPHGG3HxxRfHmDFjjjpe2/47PrpqamqKLVu2dDvv1q1b2y/Bn3322TF37tx44YUXYt26dXHJJZfE1KlTj+t4HZ9vaWmJhx56KO68885ePmsweAk8nCKuvPLKWLFixVHPt93bbm1tjVWrVsW7777bHsKXXnopVq5ceVQsFy9eHD//+c9j1apV8ZWvfOWEZ1q0aFE8/vjj3a577LHH4rrrrmtfXrJkSTz77LPx/PPPxw033HDCx4yIaGxsjGXLlsVbb731ifYDcCJcouek2rRpU5kwYUJZvnx5+fDDD8u+ffvKE088US688MIyZMiQsmbNmnLJJZd0ek1LS0s5/fTT2y/Tt/1W3L17dznrrLPKzJkz2+/rRxz/PfgdO3aU8847r3znO98pmzdvLocPHy6bN28u3/rWt8qMGTPKrl272rfdvn17OeOMM8qkSZPKe++912k/vR3vWOv90Rocevv6SPecwcMp4vzzz4+XXnopnnzyyRg/fnzMmDEjXn755Vi7dm3Mnj07Vq5cGddcc02n1zQ2NsYVV1xx1GX68ePHx+zZs2PKlCkxadKkbo/Xl++Db2pqijfeeCMOHjwY8+fPj+HDh8e8efOipaUlXn/99Zg4cWL7tpMnT44vfOEL8ZnPfCamTZt2QsfrasKECb1uA4NVpfdNqJPXIuKytoVj/C0WIL1Kx39FqarT10e65wweABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABIaGi9BwA41cyaNSsaGup3fjVlypRYt27dMdcP9PnoG4EH6GctLS3x3HPP1e348+bN63H9QJ+PvhF4gH42ffr0mDt3bt2O39jY2OP6gT4ffeMePAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQ0tN4DAJxq1q9fH3fccUfdjr9///4e1w/0+egbgQfoZy+++GI0NzfX7fi33XZbj+sH+nz0TaXeA3BMr0XEZW0LpZRSx1kA6qZSqXRtVaevj3TPPXgASEjgASAhgQeAhAQeABISeABIyLfJAcdt1qxZ0dBQv/ODKVOmxLp164653nw9622+22+/va7zTZs2Le655566HT8LgQeOW0tLSzz33HN1O/68efN6XG++nvU23y9+8Yu4//77+2maoz3yyCMCfxIIPHDcpk+fHnPnzq3b8RsbG3tcb76e9TZfU1NTXHDBBf00zdEmTpxYt2Nn4h48ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJDS03gMAg8/69evjjjvuqNvx9+/f3+N68/Wst/k2bNgQDz/8cD9Nc7SNGzfW7diZCDxw3F588cVobm6u2/Fvu+22Htebr2e9zbdx48Ze/xLwaWpqaqrbsTOp1HsAjum1iLisbaGUUuo4C0DdVCqVrq3q9PWR7rkHDwAJCTwAJOQe/CDRzSUqADgmZ/AAkJDAA0BCLtEPXPMj4r/VewiAAWhWvQcYDNzXHbgujIgz6j0EwAC0MyK21nsIAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACAU9P/B+nt1om7lArIAAAAAElFTkSuQmCC" alt="plot of chunk unnamed-chunk-7"/> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAYAAACmKP9/AAAACXBIWXMAAAsSAAALEgHS3X78AAAWLUlEQVR4nO3df4xX9b3n8deXXw4qICgM+AOVWqXU4B9kt1Y0vY2Y2Bp//6BB09K9Wqt2+2Ozu71pr2u0Se2mUdea9O66m3uv1h8ppmpbjSYVNZHW7fbe3LvaNVSqRITyo6MFis44BT77x3dmmIFh+CnfmQ+PR/KNnO8533PeMwPz5HvOGUwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIChjE7yH5K8muSDJO8l+ackf52kcRCP05ak9DwOpkaS65L87zRnfy/Jr5J8Pgd3fgAYUf5rmtH9S5qR/JfsCPF/OYjH+bAC/9/67fcPSTb0W77nIB8LAEaERpKNacbwE/2ev7TnuXdz8N4FfxiBv7BnfxuTfDrNWRtJLkvzbERJMvcgHg8ARoQx2RHd0/s9PyrJ4p7HqAyM8+VJ3kjzHf+vk8zpeU17kn9IsjbJtiSbkjyT5OM963cO/Kgkf5dkWZJ/1/Pc6UmeTvM0+8YkDyQ5doj5n+rZ31cGWXdjkv+e5K/2cba/6Tn+6z2v/U9JVvV8vP+S5N/0O8a+zgsAh8wv0wzb+0l+nuTbSS5IM3q9+gdwS5pRfqdn+Zc92yztWV6d5NE041+SLB9kH0ny3Z5fP5nmXzSOSbImSXeSB9MMcOnZz+jdzP6nnm1O3MPHuC+zrUsz7r2XLT7o+Xh7T/2/2vOa/ZkXAA6ZE5IsSdKZHZEraQbtr3q26R/A3ufmZkcEk+QfkzyUZHbP8sx+rzlip318oee/LyUZ37P9v+957rv9Znui57nPDDJ3/7MPY/o9/w9phrv3cdM+zvaxNOPdu7yg5zUn9yxv3895AaAlxiX5t0m+nuT/pRmq13vW9Q/gqJ7nxvd7rvf1f53kkST/moF/YWjbaR/dPf+9sd/xH+y3fufHrYPM20jzrENJcny/55/b6bV/u4+z9d5z0Lt85CCfg/2ZFw47o/a8CfAhGZXmNee/SfPa9P9J8670T/es/2gGvjtOdryD3flmuQeT/K80zwg8kOSSIY67pue/dySZ1PPr3nfya5P8bqdH9yD7KD3zJsk1/Z5fkGak+0d2X2bb+ePaPuhW+z4vABxS69KM2uf6PTc/O65HJ4PfAb/zc73vyqf3LJ+X3b9LnpHkpxn4o2y91+T/Y79jzEvzjvhZu5n9muy4L+Di7LiL/rNJNmfHO/h9ma1X//WDfbz7My8AHDL/OTvC9X/T/EdieoP4rZ5t9ibwb/Tbx0/SvFO9d/34Qbafk+ZZg61p3s0+K8270TuT/H2SH/XM8ac0r5kPppHmNffe/a7vefTe/Ncb+H2dLdlz4PdnXgA4ZBppXp/+pzRj1Z3mNfivZMcltL0J/Pwkv03zprt/TvNn1Dt61l+5m33c37O8tGeOT6Z5411Xkj+n+S7/Y3uYf1Sa1/L/Nc073jvS/PG7iWleCvjb/ZxtT4HPfs4LAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAsPcarR6A3fqfSba3egiAYWh1ku+0eojhTuCHr7eSXNHqIQCGoYeTzG71EMPdmFYPwG6tSvLPrR4CYBj6Y6sHGAkEfoQopZRWzwCHi0aj4ewmI96oVg8AABx8Ag8AFRJ4AKiQwI9gb775Zi688MJMmjQpkydPzqWXXpq333671WMBMAwI/Ah23XXX5YILLsg777yTtWvXZu7cuVm8eHGrx8r8+fPTaDT6Hl/+8pdbPRLAYcedosPXS0nO610Y7C76cePGpaOjIxMnTkySdHV15Ytf/GIeffTRQzflTkopOe644/LKK69k0qRJSZIxY8akra2tZTPBvnIX/bA34Psjg/MOfgT75Cc/mSuuuCI//elPs2XLlrS1tfXFfcOGDbnooosyadKknHzyyXnooYeSJB0dHbnmmmsyefLkTJkyJQsXLkxHR0ffPhuNRu6///5Mnjw5a9asycUXX5xjjz02U6dOzcKFC7N27doB2+5s/fr16e7uzmc/+9lMmTIlV155ZTZu3PghfyYA2JnAj2BPPvlkPv3pT+e73/1upk6dmksvvTRvvvlmkuSWW27Jxz/+8bz77rt56qmnctNNN2Xz5s256aabMnHixPzhD3/ImjVrMnHixNxyyy0D9vv73/8+r732WhYvXpzFixeno6Mja9asyZlnnplFixb1bffCCy/sMtO6dety9tln58c//nHWrVuX6dOn5+abb/5wPxEAMIK81H+hDOIvf/lL36/XrVtXbr755jJv3rxSSilHHXVU6ejo6Fv/zjvvlK1bt5ajjjqqrF+/fsDrJkyY0LecpGzcuLGUUsrYsWNLkgGPsWPHDjbKbq1du7YceeSR+/QaaLUW/Hln37y0503wDn4EmzBhQjZs2JAkaW9vz5133pnf/va3SZLt27dn1KgdX95169alq6srpZQBp9ZHjRqVbdu2Ddhv77XziRMn5t13300pJaWUbNmyJStWrBhypt/85jf59a9/3bc8bty4jB079sA+UAD2mcCPYBdeeGG++c1vZtOmTXn//fdz77335lOf+lSS5Pzzz8/3v//9bNu2La+99lrmz5+frq6ufOYzn8m3v/3tdHV1pbOzM9/61rdy0UUXDbr/K664Irfddlu6u7vz5z//Oddee21uv/32vvUvvvjiLq9Zt25dFi1alDVr1qS7uzt33HFHLr/88g/l4weAkWiPp+g7OjrKwoULy9FHH13Gjx9fLrzwwrJq1apSSimrV68uCxYsKBMmTCgnnXRSefjhh0sppfzxj38sV155ZZk0aVI55phjytVXXz3gVH7/Q23atKl87nOfK5MnTy5Tpkwpn//858umTZsG3bbX9u3by/e+970yY8aMMmHChHL11Vf3nfKHkaI1f+TZB07R7wU/CjJ87fHH5IAPhx+TG/b8mNxecIoeACok8ABQIYEHgAoJPABUaEyrBwAYbtzUOry46XH/eAcPABUSeACokFP0I9QNN9ww4J+iPdROPPHE3HrrrbtdP3v27JbON2PGjCxdunS36803tJE+HyDwI9bPfvaz3HHHHS07/j333DNk4Ds7O/P4448fwokGOuecc4Zcb76hjfT5AIEfsdrb23P66ae37PhTp04dcv3MmTMzb968QzTNrtra2oZcb76hjfT5ANfgAaBKAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACo0ptUDsH9effXV3H333S07/vLly4dcv2zZstx4442HaJpdbdmyZcj15hvaSJ8PEPgRa/ny5S39Jtfe3j7k+meeeSYdHR2HaJpdXX/99UOuN9/QRvp8QNJo9QDs1ktJzutdKKWUFs4C0DKNRmPnVg34/sjgXIMHgAoJPABUSOABoEICDwAVEngAqJAfkxuhZs+enVGjWvf3sxkzZmTp0qW7XT/c5+PA+PrC8CfwI1RnZ2cef/zxlh3/nHPOGXL9cJ+PA+PrC8OfwI9QM2fOzLx581p2/La2tiHXD/f5ODC+vjD8uQYPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUKExrR6A/bNs2bLceOONLTv+li1bhlw/3OfjwPj6wvAn8CPUM888k46OjpYd//rrrx9y/XCfjwPj6wvDX6PVA7BbLyU5r3ehlFJaOAtAyzQajZ1bNeD7I4NzDR4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPBwGPnVr36V8847L21tbTnqqKNy7rnn5oknnhiwzdatW9Pe3p6PfexjKaUMWNdoNLJo0aJd9nvttdem0WgM2G7MmDG7PAazYcOG3HLLLTnppJMybty4zJw5M1/5yleyYcOGJMmdd96ZSy65ZJfXLViwID/84Q/36niDrZ87d26WLVu2l585gIPnpf4LBQ7QK6+8UqZNm1Yefvjh8v7775c//elP5Uc/+lE55phjykMPPdS33XPPPVfmzJlTjjvuuPLqq68O2EeSMmPGjNLV1dX33AcffFCOP/740v+36d7+ll2/fn055ZRTyte//vWycuXKsnXr1vLmm2+Wr371q+XUU08tGzZsKG+88UY54ogjyqZNm/pe9+6775Zx48aVtWvX7tXxdl7f3d1d7r777nLaaaft1Zy01p6+PzI47+DhMPGd73wnt912WxYtWpTx48fnmGOOyXXXXZcHHnggK1as6NtuyZIl+cIXvpCrrroqjz322C77Offcc/OLX/yib/n555/P/Pnz92um22+/PZdccknuueeenHLKKRk9enROPfXU3Hvvvbnoootyxx13ZNasWZk7d25+/vOf973u6aefztlnn53p06fv13HHjh2bG264IW+99dZ+vR7gQHgHz0F13HHHlfXr1w+5TXd3d5k2bVpZvXp1eeGFF8rs2bPL9u3b+9YnKY8++mhZvHhx33Nf+tKXyiOPPLJf7+BPOOGEsmLFikHX/e53vyszZ84spZRy1113lcsuu6xv3ZVXXlnuu+++vT7ezuvfe++9cvfdd5fzzz9/r+aktfb0/RFGGoHnoBo9enT54IMP+paTDHiUUsqzzz5bFixYUEopZevWrWXGjBkDTtMnKZs3by7HH3986e7uLtu2bSszZ84sGzdu3CXwgz12NmbMmNLd3T3ovF1dXWXs2LGllFJWrVpVxo8fXzZv3lw6OzvLhAkTypo1a/b6eIOtazQa5YUXXtiPzySH2p6+PzI4p+jhMDF9+vSsXLmyb7mUklJK3n///b7nlixZkueee67vprS1a9fucpp+woQJmTdvXl588cW8/PLLOfPMMzNp0qRdjte7//6PnbW3t2fVqlWDzrt69eq+U/AnnXRS5s2bl6effjpLly7NWWedleOPP36fjtf/+c7Oztx11125+eab9/BZg5FL4OEwccEFF+SRRx7Z5fnea9vd3d154okn8tZbb/WF8Nlnn82SJUt2ieUVV1yRn/zkJ3niiSdy+eWX7/dMl1xySe67775B1/3gBz/IZZdd1re8cOHCPPbYY3nyySdz1VVX7fcxk6StrS033XRTXn/99QPaD8D+cIqeg2rFihXl2GOPLffff3957733yubNm8sDDzxQzjjjjDJ69Ojy1FNPlbPOOmvAazo7O8uRRx7Zd5q+97fiO++8U0444YTy0Y9+tO+6frLv1+DXrVtXTj311PKNb3yjrFy5smzbtq2sXLmyfO1rXyuzZs0qGzZs6Nt27dq15eijjy7Tpk0rb7/99oD97Ol4u1vvj9bIsKfvjwzOO3g4TJx22ml59tln8+CDD2bKlCmZNWtWnn/++Tz33HOZM2dOlixZkosvvnjAa9ra2nL++efvcpp+ypQpmTNnTmbMmJFp06YNery9+Tn49vb2vPzyy+nq6sq5556bI444Iuecc046Ozvzy1/+MlOnTu3bdvr06fnEJz6Rj3zkIznxxBP363g7O/bYY/e4DYxUjT1vQou8lOS83oXd/C0WoHqN/v+KUtOA748Mzjt4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUaEyrB+DQuOGGGzJq1O7/PnfiiSfm1ltvPYQTAfBharR6AHbrpSTn9S6UUsqB7OzkRiOPD7H+2jPOyPLlyw/kEAAfikajsXOrBnx/ZHDewR8mZiaZN8T6qVOnHqpRADgEXIMHgAoJPABUSOABoEICDwAVEngAqJC76A8Ty5LcOMR6PyIHUBeBP0wsX748W7Zs2e36W9vbD+E0AHzYBP4wccYZZ7R6BAAOIdfgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAh/7vYEaLRaDRaPQMAI4d38ABQIYEHgAo5RT98nZvkf7R6CIBhaHarBxgJXNcdvs5IcnSrhwAYhtYnWd3qIQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA5P/x9WXrBxAWjsiAAAAABJRU5ErkJggg==" alt="plot of chunk unnamed-chunk-7"/> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAYAAACmKP9/AAAACXBIWXMAAAsSAAALEgHS3X78AAAVwklEQVR4nO3df4zV9b3n8dfBEbEFRkBBsU7rb6VG/8Bu06rtNWBC1/qjWKTRpkv2bv3J2na7u72p99bUm2iaRndb0+5qN3vX3wmmotGm/KHVRKzbvd7cu7VrilaICEUQUfAH4xT47B/nzDCM4wwMP87Mh8cjOYHz/Z7z/b5nYM6T7/d7ZkgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIChHJLkPyR5IckHSd5L8nySv07S2If7mZCktG77UiPJ15P87zRnfy/Jb5N8I/t2fgAYU36UZnT/kmYk/zk7Q/yDfbif/RX4/9pvu39OsqHf/f+yj/cFAGNCI8nbacbws/2WX9Jatin77ih4fwR+Xmt7byc5P81ZG0kuTfNsREly5j7cHwCMCR3ZGd1T+i0fl2RR6zYuu8b5K0leSfOI/3dJZrWeMyPJPyRZl2R7ks1Jfp3k0631AwM/Lsl/S7I8yb9tLTslya/SPM3+dpK7k0wbYv7HW9tbPMi6q5P89yR/tYez/U1r/y+1nvufkqxufbz/nOQz/faxp/MCwAHzbJphez/JY0luTHJBmtHr1T+A76YZ5Tdb959tPebJ1v01SR5MM/4lyR8H2UaS3NL6/SNp/kPjiCRrk/QkuSfNAJfWdg75iNnfaj3mE8N8jHsy2+tpxr33ssUHrY+399T/C63njGReADhgjk2yJMnW7IxcSTNof9V6TP8A9i47MzsjmCT/K8l9SU5r3e/q95zDBmzj37R+fSbJ4a3H//vWslv6zba0texLg8zd/+xDR7/l/5BmuHtv1+7hbKenGe/e+3Nbz/lk6/6OEc4LAG0xPsm/SvLtJP8vzVC91FrXP4DjWssO77es9/l/neSBJP+SXf/BMGHANnpav17db//39Fs/8PZ3g8zbSPOsQ0kys9/yJwY892/3cLbe9xz03v/YIJ+DkcwLB51xwz8E2E/GpXnN+W/SvDb9f9J8V/r5rfUnZ9ej42TnEezAN8vdk+R/pHlG4O4kFw+x37WtX29O0tn6fe+R/LokKwbcegbZRmnNmySX91s+N81I94/snsw28OPaMeij9nxeADigXk8zal/rt+yc7LwenQz+DviBy3qPyo9u3T8vH32UfEySR7Prt7L1XpP/j/32MTvNd8Sf8BGzX56d7wu4KDvfRf+vk2zJziP4PZmtV//1g328I5kXAA6Y/5yd4fq/af6QmN4gfr/1mN0J/Cv9tvHLNN+p3rv+8EEePyvNswbb0nw3+wlpvht9a5L/meTe1hxvpXnNfDCNNK+59253fevW++a/3sDv6WzJ8IEfybwAcMA00rw+/XyasepJ8xr84uy8hLY7gT8nyR/SfNPdP6X5PeobW+sv+4ht3NW6/2Rrjs+l+ca77iTvpHmUf/ow849L81r+v6T5jveNaX773eQ0LwX87QhnGy7wGeG8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAALuv0e4B+Ei/SLKj3UMAjEJrkvx9u4cY7QR+9Ho1yfx2DwEwCt2f5LR2DzHadbR7AD7S6iT/1O4hAEahN9o9wFgg8GNEKaW0ewY4WDQaDWc3GfPGtXsAAGDfE3gAqJDAA0CFBH4MW7lyZebNm5fOzs5MmTIll1xySV577bV2jwXAKCDwY9jXv/71XHDBBXnzzTezbt26nHnmmVm0aFG7x0qSbN++Paed5rtYANpF4Mew559/Pt/85jfT0dGRCRMm5MYbb8z06dPbPVbuvPPOzJ07NytWrGj3KAAHLYEfwz73uc9l/vz5efTRR/Puu+9mwoQJefDBB5MkGzZsyIUXXpjOzs588pOfzH333Zck2bhxYy6//PJMmTIlU6dOzcKFC7Nx48a+bTYajdx1112ZMmVK1q5dm4suuijTpk3LUUcdlYULF2bdunW7PHYwp5xySr773e/ux48cgOEI/Bj2yCOP5Pzzz88tt9ySo446KpdccklWrlyZJLn++uvz6U9/Ops2bcrjjz+ea6+9Nlu2bMm1116byZMn589//nPWrl2byZMn5/rrr99lu3/605/y4osvZtGiRVm0aFE2btyYtWvX5owzzsgVV1zR97innnpq0LnOP//8fPnLX95/HzgAw/LDHEavZ5Kc13tnsB90s23btnR0NH9W0fr163PzzTfnd7/7XZ5//vlMnDgxr776aqZNm5Yk2bRpUzo7O9PZ2ZmVK1f2ncpfv359Tj755GzZsiVJ86j87bffTmdnZ8aPH5+//OUvu+zz0EMPTU9Pz259AI1GI34+D2ORH3Qz6u3y+sjgHMGPYZMmTcqGDRuSJDNmzMitt96aP/zhD0mSHTt2ZNy4nX+8r7/+erq7u1NK2eXU+rhx47J9+/ZdttvZ2ZkkmTx5cjZt2pRSSkopeffdd/Pyyy/v7w8LgH1A4MewefPm5Xvf+142b96c999/Pz/5yU/yxS9+MUkyZ86c/PjHP8727dvz4osv5pxzzkl3d3e+9KUv5cYbb0x3d3e2bt2a73//+7nwwgsH3f78+fNz0003paenJ++8806uvPLK/PCHP+xb//TTTx+IDxMAqvJM/ztlEBs3biwLFy4sEydOLIcffniZN29eWb16dSmllDVr1pS5c+eWSZMmleOOO67cf//9pZRS3njjjXLZZZeVzs7OcsQRR5QFCxaUjRs39m2z/642b95cvva1r5UpU6aUqVOnlm984xtl8+bNgz52MMOth9GqLV/x7Ilnhn8IrjONXsNegwf2D9fgRz3X4HeDU/QAUCGBB4AKCTwAVEjgAaBCHe0egN3jTT918ybKvePrAz7METwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFChjnYPQB1KKaXdMwCwkyN4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEId7R6AOjQajUa7ZxhKKaW0ewb2n9H+5zvavz6okyN4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4BktfpHkziFuf9e+0erx29/+Nuedd14mTJiQj3/84zn33HOzdOnSXR6zbdu2zJgxI6effnoG/i+sjUYjV1xxxYe2e+WVV6b//4jaaDTS0dHxodtgNmzYkOuvvz7HHXdcxo8fn66urixevDgbNmxIktx66625+OKLP/S8uXPn5uc///lu7W+w9WeeeWaWL1++m585gH3nmXYPcCB1JeX5IW5J/rg32y+U3//+92X69Onl/vvvL++//3556623yr333luOOOKIct999/U97oknniizZs0qRx55ZHnhhRd22UaScswxx5Tu7u6+ZR988EGZOXNm6f9p3t1P+fr168unPvWp8u1vf7usWrWqbNu2raxcubLccMMN5fjjjy8bNmwor7zySjnssMPK5s2b+563adOmMn78+LJu3brd2t/A9T09PeX2228vJ5100m7Nubf2/iuEAQ6q10fqc1D9BT43KWWIW/by83FAXsVHuQULFpSf/exnH1r+6KOPlptuuqnv/lVXXVV+9KMflWuuuab84Ac/2OWxScqCBQvKY4891rfs17/+dVmwYMGIAn/dddeVG264YdB1ixcvLosXLy6llPKZz3xml3+E3HvvveULX/jCbu9vsPXvvPNOOfTQQ3drzr21918hDHBQvT5Sn4PqL7DA739HHnlkWb9+/ZCP6enpKdOnTy9r1qwpTz31VDnttNPKjh07+tYnKQ8++GBZtGhR37KrrrqqPPDAAyMK/LHHHltefvnlQdetWLGidHV1lVJKue2228qll17at+6yyy4rd9xxx27vb+D69957r9x+++1lzpw5uzXn3tr7rxAGOKheH6nPQfUXWOD3v0MOOaR88MEHfffT/Lz23UopZdmyZWXu3LmllFK2bdtWjjnmmF1O0ycpW7ZsKTNnziw9PT1l+/btpaurq7z99tsfCvxgt4E6OjpKT0/PoPN2d3f3HWGvXr26HH744WXLli1l69atZdKkSWXt2rW7vb/B1jUajfLUU0+N4DO55/bm7y6DOqheH0fKm+zgIHH00Udn1apVffdLKSml5P333+9btmTJkjzxxBN9b0pbt25dHnrooV22M2nSpMyePTtPP/10nnvuuZxxxhnp7Oz80P56t9//NtCMGTOyevXqQedds2ZNjj766CTJcccdl9mzZ+dXv/pVnnzyyZx11lmZOXPmHu2v//KtW7fmtttuy3XXXTfMZw3GLoGHg8QFF1yQBx544EPLH3vssSRJT09Pli5dmldffbUvhMuWLcuSJUs+FMv58+fnl7/8ZZYuXZqvfOUrI57p4osvzh133DHoup/+9Ke59NJL++4vXLgwDz30UB555JF89atfHfE+k2TChAm59tpr89JLL+3VdgBG4qA6BeUU/f738ssvl2nTppW77rqrvPfee2XLli3l7rvvLqeeemo55JBDyuOPP17OOuusXZ6zdevW8rGPfazvNH3vp/LNN98sxx57bDn55JP7rutnwCnx3fH666+X448/vnznO98pq1atKtu3by+rVq0q3/rWt8oJJ5xQNmzY0PfYdevWlYkTJ5bp06eX1157bZftDLe/j1p/oP5q7PUXCAMdVK+PI+UInlFheZKrh7glOa1No1XjpJNOyrJly3LPPfdk6tSpOeGEE/Kb3/wmTzzxRGbNmpUlS5bkoosu2uU5EyZMyJw5cz50mn7q1KmZNWtWjjnmmEyfPn3Q/e3O98HPmDEjzz33XLq7u3PuuefmsMMOy+c///ls3bo1zz77bI466qi+xx599NH57Gc/mxNPPDGf+MQnRrS/gaZNmzbsY2Csagz/ENrkmSTntXuIA+jUJBOHWL8+yZqRbtxRFO3U6P9TgNgXDrbXxxEZ/p+4cGCsaPcAADVxih4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFfLfxY4R/j9zAPaEI3gAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCHe0egDo0Go1Gu2cYy0oppd0zsP/48907Xl9GxhE8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQoY52D8DuaTQajXbPwP4z2v98Syml3TMAe8YRPABUSOABoEJO0Y9df0yyo437X5dkzhDrR/t8AFUT+DGqq6vr1Icffrht+z/77LNPHGr9aJ8PoHYCP0Z1dXVl9uzZ7Ryhe6iVo30+gNq5Bg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQoY52D8DILF++PFdffXU7R5g41MrRPh9A7RrtHoCP9EyS84ZYPy/JkQdolsGsSPKPQ6wf7fOxB0oppd0zcPBqNBoDWzXc6yNxBD+WLWv3AMMY7fMBVM01eACokMADQIWcogeGNcg1UPaA9zDQDo7gAaBCAg8AFXKKfuz6RZIdbdz/miR/P8T6P6a9861LMqeN+wdoK4Efo6ZPn/7vbr755rbt/5prrlmRIQLf1dV16sMPP3wAJ9rV2WeffWLbdg4wCgj8GDVjxoyccsop7RzhjaFWdnV1Zfbs2QdqlsF0t3PnAO3mGjwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhTraPQAj88ILL+T2229v5winDbVy+fLlufrqqw/ULIOZ2M6dA7Rbo90D8JGeSXLeEOtPTXsjtj7JmiHWz0ty5AGaZTArkvxjG/cPfUoppd0zjGWNRmNgq4Z7fSSO4MeyFe0eYBjL2j0AwMHMNXgAqJDAA0CFnKIH2M8GuYYM+50jeACokMADQIWcoh+9zk1yZ7uHABiFhvw5HDS5LjR6tfv73AFGq+F+DgcAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAD7zf8HOPyL6BCMxD0AAAAASUVORK5CYII=" alt="plot of chunk unnamed-chunk-7"/> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAYAAACmKP9/AAAACXBIWXMAAAsSAAALEgHS3X78AAAW+UlEQVR4nO3df4zV9b3n8ddB0FFEBIURf4zVSkVqNC11m1JtbyMmuMYf1VKMNl2y10pV1rab7famva1pb1LTNNqtpt1d7nbv1p+ppirRRv6QaiPW/XFv7t3SNVhaSRXKj45WEB2YAp/948wMM8MAA8U5Mx8ej+QEvuf7Pef7nmE4T77f7xkmAQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA9ueoJP8+yaokO5K8neQfk/x1ksZh3E9bktJzO5waST6T5H+mOfvbSX6Z5LM5vPMDwJjynTSj++c0I/nP2RPibxzG/bxbgf9P/Z73D0k291v+3mHeFwCMCY0kb6YZww/3u//qnvveyOE7Cn43Aj+/5/neTPKJNGdtJLkmzbMRJckFh3F/ADAmjM+e6L6v3/3jkizquY3LwDh/Msnv0jzi/19JZvc8pj3JPyTZkGRXki1Jnk7y/p71gwM/Lsl/TrIyyb/tue99SX6W5mn2N5P8OMlJ+5n/qZ7nWzLEusVJ/kuSvzrI2f6mZ/+/6Xnsl5O82vPx/nOSi/rt42DnBYAR80KaYXsnyZNJvpbksjSj16t/ALelGeXXe5Zf6NlmRc/yuiQPpxn/kmT1EM+RJN/u+f0Taf5D48Qk65N0J7kvzQCXnuc5ah+z/6lnm9MP8DEezGwb04x772WLHT0fb++p/1U9jzmUeQFgxJyW5JEkXdkTuZJm0P6qZ5v+Aey974LsiWCS/I8kDySZ1bPc0e8xxwx6jn/T8+vzSY7t2f7f9dz37X6zPd5z3+VDzN3/7MP4fvf/Q5rh7r3dcpCznZdmvHuX5/U85sye5d2HOC8AtMTRSf5Vki8m+X9phuo3Pev6B3Bcz33H9ruv9/F/neShJP+Sgf9gaBv0HN09vy7ut//7+q0ffPv6EPM20jzrUJKc2u/+ZwY99m8Pcrbe9xz0Lh83xOfgUOaFI864A28CvEvGpXnN+W/SvDb9v9N8V/onetbPzMCj42TPEezgN8vdl+S/pXlG4MdJrtrPftf3/PqtJJN7ft97JL8hycuDbt1DPEfpmTdJPt3v/nlpRrp/ZA9mtsEf1+4htzr4eQFgRG1MM2rX97vvo9lzPToZ+h3wg+/rPSo/pWf5kuz7KHlGkmUZ+K1svdfk/0O/fcxJ8x3xZ+9j9k9nz/sCrsyed9H/6yRbs+cI/mBm69V//VAf76HMCwAj5j9mT7j+b5r/SUxvEL/as81wAv+7fs/x0zTfqd67/tghtp+d5lmDnWm+m/3sNN+N3pXkvye5v2eOP6V5zXwojTSvufc+76aeW++b/3oDf7CzJQcO/KHMCwAjppHm9el/TDNW3Wleg1+SPZfQhhP4jyb5dZpvuvunNL9HvbNn/XX7eI6lPcsreub4SJpvvNue5K00j/LPO8D849K8lv8vab7jvTPNb787Ic1LAX97iLMdKPA5xHkBAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAhq/R6gHYp79PsrvVQwCMQuuS/F2rhxjtBH70+n2Sa1s9BMAo9GCSWa0eYrQb3+oB2KdXk/xTq4cAGIX+2OoBxgKBHyNKKaXVM8CRotFoOLvJmDeu1QMAAIefwANAhQQeACok8GPYK6+8kvnz52fy5MmZMmVKrr766rz22mutHguAUUDgx7DPfOYzueyyy/L6669nw4YNueCCC7Jo0aJWj5Vly5Zl5syZmTJlSmbOnJmnnnqq1SMBwKjxfP+FMoQJEyaULVu29C13dXWV66+/fqhNR0x3d3c58cQTywsvvFBKKeW5554r06ZNa+lMcLBa9ZeeYXv+wJvgCH4M+8hHPpJrr702y5Yty7Zt29LW1paHH344SbJ58+ZcccUVmTx5cs4888w88MADSZLOzs58+tOfzpQpUzJ16tQsXLgwnZ2dfc/ZaDSydOnSTJkyJevXr8+VV16Zk046KdOmTcvChQuzYcOGAdsOtnv37jz00EOZO3dudu7cmZ07d2bixInv8mcCgMEEfgx74okn8olPfCLf/va3M23atFx99dV55ZVXkiS33XZb3v/+9+eNN97IU089lVtuuSVbt27NLbfckhNOOCF/+MMfsn79+pxwwgm57bbbBjzvb3/727z00ktZtGhRFi1alM7Ozqxfvz7nn39+brjhhr7tnn322b1mOuaYY3L55Zdn27ZtmTBhQubNm5d77rnn3f1EAMAYcsBT9H/+85/7fr9x48Zy6623ljlz5pRSSpk4cWLp7OzsW//666+XnTt3lokTJ5ZNmzYNeNykSZP6n5osb775ZimleQkgyYDbhAkThn2as7u7u/zoRz8q55133rAfA6NBC/6+c3Ccoh8GR/Bj2KRJk7J58+YkSXt7e+688878+te/TtI8VT5u3J4/3o0bN2b79u0ppQw4tT5u3Ljs2rVrwPNOnjw5SXLCCSfkjTfeSCklpZRs27Yta9as2e9Ma9asyZe//OUkyYQJE7Jw4cIDPgaAw0/gx7D58+fnK1/5SrZs2ZJ33nkn3//+9/Pxj388SXLppZfmu9/9bnbt2pWXXnopH/3oR7N9+/Zcfvnl+drXvpbt27enq6srX/3qV3PFFVcM+fzXXntt7rjjjnR3d+ett97KjTfemG9+85t965977rm9HjN9+vQsXbo0v/jFL1JKyU9+8pN84AMfeFc+fgAYiw54ir6zs7MsXLiwHH/88eXYY48t8+fPL6+++moppZR169aVefPmlUmTJpUzzjijPPjgg6WUUv74xz+W6667rkyePLmceOKJZcGCBQNO5fff1ZYtW8r1119fpkyZUqZOnVo++9nPDnjX/j7GKsuWLSszZ84sxx13XPngBz9YVq1a9ZedL4UR1pq/8hwEp+iHwQ9UGL2eT3JJ74IXHRg5ftjMqDfg9ZGhOUUPABUSeACokMADQIUEHgAqNL7VAwCMNkfam1q9qbBOjuABoEICDwAVcooeYITNmjVrwH8lPdJmzJiRFStW7G+T1Ul2j9A4Q9mQ5NIW7r8KAg8wwrq6uvLYY4+1bP9z587d7/qOjo5zWznfhz70ofe2bOcVEXiAEdbR0ZE5c+a0bP9tbW37Xd/q+ZJsb+XOa+EaPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFxrd6AIAjzcqVK7N48eKW7X/btm37Xd/q+ZIc38qd16LR6gHYp+eTXNK7UEopLZwFOIyWL1+ezs7Olu3/3HPPzUUXXdS33Gg0BrdgfpKTR3SogV5O8n/2s37A6yNDE/jRS+CBETFE4Ec7gR8G1+ABoEICDwAV8iY7gCPcaL8EOAYvIYwKjuABoEICDwAVcooeYITNmjUr48a17vhqxowZWbFixT7Xf+5zn2vpfKeffnq+/vWvt2z/tXBdY/TybXJQqTPPPDOPPfZYy/Y/d+7c7NixY5/r29vb861vfWsEJxroe9/7XlavXt23PMQ1eN8mNwyO4AFGWEdHR+bMmdOy/be1te13fXt7e973vveN0DR7mzZtWsv2XRPX4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqNL7VAwAcaVauXJnFixe3bP/btm3b7/pVq1bl7rvvHqFp9rZ69eqW7bsmAg8wwp5++ul0dna2bP833XTTftevXr36gP8IeDe1t7e3bN81abR6APbp+SSX9C6UUkoLZwFomUajMbhVA14fGZpr8ABQIYEHgAoJPABUSOABoEICDwAV8m1yACNs1qxZGTeudcdXM2bMyIoVK/a5frTPx/AIPMAI6+rqymOPPday/c+dO3e/60f7fAyPwAOMsI6OjsyZM6dl+29ra9vv+tE+H8PjGjwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhca3egCAI83KlSuzePHilu1/27Zt+10/2udjeAQeYIQ9/fTT6ezsbNn+b7rppv2uH+3zMTyNVg/APj2f5JLehVJKaeEsAC3TaDQGt2rA6yNDcw0eACok8ABQIdfgx4ghTlGNKi4hAIwujuABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBhyPIL3/5y1xyySVpa2vLxIkTc/HFF+fxxx8fsM3OnTvT3t6e8847L4N/CnCj0cgNN9yw1/PeeOON6f8TjRuNRsaPH7/XbSibN2/ObbfdljPOOCNHH310Ojo6smTJkmzevDlJcuedd+aqq67a63Hz5s3LD3/4w2Htb6j1F1xwQVauXDnMzxzA4fN8qwc4GIVR71e/+lWZPn16efDBB8s777xT/vSnP5X777+/nHjiieWBBx7o2+6ZZ54ps2fPLieffHJZtWrVgOdIUmbMmFG2b9/ed9+OHTvKqaeeWvp/GQz3S2LTpk3lPe95T/niF79Y1q5dW3bu3FleeeWVcvvtt5ezzjqrbN68ufzud78rxxxzTNmyZUvf4954441y9NFHlw0bNgxrf4PXd3d3l7vvvrucc845w5qT1hriJWdMvT7CYGPqC7jVLwAc2IIFC8oPfvCDve5ftmxZueOOO/qWb7755vKd73ynfP7zny/f+MY3BmybpCxYsKA8+eSTffc9/fTTZcGCBYcU+FtvvbXcfvvtQ65bsmRJWbJkSSmllIsuumjAP0Luv//+8rGPfWzY+xtq/VtvvVUmTJgwrDlprSFecsbU6yMMNqa+gFv9AsCBnXzyyWXTpk373aa7u7tMnz69rFu3rjz77LNl1qxZZffu3X3rk5SHH364LFq0qO++m2++uTz00EOHFPjTTjutrFmzZsh1L7/8cuno6CillHLXXXeVa665pm/dddddV+69995h72/w+rfffrvcfffd5dJLLx3WnLTWEC85Y+r1EQYbU1/ArX4B4MCOOuqosmPHjr7lJANupZSyfPnyMm/evFJKKTt37iwzZswYcJo+Sdm6dWs59dRTS3d3d9m1a1fp6Ogob7755l6BH+o22Pjx40t3d/eQ827fvr3vCPvVV18txx57bNm6dWvp6uoqkyZNKuvXrx/2/oZa12g0yrPPPnsIn0lG2hAvOWPq9bFVvMkOjhCnnHJK1q5d27dcSkkpJe+8807ffY888kieeeaZvjelbdiwIY8++uiA55k0aVLmzJmT5557Li+++GLOP//8TJ48ea/99T5//9tg7e3tefXVV4ecd926dTnllFOSJGeccUbmzJmTn/3sZ1mxYkUuvPDCnHrqqQe1v/73d3V15a677sqtt956gM8ajF0CD0eIyy67LA899NBe9z/55JNJku7u7jz++OP5/e9/3xfC5cuX55FHHtkrltdee21++tOf5vHHH88nP/nJQ57pqquuyr333jvkunvuuSfXXHNN3/LChQvz6KOP5oknnsinPvWpQ95nkrS1teWWW27Jb37zm7/oeQAOxZg6BTXSp+w4eGvWrCknnXRSWbp0aXn77bfL1q1by49//ONy7rnnlqOOOqo89dRT5cILLxzwmK6urnLcccf1nabv/aN+/fXXy2mnnVZmzpzZd10/Ofhr8Bs3bixnnXVW+dKXvlTWrl1bdu3aVdauXVu+8IUvlLPPPrts3ry5b9sNGzaU448/vkyfPr289tprA57nQPvb13pfumPDEC85Y+r1sVUcwcMR4pxzzsny5ctz3333ZerUqTn77LPz85//PM8880xmz56dRx55JFdeeeWAx7S1teXSSy/d6zT91KlTM3v27MyYMSPTp08fcn/D+T749vb2vPjii9m+fXsuvvjiHHPMMZk7d266urrywgsvZNq0aX3bnnLKKfnwhz+c9773vTn99NMPaX+DnXTSSQfcBsaqxoE3oUWeT3JJq4cYrn38KxvgL9bo/78oNY2p18dWcQQPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVGh8qwegDo1Go9HqGcayUkpp9Qy8e0b73w9ff3VyBA8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRofKsHoA6llNLqGWC08veDVnAEDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AK+XGxR46/T7J7P+vXJfm7EZoFgHeZwB8hOpKbHtvP+g8lL0fgAaoh8EeIjiRz9r/JH0dkEABGhGvwAFAhgQeACgk8AFRI4AGgQgIPABXyLvojxMoki/e/yawRGQSAESHwR45ZS5Pj97N+04hNAsC7TuCPHC+3egAARo5r8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokB8XCxxQo9FotHoG4OA4ggeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACo0vtUDUIdGo9Fo9QwA7OEIHgAqJPAAUCGn6Eevi5P811YPATAKzWr1AGOB66aj17lJjm/1EACj0KYk61o9BAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAwJHp/wOEZGxbhMqypQAAAABJRU5ErkJggg==" alt="plot of chunk unnamed-chunk-7"/> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAYAAACmKP9/AAAACXBIWXMAAAsSAAALEgHS3X78AAAWe0lEQVR4nO3df4yV9b3g8c8ZEcciICAMII71RxWp0T8m3aaI9TZiQtf4oyjSaNNL9qoUZW272d2aeq1pt9GYRrutaXfX3exdfyeYglYb/APUBqx3N97cu9U1VKukCAWmAwoig1Pwu3+cM8MMDDMwg5w5H16v5ESe8zzneT7Dj/Oe53kOEgEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAzkhIj4dxHxekR8HBEfRcRrEfF3EVE5isdpjohSexxNlYj4RkT8Y1Rn/ygifhcR34yjOz8ANJT7oxrdv0Y1kv8c+0P8g6N4nE8r8P+5137/HBHtvZZ/epSPBQANoRIRH0Q1hl/s9fw1tee2x9E7C/40Aj+vtr8PIuIrUZ21EhHXRvVqRImIi47i8QCgIYyK/dE9r9fzTRGxqPZoir5x/lpEvBPVM/7/HRGzaq9piYh/iIjNEbEvInZExMqI+Hxt/YGBb4qI/xIRayPi39SeOy8ifhPVy+wfRMQjETFpgPmfr+1vaT/rFkfEf42IvznC2e6sHf+t2mv/Q0RsqH29/xwRX+h1jCOdFwCOmVeiGrbdEfFcRNwVEVdENXrdegdwV1SjvK22/Eptm9W15Y0R8VRU418iYl0/+4iIuLf242ei+o3GqRGxKSK6IuLRqAa41PZzwiFmf7+2zYxBvsYjmW1LVOPefdvi49rX233p//Xaa4YyLwAcM6dHxLKI6Iz9kStRDdrf1LbpHcDu5y6K/RGMiPhfEfF4RMysLbf2es1JB+zjb2v/XRMRJ9e2/7e15+7tNduK2nNf7Wfu3lcfRvV6/h+iGu7ux5IjnO2CqMa7e3lu7TVn1pY/GeK8AFAXoyPiX0XEdyLi/0U1VG/V1vUOYFPtuZN7Pdf9+r+LiCcj4l+i7zcMzQfso6v238W9jv9or/UHPu7uZ95KVK86lIiY3uv5VQe89u+PcLbuzxx0L3+mn5+DocwLx52mwTcBPiVNUb3nfGdU703/n6h+Kv0rtfWfi75nxxH7z2AP/LDcoxHxP6J6ReCRiLh6gONuqv33RxExvvbj7jP5zRHxhwMeXf3so9TmjYi4odfzc6Ma6d6RPZLZDvy6Pul3qyOfFwCOqS1RjdrXez13Sey/Hx3R/yfgD3yu+6x8am350jj0WfK0iHg2+v5Vtu578v++1zHaovqJ+LMPMfsNsf9zAVfF/k/R/+uI2Bn7z+CPZLZuvdf39/UOZV4AOGb+Y+wP1/+N6v8kpjuI369tcziBf6fXPn4V1U+qd68/uZ/tZ0X1qsHeqH6a/eyofhq9MyL+Z0Q8Vpvj/ajeM+9PJar33Lv3u7X26P7wX3fgj3S2iMEDP5R5AeCYqUT1/vRrUY1VV1TvwS+N/bfQDifwl0TEG1H90N0/RfXvqHfU1l93iH08XFteXZvjS1H94N2eiPgwqmf5Fwwyf1NU7+X/S1Q/8d4R1b9+Ny6qtwL+foizDRb4GOK8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIevUu8BOKT/HhGf1HsIgBFoY0T8p3oPMdIJ/Mj1p4iYX+8hAEagJyJiZr2HGOlG1XsADmlDRPxTvYcAGIH+Uu8BGoHAN4hSSqn3DHC8qFQqrm7S8JrqPQAAcPQJPAAkJPAAkJDAN7B333035s2bF+PHj48JEybENddcE++99169xwJgBBD4BvaNb3wjrrjiiti2bVts3rw5Lrrooli0aFG9x+rxxhtvxJgxY+o9BsBxSeAb2GuvvRa33HJLjBo1Kpqbm+Ouu+6KKVOm1HusiIjYsWNHLFq0KHbv3l3vUQCOSwLfwL70pS/F/Pnz49lnn41du3ZFc3NzPPXUUxER0d7eHldeeWWMHz8+zjzzzHj88ccjIqKjoyNuuOGGmDBhQkycODEWLlwYHR0dPfusVCrx8MMPx4QJE2LTpk1x1VVXxaRJk2Ly5MmxcOHC2Lx5c59t+1NKiUWLFsWdd975KX71AAxE4BvYM888E1/5ylfi3nvvjcmTJ8c111wT7777bkRE3H777fH5z38+tm/fHs8//3wsWbIkdu7cGUuWLIlx48bFn//859i0aVOMGzcubr/99j77/eMf/xhvvvlmLFq0KBYtWhQdHR2xadOmuPDCC+PGG2/s2e6ll17qd677778/zjnnnLj++us/vS8eABrUmt4LpR9//etfe368ZcuWctttt5W2trZSSiljxowpHR0dPeu3bdtW9u7dW8aMGVO2bt3a53Vjx47tWY6I8sEHH5RSSjnxxBNLRPR5nHjiif2N0uPFF18sl112Wenq6urZHzSaY//HnSO0ZvBNcAbfwMaOHRvt7e0REdHS0hL33XdfvPHGGxER8cknn0RT0/5f3i1btsSePXuilNLn0npTU1Ps27evz37Hjx8fERHjxo2L7du3RyklSimxa9euePvttwecafXq1fHb3/42Ro8e3XOcSqUSa9euHf4XDMBhE/gGNm/evPje974XO3bsiN27d8fPfvazuOyyyyIi4vLLL4+f/OQnsW/fvnjzzTfjkksuiT179sRXv/rVuOuuu2LPnj3R2dkZ3//+9+PKK6/sd//z58+Pe+65J7q6uuLDDz+Mm266KX74wx/2rH/55ZcPes2Pf/zjnm8Iuk+ESikxZ86co/8TAAANaNBL9B0dHWXhwoXllFNOKSeffHKZN29e2bBhQymllI0bN5a5c+eWsWPHljPOOKM88cQTpZRS/vKXv5TrrruujB8/vpx66qllwYIFfS7l9z7Ujh07yte//vUyYcKEMnHixPLNb36z7Nixo99tB7jUeYQXR6H+6vInniPhEv1h8A8qjFxrIuLS7gVvOnDs+MdmRrw+74/0zyV6AEhI4AEgIYEHgIQEHgASGlXvAQBGGh9qHVl86HFonMEDQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACY2q9wDA8N1yyy3R1HTo79dnzJgRd9999zGcCKi3Sr0H4JDWRMSl3QullFLHWRjhzqxUYvkA6286//xYt27dMZsHjqZKpXJgq/q8P9I/Z/CQQGtEtA2wfvLkycdqFGCEcA8eABISeABISOABICGBB4CEBB4AEvIpekhgbUQsHmC9vyIHxx+BhwTWrVsXu3btOuT6u1tajuE0wEgg8JDA+eefX+8RgBHGPXgASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAh/1wswHGuUqlU6j0DR58zeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBh+PI7373u7j00kujubk5xowZE3PmzIkVK1b02Wbv3r3R0tISF1xwQZRS+qyrVCpx4403HrTfm266KXr/i6OVSiVGjRp10KM/7e3tcfvtt8cZZ5wRo0ePjtbW1li6dGm0t7dHRMR9990XV1999UGvmzt3bvzyl788rOP1t/6iiy6KtWvXHubPHMDRs6b3QoFh+v3vf1+mTJlSnnjiibJ79+7y/vvvl8cee6yceuqp5fHHH+/ZbtWqVWXWrFnltNNOK6+//nqffUREmTZtWtmzZ0/Pcx9//HGZPn166f3b9HB/y27durV89rOfLd/5znfK+vXry969e8u7775b7rjjjnLWWWeV9vb28s4775STTjqp7Nixo+d127dvL6NHjy6bN28+rOMduL6rq6s8+OCD5dxzzz2sObOrw/vbcK0ZfBMYuQSeo2rBggXlF7/4xUHPP/vss+Wee+7pWb711lvL/fffX771rW+VH/zgB322jYiyYMGC8txzz/U8t3LlyrJgwYIhBf62224rd9xxR7/rli5dWpYuXVpKKeULX/hCn29CHnvssfLlL3/5sI/X3/oPP/ywnHjiiYc1Z3Z1eH8bLoGnoQk8R9Vpp51Wtm7dOuA2XV1dZcqUKWXjxo3lpZdeKjNnziyffPJJz/qIKE899VRZtGhRz3O33nprefLJJ4cU+NNPP728/fbb/a77wx/+UFpbW0sppTzwwAPl2muv7Vl33XXXlYceeuiwj3fg+o8++qg8+OCD5fLLLz+sObOrw/vbcAk8DU3gOapOOOGE8vHHH/csR0SfRymlvPDCC2Xu3LmllFL27t1bpk2b1ucyfUSUnTt3lunTp5eurq6yb9++0traWj744IODAt/f40CjRo0qXV1d/c67Z8+enjPsDRs2lJNPPrns3LmzdHZ2lrFjx5ZNmzYd9vH6W1epVMpLL700hJ/JfD71d7OjT+APgw/ZwXFi6tSpsX79+p7lUkqUUmL37t09zy1btixWrVrV86G0zZs3x9NPP91nP2PHjo22trZ4+eWX49VXX40LL7wwxo8ff9Dxuvff+3GglpaW2LBhQ7/zbty4MaZOnRoREWeccUa0tbXFb37zm1i9enVcfPHFMX369CM6Xu/nOzs744EHHojbbrttkJ81aFwCD8eJK664Ip588smDnn/uueciIqKrqytWrFgRf/rTn3pC+MILL8SyZcsOiuX8+fPjV7/6VaxYsSK+9rWvDXmmq6++Oh566KF+1/385z+Pa6+9tmd54cKF8fTTT8czzzwT119//ZCPGRHR3NwcS5YsibfeemtY+wEYCpfoOarefvvtMmnSpPLwww+Xjz76qOzcubM88sgj5fzzzy8nnHBCef7558vFF1/c5zWdnZ3lM5/5TM9l+u7fitu2bSunn356+dznPtdzXz/iyO/Bb9mypZx11lnlu9/9blm/fn3Zt29fWb9+ffn2t79dzj777NLe3t6z7ebNm8spp5xSpkyZUt57770++xnseIda749W1bF/exs2l+gPgzN4OE6ce+658cILL8Sjjz4aEydOjLPPPjtefPHFWLVqVcyaNSuWLVsWV111VZ/XNDc3x+WXX37QZfqJEyfGrFmzYtq0aTFlypR+j3c4fw++paUlXn311dizZ0/MmTMnTjrppJg9e3Z0dnbGK6+8EpMnT+7ZdurUqfHFL34xzjnnnJgxY8aQjnegSZMmDboNNKrK4JtQJ2si4tLuhQb9LhtoAJXe/5eixtDn/ZH+OYMHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASGhUvQcAGK6R/u+Zl1JKvWcYyEifb6T/+o5UzuABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhpV7wEAhquUUuo9A4w0zuABICGBB4CEXKJvUDNnzoympvp9fzZt2rRYvXp13Y4PwMAEvkF1dnbG8uXL63b82bNn1+3YAAxO4BtUa2trtLW11e34zc3NdTs2AINzDx4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASEngASEjgASAhgQeAhAQeABISeABISOABICGBB4CEBB4AEhJ4AEhI4AEgIYEHgIQEHgASGlXvARiatWvXxuLFi+t2/F27dtXt2AAMTuAb1MqVK6Ojo6Nux7/55pvrdmwABifwDWrevHn1HgGAEcw9eABISOABICGBB4CEBB4AEhJ4AEjIp+gb1MyZM6OpqX7fn02bNi1Wr159yPXmG9hg8wEMl8A3qM7Ozli+fHndjj979uwB15tvYIPNBzBcAt+gWltbo62trW7Hb25uHnC9+QY22HwAw+UePAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQ0qt4DMDRr166NxYsX1+34u3btGnC9+QY22HwAwyXwDWrlypXR0dFRt+PffPPNA64338AGmw9guCr1HoBDWhMRl3YvlFJKHWcBqJtKpXJgq/q8P9I/9+ABICGBB4CE3IMHBtXPJdIRxS2s4Rnpv74MjTN4AEhI4AEgIZfoG9Qtt9wSTU31+/5sxowZcffddx9y/cyZM+s637Rp02L16tWHXD/S5wMYLoFvUL/+9a/jRz/6Ud2O/9Of/nTAwHd2dsby5cuP4UR9zZ49e8D1I30+gOES+AbV0tIS5513Xt2OP3ny5AHXt7a2Rltb2zGa5mDNzc0Drh/p8wEMl3vwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkJDAA0BCAg8ACQk8ACQk8ACQkMADQEICDwAJCTwAJCTwAJCQwANAQgIPAAkJPAAkJPAAkNCoeg/A0Lz++uvx4IMP1u3469atG3D92rVrY/HixcdomoPt2rVrwPUjfT6A4RL4BrVu3bq6RqKlpWXA9StXroyOjo5jNM3Bbr755gHXj/T5AIarUu8BOKQ1EXFp90IppdRxFo5zlUplRL9X+PMxPCP917cffd4f6Z978ACQkMADQELuwTeIBryEBseMPx9wMGfwAJCQwANAQi7Rj1xzIuK/1XsIgBFoZr0HaATuW41c50fEKfUeAmAE2hoRG+s9BAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAwPHp/wPfHiilvUWwCAAAAABJRU5ErkJggg==" alt="plot of chunk unnamed-chunk-7"/> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAYAAACmKP9/AAAACXBIWXMAAAsSAAALEgHS3X78AAAXHUlEQVR4nO3df4zV9b3n8dcZwA5FQFAY8cdYta3UGv1jsqtFiLdRE7rGX2ilUdNiL5Si3XS72d029bpG29Q0jbrWpLvrbvaubtUEU7SpBtM6imFadze9uXeDNbT0akQoiuMP6OjgFPzuH2dmmIEBBop8Zz48HsmJfM/3nPN5zwjznO/3ewYSAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABgfyYk+bdJ1iX5IMl7SX6b5G+TNA7jOq1Jqv7b4dRIcmOS/53m7O8l+U2SL+fwzg8A48oP04zuX9KM5D9md4j/42Fc56MK/H8a8rp/SrJ1yPa9h3ktABgXGkneTTOG5w+5/8r++97O4TsK/igCv7D/9d5N8vk0Z20kuSrNsxFVknMP43oAMC5MzO7ofnrI/S1JlvTfWjI8zlcn+ec0j/j/T5Kz+5/TluTvk2xJsivJtiSrk3y2f/+egW9J8p+TdCX5av99n07yVJqn2d9N8mCS4/cz/5P9r/eNEfYtT/JfkvzNQc72nf71/9D/3H+fZGP/x/uPSf7FkDUOdl4AOGJ+nWbY3k/yiyS3Jrk0zegNGBrAnjSj/Fb/9q/7H9PZv70pyaNpxr9Ksn6E10iSH/T/+ok0v9E4LsnmJH1JHkozwFX/60zYx+zv9D/mlAN8jAcz2+tpxn3gssUH/R/vwKn/df3POZR5AeCIOTnJyiS92R25Ks2g/U3/Y4YGcOC+c7M7gknyP5P8NMnc/u32Ic/52B6v8ZX+/65NMrn/8f+6/74fDJnt8f77vjDC3EPPPkwccv/fpxnugduKg5ztM2nGe2D7kv7nnNa//eEhzgsAtTgmyb9M8m+S/C7NUP2hf9/QALb03zd5yH0Dz//bJI8k+acM/4ahdY/X6Ov/7/Ih6z80ZP+et9tGmLeR5lmHKslJQ+5/Zo/n/t1BzjbwnoOB7Y+P8Dk4lHnhqNNy4IcAH5GWNK85fyfNa9P/N813pX++f/+nMvzoONl9BLvnm+UeSvLf0zwj8GCSK/az7ub+/96ZZHr/rweO5Lck+f0et74RXqPqnzdJrhty/yVpRnpoZA9mtj0/rg9HfNTBzwsAR9TraUbtS0PuuzC7r0cnI78Dfs/7Bo7KT+zfXpB9HyXPSfLzDP9RtoFr8v9uyBodab4j/ox9zH5ddr8v4PLsfhf9v0qyPbuP4A9mtgFD94/08R7KvABwxPyH7A7X/0vzL4kZCOJ3+x8zmsD/85DX+Fma71Qf2D95hMefneZZg51pvpv9jDTfjd6b5H8k+V/9c7yT5jXzkTTSvOY+8Lpv9N8G3vw3EPiDnS05cOAPZV4AOGIaaV6f/m2asepL8xr8N7L7EtpoAn9hkhfTfNPdP6T5M+rd/fuv2cdrPNC/3dk/x+fSfOPdjiR/TvMo/zMHmL8lzWv5/5TmO9670/zxu2lpXgr4u0Oc7UCBzyHOCwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADB6jboHYJ/+W5IP6x4CYAzalOR7dQ8x1gn82PVqkkV1DwEwBj2cZG7dQ4x1E+segH3amOQf6h4CYAx6s+4BxgOBHyeqqqrqngGOFo1Gw9lNxr2WugcAAA4/gQeAAgk8ABRI4Mexl19+OQsXLsz06dMzY8aMXHnllXnttdfqHguAMUDgx7Ebb7wxl156ad56661s2bIl5557bpYsWVL3WINefPHFTJkype4xAI5KAj+O/fa3v82yZcsyceLEtLa25tZbb83s2bPrHitJsm3btixZsiTvv/9+3aMAHJUEfhz73Oc+l0WLFuXnP/95enp60tramkcffTRJsnXr1lx22WWZPn16TjvttPz0pz9NknR3d+e6667LjBkzMnPmzCxevDjd3d2Dr9loNPLAAw9kxowZ2bx5cy6//PIcf/zxmTVrVhYvXpwtW7YMe+xIqqrKkiVL8p3vfOcj/OgB2B+BH8eeeOKJfP7zn88PfvCDzJo1K1deeWVefvnlJMktt9ySz372s3n77bfz5JNPZsWKFdm+fXtWrFiRadOm5U9/+lM2b96cadOm5ZZbbhn2un/84x/z0ksvZcmSJVmyZEm6u7uzefPmnHPOObn++usHH/fcc8+NONcPf/jDnHnmmbn22ms/ug8eAMaptUM3qhH85S9/Gfz166+/Xt18881VR0dHVVVVNWXKlKq7u3tw/1tvvVXt3LmzmjJlSvXGG28Me97UqVMHt5NU7777blVVVTVp0qQqybDbpEmTRhpl0LPPPltddNFFVV9f3+DrwXhz5P+4c5DWHvghOIIfx6ZOnZqtW7cmSdra2nLXXXflxRdfTJJ8+OGHaWnZ/b/39ddfz44dO1JV1bBT6y0tLdm1a9ew150+fXqSZNq0aXn77bdTVVWqqkpPT082bNiw35k6Ozvz/PPP55hjjhlcp9FopKur66//gAEYNYEfxxYuXJhvf/vb2bZtW95///3cd999ueiii5IkF198cX70ox9l165deemll3LhhRdmx44d+cIXvpBbb701O3bsSG9vb7773e/msssuG/H1Fy1alNtvvz19fX3585//nBtuuCF33HHH4P41a9bs9Zzvf//7g98QDBwIVVWV+fPnH/5PAACMQwc8Rd/d3V0tXry4OvbYY6vJkydXCxcurDZu3FhVVVVt2rSpuuSSS6qpU6dWp556avXwww9XVVVVb775ZnXNNddU06dPr4477rjqi1/84rBT+UOX2rZtW/WlL32pmjFjRjVz5szqy1/+crVt27YRH7ufU50HeXIU6lfLn3gOhlP0o+AfVBi71iZZMLDhiw4cOf6xmTFv2NdHRuYUPQAUSOABoEACDwAFEngAKNDEugcAGGsO95tavWmPOjiCB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgSbWPQCMB8uWLUtLy76/Hz7llFNy2223HcGJAPavUfcA7NPaJAsGNqqqqmqc5ah3WqORVfvZf8NZZ2X9+vVHbB7Gl0aj4Wvt4TXs6yMjcwQPo9CepGM/+2fNmnWkRgEYFdfgAaBAAg8ABRJ4ACiQwANAgQQeAArkXfQwCl1Jlu9nvx+RA8YagYdRWL9+fXp6eva5/7a2tiM4DcCBCTyMwllnnVX3CAAHxTV4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABTIPxcL8BGrqqqqe4bxrNFoNOqeYTxyBA8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAo0MS6BwBgbJk7d25aWuo7/pszZ046OztrW78UAg/AML29vVm1alVt68+bN6+2tUsi8AAM097eno6OjtrWb21trW3tkrgGDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeDiK/OY3v8mCBQvS2tqaKVOmZP78+Xn88ceHPWbnzp1pa2vLZz7zmVRVNWxfo9HI9ddfv9fr3nDDDWk0GsMeN3HixL1uI9m6dWtuueWWnHrqqTnmmGPS3t6eb3zjG9m6dWuS5K677soVV1yx1/MuueSS/OQnPxnVeiPtP/fcc9PV1TXKzxyMPwIPR4l169bl6quvzooVK/LOO+9k8+bN+frXv56vfvWrefjhhwcf9/zzz+eEE05Id3d3fve73+31OmvWrMkHH3wwuN3X15c1a9bs9bidO3fuddvT1q1bc/755+eYY47J2rVr09vbm+effz4TJkzIBRdckDfffDOLFy/OL3/5y2zfvn3wee+8807Wrl2bRYsWjXq9off39vbmpptuyk033XRQn0MYTwQejhLf+973cvvtt+f666/P5MmTc9xxx+XGG2/Mgw8+mA0bNgw+buXKlfnKV76Sa6+9No899therzN//vz86le/Gtx+9tlnc+GFFx7STHfccUeuuOKK3HvvvfnEJz6RCRMm5PTTT899992Xyy67LHfeeWfOOOOMnHvuufnFL34x+LynnnoqF1xwQU488cRDWnfSpElZtmxZXn311UN6PsBfY+3QjQr+SieccEL1xhtv7PcxfX191ezZs6tNmzZVzz33XDV37tzqww8/HNyfpHr00UerJUuWDN73ta99rXrkkUeqob9NR/tb9uSTT642bNgw4r7f//73VXt7e1VVVXX33XdXV1111eC+a665prr//vtHvd6e+997773qnnvuqS6++OJRzXm0mT9/fq3rT5s2bdj2gb4+wngj8BxWEyZMqD744IPB7STDblVVVU8//XR1ySWXVFVVVTt37qzmzJlTrVu3bthztm/fXp100klVX19ftWvXrqq9vb1699139wr8SLc9TZw4serr6xtx3h07dlSTJk2qqqqqNm7cWE2ePLnavn171dvbW02dOrXavHnzqNcbaV+j0aiee+65Q/hMlk/gy+AUPRwlTjzxxLzyyiuD21VVpaqqvP/++4P3rVy5Ms8888zgm9K2bNmy12n6qVOnpqOjI2vWrMkLL7yQc845J9OnT99rvYHXH3rbU1tbWzZu3DjivJs2bRo8BX/qqaemo6MjTz31VDo7O3PeeeflpJNOOqj1ht7f29ubu+++OzfffPMBPmswfgk8HCUuvfTSPPLII3vdP3Btu6+vL48//nheffXVwRA+/fTTWbly5V6xXLRoUX72s5/l8ccfz9VXX33IM11xxRW5//77R9z34x//OFddddXg9uLFi/PYY4/liSeeyLXXXnvIaybNf298xYoV+cMf/vBXvQ7AoXCKnsNqw4YN1fHHH1898MAD1XvvvVdt3769evDBB6uzzjqrmjBhQvXkk09W55133rDn9Pb2Vh//+McHT9MP/FZ86623qpNPPrn61Kc+NXhdPzn4a/Cvv/56dfrpp1ff+ta3qldeeaXatWtX9corr1Tf/OY3qzPOOKPaunXr4GO3bNlSHXvssdXs2bOr1157bdjrHGi9fe33R2tkTtGXwRE8HCU++clP5umnn85DDz2UmTNn5owzzsizzz6bZ555JmeffXZWrlyZyy+/fNhzWltbc/HFF+91mn7mzJk5++yzM2fOnMyePXvE9Ubzc/BtbW154YUXsmPHjsyfPz8f+9jHMm/evPT29ubXv/51Zs2aNfjYE088Meeff37OPPPMnHLKKYe03p6OP/74Az4GxqvGgR9CTdYmWTCwsY/vYgEOuwULFmTt2voOkqdPn55t27YNbjeG/i1KTcO+PjIyR/AAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAId+C9rBuCo0tXVleXLl9e2fk9PT21rl0TgARhm9erV6e7urm39pUuX1rZ2SQQegGEWLlxY9wgcBq7BA0CBBB4ACiTwAFAggQeAAgk8ABTIu+gBGGbu3Llpaanv+G/OnDnp7Oysbf1SCDwAw/T29mbVqlW1rT9v3rza1i6JwAMwTHt7ezo6Ompbv7W1tba1S+IaPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBJtY9AABjS1dXV5YvX17b+j09PbWtXRKBB2CY1atXp7u7u7b1ly5dWtvaJRF4AIZZuHBh3SNwGLgGDwAFEngAKJDAA0CBBB4ACiTwAFAg76IHYJi5c+empaW+4785c+aks7OztvVLIfAADNPb25tVq1bVtv68efNqW7skAg/AMO3t7eno6Kht/dbW1trWLolr8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFmlj3AACMLV1dXVm+fHlt6/f09NS2dkkEHoBhVq9ene7u7trWX7p0aW1rl0TgARhm4cKFdY/AYeAaPAAUSOABoEACDwAFEngAKJDAA0CBvIse4AibO3duWlrqO76aM2dOOjs797l/rM/H6Ag8wBHW29ubVatW1bb+vHnz9rt/rM/H6Ag8wBHW3t6ejo6O2tZvbW3d7/6xPh+j4xo8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEm1j0AwNGmq6sry5cvr239np6e/e4f6/MxOgIPcIStXr063d3dta2/dOnS/e4f6/MxOo26B2Cf1iZZMLBRVVVV4ywAtWk0Gnu2atjXR0bmGjwAFEjgAaBAAg8ABRJ4ACiQwANAgfyYHMARNnfu3LS01Hd8NWfOnHR2du5z/7Jly2qd75RTTsltt91W2/ql8GNyY5cfk4NCnXbaaVm1alVt68+bNy8ffPDBPve3tbXlzjvvPIITDXfvvfdm/fr1g9t+TO7QOIIHOMLa29vT0dFR2/qtra373d/W1pZPf/rTR2iavc2aNau2tUviGjwAFEjgAaBAAg8ABRJ4ACiQwANAgQQeAAok8ABQIIEHgAIJPAAUSOABoEACDwAFEngAKJDAA0CBBB4ACiTwAFAggQeAAgk8ABRI4AGgQAIPAAUSeAAokMADQIEEHgAKJPAAUCCBB4ACCTwAFEjgAaBAAg8ABRJ4ACiQwANAgSbWPQDA0aarqyvLly+vbf2enp797l+3bl3uueeeIzTN3tavX1/b2iUReIAjbPXq1enu7q5t/aVLl+53//r16w/4TcBHqa2trba1S9KoewD2aW2SBQMbVVVVNc4CUJtGo7Fnq4Z9fWRkrsEDQIEEHgAK5Br8ODHCKSoA2CdH8ABQIIEHgAI5RT92zU/yX+seAmAMmlv3AOOB67pj11lJjq17CIAx6I0km+oeAgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA4Oj0/wFJLDRAOOHU7AAAAABJRU5ErkJggg==" alt="plot of chunk unnamed-chunk-7"/> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAYAAACmKP9/AAAACXBIWXMAAAsSAAALEgHS3X78AAAW2UlEQVR4nO3df4zV9b3n8dfhhx2KgKAw4o+xYlupNfrHZFeLEG+jJnSNv9CWRk2LvVCKdtPtZnfb1OsabVPTNOpak+6uu9m7ulUTTNGmGkwrioHW3U1v7t1gDS29GhGK4vgDig5Owc/+cWaGGRiGH6WcmQ+PR3IC3/P9nvN5zwjznPP9nsEEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAYztgk/zbJuiQfJHkvyW+S/G2SxhFcpy1J6b0dSY0kNyb532nO/l6SXyf5Uo7s/AAwqvwgzej+Oc1I/mP2hPg/HsF1/lqB/08DnvePSbYO2L73CK8FAKNCI8m7acbwggH3X9V739s5cq+C/xqBn9/7fO8m+WyaszaSXJ3m2YiS5LwjuB4AjArjsie6nxxw/5gki3pvYzI4ztck+ec0X/H/nyTn9D6mPcnfJ9mSZHeSbUlWJvl07/69Az8myX9OsjbJV3rv+2SSp9I8zf5ukgeTnDjM/E/2Pt/Xh9i3NMl/SfI3hzjbt3vX/33vY/99ko29H+8/JvkXA9Y41HkB4Kj5VZphez/Jz5PcmuSyNKPXZ2AAd6QZ5bd6t3/Ve8yq3u1NSR5NM/4lyfohniNJvt/7+yfS/EbjhCSbk/QkeSjNAJfe5xm7n9nf6T3mtAN8jIcy2+tpxr3vssUHvR9v36n/db2POZx5AeCoOTXJ8iTd2RO5kmbQ/qb3mIEB7LvvvOyJYJL8zyQ/STK7d7tjwGM+stdzfLn31zVJJvQe/6977/v+gNke773vc0PMPfDsw7gB9/99muHuuy07xNk+lWa8+7Yv7X3MGb3bHx7mvADQEscl+ZdJ/k2S36YZqt/37hsYwDG9900YcF/f4/82ySNJ/imDv2Fo2+s5enp/XTpg/YcG7N/7dtsQ8zbSPOtQkpwy4P5n9nrs3x3ibH3vOejb/ugQn4PDmReOOWMOfAjwVzImzWvO307z2vT/TfNd6Z/t3f+JDH51nOx5Bbv3m+UeSvLf0zwj8GCSK4dZd3Pvr3cmmdL7+75X8luS/G6vW88Qz1F6502SLwy4/9I0Iz0wsocy294f14dDHnXo8wLAUfV6mlH74oD7Lsqe69HJ0O+A3/u+vlflJ/duz8v+XyXPTPKzDP5Rtr5r8v9uwBqdab4jftZ+Zv9C9rwv4IrseRf9v0qyPXtewR/KbH0G7h/q4z2ceQHgqPkP2ROu/5fmPxLTF8Tv9B5zMIH/5wHP8dM036net3/CEMefk+ZZg11pvpt9VprvRu9O8j+S/K/eOd5J85r5UBppXnPve943em99b/7rC/yhzpYcOPCHMy8AHDWNNK9P/ybNWPWkeQ3+69lzCe1gAn9RkhfTfNPdP6T5M+pdvfuv3c9zPNC7vap3js+k+ca7nUn+lOar/E8dYP4xaV7L/6c03/HeleaP301O81LA3x3mbAcKfA5zXgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAICD12j1AOzXf0vyYauHABiBNiX5bquHGOkEfuR6NcmCVg8BMAI9nGR2q4cY6ca1egD2a2OSf2j1EAAj0JutHmA0EPhRopRSWj0DHCsajYazm4x6Y1o9AABw5Ak8AFRI4AGgQgI/ir388suZP39+pkyZkqlTp+aqq67Ka6+91uqxABgBBH4Uu/HGG3PZZZflrbfeypYtW3Leeedl0aJFrR6r34svvpiJEye2egyAY5LAj2K/+c1vsmTJkowbNy5tbW259dZbM2PGjFaPlSTZtm1bFi1alPfff7/VowAckwR+FPvMZz6TBQsW5Gc/+1l27NiRtra2PProo0mSrVu35vLLL8+UKVNyxhln5Cc/+UmSpKurK1/4whcyderUTJs2LQsXLkxXV1f/czYajTzwwAOZOnVqNm/enCuuuCInnnhipk+fnoULF2bLli2Djh1KKSWLFi3Kt7/97b/iRw/AcAR+FHviiSfy2c9+Nt///vczffr0XHXVVXn55ZeTJLfccks+/elP5+23386TTz6ZZcuWZfv27Vm2bFkmT56cP/7xj9m8eXMmT56cW265ZdDz/uEPf8hLL72URYsWZdGiRenq6srmzZtz7rnn5vrrr+8/7rnnnhtyrh/84Ac566yzct111/31PngAGKXWDNwoQ/jzn//c//vXX3+93HzzzaWzs7OUUsrEiRNLV1dX//633nqr7Nq1q0ycOLG88cYbgx43adKk/u0k5d133y2llDJ+/PiSZNBt/PjxQ43S79lnny0XX3xx6enp6X8+GG2O/l93DtGaAx+CV/Cj2KRJk7J169YkSXt7e+666668+OKLSZIPP/wwY8bs+c/7+uuvZ+fOnSmlDDq1PmbMmOzevXvQ806ZMiVJMnny5Lz99tsppaSUkh07dmTDhg3DzrRq1ao8//zzOe644/rXaTQaWbt27V/+AQNw0AR+FJs/f36+9a1vZdu2bXn//fdz33335eKLL06SXHLJJfnhD3+Y3bt356WXXspFF12UnTt35nOf+1xuvfXW7Ny5M93d3fnOd76Tyy+/fMjnX7BgQW6//fb09PTkT3/6U2644Ybccccd/ftXr169z2O+973v9X9D0PdCqJSSuXPnHvlPAACMQgc8Rd/V1VUWLlxYjj/++DJhwoQyf/78snHjxlJKKZs2bSqXXnppmTRpUjn99NPLww8/XEop5c033yzXXnttmTJlSjnhhBPK5z//+UGn8gcutW3btvLFL36xTJ06tUybNq186UtfKtu2bRvy2GFOdR7iyVFovZb8jedQOEV/EPwPFUauNUnm9W34ogNHj//ZzIg36OsjQ3OKHgAqJPAAUCGBB4AKCTwAVGhcqwcAGGm8qXVk8abHw+MVPABUSOABoEJO0Y9Ss2fPHvRP0R5tM2fOzKpVq/a7f6TPB1A7gR+luru7s2LFipatP2fOnGH3j/T5AGon8KNUR0dHOjs7W7Z+W1vbsPtH+nwAtXMNHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBC41o9AIdn7dq1Wbp0acvW37Fjx7D7R/p8ALUT+FFq5cqV6erqatn6ixcvHnb/SJ8PoHaNVg/Afq1JMq9vo5RSWjgLQMs0Go29WzXo6yNDcw0eACok8ABQIYEHgAoJPABUSOABoEJ+TG6UWrJkScaMad33Z6eddlpuu+22/e6fPXt2S+ebOXNmVq1a1bL1AVrNj8mNXMP+mFx7e3vuvPPOoz5Un3vvvTfr16/f7/4zzjgjK1asOIoTDTZnzpx88MEHLVsfOHL8mNzh8Qp+lGpvb88nP/nJlq0/ffr0Yfd3dHSks7PzKE2zr7a2tpatDTASuAYPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUKFxrR6Aw7Nu3brcc889LVt//fr1w+5fu3Ztli5depSm2deOHTtatjbASCDwo9T69etbGrH29vZh969cuTJdXV1HaZp9LV68uGVrA4wEjVYPwH6tSTKvb6OUUlo4C0DLNBqNvVs16OsjQ3MNHgAqJPAAUCGBB4AKCTwAVEjgAaBCfkxulJo9e3bGjGnd92czZ87MqlWr9rt/pM8HUDuBH6W6u7uzYsWKlq0/Z86cYfeP9PkAaifwo1RHR0c6Oztbtn5bW9uw+0f6fAC1cw0eACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoELjWj0Ah2ft2rVZunRpy9bfsWPHsPtH+nwAtRP4UWrlypXp6upq2fqLFy8edv9Inw+gdo1WD8B+rUkyr2+jlFJaOAtAyzQajb1bNejrI0NzDR4AKiTwAFAhgQeACgk8AFRI4AGgQn5MbpSaPXt2xoxp3fdnM2fOzKpVq1q2PgDDE/hRqru7OytWrGjZ+nPmzGnZ2gAcmMCPUh0dHens7GzZ+m1tbS1bG4ADcw0eACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwcAz59a9/nXnz5qWtrS0TJ07M3Llz8/jjjw86ZteuXWlvb8+nPvWplFIG7Ws0Grn++uv3ed4bbrghjUZj0HHjxo3b5zaUrVu35pZbbsnpp5+e4447Lh0dHfn617+erVu3JknuuuuuXHnllfs87tJLL82Pf/zjg1pvqP3nnXde1q5de5CfORh9BB6OEevWrcs111yTZcuW5Z133snmzZvzta99LV/5ylfy8MMP9x/3/PPP56STTkpXV1d++9vf7vM8q1evzgcffNC/3dPTk9WrV+9z3K5du/a57W3r1q254IILctxxx2XNmjXp7u7O888/n7Fjx+bCCy/Mm2++mYULF+YXv/hFtm/f3v+4d955J2vWrMmCBQsOer2B93d3d+emm27KTTfddEifQxhNBB6OEd/97ndz++235/rrr8+ECRNywgkn5MYbb8yDDz6YDRs29B+3fPnyfPnLX851112Xxx57bJ/nmTt3bn75y1/2bz/77LO56KKLDmumO+64I1deeWXuvffefOxjH8vYsWNz5pln5r777svll1+eO++8M7Nmzcp5552Xn//85/2Pe+qpp3LhhRfm5JNPPqx1x48fnyVLluTVV189rMcD/CXWDNwoe5k7d+7edx1VkydPbun6HLqTTjqpvPHGG8Me09PTU2bMmFE2bdpUnnvuuTJ79uzy4Ycf9u9PUh599NGyaNGi/vu++tWvlkceeaQM/GM6xB/ZIZ166qllw4YNQ+773e9+Vzo6Okoppdx9993l6quv7t937bXXlvvvv/+g19t7/3vvvVfuueeecskllxzUnLTWgb4+wmgj8BxRY8eOLR988EH/dpJBt1JKefrpp8ull15aSill165dZebMmWXdunWDHrN9+/ZyyimnlJ6enrJ79+7S0dFR3n333X0CP9Rtb+PGjSs9PT1Dzrtz584yfvz4UkopGzduLBMmTCjbt28v3d3dZdKkSWXz5s0Hvd5Q+xqNRnnuuecO4zPJ0Xagr48MzSl6OEacfPLJeeWVV/q3SykppeT999/vv2/58uV55pln+t+UtmXLln1O00+aNCmdnZ1ZvXp1XnjhhZx77rmZMmXKPuv1Pf/A297a29uzcePGIefdtGlT/yn4008/PZ2dnXnqqaeyatWqnH/++TnllFMOab2B93d3d+fuu+/OzTfffIDPGoxeAg/HiMsuuyyPPPLIPvf3Xdvu6enJ448/nldffbU/hE8//XSWL1++TywXLFiQn/70p3n88cdzzTXXHPZMV155Ze6///4h9/3oRz/K1Vdf3b+9cOHCPPbYY3niiSdy3XXXHfaaSdLW1pZly5bl97///V/0PACHwyl6jqgNGzaUE088sTzwwAPlvffeK9u3by8PPvhgOfvss8vYsWPLk08+Wc4///xBj+nu7i4f/ehH+0/T9/1RfOutt8qpp55aPvGJT/Rf108O/Rr866+/Xs4888zyzW9+s7zyyitl9+7d5ZVXXinf+MY3yqxZs8rWrVv7j92yZUs5/vjjy4wZM8prr7026HkOtN7+9h/snLTWgb4+MjSv4OEY8fGPfzxPP/10HnrooUybNi2zZs3Ks88+m2eeeSbnnHNOli9fniuuuGLQY9ra2nLJJZfsc5p+2rRpOeecczJz5szMmDFjyPUO5ufg29vb88ILL2Tnzp2ZO3duPvKRj2TOnDnp7u7Or371q0yfPr3/2JNPPjkXXHBBzjrrrJx22mmHtd7eTjzxxAMeA6NV48CH0CJrkszr29j7u9h58+ZlzZrWfRM7ZcqUbNu2rWXrA8eOxsB/Ralp0NdHhuYVPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgf+x5oZkdauXZulS5e2bP0dO3a0bG0ADkzgR6mVK1emq6urZesvXry4ZWsDcGACP0rNnz+/1SMAMIK5Bg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUa1+oBgL/ckiVLMmbM/r9fP+2003LbbbcdxYmAVmu0egD2a02SeX0bpZTSwlkY4c5oNLJimP03nH121q9ff9TmgSOp0Wjs3apBXx8ZmlfwUIGOJJ3D7J8+ffrRGgUYIVyDB4AKCTwAVEjgAaBCAg8AFRJ4AKiQd9FDBdYmWTrMfj8iB8cegYcKrF+/Pjt27Njv/tva24/iNMBIIPBQgbPPPrvVIwAjjGvwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKjQuFYPwMFpNBqNVs8AwOjhFTwAVEjgAaBCTtGPXHOT/NdWDwEwAs1u9QCjgeu6I9fZSY5v9RAAI9AbSTa1eggAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIBj0/8HeWg0f67vYjEAAAAASUVORK5CYII=" alt="plot of chunk unnamed-chunk-7"/> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAYAAACmKP9/AAAACXBIWXMAAAsSAAALEgHS3X78AAAWaElEQVR4nO3df4zV9b3n8dfhhw4VRFBA/DEqtZZSg3+wu21R09uIiV0jWrXSqGnpXi1Vu7Xd7G5v6m390bRm06hrTbq7utm7Wn8kmKptNZJU1ERa72288W7dNSqtRITyo6MVBAen4Gf/ODPDDAwzA6VzZj48HskJfM/3e77f95zAec75fs9AAgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADGZ8kv+Q5KUk7yfZnuSFJH+bpHEQj9OWpHTfDqZGkiuT/GOas29P8uskX8zBnR8AxpT/kmZ0/5xmJF/M7hB/9yAe568V+P/aZ79/SLK5z/IdB/lYADAmNJK8k2YMP9Hn/gu773s7B+9d8F8j8Od17++dJJ9Jc9ZGkovSPBtRksw/iMcDgDFhQnZH97Q+949LsrT7Ni794/y5JL9P8x3/PyWZ1/2YWUn+IcmGJLuSbEnyZJKPd6/fM/Djkvy3JKuS/Lvu+05L8kSap9nfSXJvkqMHmf/x7v19bYB1y5L89yR/s5+z/V338V/rfux/SrK2++t9Mcm/7nOM/Z0XAEbMr9IM23tJfpHkhiTnphm9Hn0DuC3NKL/Vvfyr7m1Wdi+vS/JQmvEvSV4ZYB9J8oPu3z+W5jcaRyVZn6QryX1pBrh072f8Pmb/U/c2JwzxNe7PbBvTjHvPZYv3u7/enlP/L3U/5kDmBYARc3yS5Uk6sztyJc2g/U33Nn0D2HPf/OyOYJL87yT3J5nbvdze5zGH77GPL3X/+lySSd3b//vu+37QZ7ZHu+/77ABz9z37MKHP/f+QZrh7btfs52wfSzPePcuLuh9zUvfyBwc4LwC0xGFJ/k2SbyT5f2mG6rXudX0DOK77vkl97ut5/N8meTDJv6T/Nwxte+yjq/vXZX2Of1+f9XvevjPAvI00zzqUJMf1uf+pPR779/s5W89nDnqWPzTAc3Ag88IhZ9zQmwB/JePSvOb8d2lem/5Nmp9K/0z3+o+k/7vjZPc72D0/LHdfkv+Z5hmBe5MsHuS467t/vSXJ1O7f97yT35Dk1T1uXQPso3TPmySX9bl/UZqR7hvZ/Zltz6/rgwG32v95AWBEbUwzal/oc9+Z2X09Ohn4E/B73tfzrvzY7uWzs+93ybOT/Cz9f5St55r8f+xzjAVpfiJ+zj5mvyy7PxdwQXZ/iv7fJtma3e/g92e2Hn3XD/T1Hsi8ADBi/nN2h+v/pPmPxPQE8dvd2wwn8L/vs4+fpvlJ9Z71kwbYfl6aZw12pvlp9jlpfhq9M8n/SvKT7jn+lOY184E00rzm3rPfTd23ng//9QR+f2dLhg78gcwLACOmkeb16RfSjFVXmtfgv5bdl9CGE/gzk/zfND90989p/ox6R/f6S/axj7u7l1d2z/GpND94tyPJu2m+y//YEPOPS/Na/r+k+Yn3jjR//O7INC8F/P0BzjZU4HOA8wIAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAMX6PVA7BP9yT5oNVDAIxC65J8r9VDjHYCP3q9keTiVg8BMAo9kGRuq4cY7Sa0egD2aW2Sf271EACj0B9bPcBYIPBjRCmltHoGOFQ0Gg1nNxnzxrV6AADg4BN4AKiQwANAhQR+DHv99ddz3nnnZerUqZk2bVouvPDCvPnmm60eC4BRQODHsCuvvDLnnntu3nrrrWzYsCHz58/P0qVLWz1WGo3GgDcARo5P0Y9hL7zwQlasWJEJEyZkwoQJueGGG/LlL3+51WPl3Xff7bd80003ZdKkSS2aBuDQ5B38GPapT30qF198cX72s59l27ZtaWtry0MPPZQk2bx5c84///xMnTo1J510Uu6///4kSUdHRy677LJMmzYt06dPz5IlS9LR0dG7z0ajkbvvvjvTpk3L+vXrc8EFF+Too4/OjBkzsmTJkmzYsKHftgOZPHly7+03v/lNXnjhhdx0001/vScCgL0I/Bj22GOP5TOf+Ux+8IMfZMaMGbnwwgvz+uuvJ0muu+66fPzjH8/bb7+dxx9/PNdcc022bt2aa665JkceeWT+8Ic/ZP369TnyyCNz3XXX9dvv7373u7z88stZunRpli5dmo6Ojqxfvz6nn356Lr/88t7tnnnmmUHn27ZtW7761a/mnnvuyfjx4w/+EwAAY9BzfRfKAP785z/3/n7jxo3l2muvLQsWLCillHLEEUeUjo6O3vVvvfVW2blzZzniiCPKpk2b+j1uypQpvctJyjvvvFNKKWXixIklSb/bxIkTBxplQN/97nfL9ddfP+ztYbRowd939s9zQ2+Ca/Bj2JQpU/LGG29k5syZmTVrVm699dbMnDkzSfLBBx9k3LjdJ2g2btyYww8/PKWUfqfWx40bl127dvXb79SpU5MkRx55ZFavXp1p06YlSbZv397vdP5gtm3blrvuuisvvvjiX/Q1AnBgnKIfw84777x861vfypYtW/Lee+/lzjvvzKc//ekkyTnnnJMf/vCH2bVrV15++eWceeaZ2bFjRz772c/mhhtuyI4dO9LZ2Zlvf/vbOf/88wfc/8UXX5wbb7wxXV1deffdd3PFFVfk5ptv7l3/7LPP7nO2hx9+OPPnz89JJ510UL9mABjrhjxF39HRUZYsWVImT55cJk2aVM4777yydu3aUkop69atK4sWLSpTpkwpJ554YnnggQdKKaX88Y9/LJdcckmZOnVqOeqoo8rnP//5fqfy+x5qy5Yt5Qtf+EKZNm1amT59evniF79YtmzZMuC2e1q8eHH5/ve//xedJoVWac1fefaDU/TD4IeTR6/nkpzds+BFB0aO/2xm1Ov3+sjAnKIHgAoJPABUSOABoEICDwAV8nPwAHs41D7U6kOFdfIOHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKjQhFYPACOh0Wg0Wj3DYEoppdUzcOga7X/+Rvvf39HKO3gAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVMh/FwvDc0+SDwZZvy7J90ZoFoAhCTwMQ3ty1SODrP9XyasReGAUEXgYhvYkCwbf5I8jMgjAMLkGDwAVEngAqJDAA0CFBB4AKiTwAFAhn6KHYViVZNngm8wdkUEAhkngYXjm3p1MHmT9phGbBGAYBB6G59VWDwCwP1yDB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CF/HexHBJKKaXVMwCMJO/gAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgYdDyK9//eucffbZaWtryxFHHJGzzjorjz76aL9tdu7cmVmzZuVjH/tY9vxfdhuNRi6//PK99nvFFVek0Wj0227ChAl73QayefPmXHfddTnxxBNz2GGHpb29PV/72teyefPmJMmtt96axYsX7/W4RYsW5cc//vGwjjfQ+vnz52fVqlXDfOYADp7n+i4U+Av99re/LTNnziwPPPBAee+998qf/vSn8pOf/KQcddRR5f777+/d7qmnnirz5s0rxxxzTHnppZf67SNJmT17dtmxY0fvfe+//3457rjjSt8/psP9I7tp06Zy8sknl2984xtlzZo1ZefOneX1118vX//618spp5xSNm/eXH7/+9+Xww8/vGzZsqX3cW+//XY57LDDyoYNG4Z1vD3Xd3V1ldtvv72ceuqpw5qT1hrq9ZGBeQcPh4jvfe97ufHGG3P55Zdn0qRJOeqoo3LllVfm3nvvzerVq3u3W758eb70pS/l0ksvzcMPP7zXfs4666z88pe/7F1++umnc+aZZx7QTDfffHMWL16cO+64IyeffHLGjx+fU045JXfeeWfOP//83HLLLZkzZ07mz5+fX/ziF72Pe+KJJ/LJT34yxx577AEdd+LEibn66qvzxhtvHNDjAf4S3sFzUB1zzDFl06ZNg27T1dVVZs6cWdatW1eeeeaZMnfu3PLBBx/0rk9SHnroobJ06dLe+77yla+UBx988IDewR9//PFl9erVA6579dVXS3t7eymllNtuu61cdNFFvesuueSSctdddw37eHuu3759e7n99tvLOeecM6w5aa2hXh9hrBF4Dqrx48eX999/v3c5Sb9bKaWsWLGiLFq0qJRSys6dO8vs2bP7naZPUrZu3VqOO+640tXVVXbt2lXa29vLO++8s1fgB7rtacKECaWrq2vAeXfs2FEmTpxYSill7dq1ZdKkSWXr1q2ls7OzTJkypaxfv37YxxtoXaPRKM8888wBPJOMtKFeHxmYU/RwiDj22GOzZs2a3uVSSkopee+993rvW758eZ566qneD6Vt2LBhr9P0U6ZMyYIFC/Lss8/m+eefz+mnn56pU6fudbye/fe97WnWrFlZu3btgPOuW7eu9xT8iSeemAULFuSJJ57IypUrc8YZZ+S4447br+P1vb+zszO33XZbrr322iGeNRi7BB4OEeeee24efPDBve7vubbd1dWVRx99NG+88UZvCFesWJHly5fvFcuLL744P/3pT/Poo4/mc5/73AHPtHjx4tx1110DrvvRj36Uiy66qHd5yZIlefjhh/PYY4/l0ksvPeBjJklbW1uuueaavPbaa3/RfgAOhFP0HFSrV68uRx99dLn77rvL9u3by9atW8u9995bPvrRj5bx48eXxx9/vJxxxhn9HtPZ2Vk+9KEP9Z6m7/mj+NZbb5Xjjz++fOQjH+m9rp/s/zX4jRs3llNOOaV885vfLGvWrCm7du0qa9asKddff32ZM2dO2bx5c++2GzZsKJMnTy4zZ84sb775Zr/9DHW8fa33V2tsGOr1kYF5Bw+HiFNPPTUrVqzIfffdl+nTp2fOnDl5+umn89RTT2XevHlZvnx5Lrjggn6PaWtryznnnLPXafrp06dn3rx5mT17dmbOnDng8Ybzc/CzZs3K888/nx07duSss87K4YcfnoULF6azszO/+tWvMmPGjN5tjz322HziE5/Ihz/84ZxwwgkHdLw9HX300UNuA2NVY+hNaJHnkpzds7CP72IBqtfo+68oNfV7fWRg3sEDQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACo0odUDAHubO3duxo1r3fffs2fPzsqVK/e53nyDG2o+GAkCD6NQZ2dnHnnkkZYdf+HChYOuN9/ghpoPRoLAwyjU3t6eBQsWtOz4bW1tg6433+CGmg9GgmvwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABWa0OoBgL2tWrUqy5Yta9nxt23bNuh68w1uqPlgJAg8jEJPPvlkOjo6Wnb8q666atD15hvcUPPBSGi0egD26bkkZ/cslFJKC2cBaJlGo7Fnq/q9PjIw1+ABoEICDwAVEngAqJDAA0CFBB4AKuTH5DgkzZ07N+PGte7729mzZ2flypX7XD/a57v66qtbOt8JJ5yQ73znO/tcP9qfPxgJAs8hqbOzM4888kjLjr9w4cJB14/2+X7+85/nlltuGaFp9nbHHXcMGvjR/vzBSBB4Dknt7e1ZsGBBy47f1tY26PrRPt+sWbNy2mmnjdA0e5sxY8ag60f78wcjwTV4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQIYEHgApNaPUA0AqrVq3KsmXLWnb8bdu2Dbp+tM/30ksv5fbbbx+hafb2yiuvDLp+tD9/MBIEnkPSk08+mY6OjpYd/6qrrhp0/Wif75VXXmlpxGbNmjXo+tH+/MFIaLR6APbpuSRn9yyUUkoLZwFomUajsWer+r0+MjDX4AGgQgIPABUSeACokMADQIUEHgAq5MfkOCTNnTs348a17vvb2bNnZ+XKlftcb77BDTUfIPAcojo7O/PII4+07PgLFy4cdL35BjfUfIDAc4hqb2/PggULWnb8tra2Qdebb3BDzQe4Bg8AVRJ4AKiQwANAhQQeACok8ABQIYEHgAoJPABUSOABoEICDwAVEngAqJDAA0CFBB4AKiTwAFAhgQeACgk8AFRI4AGgQgIPABUSeACokMADQIUEHgAqJPAAUCGBB4AKCTwAVEjgAaBCAg8AFRJ4AKiQwANAhQQeACok8ABQoQmtHgBaYdWqVVm2bFnLjr9t27ZB15tvcEPNBwg8h6gnn3wyHR0dLTv+VVddNeh68w1uqPmApNHqAdin55Kc3bNQSiktnAWgZRqNxp6t6vf6yMBcgweACgk8AFTINfgxYoBTVACwT97BA0CFBB4AKuQU/eh1VpL/0eohAEahua0eYCxwXXf0+miSya0eAmAU2pRkXauHAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAODT9f9ZzwyY4M/XRAAAAAElFTkSuQmCC" alt="plot of chunk unnamed-chunk-7"/> </p>

<pre><code>##  [1] 5 5 3 5 1 3 4 4 4 7
</code></pre>

<p>ヘビが動き回るのが見えて、けっこう面白いのですが、多くの回数繰り返すには動きが見える必要はありません。そこで、plot=FALSE パラメータを追加して、100回繰り返してヒストグラムを書いてみます。</p>

<pre><code class="r">r &lt;- unlist(sapply(1:100, function(x) {
    rsnake(auto = random, plot = FALSE)
}))
hist(r)
</code></pre>

<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAMAAACR9g9NAAAAdVBMVEX9/v0AAAAAADkAAGUAOTkAOY8AZo8AZrU5AAA5ADk5AGU5OWU5OY85j9plAABlADllAGVlOQBlOY9lZjlltf2POQCPOTmPZgCPtY+P29qP2/21ZgC1/rW1/v3ajzna/rXa/tra/v39tWX924/9/rX9/tr9/v3yO58FAAAAJ3RSTlP//////////////////////////////////////////////////wCDVpfZAAAACXBIWXMAAAsSAAALEgHS3X78AAAKUUlEQVR4nO3dbUPiVh5A8UbryrbqTLfS3ZWFrgTy/T9ik5A4OAMmwZsmf845L6rFcC/OjzwhDz8VhuynqW+ATZPw0ISHJjw04aEJD014aMJDEx6a8NCEhyY8NOGhCQ9NeGjCQxMemvDQhIcmPDThoQkPTXhowkMTHtrVwu+ebtdFsV/erg/fta2ePzFovshuXj5902YRAP7dxavsM/Cfu/asAsDX35Xrana73i+zLLsvik355bGofpzd/u/p7nW/vPm9XJe35cXlGr3N/lFevi4XapkPy9fXvnutr1cvP+Gv9+kg8LunrEL7/wF+Vf1P+bWW/HlRwVc//fOw0Os2a2s2F83yx/Dtd2G7YvjWroLPFw1TtbHOF+XKunu6eckXpewmq+Hv26vdrrdZdYXsuVqkuuxt+bdN/bflwwaBr/+vWn8rum2ttsme6292h039S3ulCv6+vgO0l74tfwQfeztfXDX8D/v4rFlnz8HXy7zBf7s7CB+pE6dzq4bu1Ka+gtyUB3DVRd/Dn9rUCz/X3sO3K3N1nHby4K6C3Gan1/i35YUP0IlNfbM5Lw/z2tO5cp/ens5VkCXwY7MHeAf/trzwV9Q2+HnZhaHhmyP/x6lvxxSh4Q979ehn5JfFhgcnPDThoQkPTXhowkMTHprw0ISHJjw04aEJD014aMJDEx6a8NCEhyY8NOGhCQ9NeGjCQxMemvDQhIcmPDThoQkPTXhowkMTHprw0ISHlhA+G1C6We2yUsKPsqiNk/DQhIcmPDThoQkPTXhowkMTHprw0ISHJjw04aEJD014aMJDEx6a8NCEhyY8NOGhCQ9NeGjCQxMemvDQhIcmPDThoQkPTXhowkMTHprw0ISHJjw04aEJD014aMJDEx6a8NCEhyY8NOGhCQ9NeGjCQxMemvDQhIcmPDThoQkPTXhowkMTHprw0ISHJjw04aEJD014aMJDEx6a8NC64fNFVnW77lpQ+Eh1wu+Xz/XX7d1rx5LCR6oTfvd1/e7r+YSPlGs8tO59/O7JffwV5lE9tEvgs7bvLh4wwgWzWtJ6nc49rk5u6oUPXK+Du1Wpnj/8cHAnfOB6nc5tH0+ezgkfuN6nc67x11Wf07lKfuM+/rr6zOmc8IETHprw0ISHJjw04aEJD014aMJDEx6a8NCEhyY8NOGhCQ9NeGjCQxMemvDQhIcmPDThoQkPTXhowkMTHprw0ISHJjw04aEJD014aMJDEx6a8NCEhyY8NOGhCQ9NeGjCQxMemvDQhIcmPDThoQkPTXhowkMTHprw0ISHJjw04aEJD014aMJDEx6a8NCEhyY8NOGhCQ9NeGjCQzvA757uL7iu8IFr1/htlt28DLyu8IE72tTvl1n2POS6wgeuhc8X1Rq/+7oecF3hA9fu4+9eh19X+MB5VA+tgd+We/fN0KM74QPXbOq/VOb5L0P28MKH7gC/X1aH89uBO3rhA9ds6ndPWZbdDlvhhY/cRAd3/fvE7bMPmvtRvfAj9XZUnw3f1gsfuPYBnPOP1eaLc3cK4QPXwJ9/qPZwwH/ykF/4wDWb+tXjuQXa+8SJ+4bwgWs39Wf38a7x11n3Uf35O4XwgfN0DloDv19md39+6flHmjMPrggfqfax+sf84fX0Y/WHp2h4cHdlvZ3OlfAnT+qqg7vyfiH8lXW8xm9OrfEH8NW98NfVt318lp3c0jenc5uff/xjvfCB63M6Vz+4s/nxfE74wHk6B63zkbsPEj5wx2v85uwD9qcTPnDH8MNeTiF86I7ht27qOb3bxw965ZzwofOoHprw0N5t6gee0AkfuGaN39y3/xmQ8IE7frKlp3Og3v46V7jGozr+69zQN0ASPnAe1UMTHtolT7ZsEz5w3U+2PJ/wget8suUHCR+4zidbfpDwget8suUHCR84j+qhdb4+/oOED1yzj/9j6DtXVwkfOJ9lC819PDThoVXwlx3aCR+6Fj5/GHgSXwgfOuGhCQ+thr/oObbCh86jemjCQxMemvDQhIcmPDThoQkPTXhowkMTHprw0ISHJjw04aEJD014aMJDuyL4AR9K793pquBHWfRaEx6a8NCEhyY8NOGhCQ9NeGjCQxMemvDQhIcmPDThoQkPTXhowkMTHprw0ISHJjw04aEJD014aMJDEx6a8NCEhyY8tG74fHHuQwyED1wn/H75XH898eHywgeuE779bLITn1EmfOBc46F17+PPf+6s8IHzqB7aJfBn3klG+Ei5xkMTHlr36dz5j6EUPnDda/x++XjmJ8IHrsemfvfl5fQPhA+c+3howkMTHtrs4Ud5Z2Lh5w8/xqDCC49NeGjCQxMemvDQhIcmPDThoQkPTXhowkMTHprw0ISHJjw04aFB4f3AcSj8GIPGSvhkg8ZK+GSDxkr4ZIPGSvhkg8ZK+GSDxkr4ZIPGSvhkg8ZK+GSDxkr4ZIPGSvhkg8ZK+GSDxkr4ZIPGSvhkg8ZK+GSDxkr4ZIPGSvhkg8ZK+GSDxkr4ZIPGSvhkg8ZK+GSDxkr4ZIPGSvhkg8ZK+GSDxkr4ZIPGSvhkg8ZK+I4lr/XVVsJPMf8MEn6K+WeQ8FPMP4OEn2L+GST8FPPPIOGnmH8GCT/F/DNI+Cnmn0HCTzL/KPWfvxCeOX8hPHP+Qnjm/IXwzPkL4ZnzF8Iz5y+EZ85fCM+cvxCeOX8hPHP+Qnjm/IXwzPkL4ZnzF8Iz5y+EZ85fCM+cvxCeOX8hfMpBx3iyjPATLXqd8xfCM+cvhGfOXwjPnL8Qnjl/ITxz/qIPfL6oT0Bu110zBfrF6fMXPeD3y+f66/butWOmQL84ff6iB/zu6/rd1+LokYrvJrZpSwv/wRpvgevex++e6rvTiX28Be4zR/UWOOGhCQ9NeGjCQxMemvDQhIcmPDThoQkPbRr4if+MdaUNIpgInj3oHG6q8BMMOoebKvwEg87hpgo/waBzuKnCTzDoHG6q8BMMOoebKvwEg87hpvoADjThoQkPTXhowkMTHprw0ISHJjw04aFNAF+9xcZz+mHb13MnHTK7eUk9aPX7J3/pcf7Lun5dc//Xsv/98LsvL0X+z/T/npv096bV8whvC1D9/pvEo26ru1J1z9/c973K3w+/rW7bKjlS/utvqcc8ehOQhOUPr6lHXt38p1zjq0HrNb9X0+zjq3t92vZ//Df5pj5/+PcIm/ox1vgavL5H9f6HnQR+v3xMPeTmMf0+Pl881/+aiRu0K+5ZBV/tleYNv3tK7l76jAA/aBXqPWp5fLNNfXQXYo2vVqTUbeonlie+P+3+NQb8sBWzZ3mAffwY7lUjnM6txtjUj7bGVzvQOR/VH1bO9PgjwJd74xHe7GubpT9kjHAeb7NIeGjCQxMemvDQhIcmPDThoQkPTXhowkMTHprw0ISHJjw04aEJD014aMJDEx6a8NCEhyZ8U/7rb6jPTRa+aazXecw14Zv6v/joOhK+SXhowkMTHprwhkh4aMJDEx6a8NCEhyY8NOGhCQ9NeGjCQxMemvDQhIcmPDThoQkPTXhofwESvVRIvv5YwgAAAABJRU5ErkJggg==" alt="plot of chunk unnamed-chunk-8"/> </p>

<p>さすがにこのアルゴリズムではあまり良い点は出ないようですね。では、せっかくエサの位置の情報もあるので、エサに向かって一直線という本能に忠実なヘビを作ってみましょう。</p>

<pre><code class="r">towardsFood2 &lt;- function(snake, food, x, y) {
    bodycheck &lt;- function(newhead) {
        if (newhead[1] &lt; 1 || newhead[1] &gt; x || newhead[2] &lt; 1 || newhead[2] &gt;
            y)
            return(FALSE)
        return(sum(sapply(1:length(snake), function(x) {
            identical(snake[[x]], newhead)
        })) == 0)
    }
    rad &lt;- atan2(food[2] - snake[[1]][2], food[1] - snake[[1]][1])/pi + 0.25
    if (rad &gt; 1)
        rad &lt;- rad - 2
    v &lt;- as.numeric(cut(rad, breaks = c(-1, -0.5, 0, 0.5, 1)))
    nh &lt;- snake[[1]] + list(c(-1, 0), c(0, -1), c(1, 0), c(0, 1))[[v]]
    while (!bodycheck(nh)) {
        v &lt;- v + 1
        if (v &gt; 4)
            v &lt;- 1
        nh &lt;- snake[[1]] + list(c(-1, 0), c(0, -1), c(1, 0), c(0, 1))[[v]]
    }
    return(c(&quot;h&quot;, &quot;j&quot;, &quot;l&quot;, &quot;k&quot;)[v])
}
</code></pre>

<p>単にエサに向かっていくだけだと、その方角に自分の体があったときにそこでスタックしてしまいます。そこで、rsnake()本体でも使っているbodycheck()関数をいれておき、エサの方向が通れないようであったら方向をずらすようにしてあります。</p>

<pre><code class="r">r &lt;- unlist(sapply(1:100, function(x) {
    rsnake(auto = towardsFood2, plot = FALSE)
}))
hist(r)
</code></pre>

<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAMAAACR9g9NAAAAdVBMVEX9/v0AAAAAADkAAGUAOTkAOY8AZo8AZrU5AAA5ADk5AGU5OWU5OY85ZrU5j9plAABlADllAGVlOQBlZjlltf2POQCPOTmPOWWPtY+P29qP2/21ZgC1/rW1/v3ajzna/rXa/tra/v39tWX924/9/rX9/tr9/v3/OhVWAAAAJ3RSTlP//////////////////////////////////////////////////wCDVpfZAAAACXBIWXMAAAsSAAALEgHS3X78AAAN60lEQVR4nO2df0PjxhFAI0gKTTm4tjjXg+IEG/v7f8RqZdln1z52R4LZ2Zn3/shdwO+k1UM/bGztL1sIyS+1VwDqQPigED4ohA8K4YNC+KAQPiiEDwrhg0L4oBA+KIQPCuGDQvigED4ohA8K4YNC+KAQPiiEDwrhg0L4oBA+KIQPCuGD4jb828P1y3a7WVy/7P625/lxxj+6vu2unmavmgkChD/58nM3J/w82xQBwg9/6/fV7vpls+i67ma7XfZ/3G/Tt7vr/z789rpZXP2735dX/Zf7PXrV/a3/+kv/oH3m3eMH+7fXwRseX3F4swkS/u2hS9H+2oV/Tv/T/zmU/PU2hU/f/XP3oNdVt2c8XIyPPw6//1uzOA6/b5fCr2/HTOlgvb7td9a3h6un9W1fdtkN4W/22vXLqktC95gekr52ePzhUP/j8c0SJPzwf2n/TelWQ7Vl9zj85W13qH/aSyn8zfADsP/q4fFH4ds+zm9dhz87x3fjPvuz8MNjDuF//DgQviUuPJ17HtNdOtSnkMv+Ai596f/DXzrUE94qp+H3O3O6Trt4cZdCrrrLe/zh8YRvgAuH+vFw3l/m7Z/O9ef0/dO5FLIPfD+eAU7CHx5PeEesGn9eNpHQ4ccr//va61GD0OF3Z/XWn5FPI3b4wBA+KIQPCuGDQvigED4ohA8K4YNC+KAQPiiEDwrhg0L4oBA+KIQPCuGDQvigZMO//fNleItSyDemOaYkfGq/XX/RWB3QoiT8+u51t+eDH/LhH66+f0t7/B3Hek8UXNwNHw1dXbPDu4Kr+qBMCX+4b8CHrw2okQ+/vu3unw83hziG8A2TDb9ZPG6f06cOzy/uCN8wRS/grO4vPp0jfMMU7fEJ9nhf5M/xbw+p/JJzvC/mPJ0jfMMQPiiEDwrhg0L4A91cag9ABOEPzB1OW5uD8AcIX0pbI81C+FLaGmkWwpfS1kizEL6UtkaahfCltDXSLIQvpa2RZiF8KW2NNAvhS2lrpFkIX0pbI81C+FLaGmkWwpfS1kizEL6UtkaahfCltDXSLIQvpa2RZiF8KW2NNAvhS2lrpFkIX0pbI81C+FLaGmkWwpfS1kizEL6UtkaahfCltDXSLIQvpa2RZiF8KW2NNAvhS2lrpFkIX0pbI81C+FLaGmkWwpfS1kizEL6UtkaahfCltDXSLIQvpa2RZiF8KW2NNAvhS2lrpFkIf8r6drjPh//73BH+hP2dLVfnk9K0NdIshD9hfw9b//eyJfwJ7PFavi4l97LlHK/i68JV/QHC52jzjn5ZCH9KmKlJCH9CnKlJCH9CnKlJCH9CnKlJCH9KmKlJCF9KWyPNQvhS2hppFsKX0tZIsxC+lLZGmoXwpRgbae2pRYxtjgyewjfu60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj60J4M74uhDfj61J0g8N0A+vzu5ZbG2ntcMY2R4aS8MNNLddfzr5lbKS1wxnbHBlKwg83tbR/S9Pa4Yxtjgz58A9X37+9tHBL09rhjG2ODAUXd5tFd7Nd2b+lae1wxjZHBq7qzfi6OJqhonY4Y5sjQ9EMFVdPXNx9vq9L0f3qN4t7wn+6r0vhhIPPN4T/bF+X0hkqlr/+TvjP9XUpmaHiPv1xYYoKYyOtHc7Y5sjA0zkzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi6EN+PrQngzvi5F97JNcNvyz/Z1Kb2z5XZ1PimNsZHWDmdsc2QovJctd6/+fF8X9ngzvi4l97LlHK/i68JVvRlfF6YmMePrwtQkZnxdduHfHm5+9gCmJtHyddnv8atu2K/PYWoSLV+Xo0P9ZtF1j2cPYGoSLV+XffjdmfzCbs3UJFq+Lvtz/IVZwrMYG2ntcMY2Rwaex5vxdRnDr/qz+/Ly1d3PMTbS2uGMbY4M46H+a2q+Pr9+exdjI60dztjmyLALv7t0v/B7mHcxNtLa4YxtjgzjoX74TcyF38O8i7GR1g5nbHNk4OLOjK8L4c34uhyu6n/yO/f3MDbS2uGMbY4M+xdwzl+rzWNspLXDGdscGcbwF16qzWNspLXDGdscGcZD/fP9BNfYSGuHM7Y5MuwP9Zzjq/u6cFVvxteF8GZ8Xcbwm0X3259f+SVNTV+X/Wv19+u7V16rr+rrcng614eXPqkzNtLa4YxtjgzHe/ySPb6mr8uPc3zXSd9+ZWyktcMZ2xwZuKo34+tCeDO+LrxyZ8bX5XiPXwpfsDc20trhjG2ODMfheTpX1dflOPyKQ31NX5eTc7zw3RjGRlo7nLHNkYGrejO+LoQ34+tycqgXPqEzNtLa4YxtjgzjHr+82f9HgLGR1g5nbHNkOH6zJU/nqvq6HH47t2WPr+zrcvzbOWF3ayOtHc7Y5sjAVb0ZXxfCm/F14c2WZnxdeLOlGV8X3mxpxtcl+2bL9NOQXti78D1jI60dztjmyJB9s2UffjgSrL+cfcvYSGuHM7Y5MpRMTdIfDLh79ef7umQ/H//2cPX9W9rj786OB8ZGWjucsc2RYTzH//HOM7n+NHBz8d05xkZaO5yxzZGBd9ma8XUxNDVJN5e5y298/WUYesm2erjGfRkp/LRbHxHemC9jH/7CRfuOd96VRXhTvoxs+PFNGpcgvClfRj78eE/zCxDelC9jCD/pPbaEN+bL4KrejS+D8G58GYR348sgvBtfBuHd+DII78aXQXg3vgzCu/FlEN6NL4PwbnwZhHfjyyC8G18G4d34MgjvxpdBeDe+DMK78WUQ3o0vg/BufBmEd+PLILwbXwbh3fgyCO/Gl0F4N74MwrvxZRDejS+D8G58GYR348sgvBtfBuHd+DII78aXQXg3vgzCu/FlEN6NL4PwbnwZhHfjyyC8G18G4d34MvLh17c/uyMW4U35MrLhN4vdrPIXZioivClfRskMFSd/HkF4U74M9ng3voz8Of7nkxgQ3pQvg6t6N74MSzNU4CvCHu/Gl0F4N76M/NM5ZqhoxJeR3+OZoaIRX0bBoZ4ZKtrwZXCOd+PLILwbXwbh3fgyCO/Gl0F4N74MwrvxZRDejS+D8G58GYR348sgvBtfBuHd+DII78aXQXg3vgzCu/FlEN6NL4PwbnwZhHfjyyC8G18G4d34MgjvxpdBeDe+DMK78WUQ3o0vg/BufBmEd+PLILwbXwbh3fgyCO/Gl0F4N74MwrvxZRDejS+D8G58GYR348sgvBtfBuHd+DII78aXQXg3vgzCu/FlEN6NL4PwbnwZhHfjyyC8G18G4d34MoomI0p3sD6fkobwtnwZJeGHCajWX86+RXhTvoyS8Ou715Ppx5iaxKQvo2Cigqvv39Ief3d2rCe8KV9GwcXdZtHdbFdMP2bdl8FVvRtfBuHd+DII78aXQXg3vgzCu/FlEN6NL4PwbnwZhLfjz0W0NMIH9Qkf1Cd8UJ/wQX3CB/UJH9QnfFCf8EF9wgf1PzC87itP2dXBf5ePDD/jn8InPL6KT/igPuGD+oQP6hM+qE/4oD7hg/qED+oTPqhP+KA+4YP6hA/qEz6oT/igPuGD+oQP6hM+qE/4oD7hg/qED+oTPqhP+KA+4YP6hA/qEz6oT/igfj78+nb4YFv+7tVtDTy6nw2/WTwOf67O5yYhfMN+0Zw0x3+mJVz+eOvsT8vCPD42/Dt7PDRM/hyfpqDqLp7joWHmXNVDwxA+KIQPCuGDQvigED4ohA8K4YNC+KAQPiiED8oHhq/8uymoFh6/IZ/wQX3CB/UJH9QnfFCf8EF9wgf1eQEnKIQPCuGDQvigED4ohA8K4YNC+KAQPiiED8pHhX976GZ9jnrZzfpA7vr3l1nrMPiT1yHdNORxxvJHf/LyVztPtPwPCp8+Rb+8mfEPPD/OWfwqDXzGOgz+5HV4+/q0Xf/9afLyR3/y8tNPbb9g2fI/KHy6X8aw10xk88fTjKU/X/2nX/b0ddj5k9dhlbb28+Pk5Y/+rG3QL1i2/A8Kv757HX5wpzLcfWHGTp8GPGcdkj9rHfoFz9oGvThr+f2uLlv+B4VPN0qZEz4d6eb8xKdwc9Zh+MGZsQ6bxf2s5Sd/xvLXt1dPwuUb2eMHZpznP2KPn7EObw/3s7bB4M9Y/oQjjpFz/MDM8HPWYV749W2ypi9/509f/k6sco5Ph6o5V/XpMLX5Nu/p3Jx12J8qJq3D2G3y8kd/8vLHY7xs+Yaex1/Nu6j9kOfxk9ZhOXyO5XHy8vf+5G0wijWex0NrED4ohA8K4YNC+KAQPiiEDwrhg0L4oBA+KIQPCuGDQvigED4ohA8K4YNC+KAQPiiEDwrhg0L4oBA+KIQfWf/jX6HmTSb8yNGnWUJA+JHZHwBrDMKPED4ohA8K4YNCeAgB4YNC+KAQPiiEDwrhg0L4oBA+KIQPCuGDQvigED4ohA8K4YNC+KAQPiiEDwrhg/I/GboWP5vYCuUAAAAASUVORK5CYII=" alt="plot of chunk unnamed-chunk-10"/> </p>

<p>だいぶ良いスコアになりました。
さらに、移動可能なスペースを稼ぐためになるべく直線的に移動するとか、高スコア化するアイデアも考えられますね。皆さんもぜひより良いヘビ駆動アルゴリズムを考えてみてください！</p>

</body>

</html>