stat 545a class meeting #3 wednesday, september 11, 2013jenny/stat545a/2013-lectures/cm...stat 545a...

28
STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department of Statistics and Michael Smith Laboratories Wednesday, 11 September, 13

Upload: others

Post on 13-Sep-2020

0 views

Category:

Documents


0 download

TRANSCRIPT

Page 2: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

any questions from tutorial re: Basic care and feeding of data in R?

HW #1 “Dry Run”. Due as soon as possible. Publish 1 page of something -- practically anything -- to Rpubs, then share that link w/ rest of the class. Instructions/”assignment” is posted but will get some editing. Content not important.

HW #2 “Dry Run 2 + Hello Gapminder”. Due 9:30am Monday Sept. 16. Perform basic data intake and smell testing for the Gapminder data. Write it up via R Markdown. Compile, publish, share. Content can be drawn from tutorials, your job is to pull together. Instructions coming.

Wednesday, 11 September, 13

Page 3: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

today

data analysis vs programming vs developing

combo lecture and hands on computing:R objects, beyond the data.framestudents will complete on their own

R Markdown

Wednesday, 11 September, 13

Page 4: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

three different modes of working with R:

data analysis

programming

developing

Wednesday, 11 September, 13

Page 5: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

three different modes of working with R:

data analysis

programming

developing

I am mostly trying to prepare you for this.

Wednesday, 11 September, 13

Page 6: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

three different modes of working with R:

data analysis

programming

developing

But we need to talk .... about some things more related to this.

Wednesday, 11 September, 13

Page 7: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

hands-on computing phase here

Wednesday, 11 September, 13

Page 8: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

Simple view

Technically correct R viewTechnically correct R viewTechnically correct R viewSimple view

mode class typeof

character character character character

logical logical logical logical

numeric numeric integer or numeric

integer or double

factor numeric factor integer

JB’s typology of R objects, i.e. “flavors”a simple view of simple R objects that will get you pretty far

Wednesday, 11 September, 13

Page 9: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

Simple view

Technically correct R viewTechnically correct R viewTechnically correct R viewSimple view

mode class typeof

character character character character

logical logical logical logical

numeric numeric integer or numeric

integer or double

factor numeric factor integer

JB’s typology of R objects, i.e. “flavors”a simple view of simple R objects that will get you pretty far

Wednesday, 11 September, 13

Page 10: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

JB’s typology of R objects, i.e. “flavors”a simple view of simple R objects that will get you pretty far

Simple view

Technically correct R viewTechnically correct R viewTechnically correct R viewSimple view

mode class typeof

character character character character

logical logical logical logical

numeric numeric integer or numeric

integer or double

factor numeric factor integer

We will spend some time talking about factors in a future class.Wednesday, 11 September, 13

Page 11: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

“Simple view” of data collections

vector

list

Collecting scalars of the same mode?

yes

no

matrix

yes

Collecting vectors of the same length?

no

Collecting vectors of the same mode?data.frame

yes

no

Wednesday, 11 September, 13

Page 12: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

weak links in the chain:process, packaging and presentation

Wednesday, 11 September, 13

Page 13: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

Excerpt from “Real Advice from Real People” by Tom Loughin, Statistical Society of Canada Liaison, Vol. 22 No. 4 (November 2008).

Wednesday, 11 September, 13

Page 14: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

http://daringfireball.net/projects/markdown/

Wednesday, 11 September, 13

Page 15: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

http://www.rstudio.com/ide/docs/r_markdown

http://www.rstudio.com/ide/docs/authoring/using_markdown

http://yihui.name/knitr/

Wednesday, 11 September, 13

Page 16: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

File > New > R MarkdownSave As ... test.rmd

Wednesday, 11 September, 13

Page 17: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

Displays Markdown Quick Reference in the Help tab/pane

Wednesday, 11 September, 13

Page 18: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

File > New > R MarkdownSave As ... test.rmd

Knit HTML.Yes this can be accomplished outside of RStudio, using knitr functions at the command line, so we are not creating unhealthy dependency on RStudio.

Wednesday, 11 September, 13

Page 19: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

This is where you’ll need that RPubs account.Wednesday, 11 September, 13

Page 20: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

You author this. You’ve seen it already.

What’s actually happening?

Wednesday, 11 September, 13

Page 21: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

knit() function from knitr package does R Markdown to Markdown conversion

R chunks are run and replaced by their output.Yes that’s oversimplified but you get the idea.

Wednesday, 11 September, 13

Page 22: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

Title========================================================

This is an R Markdown document. Markdown is a simple formatting syntax for authoring web pages (click the **MD** toolbar button for help on Markdown).

When you click the **Knit HTML** button a web page will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

```rsummary(cars)```

```## speed dist ## Min. : 4.0 Min. : 2 ## 1st Qu.:12.0 1st Qu.: 26 ## Median :15.0 Median : 36 ## Mean :15.4 Mean : 43 ## 3rd Qu.:19.0 3rd Qu.: 56 ## Max. :25.0 Max. :120```

You can also embed plots, for example:

```rplot(cars)```

![plot of chunk unnamed-chunk-2](figure/unnamed-chunk-2.png)

Title========================================================

This is an R Markdown document. Markdown is a simple formatting syntax for authoring web pages (click the **MD** toolbar button for help on Markdown).

When you click the **Knit HTML** button a web page will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

```{r}summary(cars)```

You can also embed plots, for example:

```{r fig.width=7, fig.height=6}plot(cars)```

test.rmd test.md

Wednesday, 11 September, 13

Page 23: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

markdownToHTML() function from markdown package does Markdown to HTML conversion

Pandoc, the Swiss army knife of document conversion, is another option for Markdown to HTML conversion and, importantly, for Markdown to PDF conversion.

You can set things up so that RStudio and knitr use Pandoc for the second conversion step.

Wednesday, 11 September, 13

Page 24: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

<!DOCTYPE html><html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"/>

<title>Title</title>

<!-- 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>

<!-- MathJax scripts --><script type="text/javascript" src="https://c328740.ssl.cf1.rackcdn.com/mathjax/2.0-latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>

<style type="text/css">body { font-family: Helvetica, arial, sans-serif; font-size: 14px; line-height: 1.6; padding-top: 10px; padding-bottom: 10px; background-color: white; padding: 30px; }

body > *:first-child { margin-top: 0 !important; }body > *:last-child { margin-bottom: 0 !important; }

a { color: #4183C4; }a.absent { color: #cc0000; }a.anchor { display: block; padding-left: 30px; margin-left: -30px; cursor: pointer; position: absolute; top: 0; left: 0; bottom: 0; }

h1, h2, h3, h4, h5, h6 { margin: 20px 0 10px; padding: 0; font-weight: bold; -webkit-font-smoothing: antialiased; cursor: text; position: relative; }

h1:hover a.anchor, h2:hover a.anchor, h3:hover a.anchor, h4:hover a.anchor, h5:hover a.anchor, h6:hover a.anchor { text-decoration: none; }

h1 tt, h1 code { font-size: inherit; }

h2 tt, h2 code { font-size: inherit; }

h3 tt, h3 code { font-size: inherit; }

h4 tt, h4 code { font-size: inherit; }

h5 tt, h5 code { font-size: inherit; }

h6 tt, h6 code { font-size: inherit; }

h1 { font-size: 28px; color: black; }

h2 { font-size: 24px; border-bottom: 1px solid #cccccc; color: black; }

h3 { font-size: 18px; }

h4 { font-size: 16px; }

h5 { font-size: 14px; }

h6 { color: #777777; font-size: 14px; }

p, blockquote, ul, ol, dl, li, table, pre { margin: 15px 0; }

hr { border: 0 none; color: #cccccc; height: 4px; padding: 0;}

body > h2:first-child { margin-top: 0; padding-top: 0; }body > h1:first-child { margin-top: 0; padding-top: 0; }body > h1:first-child + h2 { margin-top: 0; padding-top: 0; }body > h3:first-child, body > h4:first-child, body > h5:first-child, body > h6:first-child { margin-top: 0; padding-top: 0; }

a:first-child h1, a:first-child h2, a:first-child h3, a:first-child h4, a:first-child h5, a:first-child h6 { margin-top: 0; padding-top: 0; }

h1 p, h2 p, h3 p, h4 p, h5 p, h6 p { margin-top: 0; }

li p.first { display: inline-block; }li { margin: 0; }ul, ol { padding-left: 30px; }

ul :first-child, ol :first-child { margin-top: 0; }

dl { padding: 0; }dl dt { font-size: 14px; font-weight: bold; font-style: italic; padding: 0; margin: 15px 0 5px; }dl dt:first-child { padding: 0; }dl dt > :first-child { margin-top: 0; }dl dt > :last-child { margin-bottom: 0; }dl dd { margin: 0 0 15px; padding: 0 15px; }dl dd > :first-child { margin-top: 0; }dl dd > :last-child { margin-bottom: 0; }

blockquote { border-left: 4px solid #dddddd; padding: 0 15px; color: #777777; }blockquote > :first-child { margin-top: 0; }blockquote > :last-child { margin-bottom: 0; }

table { padding: 0;border-collapse: collapse; }table tr { border-top: 1px solid #cccccc; background-color: white; margin: 0; padding: 0; }table tr:nth-child(2n) { background-color: #f8f8f8; }table tr th { font-weight: bold; border: 1px solid #cccccc; margin: 0; padding: 6px 13px; }table tr td { border: 1px solid #cccccc; margin: 0; padding: 6px 13px; }table tr th :first-child, table tr td :first-child { margin-top: 0; }table tr th :last-child, table tr td :last-child { margin-bottom: 0; }

img { max-width: 100%; }

span.frame { display: block; overflow: hidden; }span.frame > span { border: 1px solid #dddddd; display: block; float: left; overflow: hidden; margin: 13px 0 0; padding: 7px; width: auto; }span.frame span img { display: block; float: left; }span.frame span span { clear: both; color: #333333; display: block; padding: 5px 0 0; }span.align-center { display: block; overflow: hidden; clear: both; }span.align-center > span { display: block; overflow: hidden; margin: 13px auto 0; text-align: center; }span.align-center span img { margin: 0 auto; text-align: center; }span.align-right { display: block; overflow: hidden; clear: both; }span.align-right > span { display: block; overflow: hidden; margin: 13px 0 0; text-align: right; }span.align-right span img { margin: 0; text-align: right; }span.float-left { display: block; margin-right: 13px; overflow: hidden; float: left; }span.float-left span { margin: 13px 0 0; }span.float-right { display: block; margin-left: 13px; overflow: hidden; float: right; }span.float-right > span { display: block; overflow: hidden; margin: 13px auto 0; text-align: right; }

code, tt { margin: 0 2px; padding: 0 5px; white-space: nowrap; border: 1px solid #eaeaea; background-color: #f8f8f8; border-radius: 3px; }

pre code { margin: 0; padding: 0; white-space: pre; border: none; background: transparent; }

.highlight pre { background-color: #f8f8f8; border: 1px solid #cccccc; font-size: 13px; line-height: 19px; overflow: auto; padding: 6px 10px; border-radius: 3px; }

pre { background-color: #f8f8f8; border: 1px solid #cccccc; font-size: 13px; line-height: 19px; overflow: auto; padding: 6px 10px; border-radius: 3px; }pre code, pre tt { background-color: transparent; border: none; }

sup { font-size: 0.83em; vertical-align: super; line-height: 0;}* {! -webkit-print-color-adjust: exact;}@media screen and (min-width: 914px) { body { width: 854px; margin:0 auto; }}@media print {! table, pre {! ! page-break-inside: avoid;! }! pre {! ! word-wrap: break-word;! }}</style>

</head>

<body><h1>Title</h1>

<p>This is an R Markdown document. Markdown is a simple formatting syntax for authoring web pages (click the <strong>MD</strong> toolbar button for help on Markdown).</p>

<p>When you click the <strong>Knit HTML</strong> button a web page will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:</p>

<pre><code class="r">summary(cars)</code></pre>

<pre><code>## speed dist ## Min. : 4.0 Min. : 2 ## 1st Qu.:12.0 1st Qu.: 26 ## Median :15.0 Median : 36 ## Mean :15.4 Mean : 43 ## 3rd Qu.:19.0 3rd Qu.: 56 ## Max. :25.0 Max. :120</code></pre>

<p>You can also embed plots, for example:</p>

<pre><code class="r">plot(cars)</code></pre>

<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAGwCAYAAABFI3d+AAAEJGlDQ1BJQ0MgUHJvZmlsZQAAOBGFVd9v21QUPolvUqQWPyBYR4eKxa9VU1u5GxqtxgZJk6XtShal6dgqJOQ6N4mpGwfb6baqT3uBNwb8AUDZAw9IPCENBmJ72fbAtElThyqqSUh76MQPISbtBVXhu3ZiJ1PEXPX6yznfOec7517bRD1fabWaGVWIlquunc8klZOnFpSeTYrSs9RLA9Sr6U4tkcvNEi7BFffO6+EdigjL7ZHu/k72I796i9zRiSJPwG4VHX0Z+AxRzNRrtksUvwf7+Gm3BtzzHPDTNgQCqwKXfZwSeNHHJz1OIT8JjtAq6xWtCLwGPLzYZi+3YV8DGMiT4VVuG7oiZpGzrZJhcs/hL49xtzH/Dy6bdfTsXYNY+5yluWO4D4neK/ZUvok/17X0HPBLsF+vuUlhfwX4j/rSfAJ4H1H0qZJ9dN7nR19frRTeBt4Fe9FwpwtN+2p1MXscGLHR9SXrmMgjONd1ZxKzpBeA71b4tNhj6JGoyFNp4GHgwUp9qplfmnFW5oTdy7NamcwCI49kv6fN5IAHgD+0rbyoBc3SOjczohbyS1drbq6pQdqumllRC/0ymTtej8gpbbuVwpQfyw66dqEZyxZKxtHpJn+tZnpnEdrYBbueF9qQn93S7HQGGHnYP7w6L+YGHNtd1FJitqPAR+hERCNOFi1i1alKO6RQnjKUxL1GNjwlMsiEhcPLYTEiT9ISbN15OY/jx4SMshe9LaJRpTvHr3C/ybFYP1PZAfwfYrPsMBtnE6SwN9ib7AhLwTrBDgUKcm06FSrTfSj187xPdVQWOk5Q8vxAfSiIUc7Z7xr6zY/+hpqwSyv0I0/QMTRb7RMgBxNodTfSPqdraz/sDjzKBrv4zu2+a2t0/HHzjd2Lbcc2sG7GtsL42K+xLfxtUgI7YHqKlqHK8HbCCXgjHT1cAdMlDetv4FnQ2lLasaOl6vmB0CMmwT/IPszSueHQqv6i/qluqF+oF9TfO2qEGTumJH0qfSv9KH0nfS/9TIp0Wboi/SRdlb6RLgU5u++9nyXYe69fYRPdil1o1WufNSdTTsp75BfllPy8/LI8G7AUuV8ek6fkvfDsCfbNDP0dvRh0CrNqTbV7LfEEGDQPJQadBtfGVMWEq3QWWdufk6ZSNsjG2PQjp3ZcnOWWing6noonSInvi0/Ex+IzAreevPhe+CawpgP1/pMTMDo64G0sTCXIM+KdOnFWRfQKdJvQzV1+Bt8OokmrdtY2yhVX2a+qrykJfMq4Ml3VR4cVzTQVz+UoNne4vcKLoyS+gyKO6EHe+75Fdt0Mbe5bRIf/wjvrVmhbqBN97RD1vxrahvBOfOYzoosH9bq94uejSOQGkVM6sN/7HelL4t10t9F4gPdVzydEOx83Gv+uNxo7XyL/FtFl8z9ZAHF4bBsrEwAAQABJREFUeAHtnQn8ZXP9/5swhgyRPcsYWcv2ozRKSNTfUpYoRMlaaZFKi5TK1oKSRKVIo0RRJBRjiKQMmlQY65Bdg7Ez/+dr3Ks797vd73fucu65z/fj8fp+lnPOZ3l+zj3v7+eso2bNmvUyTQISkIAEJCCBchF4ebm6Y28kIAEJSEACEggBHbz7gQQkIAEJSKCEBHTwJRxUuyQBCUhAAhLQwbsPSEACEpCABEpIQAdfwkG1SxKQgAQkIAEdvPuABCQgAQlIoIQEdPAlHFS7JAEJSEACEtDBuw9IQAISkIAESkhAB1/CQbVLEpCABCQgAR28+4AEJCABCUighAR08CUcVLskAQlIQAIS0MG7D0hAAhKQgARKSEAHX8JBtUsSkIAEJCABHbz7gAQkIAEJSKCEBHTwJRxUuyQBCUhAAhLQwbsPSEACEpCABEpIQAdfwkG1SxKQgAQkIAEdvPuABCQgAQlIoIQEdPAlHFS7JAEJSEACEtDBuw9IQAISkIAESkhAB1/CQbVLEpCABCQgAR28+4AEJCABCUighAR08CUcVLskAQlIQAIS0MG7D0hAAhKQgARKSEAHX8JBtUsSkIAEJCABHbz7gAQkIAEJSKCEBHTwJRxUuyQBCUhAAhLQwbsPSEACEpCABEpIQAdfwkG1SxKQgAQkIAEdvPuABCQgAQlIoIQEdPAlHFS7JAEJSEACEtDBuw9IQAISkIAESkhAB1/CQbVLEpCABCQgAR28+4AEJCABCUighAR08CUcVLskAQlIQAIS0MG7D0hAAhKQgARKSEAHX8JBtUsSkIAEJCABHbz7gAQkIAEJSKCEBHTwJRxUuyQBCUhAAhLQwbsPSEACEpCABEpIQAdfwkG1SxKQgAQkIAEdvPuABCQgAQlIoIQEdPAlHFS7JAEJSEACEtDBuw9IQAISkIAESkhAB1/CQbVLEpCABCQgAR28+4AEJCABCUighAR08CUcVLskAQlIQAIS0MG7D0hAAhKQgARKSEAHX8JBtUsSkIAEJCABHbz7gAQkIAEJSKCEBHTwJRxUuyQBCUhAAhLQwbsPSEACEpCABEpIQAdfwkG1SxKQgAQkIAEdvPuABCQgAQlIoIQEdPAlHFS7JAEJSEACEtDBuw9IQAISkIAESkhAB1/CQbVLEpCABCQggXl7CcGoUaPeTX97qs+9NL72VQISkEABCdw/a9asSzrRrlFU3Il6214nzn1HKj0Indr2yq1QAhKQgAR6lcDH6Phu+Nrr2g2gcLNZHHHaNBYYjzQZRso9jXJPanK5FicBCUhAAhLolwA+bVUWdORyeEcqracAgNHoCHQXy55BDxOfiaaiPevXNy0BCUhAAhKQwOAEijKDP55mLo22RreimWhhtCY6Dic/hpn3icQ1CUhAAhKQgAQaIFCIGTzt3BLthxO/AT2OYjPQVeR/HG2HNAlIQAISkIAEGiRQFAc/lfZuNkCbtyH/gQGWmS0BCUhAAhKQQD8EinKK/lDaNpFT8QcSTkOPokXQGiht3AppEpCABCQgAQk0SKAQDp5T8VNw7uvR5gloHMr1+Mzac919cs7XE2oSkIAEJCABCTRIoBAOPm3Fhz9FcGl9u3H886B5Wf50/TLTEpCABCQggVoC+ItRpD+KtkXzoefQNhUfQ7R3rBDX4BmP5dFp6HF0MXpNzRDsRPynNWmjEpCABCQggYEIHM6CldBWOPVNCc9E38GvzEPYU1aUGXyuvf8HbYB2RZMZjE0ZnJuIN2xsM46Vlx9gg/8jP8/YaxKQgAQkUF4Cm9O1N+E/MnPP2eGT8Q3rEF0X/S15vWJFcfC5iW49BuJJwkMZjBsJLyR88zAHYhXW33iAbd5Cfm7e0yQgAQlIoLwEHqo695ou5hLwmJp0T0SL4uDj0DN7vzzUGZyf49yXJXoBavjVsmx3MetHfYzyxpK5VJ8FZkhAAhKQQJkI3MTx/iv4gzyd9TLimdF/En0Z9ZQVxcF/H+q/ZCCOZVCOzggQHlNxyseSPCd5mgQkIAEJSGAIAl9k+bX4jxUI70WvRePxKY8R9pQVwsED/iIGY2XIj6+lT/5h5F9GXpZpEpCABCQggUEJxJHjN1ZnpZwVHo0O70XnHkiFcPBpCAOQ98//PfFaI38S6UiTgAQkIAEJDEkAv/E8K1095IolX6EQj8mVnLHdk4AEJCABCbSdQCFm8JxOOYie54UEA9m/+I/M6/AD0TFfAhKQgAQkUEegEA6eNo1DB6BTUU7V15sfm6knYloCEpCABCQwCIFCOHhm5x9lFp/LBS8n/pFB2usiCUhAAhKQgAQaIFCka/AH096FcfQLNdBuV5GABCQgAQlIYBAChZjBp33M3B8n2G2QtrpIAhKQgAQkIIEGCRRpBt9gk11NAhKQgAQkIIGhCOjghyLkcglIQAISkEAXEtDBd+Gg2WQJSEACEpDAUAR08EMRcrkEJCABCUigCwno4Ltw0GyyBCQgAQlIYCgCOvihCLlcAhKQgAQk0IUEdPBdOGg2WQISkIAEJDAUAR38UIRcLgEJSEACEuhCAjr4Lhw0mywBCUhAAhIYioAOfihCLpeABCQgAQl0IQEdfBcOmk2WgAQkIAEJDEVABz8UIZdLQAISkIAEupCADr4LB80mS0ACEpCABIYioIMfipDLJSABCUhAAl1IQAffhYNmkyUgAQlIQAJDEdDBD0XI5RKQgAQkIIEuJKCD78JBs8kSkIAEJCCBoQjo4Ici5HIJSEACEpBAFxLQwXfhoNlkCUhAAhKQwFAEdPBDEXK5BCQgAQlIoAsJ6OC7cNBssgQkIAEJSGAoAjr4oQi5XAISkIAEJNCFBHTwXThoNlkCEpCABCQwFAEd/FCEXC4BCUhAAhLoQgI6+C4cNJssAQlIQAISGIqADn4oQi6XgAQkIAEJdCGBwjn4UaNGzYsW7UKWNlkCEpCABCRQGAKFcPA49NHoCHQXZJ5BDxOfiaaiPQtDy4ZIQAISkIAEuoTAvAVp5/G0Y2m0NboVzUQLozXRcTj5MbNmzTqRuCYBCUhAAhKQQAMECjGDp51bov1w4jegx1FsBrqK/I+j7ZAmAQlIQAISkECDBIri4KfS3s0GaPM25D8wwDKzJSABCUhAAhLoh0BRTtEfStsmcir+QMJp6FG0CFoDpY1bIU0CEpCABCQggQYJFMLBcyp+Cs59Pdo8AY1DuR6fWXuuu0/O+XpCTQISkIAEJCCBBgkUwsGnrfjwpwguxdGnTWNJP5J8TQISkIAEJCCB4RMoxDV4nLqPyQ1/7NxCAhKQgAQkMCCBoszgm/KYHP8ojKOnyw3Q22XJHz3AMrMlIAEJSEACpSJQFAefx+QmcFr+3hq6M4hfhdPOY3KHoUaeg1+F9d6C+rMVyUyZmgQkIAEJSKD0BIri4KuPyZ3RD/GGH5PjH4SL2T7qY/yjsBCZS/VZYIYEJCABCUighASK4uB9TK6EO5ddkoAEJCCBzhEohINn5u1jcp3bB6xZAhKQgARKSKAQDj5ccfKzH5MrIWO7JAEJSEACEmg7gUI8Jtf2XluhBCQgAQlIoOQECjGD5wa4g+A83yCs/8UM/5xBlrtIAhKQgAQkIIEaAoVw8LRnHDoAnYryqdh682Mz9URMS0ACEpCABAYhUAgHz+z8o8zic7ng5cQ/Mkh7XSQBCUhAAhKQQAMEinQN/mDau3DlefUGmu4qEpCABCQgAQkMRKAQM/g0jpn74wS7DdRQ8yUgAQlIQAISaJxAkWbwjbfaNSUgAQlIQAISGJSADn5QPC6UgAQkIAEJdCeBwpyi7058tloCEpCABIpGgHu5lqdN26N8QXQSl4D/WrQ2tqM9zuDbQdk6JCABCUigLQRw7qtTUR65fghNQ9eQtx1hz5kz+J4bcjssAQlIoNQEzqR3uWH7WTQGLYl+gJP/KzP56cR7xnTwPTPUdlQCEpBATxC4nV7ug8aiOPn3oe+jcainHLyn6BlxTQISkIAESkPgjfTkbmbre6J9icfZ74X6e0sq2eU1Z/DlHVt7JgEJSKAXCdxDp9/KKflJhE+gOPd/otxw11Omg++p4bazEpCABEpP4G56+Hm0J1oAnYTehp5DPWWeou+p4bazEpCABEpP4Ap6uCM6kFP0+xGugt6BrkM9ZTr4nhpuOysBCUig9AS+Tg+XQ5dxmv7sSvx1OPvnS9/zug56ir4OiEkJSEACEuheAhVH/kGc+2L0Ij7uAfJmdW+PRt5yHfzI2bmlBCQgAQkUlAA+/eGCNq1tzfIUfdtQW5EEJCABCUigfQR08O1jbU0SkIAEJCCBthHQwbcNtRVJQAISkIAE2kdAB98+1tYkAQlIQAISaBsBHXzbUFuRBCQgAQlIoH0EdPDtY21NEpCABCQggbYR0MG3DbUVSUACEpCABNpHQAffPtbWJAEJSEACEmgbAV900zbUViQBCUhAAt1EgLfhrUB7P4oWRbeioytvyiNafHMGX/wxsoUSkIAEJNBmAjj3RajyTPQXdASagX5Dfr5Q1xXmDL4rhslGSkACnSDAwXwU9W6JMoObxuztmk60wzo7QuCL1HoUY35OpfYT2B1eRXw39MNKXqEDZ/CFHh4bJwEJdJhAviW+Pcps7pcc4A/pcHusvn0EFqaqf9RVdy3pV9blFTapgy/s0NgwCUigkwRw5vtT/4PM4PZHcfSrognkb97Jdll32whMo6b31NX2BdK31OUVNukp+sIOjQ2TgAQ6TGBD6j+q2gac/DM49x+Q3gD9sZpvWFoCx9OzqYx5vi3/O7QduqXmlD3JYpsz+GKPj62TgAQ6R+Bxql6+rvrVSc+syzNZQgI48ifo1ioo912sjM4lb3fCrjFn8F0zVDZUAhJoM4ETqe/7zOAeIvwX2hEdiRZCWg8QwKE/Tzd/1K1d1cF368jZbglIoKUEOLjfiHPfl0qOQS+gaWgZ8p3BA0IrPoHCOXh+UGnTWH5EjxQfny2UgATKTIDjUGbuW5W5j/atvAQKcQ0epz4aHYHuAvUz6GHiM1FucNizvPjtmQQkIAEJSKA1BIoyg8/dikujrVFeB5hTYHkGcU10HE5+DP9J53qYJgEJSEACEpBAAwQKMYOnnXlT1H448RvQ4yg2A11F/sdRHk/QJCCBJhLgH+dd0R/QJehytFgTi7coCUigwwSK4uCnwmGzAVhsQ/4DAywzWwISGAEBnPm72exdaDv+kX4r4dHox+S/glCTgARKQKAop+gPheVEDi4HEuZO1UdRXg25BkobG7rJhe1XYN28lKA/W4bM+ftbYJ4EepDA3vR5L5x7nvV+GeF5/H7WJboFOid5mgQk0N0ECuHgObhM4eCyHignoHEo1+MfRD9BF7J8FmEjlpdQbDLAiuPInzHAMrMl0IsEHqvrdP6xHlOXZ1ICEuhSAoVw8Dj33FCX503jnPPMaV4skZcLZMb9c5Z/Eh8/e6ZBekBjnYtYGPUxyliQzKX6LDBDAr1JYBLd/hbaJ93n9/Eagm+jcUiTgARKQKAQDh6On0U5wOR9vznIpF3vRP9Gcfg7oR8jTQISaA6B/K5yWv48wiloffQm/km+g1CTgARKQKAoDj43+7yBg0uefc8se3HiuYM+M4sjCXIw0sEHiCaBJhDg95X3TWzJ72tDwoXRd8m7rwlFW4QEJFAQAkVx8P+ExxYcbC4lfAtaoIbP2sSvrUkblYAEmkQAp351k4qyGAlIoGAEiuLgcy3wFDQefQeNxdnH6V+P3ow2RZoEJCABCUhAAg0SKISDz+l4HHreWrcY8YeI5+a6t6NXoj3Je5JQk4AEJNBWAhyLxlLhB1GORdM4Fp3e1gZYmQTmgsDL52Lbpm7KDyeWu+fzTO7T6DfoNKRzbyppC5OABBohgHMfzXq/RvOhy9Ae5OXGxMIcN2mTJoEBCbijDojGBRKQQI8T+Az9P4NJxjfRJJRXat+G8l14TQKFJ6CDL/wQ2UAJSKBDBJal3sl1dV9AeqC3ZdatalICnSWgg+8sf2uXgASKSyBvvnxjXfPeQ/rhujyTEigkgULcZFdIMjZKAhLodQLHAeBerrkvTngJ2hmtgz6INAkUnoAz+MIPkQ2UgAQ6QYBr7nnxz0JoEfQB9B+0IfnPE2oSKDwBZ/CFHyIbKAEJdIoAznwmdX+5U/VbrwTmhoAz+Lmh57YSkIAEJCCBghLQwRd0YGyWBCQgAQlIYG4I6ODnhp7bSkACEpCABApKQAdf0IGxWRKQgAQkIIG5IaCDnxt6bisBCUhAAhIoKAEdfEEHxmZJQAISkIAE5oaAj8nNDT23lYAEJCCBwhHg5URL0qhtUL5MOpnHHf9RuEa2oUHO4NsA2SokIAEJSKA9BHDu46jpZyhfAXwMTSVvC8KeM2fwPTfkdlgCEpBAqQmcS+/2Ydb+l/QS534xwQ8Ip5KXtxH2jDmD75mhtqMSkIAEeoLAHVXnnt4SzyuHp6KVk+4l08H30mjbVwlIQALlJ7AAs/VXVbtJfBTxnKL/bzWvV0IdfK+MtP2UgAQ6RgAfE6fzWfQLdCoa17HGlL/i79PFU2D8OjSe+Pnor8zkM4vvKfMafE8Nt52VgATaTQAnk4nU39DZ6CC0IvoZ+R/C6dxAXGsiAZieDdv7KfKTKDfanYV+jHrOdPA9N+R2WAISaDOBd1PfX3A8X6zUOx0H9BniB6I9K3kGTSQA68spLupp8xR9Tw+/nZeABNpAIN+Tv6Sunn+TTr4mgZYR0MG3DK0FS0ACEphNIC9Z2Z5Z++gaHu8kfndNuquj9C33GCyHxnZ1R0rWeE/Rl2xA7Y4EJFAsApwuvhLH91dadSvhpwlXQ+ujnVDXG33amE7k3oIZaBPSH6PPv+n6jpWgAzr4EgyiXZCABIpNAId3OI7vClq5FroVHU3eU8Vu9dCto095tnwyWpn+5B+YpYhfR3g/6T8PXYJrtJKADr6VdC1bAhKQQIUADu8yolGZbHs6s32cezpFeB/OfWei2yIdfKB00LwG30H4Vi0BCUigywnkMbSn6/rwLGknj3VQOpF0EDpB3TolIAEJjJAAM+TF2fTN6AU0iVnzoyMsqhmb/YFCfkKb8hjgQ4QLkL4Kbd6Mwi1j7ggMOINnoMai2rs+89L+MWiluavSrSUgAQlIYCQEOP6uznbnoRyHX4dmkDeOsCOGU7+Gir+Ect39WMIz0Z7k1z8WSLbWbgJ9ZvAMUpx6HP/X0O9JX1rTqO2I7462rskzKgEJSEACLSbAsXhBqrgavQsHOinVkXcdwVGEe5D3TPLabdR7FvXn3oJl0EOkS/P4X7tZNru+/mbwH6SSJ9HH0O8q8aSjk1He66tJQAISkEB7CeSO9TOrzj1VE88x+nm0XNKdMtrxALoB6dw7NQj91NvHwTNAeVF/bpzIc40TKvGk52XZwuh7xDUJSEACEmgvgcepbonaKpk55xieU/VZpklgDgJ9HHyW4sSfQ8cQvQuNQrlx4gB2pjwSoUlAAhKQQJsJcEy+jSqv4Tj8DbQoirPPC3TOY1k+rqJJYA4Cfa7BV5ey82xE/GKUty4dijZAo8lfjJ3pR8Q1CUhAAhJoL4EjqO6L6DSUx9O+i3ryS2n0WxuCwIAOnu3eh/ZG96G8uCCPZWT9I1HLHDz/QKSOsfwT8QihJgEJSEACFQIcF2cR/YpAJNAIgX5P0Vc2zJeOHkAbo7x2cCrh/Kjpz1zi1HNm4AiUSwK5E/Rh4jPRVLQnaU0CEpCABCQggWEQGGwGn7vlj0PzoLzIYE3CU9HhqNl2PAUujfL4XV55OBMtjFLncdQ9hn8wTiSuSUACEpCABCTQAIEBHTwOdSKO9UHKeCU6C41HHyK/9rl4sppiW1LKBMq+t6a0fJnoKtrwccLDkA6+Bo5RCcwtAX5b76WMfVD+ic/ZvHfyG/wvoSYBCZSAwIAOPn3jx35RTR9vIR61wnL6fzN0Rj+Fb0NeLhVoEpBAkwjg3HekqDwVsx2/88dI5+MgpxLuQvqJJlVjMRKQQAcJ9HHw/MDzpqTPoNxF39/17ws4AGRW3UzLXfo5Y3Ag4TSU6/y5B2ANlDZuhYY0tl+BlV49wIp5y9Icr94dYD2zJdALBDJz3yvOPZ0l/C2/n3WJ5mzaOcnTJCCB7ibQx8HTnX3R7SiOtnYGT3K2PVKNNCvk4DKFg8t6lJcX64xDuR6fek5Gl7I8d442YvmHYJMBVlyJfE8/DgDH7J4jkN9U/Q2zuSw2pudI2GEJlJRAHwePL72+0tf82Ke3o9+VmffXqOs76Dx0FMopxLxB7xcs3592PUN8UGOdC1kh6mOUcSyZS/VZYIYEepPAZXQ7L7PKTD7vNF+Z4NtoHNIkIIESEOjj4PmhX0K/XjlI3y7GkR48yPKRLMpznXeif6DPobTrdSiP5cXZf7EiAk0CEmgCgTj3C/i9/5bwWrQh2pjf9h2EmgQkUAICfRw8ffoCSv4q6BCUd89fifLI2ofRFNRsewsFrp5ZOgec7Ynnxp/pqYR0nHvej69JQAJNIpDfGkVtzu/rjYRj0ffIy0utNAlIoCQE8mjMHMaP/Cp0OZkboC8R/ya6Ev2Q9AFod9Rsu4kC96gUOomw9qa63EV/c2WZgQQk0EQC/K7/jHJWTufeRK4WJYEiEOhvBl9tV+6uHVdNVMKcNn+wLq8ZyY9QyHnMJvYizKN43yT+QcIXUF54kxm+JgEJSEACEpBAgwQGc/CnUMbvcbTvIPwLWh+NR5lRN9WYPUyjnlwC2AKthnI9PnfRZ+Z+PsufI9QkIAEJtJUAx6WFqDCTjdyXdAvHooltbUA/ldGm9ch+G8oE6FzalEmRJoE+BAZ08Ow0/2ZHyo0370K5Hn86ymcJ7yFsulFuHtu5qKKml2+BEpCABIZDgOPfaNb/FfoDugJ9lrxdCN/F4SrOte1G/e+k0m+gg1CO3zeTtyntyVMRmgTmIDCgg89a7DT3E/xgji1MSEACEugNAp+mm7/kOFg9Bl6CMz2BvB2T324E1P0q6swNx/9Hm+5N/eStRPBtwmvJeyx5mgSqBF5ejRhKQAISkMAcBJYjNWmOnJe97Hekk98JW5JKz6o69zSA+O0EeZX3YkiTwBwEdPBz4DAhAQlI4CUCefPlhJdSL0beTZD7gzphqXc8s/UFq5VXZvWbk/YtnVUohi8RGPQU/UtrGZGABCQwAAGczAosysuvxlVWOYCZ5W2VeDcHx9H4e+nfooSXop3Q+mhv1HbLzJ22/IKKZxLm/qjcIzARHcSyGYRdbfQpLzb7FNoYZfKZdzOcQ6iNkIAOfoTg3EwCEph9DXgRONyBtkO5SXYddBoH6/dzcL6VeNca7b+PfuQu+nx8q/oI7+vJf75TnaLun9Kmu6h/B5Qbk3cl74pOtafJ9V5CeVNQzpJkvzqOvs5H/9p+vwN1l8JGAa8UHRmqE+wox7LOUvR316HWdbkEJNAYAX5Xn2XNx/ldfbe6BXnvIv5G8vLaaU0CQxJgn8ms/RPsM7mBcbaRlzMnvyBvy0pWVwb041s0/Gf0I6+Ebqt5Db6tuK1MAqUjkOvB/6rr1XTSmflqEmiUQPaX62tXxiE+QnpMbZ7x4RHQwQ+Pl2tLQAJzEriG5AHMUkbVZJ9CPPlaCwiAeh60OloN1XJvQW1tK/IGatqI7qxYrZF4TtU/XU0bDp+A1+CHz8wtJCCBCgFmWb/lQPwWktcT5sNUb0DXkH9aZRWDJhKAca5N53Ljsyiz3vXIewO8HyfetUb776YfR9KB2wkPJMzp+bVRnLw2QgLO4EcIzs0kIIEXCXBwzgth9kUPodNJ7/3iEv82kwCObx7Km4puhvF+aDfiv0VHo643+nMZnRiH7kF/Q7mBcAahNkICzuBHCM7NJCCB/xHgQPxnUpHWOgLjKXoSrDPTnW3ED8bxn4GWJR7H2NVGH+6gA5HWBALO4JsA0SIkIAEJtIFAHs97rp96Xklexx7d66c9ZhWEgA6+IANhMyQggeYRYEa7BPoA2h+t0bySO1cSs9u8VyAvuTm52griZxDP4873VfMMJVAloIOvkjCUgARKQQCnN46OTEQLoJnoRvLeRlgGy/0Oa9Gf3Nx4NvG89CbvHdAk0IeA1+D7IDFDAhLocgLn0v59mNX+Jf3AEeZzrycRTiVv9lfYkt+NRvufpN0T6MsShCRnPdiN/bDN7SHgDL49nK1FAhJoH4E7cXyznXuqJP4fgtx9vkrSZTD69IDOvQwj2do+6OBby9fSJSCB9hMYwww3z1HPNuKjiOQU/X9fzPGvBHqDgA6+N8bZXkqglwicRGd/gl9/LRpH/Dx0LTPevxN2vdGnndAPK9qm6ztkB1pGwGvwLUNrwRKQQCcI4MjPwvndT92fQqNRbkb7Cep6o19fpBP5ct97UY7fnycvH9H6EXFNAnMQ0MHPgcOEBCRQBgI4vMn0IyqN4ciXozN5e91r6d/s597J25v0bwjPJc8b7koz2s3piKfom8PRUiQgAQm0mkBeaHNp1bmnMuL5GEuegR+btCaBWgI6+FoaxiXQJQSYsb0q6pLm2kwIMF6j0fIojnoklmfel2H71aobE8/HfTZHeVJAk8AcBDxFPwcOExIoNgEO6Pn++uEop2sXJp2viu3ATO4ZQq2gBBin19C0I1BOrecLcCcwZscPp7msP4PtUsa/CN9PmOP3wWgrlj1FqElgDgLO4OfAYUIChSdwGS2cnwP6TujtxPMSl68VvtU93ECc8WJ0/2Z0ImO2C2E+g5qX1exIOCxj+zzfvxLKzYOz0ObkXU+oSaAPAR18HyRmSKCYBHAI42jZHRzQP1xtIfHjiK/IsiWreYaFI7AFLfoEY3VpWkaYsy1fQCN6xI3tb0c/RD9G01OmJoH+COjg+6NingSKSSCztv/207Scqs23wrViEpifZuWd+LWWSyuvqM0wLoFmE9DBN5uo5UmgdQSmUfTLma1/ploF8W8QX4WZnDdZVaEUL5xEkz7CWL02TSMcRZAZfC63aBJoGYH8569JQAJdQAAn/jy+4SCa+g/CNQnziFScxUZIKygBxu1Oxmt/mncl4amEy6OryT+hoE22WSUhoIMvyUDajd4ggFN4BCcRB7Eqyk1WN5H3Qm/0vnt7yRhdzbitQA/Go0dJ52yMJoGWEtDBtxSvhUug+QRwDs9T6j+bX7IltpIA4zaD8qe0sg7LlkAtAa/B19IwLgEJSEACEigJAWfwJRlIuyEBCTSfAKfV8/hhHktcBN2Kvlc5g0JUk0CxCTiDL/b42DoJSKBDBHDuC1L1r9Cd6Mcor4W9hvz5CDUJFJ5A4Rw8P5550aKFJ2cDJSCBshM4hA4ez4z9FHQD2p30BWjXsnfc/pWDQCEcPA49H2E4AuVjCnnL08PEZ6KpaM9yoLYXEpBAlxHIRKP+prg/kZdXz2oSKDyBQjh4KOWjC3kJxNZoYZR2LYv2Qfvj5D9EqElAAhJoJ4G7qewddRUeQHp6XZ5JCRSSQFFustsSOhM4BXZvDaU8UnIVzv3jhIehE2uWGZWABEpGgN/6QnTpK2hdlH/0f8ox4duEwzLKyct/Mil4F1oA5ViyI2Xl7OBw7DusfDvFZbLxe7QTophZvyTUJFB4AkWZwU+F1GYD0MoHGR4YYJnZEpBACQjgRPMu/evRDBzoWwnfgtYmfySX6L7Etmugd6JNUJzztylrWMc72vEo2y2N7kEp50qUMjUJdAWBUezEHW8oP7z1aMRE9BjKG57yw8pjKfmR5ixDvnd8B+GgRjnLs8KrB1gpZwLymc0dBlhutgQk0CEC/Hbj0Pfg97l3tQnkZUZ/Nnn5LG7DxnbXsvIb2O656kbk5Qzgj8j7azXPUALtIMC+9y3q+Rn7XvbLtlohTtHT8SlAiJOfgMah/Nf8CDoV/Z7ljf4Xkuv4+U+7P1uZzJSpSUACxSMwhiZNr2tWvsCWU+zDtfs5ZLzk3Csbp6x81U2TQM8QKISDD21+kE8RXFolj8PfjXium+WxlIaMdXMqLupjlHcsmUv1WWCGBCRQBAKZ3Xw6/+jzO55SaVBusu3v87hDtfdmyjmEcr6WFYnnn/58pOewpDUJ9AqBQjh4foA3A3zxOuijSeeZ+B0Jz+HHumfdcpMSkEBJCPD7fpDf+qfoznWEccwLo7EoN7YN1/L8+t8pZ0XCfEb39WhV6sglQE0CPUOgEA4e2nHep6DTUU7Lx7ZDOWV/MMrpNU0CEigxARzw9Tjl/KP/f+hJdBV5+bDOsIxtZlDOSmyU40cmCseQN5IzAcOq15UlUDQChXDw/Piu4Ae5AXC+i45B+6EH0eMsu4NQk4AEeoAAv/eH6ObFc9tVysk/BlfMbTluL4FuJlAIBx+A/CBz5/weOPqdCSejq9Gw/3tnG00CEpCABCTQ8wQK4+CrI4GjPxMnn+dNT0R5LlaTgARaQIDf2RIUuwfKte5r+O2d34JqLLJCAN7bE10b5V6AE+GdyxCaBFpGYFgvfmhZK+oKZsefjrZFuVlGk4AEmkyg4tzPpNgZKO9XP4m8rzS5GourEIDtSUTfja5CuS/g8coYENUk0BoChZvBt6ablioBCdQROIH0UfwTfWHycTa54/x0wjeS9+fkac0hANO3UdL6cM19RrGLyMtruT+NPpMMTQKtIFDIGXwrOmqZEpDAHATyYpnc6zLbcD653yXvoVjhxRz/NpFA3q6Z93DUWi6HyLqWiPGmE9DBNx2pBUqgKwg8TCurM8pqg99LJHexa80lkDdo5pG9WluHxHA/flO7vXEJDElABz8kIlfoRgKcAt0ZnYV+j45B+ZhJ11v6gSai6ehelOvoI7Gvs9Fktt8RrYvy/okHmcn/cSSFuc2gBH7L0vEwvhitg3It/jCUU/SaBFpGQAffMrQW3CkCHED3p+68KOmj6D3ocXQ0KoPlNa7vQJtWtBr9/RXxYRmO/EY2WBptiPZFk1Bm8FqTCcB6FkVujfKZ2bzj4/VoF7LvI9Qk0DIC3mTXMrQW3EECcewbcACtPoZ0KE7wx6j2PecdbN7Iqqb9G7PlePqVr6zNNvLeQCSz+TXJj9Nu2CoOxpu8GiY28hVhHSd/8shLcEsJDJ+AM/jhM3OL4hO4g+Np1blXW5u7lvO8dzfbq2j8nbUdoJ9Pk84ZiiVq841LQAIS0MG7D5SRwH+Y0eYFLrMts1siud459cWcrv17OS1fhP7kEsRsI/41IsujUnznnP6MQWuh17zYw5H/pYxV0OtQnjsvjdGf3IexOlqlNJ2yIy0h4Cn6lmC10A4T+Dz1380BcFXC+1GuLb+Z2W7uHO9ao/0P0ae96MAFhFsRjkKbok1YNpOwq40+5R+VY1DuOs9NaTkzsRN9e5awYWO7OPTcc7E4yra5kfA1lPMA8a42+rEIHQijfO8+/+y9jjDvLggrTQJzEHAGPwcOE2UgwMEuNy/ldPw1KJ8L3Zm8PxN2vdGP39OJzNxyNiKvcl6LvLyJrqsNR5X7CnL54TT6sy96G/E8RjaSewTOYrtbKWN39EHin0AnUMcYwq412p8nQf6O0rf9UP5xvQCV5QZSuqI1k4AOvpk0LaswBDj4PYnORb9E0wvTsCY0hP7cgj6PDkG3N6HIIhSRO8uPpD95pKxquxLZpJoYRph/7r5bXZ8yf0w8X6Vcu5rXpeFKtHsy/Tm82n7iufS0GM5/2WqeoQSqBHTwVRKGEpBAJwk8T+Uv1DUglxBHMut+iu1y+aLWcmo7dXSzhU9OzddbGfpW3yfTTSCgg28CRIuQgATmmkAup+QlMB9JSYRx0F9CI7n8kJsRcyljtlFUTvNvg657Mac7/zJbv5WWP0Z/8uGa2Ub8Z0TmZZnP1FeYGPyPgDfZ/Y+FMQlIoEMEcFBP4qxyvfxqws0IM3O/AuWNe8O1bLMm5UwizE2Wj6LXUke3z+Dpxux7EibRt3OJ5wbCOP1w0yTQh4AOvg8SMyQggU4QwAE/gONahbqXQU8nPZJ2sF1OY7+PsnIXfe6ov5e8+tP/Iym649vQj7zfYUP6tiThC6Qf7HijbEBhCejgCzs0NkwCvUcAh5VZ9vRm9LzMzo++5cyEJoFBCXgNflA8LpSABCQgAQl0JwFn8N05bra6TQQ4FTqeqj6JcqfynegrzJ7yethhGeUswAYfQ+uinGY9lHKaMlOlLE0CEpBAHwLO4PsgMUMCLxLAKS9KbCK6EOXteLehc8l/BWHDxvr5nV2N8k/Cwegn6Ezy1yLUJCABCbSEgA6+JVgttCQEDqEfmbH/Ft2Ffkj6MrT7MPu3I+tfx/Z5Oc2daDLpOPoDh1mOq0tAAhJomIAOvmFUrtiDBF5Jn2+p63deEZuZ+HAs6/+hboMbSad8TQISkEBLCOjgW4LVQotAgFPgeYXnUqj+rWaNNi9OuH62/lXykj8c+zcrv5NmzFez0bbE765JDytKWa9Ciw1ro35WpoxF0bJoru7HYfsFKuXknfKaBCRQAAJz9aMuQPttggT6EMDR5KMcn0NvqCzMl8kmcGr8sT4rD57xPRb/hW2XIzwfvRP9mXJq35dO1uDG+pdTRl7eMo3wU4SroTeid6NhGdvnZr3D0QpoYdJ5Leu7qSMfZhmWse2H2WB7NAOtR/otlDPsfzrY7k1sn37li3Ybkf4Y5ZxHXJOABDpIQAffQfhW3TICJ1ByXt35LhzNLBzOzsSPJ9yHZMOfHmXdvF0td73vgpZGE8l76RWopBs2tvsKZeUVqrmxLnfPxynnbvrh2iQ2mMK2s/85oMxcx/8a+gxq2NhuX1ZeB21FWc+SnkD8ZMJdScfhN2SsP54Vr0Crst3NpBcnfjrhA6RzY6EmAQl0iIAOvkPgrbalBNag9A/jYGalFoLcsR4HFud6bfIaNbbNi1dOb3T9wdajrEtZHo3I6MNKbJib9PavFkD8WPJ/gZYkPpyXn2xNGS/9w8O2V1FG2rYxGs7sO2cAtmf7mwnD+kHKyT8cefe7Dj5QNAl0iIDX4DsE3mpbSmAmjqb+1aS5Dl97DbylDWhR4Wl/f7Pr5A/3n/Xwea6unfmHaLhfb5ufbeovD6Tc4ZZT1xSTEpDA3BLQwc8tQbcvIoEbmEWeVG0Y8fcSPwBdX83rREg7Xoveg/4fGsmNf9No9zxselC1/cS/Tnw1/qG5p5rXYPhH1jutui7lvJ74N1EuIwzH8nTA0Wy/aDYijGO/Cl2UtCYBCXSOwHD/6+9cS61ZAo0TOIxV/4Cz+TXhf1DuNl8eJ5gb0jpitGUPKs4/GhegPBf/ZfJyU9vTxBsy1n2ebfJWvX8SvpYw9xPkH4UN0XDtRDZ4PeVcSZjLFq9Ga1FH7l1o2Fg/NyHmlPzfCX9BuCr6IPkXN1yIK0pAAi0hMIofYksKLlqhHHyOpU1L0d9di9Y229MaAoz56pScGeWNjHv9aeTWVNpPqbRjPbL/ghajHbPv5K/sj/nK2dH9bDJoFtvmKYH0LafZb6KM3CcwIqOsOOSxKOXMbttICqKcJdguTxvcTznDvhN/JHW6jQS6gQC/jW/Rzp/xuxjW/T/N6Jsz+GZQtIxCEuAH9a+CNGwD2rFfnQM9grwfjKR9lBOH/o+RbFu/DWXdVJ83kjTl5NOukSYBCRSEgNfgCzIQNqPUBJ6gdzkFXmvLkBjxzLu2IOMSkIAE+iOgg++PSo/kcepoJXQUyjPiu/dItzvRzdwL8AYY5wa70SgvujkGHd6JxlinBCTQGwQ8Rd8b49ynlziZFcnMXdS5PjQd/Yi8TTjVujdxrYkEYPoEbHPvR07J74seQV8lv+3X5KhXk4AEeoRA4Rw8B8K0aSwHvxwEtdYROJ6iD4bzlZUq1oH9Keit5F3Sump7s2SY5ga23EWvSUACEmgLgUKcosep5LTlEeguep27nR8mPhNNRXu2hUTvVZKxv66u23nzWO6G1iQgAQlIoMsJFMLBwzCzyTzXuzVaGKVdy6J90P44+Q8Ras0lcDvFbVEtEsZ5njqc76zmGUpAAhKQQPcSKIqD3xKEeYzoBvQ4is1AeSPWx9F23Yu4sC3PG9DOwa9/GL2N+NnojxXmhW10tzYMxkujH6C8gOd69M5u7Ut9u+nLmuin6Pfob2j9+nVMS0AC7SdQFAc/la5vNkD389EKn68dAM5Is3Hkmam/EuUlJ5uin5B3EKHWZAI4vAUpMm/UuxTG+Wdqc7QP+S+dQSHdlUYflqLheSb/+/TtHYS5z+Ab5L+OUJOABDpIoBBvsuNgkDd9TUS5EWkaehQtgtZAuekun7S8g3BQo5y8SSun9vuzT5A5P+Xs2N9C8yTQKgLsl7tR9nj2va9W6yBvVeKHkJdX2Hat0Y98pjZv5Dut2gnyckZuM/I+V80zlECvEuD30NtvsuNAMAUIcfIT0Di0NMqsPe/LvgLFyTdia7HSpgOs+BryvTN/ADhmt5RAXpd7d10NM0jnfpNut5yduLeuE/mdvaIuz6QEJNBmAo06zpY2C+e+PBUcjnZAue7+IZz+LYT5OlVO+SV/56QHM7a5gOVRH6OcY8nM6UStBwgw3tm384rY+dB17Bs5O9Qpu5KK856B82nHfZVGfJRwSqca1MR6L6esQ+nbZPr2FOHLSV+E9m5iHRYlAQmMgEAhHDztPhDlGmUOyLuiyRwoNuWAcRNxTQLDIsC+k1nl91Bmks+i7E85RX4b8bYb9ebrb0dR8W2ERxKuhHIZKvt9Vxt9+yN9WodOpI8nE+ZM3DHkn93VHbPxEigBgaI4+K1guR4HhScJMxu4kfBCwjeXgLFdaD+BzCBz09fpqZr96EKCYwl3J68jM3nq/Q31r0Y71kR5WiGz+lIYfTmGvp1PZ3IPzO9IX1+KjtkJCXQ5gaI4+Dj0zN5zuu9lHCB+zgFjWaI53X5S8jQJDIPAU+xDs517tiGeWebbia6N/pS8Ro3tcl/IQWhR9DTak/Kqp9lJNm5sdxdrR6Uz+vZvOhVpEpBAQQjkelkR7Ps04pccTA+uNoYDxjHEc5ov1841CQyHwCz2pXnqNsh9HnHQDRtljGPlr6OvoLeh7JNnkL84oSYBCUig0AQK4eBx5jmlujL6XS0t8g8jnUduLqzNN14sAji8edFu6PMoz3fXO9d2N/g8KvwW7ZgPxb5JOqfG/zbMhnyS9T/HfpgXMMX+QPqXKDd+ahKQgAQKTaAQDj6EOHjORH+vp0XeJPSj+nzTxSAQ70lLLka5ueoatD66kewFCDtl36Hi3MR2Bco/j7PQBPajhMOxsaz8YN0GeSTMR8DqoJiUgASKR6AwDr54aGxRgwTyopZH8J2fQhej/Un/DO3b4PZNX402xA5FG6It0KfREyOoKNfrD6luxz8tY4hPRKW5Qa7aN0MJSKB8BHTw5RvTdvco17a/XVfpuaRXqMvrxuQpNDqn+S9DBxL/OTqIfxYu78bO2GYJSKC3CBTlLvreol6u3v6X7myELqvp1tbEH65Jd2UUR/4CDX8Pzn1jwsVQzlBM7crO2GgJSKDnCOjge27Im97hzHL/hRNckjDxOPsvojjEUlgzZuzwWQQYn0PrVqB8nnKvLQUgOyEBCRSSgKfoCzks3dMonFSubY9H09FeKN8ReDX5eWmRBgGce16XezPK8/M7oE+ifHFtAqEmAQlIoCUEnMG3BGtvFYozf44ef6u3ej2s3u7M2mfAqfpOhzxlkJfnxNFfNaySXFkCEpBAgwScwTcIytUkMBcEFmLbKXXb30M6j+FpEpCABFpCQAffEqwWKoE5CPyV1C7M2uPoq5ZvxOe0vSYBCUigJQQ8Rd8SrBYqgf8R4NT833DuZ5NzF+HnCVdFS6A9kSYBCUigJQScwbcEq4VKYE4COPmTydkEPYouQe8n71lCTQISkEBLCDiDbwlWC5VAXwI49BvIjTQJSEACLSfgDL7liK1AAhKQgAQk0H4CzuCHYM410zDaES2KbmEWli+KaS0iAO+1KfpNKM/RnwXvx1tUlcVKQAISKDUBZ/CDDC/OJnzy4ZTXo0fQMeR9h1BrAQHY7kSx30Nx6suix8hLqElAAhKQwDAJOIMfHNgnWDyFWeRRWQ1ncybBGYRbk3d+8rTmEIDpSpSUV92Og+1DKZW8aQRfRvsiTQISkIAEhkHAGfzgsF7L4l9VV8Hx5HvimdG/rppn2DQCq1HSUVXnnlKJ/4Jg5cQ1CUhAAhIYHgEd/OC8ZrB49bpVNiKdfK25BPL42BzOnBn8MuQt1txqLE0CEpBAbxDwFP3g4/xdFp+Oo8k14RtRrhHvj5ZCWhMJMFu/Es47oyMp9gSUmxq/gT6LNAlIQAISGCYBHfwgwHA6t+JwdmCVOJ35Ua4J5xrxM4Ra8wkcSJEHo3yU5TF0JKxrvzNPliYBCUhAAo0Q0MEPQQkHcy+r+ErRITg1YzGsc4/D7Bsam1GeZUhAAhLoZQJeg+/l0bfvEpCABCRQWgI6+NIOrR2TgAQkIIFeJqCD7+XRt+8SkIAEJFBaAjr40g6tHZOABCQggV4moIPv5dG37xKQgAQkUFoCOvjSDq0dk4AEJCCBXibgY3I9PPo845/xfztaEF3HU2o39zAOuy4BCUigVAScwZdqOBvvDM49L+75CdocvRr9k7x3EWoSkIAEJFACAs7gSzCII+zC8Wz3J2btJ2Z7nPtEgnwp71/k/Tt5mgQkIAEJdC8BZ/DdO3Zz2/J8nvW0aiE49fuJ58t561fzDCUgAQlIoHsJ6OC7d+zmtuX5elv9R3PyedyZc1uw20tAAhKQQOcJ6OCHGANOWS+BvoSORR9GZWF2Ml0/nv6siOZH+Wrbtuj8IZC4WAISkIAEuoBAWZxVS1Dj9Bag4F+je9Bp6M3oavLnI+xq45T8hXQgX207FcWp5/Osa5L/HKEmAQlIQAJdTsCb7AYfwC+w+Ls4vZ9XVtsV555Px74HnV7J69qAfv2BxkeaBCQgAQmUjIAz+MEHdHEWX1u3yhWkl6zLMykBCUhAAhIoFAEd/ODDMZ3FW9Stsj/pu+ryTEpAAhKQgAQKRaBwp+g5BZ42jeX08SMFIJVnxe+gTcsQXoR2RqNp2y8JNQlIQAISkEBhCRRiBo8DHY2OQJkZP4MeJj4TTUV7dooejnwGded0/IPorehqtBXSJCABCUhAAoUmUJQZfGbKS6Ot0a0oz2IvjNZEx+Hkx+BsZ79xjfSAxnrLsTCz7f4sz3zn9azDMurNPxzHDWsjV5aABCQgAQl0mEBRHPyWcJiAM723hkdmz1fhtD9OeBga0sGzzlpoM9SfxfFP72+BeRKQgAQkIIGyESiKg58K2DjmM/oBvA15D/ST3yeLfxAuIDPqY/yj8B4y86y3JgEJSEACEig9gaI4+EMhPREnfCDhNJTXqC6C1kBpo9e9gaBJQAISkIAEGiVQCAfPzHsKzn09Gj0BjUO5Hp9Ze07LT2b5LEJNAhKQgAQkIIEGCRTCwaet+PCnCC5tsN2uJgEJSEACEpDAIAQK8ZjcIO1zkQQkIAEJSEACIyAwqlfOfnMJYF345KMqU0bAqZmbbERhueTgZYdmUu2/rJyhmgc93f9ic5tM4BWUl0dctdYTGE0VOYY82/qqrAEC4T15hCTGs90W+Nq7R7j9iDfrGQc/YkJN3pB/NM6hyPcz2HkMUGshAVjn5URvgvVXW1iNRVcIwHsSrDcVSOsJwHo/apkB7+qHsFpfaQ/X0K37tqfoe3intesSkIAEJFBeAjr48o6tPZOABCQggR4moIPv4cG36xKQgAQkUF4COvjyjq09k4AEJCCBHiagg+/hwbfrEpCABCRQXgI6+PKOrT2TgAQkIIEeJuBjcm0efB63WIoqH+DxlhfaXHXPVQfrMXQ6nxr+b891vgMdhvcysP5PB6ruuSphPZZOvwBv3zvQhtHv1n1bB9+GncMqJCABCUhAAu0m4Cn6dhO3PglIQAISkEAbCOjg2wDZKiQgAQlIQALtJqCDbzdx65OABCQgAQm0gYAOvg2QrUICEpCABCTQbgI6+HYTtz4JSEACEpBAGwjo4NsA2SokIAEJSEAC7Sagg283ceuTgAQkIAEJtIGADr4NkK2ivQR4KcV87a2xd2uD9bxoVO8SsOcSKC4BHXybxoZj4Abozjq9uk3V90w18N2Fzl5V22HyNkVXoNvQr9GitcuNj4wAHJdnyzvQ+GoJ5LmfV2E0KYTpLugSdD06Ha1RLZr459ANKPv256r5hiMjAMP50DfQXys6knB0SiPsun1bBz+y/WAkW23ARhej1Wt0z0gKcpu+BPjxLYq+y5Jvo5dmlOQtTnoi+jBaFd2GvoW0uSAA173Y/FK0RF0x7ud1QOYmCeel2f5YtAuvpV2H8I/oOBSHsxPB1mhjNAG9h7z/R6iNnMD72XRlFJ7RmmgPFOu6fVsH/+LAtePvulRyNVoSLcGP9Qk0qx0V90gdm9PPJ1B+oLWWH+U/QX0Depb48WiH2hWMD48ATiQzmp3RVqj+Pf/u58PDOdTaOUbvzL57X2XF6wk3qsTfQXg6y2age4mfgbavLDMYGYHw/XSOFRHxG9GbKkV13b6tg6+MXBuC7ByfQheh2zlIHt2GOnumCn6MZ6HP0OEn6zq9AunaD6DkQLkI/OevW89kgwTg/Ax6O7qpn03cz/uBMtIsGN+DJtdsvy/x8yvp+n07Tj4fs9JGSADW16Bp2ZxjxCsIdkXnJY113b6tg39x4Nrx929Ushc7T04T/x/6KDtQ/enNdrSj1+p4FR2eWdPp6j8AC9bkGW0eAffz5rGcoySOF3uTsS3KRCFWv2/nDFackjaXBGCds1Q/R3H4Z1eK67p9Wwc/lztCo5uzk3wEXZ71CacQ/Al5qjhAWmsPUvzCNVWMJf4UY/BITZ7RJhFwP28SyLpicDj7kfU19DYYT68srt+3s5/fU1lmMEICFef+KzafB2UGP9u6cd/WwVdHr4UhO8wY9KWENdVkBvlATdpoawjkYDiupujE76pJG20SAffzJoGsKwauua/kyyjO/Z81i7Nvr1iTHkfcfbsGyHCjsJ6XbTJzj3PfAd7PpIxu3bd18Bm9Fhs7yVNU8Va0V6piZ9mQYD10UdJaSwlcQunjYb45ynX3g1D1lFtLK+61wt3Pmz/i7LMrUeoJ6L3oHtKLRZWaziT8AOll0TjiWefXSBs5gQPYdGX0QbRghfdC3bpv6+BHviMMd8s8o7orO8xNhL9HuR7/+HALcf3hEYDx02yRH+05KDfPLI8OR1prCLifN5frhyku19UnoYeq4jiSM4AXolwX/ge6Cv2M/f2vhNrICXyCTddGudRR5f2LSnFdt2+PYoeotN2gHQTyHyH1/BfuL7SjPut4kQDcc+ptLNy99t6GncL9vA2QK1XAOtfen2bfzj+zWosJdNO+rYNv8c5g8RKQgAQkIIFOEPAUfSeoW6cEJCABCUigxQR08C0GbPESkIAEJCCBThDQwXeCunVKQAISkIAEWkxAB99iwBYvAQlIQAIS6AQBHXwnqFunBCQgAQlIoMUEdPAtBmzxEpCABCQggU4Q0MF3grp1SkACEpCABFpMQAffYsAWLwEJSEACEugEAR18J6hbpwQkIAEJSKDFBHTwLQZs8RKQgAQkIIFOENDBd4K6dUpAAhKQgARaTEAH32LAFi8BCUhAAhLoBAEdfCeoW6cEJCABCUigxQR08C0GbPESkIAEJCCBThDQwXeCunVKQAISkIAEWkxAB99iwBYvAQlIQAIS6AQBHXwnqFunBCQgAQlIoMUEdPAtBmzxEpCABCQggU4Q0MF3grp1SqBHCIwaNWpjdG2PdNduSqBQBHTwhRoOGyMBCUhAAhJoDgEdfHM4WooECk2AWfS86BT0X3QHOjgNJlwLnYF+gu5H56IlK8uIjjoETUd3oy+gUZVlSxD9FUp516O3JD9GfAf0d3Q7ye2Tp0lAAu0noINvP3NrlEAnCBFKjJgAAANPSURBVOxIpa9BK6N3oDjrpMeg96L70OroXvQDFNsdvQ9ti7ZDu6A3oNgpaAbKNsehH6M49/EEJ6Mvomy3OdIkIIEOENDBdwC6VUqgAwSepc7l0UboVrTErFmzbiGMPYe+Rvphwm+grXHU8xO+H8VxT0P/RnHq27JsMcKt0NfRE+hsdA/5axNugf5BWeegvxOf7fgJNQlIoM0E5m1zfVYnAQl0hsA5VLs++hEajU7DIX+aMBaH/FgihLeQH6e9Ino1yjqfQFXLDXPLoVnokmpmJcw/DzlD8Lea/KuJ71GTNioBCbSJgA6+TaCtRgIdJpAZ+bfQl9CW6AT0DxSHvRJOfR6c+/OEK5BeEN2FrkGTyZ99yp5lC5GeBz2Jcnp+LZY9SJhT80tU8vYj3Dp5FVupGjGUgATaS8BT9O3lbW0S6BSBXGc/E2XmfQHKKfeqLUxkO5x0bqD7APoTjjtO/Fy0J9mLVpadTvpAlj1D+Ef0EfJfjpYmfiPK9fjL0BvJWxWNIb4z0iQggQ4Q0MF3ALpVSqADBH5KnTNRrqffiV5AE1EsN9hlZp/8nVBurIv9DuWmu9vRzSiz96NR7Ei0G7oNXYm+ieO/ISJ+SCXvFsKnkCYBCXSAwCh+kB2o1iolIIFOEKjMqkfzu3809ZN+PcFPSa9OfDHCh+vbRf4rksey/IMwh7FsCTIeZNkcBxLy5yN/DNmzr+3PsZEJCUigLQS8Bt8WzFYigWIQwOFmRt3vrJplfZx7Wk1+H8de7Q3LHqjGa0PynyUdaRKQQIcIeIq+Q+CtVgIFIZBT7F8oSFtshgQk0EQCnqJvIkyLkoAEJCABCRSFgDP4ooyE7ZCABCQgAQk0kYAOvokwLUoCEpCABCRQFAI6+KKMhO2QgAQkIAEJNJGADr6JMC1KAhKQgAQkUBQCOviijITtkIAEJCABCTSRgA6+iTAtSgISkIAEJFAUAjr4ooyE7ZCABCQgAQk0kYAOvokwLUoCEpCABCRQFAI6+KKMhO2QgAQkIAEJNJGADr6JMC1KAhKQgAQkUBQCOviijITtkIAEJCABCTSRgA6+iTAtSgISkIAEJFAUAjr4ooyE7ZCABCQgAQk0kYAOvokwLUoCEpCABCRQFAL/HyByKW1tHRa9AAAAAElFTkSuQmCC" alt="plot of chunk unnamed-chunk-2"/> </p>

</body>

</html>

test.htmlTitle========================================================

This is an R Markdown document. Markdown is a simple formatting syntax for authoring web pages (click the **MD** toolbar button for help on Markdown).

When you click the **Knit HTML** button a web page will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

```rsummary(cars)```

```## speed dist ## Min. : 4.0 Min. : 2 ## 1st Qu.:12.0 1st Qu.: 26 ## Median :15.0 Median : 36 ## Mean :15.4 Mean : 43 ## 3rd Qu.:19.0 3rd Qu.: 56 ## Max. :25.0 Max. :120```

You can also embed plots, for example:

```rplot(cars)```

![plot of chunk unnamed-chunk-2](figure/unnamed-chunk-2.png)

test.md

Wednesday, 11 September, 13

Page 25: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

Title========================================================

This is an R Markdown document. Markdown is a simple formatting syntax for authoring web pages (click the **MD** toolbar button for help on Markdown).

When you click the **Knit HTML** button a web page will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

```{r}summary(cars)```

You can also embed plots, for example:

```{r fig.width=7, fig.height=6}plot(cars)```

test.rmd

<!DOCTYPE html><html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"/>

<title>Title</title>

<!-- 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>

<!-- MathJax scripts --><script type="text/javascript" src="https://c328740.ssl.cf1.rackcdn.com/mathjax/2.0-latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>

<style type="text/css">body { font-family: Helvetica, arial, sans-serif; font-size: 14px; line-height: 1.6; padding-top: 10px; padding-bottom: 10px; background-color: white; padding: 30px; }

body > *:first-child { margin-top: 0 !important; }body > *:last-child { margin-bottom: 0 !important; }

a { color: #4183C4; }a.absent { color: #cc0000; }a.anchor { display: block; padding-left: 30px; margin-left: -30px; cursor: pointer; position: absolute; top: 0; left: 0; bottom: 0; }

h1, h2, h3, h4, h5, h6 { margin: 20px 0 10px; padding: 0; font-weight: bold; -webkit-font-smoothing: antialiased; cursor: text; position: relative; }

h1:hover a.anchor, h2:hover a.anchor, h3:hover a.anchor, h4:hover a.anchor, h5:hover a.anchor, h6:hover a.anchor { text-decoration: none; }

h1 tt, h1 code { font-size: inherit; }

h2 tt, h2 code { font-size: inherit; }

h3 tt, h3 code { font-size: inherit; }

h4 tt, h4 code { font-size: inherit; }

h5 tt, h5 code { font-size: inherit; }

h6 tt, h6 code { font-size: inherit; }

h1 { font-size: 28px; color: black; }

h2 { font-size: 24px; border-bottom: 1px solid #cccccc; color: black; }

h3 { font-size: 18px; }

h4 { font-size: 16px; }

h5 { font-size: 14px; }

h6 { color: #777777; font-size: 14px; }

p, blockquote, ul, ol, dl, li, table, pre { margin: 15px 0; }

hr { border: 0 none; color: #cccccc; height: 4px; padding: 0;}

body > h2:first-child { margin-top: 0; padding-top: 0; }body > h1:first-child { margin-top: 0; padding-top: 0; }body > h1:first-child + h2 { margin-top: 0; padding-top: 0; }body > h3:first-child, body > h4:first-child, body > h5:first-child, body > h6:first-child { margin-top: 0; padding-top: 0; }

a:first-child h1, a:first-child h2, a:first-child h3, a:first-child h4, a:first-child h5, a:first-child h6 { margin-top: 0; padding-top: 0; }

h1 p, h2 p, h3 p, h4 p, h5 p, h6 p { margin-top: 0; }

li p.first { display: inline-block; }li { margin: 0; }ul, ol { padding-left: 30px; }

ul :first-child, ol :first-child { margin-top: 0; }

dl { padding: 0; }dl dt { font-size: 14px; font-weight: bold; font-style: italic; padding: 0; margin: 15px 0 5px; }dl dt:first-child { padding: 0; }dl dt > :first-child { margin-top: 0; }dl dt > :last-child { margin-bottom: 0; }dl dd { margin: 0 0 15px; padding: 0 15px; }dl dd > :first-child { margin-top: 0; }dl dd > :last-child { margin-bottom: 0; }

blockquote { border-left: 4px solid #dddddd; padding: 0 15px; color: #777777; }blockquote > :first-child { margin-top: 0; }blockquote > :last-child { margin-bottom: 0; }

table { padding: 0;border-collapse: collapse; }table tr { border-top: 1px solid #cccccc; background-color: white; margin: 0; padding: 0; }table tr:nth-child(2n) { background-color: #f8f8f8; }table tr th { font-weight: bold; border: 1px solid #cccccc; margin: 0; padding: 6px 13px; }table tr td { border: 1px solid #cccccc; margin: 0; padding: 6px 13px; }table tr th :first-child, table tr td :first-child { margin-top: 0; }table tr th :last-child, table tr td :last-child { margin-bottom: 0; }

img { max-width: 100%; }

span.frame { display: block; overflow: hidden; }span.frame > span { border: 1px solid #dddddd; display: block; float: left; overflow: hidden; margin: 13px 0 0; padding: 7px; width: auto; }span.frame span img { display: block; float: left; }span.frame span span { clear: both; color: #333333; display: block; padding: 5px 0 0; }span.align-center { display: block; overflow: hidden; clear: both; }span.align-center > span { display: block; overflow: hidden; margin: 13px auto 0; text-align: center; }span.align-center span img { margin: 0 auto; text-align: center; }span.align-right { display: block; overflow: hidden; clear: both; }span.align-right > span { display: block; overflow: hidden; margin: 13px 0 0; text-align: right; }span.align-right span img { margin: 0; text-align: right; }span.float-left { display: block; margin-right: 13px; overflow: hidden; float: left; }span.float-left span { margin: 13px 0 0; }span.float-right { display: block; margin-left: 13px; overflow: hidden; float: right; }span.float-right > span { display: block; overflow: hidden; margin: 13px auto 0; text-align: right; }

code, tt { margin: 0 2px; padding: 0 5px; white-space: nowrap; border: 1px solid #eaeaea; background-color: #f8f8f8; border-radius: 3px; }

pre code { margin: 0; padding: 0; white-space: pre; border: none; background: transparent; }

.highlight pre { background-color: #f8f8f8; border: 1px solid #cccccc; font-size: 13px; line-height: 19px; overflow: auto; padding: 6px 10px; border-radius: 3px; }

pre { background-color: #f8f8f8; border: 1px solid #cccccc; font-size: 13px; line-height: 19px; overflow: auto; padding: 6px 10px; border-radius: 3px; }pre code, pre tt { background-color: transparent; border: none; }

sup { font-size: 0.83em; vertical-align: super; line-height: 0;}* {! -webkit-print-color-adjust: exact;}@media screen and (min-width: 914px) { body { width: 854px; margin:0 auto; }}@media print {! table, pre {! ! page-break-inside: avoid;! }! pre {! ! word-wrap: break-word;! }}</style>

</head>

<body><h1>Title</h1>

<p>This is an R Markdown document. Markdown is a simple formatting syntax for authoring web pages (click the <strong>MD</strong> toolbar button for help on Markdown).</p>

<p>When you click the <strong>Knit HTML</strong> button a web page will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:</p>

<pre><code class="r">summary(cars)</code></pre>

<pre><code>## speed dist ## Min. : 4.0 Min. : 2 ## 1st Qu.:12.0 1st Qu.: 26 ## Median :15.0 Median : 36 ## Mean :15.4 Mean : 43 ## 3rd Qu.:19.0 3rd Qu.: 56 ## Max. :25.0 Max. :120</code></pre>

<p>You can also embed plots, for example:</p>

<pre><code class="r">plot(cars)</code></pre>

<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAGwCAYAAABFI3d+AAAEJGlDQ1BJQ0MgUHJvZmlsZQAAOBGFVd9v21QUPolvUqQWPyBYR4eKxa9VU1u5GxqtxgZJk6XtShal6dgqJOQ6N4mpGwfb6baqT3uBNwb8AUDZAw9IPCENBmJ72fbAtElThyqqSUh76MQPISbtBVXhu3ZiJ1PEXPX6yznfOec7517bRD1fabWaGVWIlquunc8klZOnFpSeTYrSs9RLA9Sr6U4tkcvNEi7BFffO6+EdigjL7ZHu/k72I796i9zRiSJPwG4VHX0Z+AxRzNRrtksUvwf7+Gm3BtzzHPDTNgQCqwKXfZwSeNHHJz1OIT8JjtAq6xWtCLwGPLzYZi+3YV8DGMiT4VVuG7oiZpGzrZJhcs/hL49xtzH/Dy6bdfTsXYNY+5yluWO4D4neK/ZUvok/17X0HPBLsF+vuUlhfwX4j/rSfAJ4H1H0qZJ9dN7nR19frRTeBt4Fe9FwpwtN+2p1MXscGLHR9SXrmMgjONd1ZxKzpBeA71b4tNhj6JGoyFNp4GHgwUp9qplfmnFW5oTdy7NamcwCI49kv6fN5IAHgD+0rbyoBc3SOjczohbyS1drbq6pQdqumllRC/0ymTtej8gpbbuVwpQfyw66dqEZyxZKxtHpJn+tZnpnEdrYBbueF9qQn93S7HQGGHnYP7w6L+YGHNtd1FJitqPAR+hERCNOFi1i1alKO6RQnjKUxL1GNjwlMsiEhcPLYTEiT9ISbN15OY/jx4SMshe9LaJRpTvHr3C/ybFYP1PZAfwfYrPsMBtnE6SwN9ib7AhLwTrBDgUKcm06FSrTfSj187xPdVQWOk5Q8vxAfSiIUc7Z7xr6zY/+hpqwSyv0I0/QMTRb7RMgBxNodTfSPqdraz/sDjzKBrv4zu2+a2t0/HHzjd2Lbcc2sG7GtsL42K+xLfxtUgI7YHqKlqHK8HbCCXgjHT1cAdMlDetv4FnQ2lLasaOl6vmB0CMmwT/IPszSueHQqv6i/qluqF+oF9TfO2qEGTumJH0qfSv9KH0nfS/9TIp0Wboi/SRdlb6RLgU5u++9nyXYe69fYRPdil1o1WufNSdTTsp75BfllPy8/LI8G7AUuV8ek6fkvfDsCfbNDP0dvRh0CrNqTbV7LfEEGDQPJQadBtfGVMWEq3QWWdufk6ZSNsjG2PQjp3ZcnOWWing6noonSInvi0/Ex+IzAreevPhe+CawpgP1/pMTMDo64G0sTCXIM+KdOnFWRfQKdJvQzV1+Bt8OokmrdtY2yhVX2a+qrykJfMq4Ml3VR4cVzTQVz+UoNne4vcKLoyS+gyKO6EHe+75Fdt0Mbe5bRIf/wjvrVmhbqBN97RD1vxrahvBOfOYzoosH9bq94uejSOQGkVM6sN/7HelL4t10t9F4gPdVzydEOx83Gv+uNxo7XyL/FtFl8z9ZAHF4bBsrEwAAQABJREFUeAHtnQn8ZXP9/5swhgyRPcsYWcv2ozRKSNTfUpYoRMlaaZFKi5TK1oKSRKVIo0RRJBRjiKQMmlQY65Bdg7Ez/+dr3Ks797vd73fucu65z/fj8fp+lnPOZ3l+zj3v7+eso2bNmvUyTQISkIAEJCCBchF4ebm6Y28kIAEJSEACEggBHbz7gQQkIAEJSKCEBHTwJRxUuyQBCUhAAhLQwbsPSEACEpCABEpIQAdfwkG1SxKQgAQkIAEdvPuABCQgAQlIoIQEdPAlHFS7JAEJSEACEtDBuw9IQAISkIAESkhAB1/CQbVLEpCABCQgAR28+4AEJCABCUighAR08CUcVLskAQlIQAIS0MG7D0hAAhKQgARKSEAHX8JBtUsSkIAEJCABHbz7gAQkIAEJSKCEBHTwJRxUuyQBCUhAAhLQwbsPSEACEpCABEpIQAdfwkG1SxKQgAQkIAEdvPuABCQgAQlIoIQEdPAlHFS7JAEJSEACEtDBuw9IQAISkIAESkhAB1/CQbVLEpCABCQgAR28+4AEJCABCUighAR08CUcVLskAQlIQAIS0MG7D0hAAhKQgARKSEAHX8JBtUsSkIAEJCABHbz7gAQkIAEJSKCEBHTwJRxUuyQBCUhAAhLQwbsPSEACEpCABEpIQAdfwkG1SxKQgAQkIAEdvPuABCQgAQlIoIQEdPAlHFS7JAEJSEACEtDBuw9IQAISkIAESkhAB1/CQbVLEpCABCQgAR28+4AEJCABCUighAR08CUcVLskAQlIQAIS0MG7D0hAAhKQgARKSEAHX8JBtUsSkIAEJCABHbz7gAQkIAEJSKCEBHTwJRxUuyQBCUhAAhLQwbsPSEACEpCABEpIQAdfwkG1SxKQgAQkIAEdvPuABCQgAQlIoIQEdPAlHFS7JAEJSEACEtDBuw9IQAISkIAESkhAB1/CQbVLEpCABCQgAR28+4AEJCABCUighAR08CUcVLskAQlIQAIS0MG7D0hAAhKQgARKSEAHX8JBtUsSkIAEJCABHbz7gAQkIAEJSKCEBHTwJRxUuyQBCUhAAhLQwbsPSEACEpCABEpIQAdfwkG1SxKQgAQkIAEdvPuABCQgAQlIoIQEdPAlHFS7JAEJSEACEtDBuw9IQAISkIAESkhAB1/CQbVLEpCABCQggXl7CcGoUaPeTX97qs+9NL72VQISkEABCdw/a9asSzrRrlFU3Il6214nzn1HKj0Indr2yq1QAhKQgAR6lcDH6Phu+Nrr2g2gcLNZHHHaNBYYjzQZRso9jXJPanK5FicBCUhAAhLolwA+bVUWdORyeEcqracAgNHoCHQXy55BDxOfiaaiPevXNy0BCUhAAhKQwOAEijKDP55mLo22RreimWhhtCY6Dic/hpn3icQ1CUhAAhKQgAQaIFCIGTzt3BLthxO/AT2OYjPQVeR/HG2HNAlIQAISkIAEGiRQFAc/lfZuNkCbtyH/gQGWmS0BCUhAAhKQQD8EinKK/lDaNpFT8QcSTkOPokXQGiht3AppEpCABCQgAQk0SKAQDp5T8VNw7uvR5gloHMr1+Mzac919cs7XE2oSkIAEJCABCTRIoBAOPm3Fhz9FcGl9u3H886B5Wf50/TLTEpCABCQggVoC+ItRpD+KtkXzoefQNhUfQ7R3rBDX4BmP5dFp6HF0MXpNzRDsRPynNWmjEpCABCQggYEIHM6CldBWOPVNCc9E38GvzEPYU1aUGXyuvf8HbYB2RZMZjE0ZnJuIN2xsM46Vlx9gg/8jP8/YaxKQgAQkUF4Cm9O1N+E/MnPP2eGT8Q3rEF0X/S15vWJFcfC5iW49BuJJwkMZjBsJLyR88zAHYhXW33iAbd5Cfm7e0yQgAQlIoLwEHqo695ou5hLwmJp0T0SL4uDj0DN7vzzUGZyf49yXJXoBavjVsmx3MetHfYzyxpK5VJ8FZkhAAhKQQJkI3MTx/iv4gzyd9TLimdF/En0Z9ZQVxcF/H+q/ZCCOZVCOzggQHlNxyseSPCd5mgQkIAEJSGAIAl9k+bX4jxUI70WvRePxKY8R9pQVwsED/iIGY2XIj6+lT/5h5F9GXpZpEpCABCQggUEJxJHjN1ZnpZwVHo0O70XnHkiFcPBpCAOQ98//PfFaI38S6UiTgAQkIAEJDEkAv/E8K1095IolX6EQj8mVnLHdk4AEJCABCbSdQCFm8JxOOYie54UEA9m/+I/M6/AD0TFfAhKQgAQkUEegEA6eNo1DB6BTUU7V15sfm6knYloCEpCABCQwCIFCOHhm5x9lFp/LBS8n/pFB2usiCUhAAhKQgAQaIFCka/AH096FcfQLNdBuV5GABCQgAQlIYBAChZjBp33M3B8n2G2QtrpIAhKQgAQkIIEGCRRpBt9gk11NAhKQgAQkIIGhCOjghyLkcglIQAISkEAXEtDBd+Gg2WQJSEACEpDAUAR08EMRcrkEJCABCUigCwno4Ltw0GyyBCQgAQlIYCgCOvihCLlcAhKQgAQk0IUEdPBdOGg2WQISkIAEJDAUAR38UIRcLgEJSEACEuhCAjr4Lhw0mywBCUhAAhIYioAOfihCLpeABCQgAQl0IQEdfBcOmk2WgAQkIAEJDEVABz8UIZdLQAISkIAEupCADr4LB80mS0ACEpCABIYioIMfipDLJSABCUhAAl1IQAffhYNmkyUgAQlIQAJDEdDBD0XI5RKQgAQkIIEuJKCD78JBs8kSkIAEJCCBoQjo4Ici5HIJSEACEpBAFxLQwXfhoNlkCUhAAhKQwFAEdPBDEXK5BCQgAQlIoAsJ6OC7cNBssgQkIAEJSGAoAjr4oQi5XAISkIAEJNCFBHTwXThoNlkCEpCABCQwFAEd/FCEXC4BCUhAAhLoQgI6+C4cNJssAQlIQAISGIqADn4oQi6XgAQkIAEJdCGBwjn4UaNGzYsW7UKWNlkCEpCABCRQGAKFcPA49NHoCHQXZJ5BDxOfiaaiPQtDy4ZIQAISkIAEuoTAvAVp5/G0Y2m0NboVzUQLozXRcTj5MbNmzTqRuCYBCUhAAhKQQAMECjGDp51bov1w4jegx1FsBrqK/I+j7ZAmAQlIQAISkECDBIri4KfS3s0GaPM25D8wwDKzJSABCUhAAhLoh0BRTtEfStsmcir+QMJp6FG0CFoDpY1bIU0CEpCABCQggQYJFMLBcyp+Cs59Pdo8AY1DuR6fWXuuu0/O+XpCTQISkIAEJCCBBgkUwsGnrfjwpwguxdGnTWNJP5J8TQISkIAEJCCB4RMoxDV4nLqPyQ1/7NxCAhKQgAQkMCCBoszgm/KYHP8ojKOnyw3Q22XJHz3AMrMlIAEJSEACpSJQFAefx+QmcFr+3hq6M4hfhdPOY3KHoUaeg1+F9d6C+rMVyUyZmgQkIAEJSKD0BIri4KuPyZ3RD/GGH5PjH4SL2T7qY/yjsBCZS/VZYIYEJCABCUighASK4uB9TK6EO5ddkoAEJCCBzhEohINn5u1jcp3bB6xZAhKQgARKSKAQDj5ccfKzH5MrIWO7JAEJSEACEmg7gUI8Jtf2XluhBCQgAQlIoOQECjGD5wa4g+A83yCs/8UM/5xBlrtIAhKQgAQkIIEaAoVw8LRnHDoAnYryqdh682Mz9URMS0ACEpCABAYhUAgHz+z8o8zic7ng5cQ/Mkh7XSQBCUhAAhKQQAMEinQN/mDau3DlefUGmu4qEpCABCQgAQkMRKAQM/g0jpn74wS7DdRQ8yUgAQlIQAISaJxAkWbwjbfaNSUgAQlIQAISGJSADn5QPC6UgAQkIAEJdCeBwpyi7058tloCEpCABIpGgHu5lqdN26N8QXQSl4D/WrQ2tqM9zuDbQdk6JCABCUigLQRw7qtTUR65fghNQ9eQtx1hz5kz+J4bcjssAQlIoNQEzqR3uWH7WTQGLYl+gJP/KzP56cR7xnTwPTPUdlQCEpBATxC4nV7ug8aiOPn3oe+jcainHLyn6BlxTQISkIAESkPgjfTkbmbre6J9icfZ74X6e0sq2eU1Z/DlHVt7JgEJSKAXCdxDp9/KKflJhE+gOPd/otxw11Omg++p4bazEpCABEpP4G56+Hm0J1oAnYTehp5DPWWeou+p4bazEpCABEpP4Ap6uCM6kFP0+xGugt6BrkM9ZTr4nhpuOysBCUig9AS+Tg+XQ5dxmv7sSvx1OPvnS9/zug56ir4OiEkJSEACEuheAhVH/kGc+2L0Ij7uAfJmdW+PRt5yHfzI2bmlBCQgAQkUlAA+/eGCNq1tzfIUfdtQW5EEJCABCUigfQR08O1jbU0SkIAEJCCBthHQwbcNtRVJQAISkIAE2kdAB98+1tYkAQlIQAISaBsBHXzbUFuRBCQgAQlIoH0EdPDtY21NEpCABCQggbYR0MG3DbUVSUACEpCABNpHQAffPtbWJAEJSEACEmgbAV900zbUViQBCUhAAt1EgLfhrUB7P4oWRbeioytvyiNafHMGX/wxsoUSkIAEJNBmAjj3RajyTPQXdASagX5Dfr5Q1xXmDL4rhslGSkACnSDAwXwU9W6JMoObxuztmk60wzo7QuCL1HoUY35OpfYT2B1eRXw39MNKXqEDZ/CFHh4bJwEJdJhAviW+Pcps7pcc4A/pcHusvn0EFqaqf9RVdy3pV9blFTapgy/s0NgwCUigkwRw5vtT/4PM4PZHcfSrognkb97Jdll32whMo6b31NX2BdK31OUVNukp+sIOjQ2TgAQ6TGBD6j+q2gac/DM49x+Q3gD9sZpvWFoCx9OzqYx5vi3/O7QduqXmlD3JYpsz+GKPj62TgAQ6R+Bxql6+rvrVSc+syzNZQgI48ifo1ioo912sjM4lb3fCrjFn8F0zVDZUAhJoM4ETqe/7zOAeIvwX2hEdiRZCWg8QwKE/Tzd/1K1d1cF368jZbglIoKUEOLjfiHPfl0qOQS+gaWgZ8p3BA0IrPoHCOXh+UGnTWH5EjxQfny2UgATKTIDjUGbuW5W5j/atvAQKcQ0epz4aHYHuAvUz6GHiM1FucNizvPjtmQQkIAEJSKA1BIoyg8/dikujrVFeB5hTYHkGcU10HE5+DP9J53qYJgEJSEACEpBAAwQKMYOnnXlT1H448RvQ4yg2A11F/sdRHk/QJCCBJhLgH+dd0R/QJehytFgTi7coCUigwwSK4uCnwmGzAVhsQ/4DAywzWwISGAEBnPm72exdaDv+kX4r4dHox+S/glCTgARKQKAop+gPheVEDi4HEuZO1UdRXg25BkobG7rJhe1XYN28lKA/W4bM+ftbYJ4EepDA3vR5L5x7nvV+GeF5/H7WJboFOid5mgQk0N0ECuHgObhM4eCyHignoHEo1+MfRD9BF7J8FmEjlpdQbDLAiuPInzHAMrMl0IsEHqvrdP6xHlOXZ1ICEuhSAoVw8Dj33FCX503jnPPMaV4skZcLZMb9c5Z/Eh8/e6ZBekBjnYtYGPUxyliQzKX6LDBDAr1JYBLd/hbaJ93n9/Eagm+jcUiTgARKQKAQDh6On0U5wOR9vznIpF3vRP9Gcfg7oR8jTQISaA6B/K5yWv48wiloffQm/km+g1CTgARKQKAoDj43+7yBg0uefc8se3HiuYM+M4sjCXIw0sEHiCaBJhDg95X3TWzJ72tDwoXRd8m7rwlFW4QEJFAQAkVx8P+ExxYcbC4lfAtaoIbP2sSvrUkblYAEmkQAp351k4qyGAlIoGAEiuLgcy3wFDQefQeNxdnH6V+P3ow2RZoEJCABCUhAAg0SKISDz+l4HHreWrcY8YeI5+a6t6NXoj3Je5JQk4AEJNBWAhyLxlLhB1GORdM4Fp3e1gZYmQTmgsDL52Lbpm7KDyeWu+fzTO7T6DfoNKRzbyppC5OABBohgHMfzXq/RvOhy9Ae5OXGxMIcN2mTJoEBCbijDojGBRKQQI8T+Az9P4NJxjfRJJRXat+G8l14TQKFJ6CDL/wQ2UAJSKBDBJal3sl1dV9AeqC3ZdatalICnSWgg+8sf2uXgASKSyBvvnxjXfPeQ/rhujyTEigkgULcZFdIMjZKAhLodQLHAeBerrkvTngJ2hmtgz6INAkUnoAz+MIPkQ2UgAQ6QYBr7nnxz0JoEfQB9B+0IfnPE2oSKDwBZ/CFHyIbKAEJdIoAznwmdX+5U/VbrwTmhoAz+Lmh57YSkIAEJCCBghLQwRd0YGyWBCQgAQlIYG4I6ODnhp7bSkACEpCABApKQAdf0IGxWRKQgAQkIIG5IaCDnxt6bisBCUhAAhIoKAEdfEEHxmZJQAISkIAE5oaAj8nNDT23lYAEJCCBwhHg5URL0qhtUL5MOpnHHf9RuEa2oUHO4NsA2SokIAEJSKA9BHDu46jpZyhfAXwMTSVvC8KeM2fwPTfkdlgCEpBAqQmcS+/2Ydb+l/QS534xwQ8Ip5KXtxH2jDmD75mhtqMSkIAEeoLAHVXnnt4SzyuHp6KVk+4l08H30mjbVwlIQALlJ7AAs/VXVbtJfBTxnKL/bzWvV0IdfK+MtP2UgAQ6RgAfE6fzWfQLdCoa17HGlL/i79PFU2D8OjSe+Pnor8zkM4vvKfMafE8Nt52VgATaTQAnk4nU39DZ6CC0IvoZ+R/C6dxAXGsiAZieDdv7KfKTKDfanYV+jHrOdPA9N+R2WAISaDOBd1PfX3A8X6zUOx0H9BniB6I9K3kGTSQA68spLupp8xR9Tw+/nZeABNpAIN+Tv6Sunn+TTr4mgZYR0MG3DK0FS0ACEphNIC9Z2Z5Z++gaHu8kfndNuquj9C33GCyHxnZ1R0rWeE/Rl2xA7Y4EJFAsApwuvhLH91dadSvhpwlXQ+ujnVDXG33amE7k3oIZaBPSH6PPv+n6jpWgAzr4EgyiXZCABIpNAId3OI7vClq5FroVHU3eU8Vu9dCto095tnwyWpn+5B+YpYhfR3g/6T8PXYJrtJKADr6VdC1bAhKQQIUADu8yolGZbHs6s32cezpFeB/OfWei2yIdfKB00LwG30H4Vi0BCUigywnkMbSn6/rwLGknj3VQOpF0EDpB3TolIAEJjJAAM+TF2fTN6AU0iVnzoyMsqhmb/YFCfkKb8hjgQ4QLkL4Kbd6Mwi1j7ggMOINnoMai2rs+89L+MWiluavSrSUgAQlIYCQEOP6uznbnoRyHX4dmkDeOsCOGU7+Gir+Ect39WMIz0Z7k1z8WSLbWbgJ9ZvAMUpx6HP/X0O9JX1rTqO2I7462rskzKgEJSEACLSbAsXhBqrgavQsHOinVkXcdwVGEe5D3TPLabdR7FvXn3oJl0EOkS/P4X7tZNru+/mbwH6SSJ9HH0O8q8aSjk1He66tJQAISkEB7CeSO9TOrzj1VE88x+nm0XNKdMtrxALoB6dw7NQj91NvHwTNAeVF/bpzIc40TKvGk52XZwuh7xDUJSEACEmgvgcepbonaKpk55xieU/VZpklgDgJ9HHyW4sSfQ8cQvQuNQrlx4gB2pjwSoUlAAhKQQJsJcEy+jSqv4Tj8DbQoirPPC3TOY1k+rqJJYA4Cfa7BV5ey82xE/GKUty4dijZAo8lfjJ3pR8Q1CUhAAhJoL4EjqO6L6DSUx9O+i3ryS2n0WxuCwIAOnu3eh/ZG96G8uCCPZWT9I1HLHDz/QKSOsfwT8QihJgEJSEACFQIcF2cR/YpAJNAIgX5P0Vc2zJeOHkAbo7x2cCrh/Kjpz1zi1HNm4AiUSwK5E/Rh4jPRVLQnaU0CEpCABCQggWEQGGwGn7vlj0PzoLzIYE3CU9HhqNl2PAUujfL4XV55OBMtjFLncdQ9hn8wTiSuSUACEpCABCTQAIEBHTwOdSKO9UHKeCU6C41HHyK/9rl4sppiW1LKBMq+t6a0fJnoKtrwccLDkA6+Bo5RCcwtAX5b76WMfVD+ic/ZvHfyG/wvoSYBCZSAwIAOPn3jx35RTR9vIR61wnL6fzN0Rj+Fb0NeLhVoEpBAkwjg3HekqDwVsx2/88dI5+MgpxLuQvqJJlVjMRKQQAcJ9HHw/MDzpqTPoNxF39/17ws4AGRW3UzLXfo5Y3Ag4TSU6/y5B2ANlDZuhYY0tl+BlV49wIp5y9Icr94dYD2zJdALBDJz3yvOPZ0l/C2/n3WJ5mzaOcnTJCCB7ibQx8HTnX3R7SiOtnYGT3K2PVKNNCvk4DKFg8t6lJcX64xDuR6fek5Gl7I8d442YvmHYJMBVlyJfE8/DgDH7J4jkN9U/Q2zuSw2pudI2GEJlJRAHwePL72+0tf82Ke3o9+VmffXqOs76Dx0FMopxLxB7xcs3592PUN8UGOdC1kh6mOUcSyZS/VZYIYEepPAZXQ7L7PKTD7vNF+Z4NtoHNIkIIESEOjj4PmhX0K/XjlI3y7GkR48yPKRLMpznXeif6DPobTrdSiP5cXZf7EiAk0CEmgCgTj3C/i9/5bwWrQh2pjf9h2EmgQkUAICfRw8ffoCSv4q6BCUd89fifLI2ofRFNRsewsFrp5ZOgec7Ynnxp/pqYR0nHvej69JQAJNIpDfGkVtzu/rjYRj0ffIy0utNAlIoCQE8mjMHMaP/Cp0OZkboC8R/ya6Ev2Q9AFod9Rsu4kC96gUOomw9qa63EV/c2WZgQQk0EQC/K7/jHJWTufeRK4WJYEiEOhvBl9tV+6uHVdNVMKcNn+wLq8ZyY9QyHnMJvYizKN43yT+QcIXUF54kxm+JgEJSEACEpBAgwQGc/CnUMbvcbTvIPwLWh+NR5lRN9WYPUyjnlwC2AKthnI9PnfRZ+Z+PsufI9QkIAEJtJUAx6WFqDCTjdyXdAvHooltbUA/ldGm9ch+G8oE6FzalEmRJoE+BAZ08Ow0/2ZHyo0370K5Hn86ymcJ7yFsulFuHtu5qKKml2+BEpCABIZDgOPfaNb/FfoDugJ9lrxdCN/F4SrOte1G/e+k0m+gg1CO3zeTtyntyVMRmgTmIDCgg89a7DT3E/xgji1MSEACEugNAp+mm7/kOFg9Bl6CMz2BvB2T324E1P0q6swNx/9Hm+5N/eStRPBtwmvJeyx5mgSqBF5ejRhKQAISkMAcBJYjNWmOnJe97Hekk98JW5JKz6o69zSA+O0EeZX3YkiTwBwEdPBz4DAhAQlI4CUCefPlhJdSL0beTZD7gzphqXc8s/UFq5VXZvWbk/YtnVUohi8RGPQU/UtrGZGABCQwAAGczAosysuvxlVWOYCZ5W2VeDcHx9H4e+nfooSXop3Q+mhv1HbLzJ22/IKKZxLm/qjcIzARHcSyGYRdbfQpLzb7FNoYZfKZdzOcQ6iNkIAOfoTg3EwCEph9DXgRONyBtkO5SXYddBoH6/dzcL6VeNca7b+PfuQu+nx8q/oI7+vJf75TnaLun9Kmu6h/B5Qbk3cl74pOtafJ9V5CeVNQzpJkvzqOvs5H/9p+vwN1l8JGAa8UHRmqE+wox7LOUvR316HWdbkEJNAYAX5Xn2XNx/ldfbe6BXnvIv5G8vLaaU0CQxJgn8ms/RPsM7mBcbaRlzMnvyBvy0pWVwb041s0/Gf0I6+Ebqt5Db6tuK1MAqUjkOvB/6rr1XTSmflqEmiUQPaX62tXxiE+QnpMbZ7x4RHQwQ+Pl2tLQAJzEriG5AHMUkbVZJ9CPPlaCwiAeh60OloN1XJvQW1tK/IGatqI7qxYrZF4TtU/XU0bDp+A1+CHz8wtJCCBCgFmWb/lQPwWktcT5sNUb0DXkH9aZRWDJhKAca5N53Ljsyiz3vXIewO8HyfetUb776YfR9KB2wkPJMzp+bVRnLw2QgLO4EcIzs0kIIEXCXBwzgth9kUPodNJ7/3iEv82kwCObx7Km4puhvF+aDfiv0VHo643+nMZnRiH7kF/Q7mBcAahNkICzuBHCM7NJCCB/xHgQPxnUpHWOgLjKXoSrDPTnW3ED8bxn4GWJR7H2NVGH+6gA5HWBALO4JsA0SIkIAEJtIFAHs97rp96Xklexx7d66c9ZhWEgA6+IANhMyQggeYRYEa7BPoA2h+t0bySO1cSs9u8VyAvuTm52griZxDP4873VfMMJVAloIOvkjCUgARKQQCnN46OTEQLoJnoRvLeRlgGy/0Oa9Gf3Nx4NvG89CbvHdAk0IeA1+D7IDFDAhLocgLn0v59mNX+Jf3AEeZzrycRTiVv9lfYkt+NRvufpN0T6MsShCRnPdiN/bDN7SHgDL49nK1FAhJoH4E7cXyznXuqJP4fgtx9vkrSZTD69IDOvQwj2do+6OBby9fSJSCB9hMYwww3z1HPNuKjiOQU/X9fzPGvBHqDgA6+N8bZXkqglwicRGd/gl9/LRpH/Dx0LTPevxN2vdGnndAPK9qm6ztkB1pGwGvwLUNrwRKQQCcI4MjPwvndT92fQqNRbkb7Cep6o19fpBP5ct97UY7fnycvH9H6EXFNAnMQ0MHPgcOEBCRQBgI4vMn0IyqN4ciXozN5e91r6d/s597J25v0bwjPJc8b7koz2s3piKfom8PRUiQgAQm0mkBeaHNp1bmnMuL5GEuegR+btCaBWgI6+FoaxiXQJQSYsb0q6pLm2kwIMF6j0fIojnoklmfel2H71aobE8/HfTZHeVJAk8AcBDxFPwcOExIoNgEO6Pn++uEop2sXJp2viu3ATO4ZQq2gBBin19C0I1BOrecLcCcwZscPp7msP4PtUsa/CN9PmOP3wWgrlj1FqElgDgLO4OfAYUIChSdwGS2cnwP6TujtxPMSl68VvtU93ECc8WJ0/2Z0ImO2C2E+g5qX1exIOCxj+zzfvxLKzYOz0ObkXU+oSaAPAR18HyRmSKCYBHAI42jZHRzQP1xtIfHjiK/IsiWreYaFI7AFLfoEY3VpWkaYsy1fQCN6xI3tb0c/RD9G01OmJoH+COjg+6NingSKSSCztv/207Scqs23wrViEpifZuWd+LWWSyuvqM0wLoFmE9DBN5uo5UmgdQSmUfTLma1/ploF8W8QX4WZnDdZVaEUL5xEkz7CWL02TSMcRZAZfC63aBJoGYH8569JQAJdQAAn/jy+4SCa+g/CNQnziFScxUZIKygBxu1Oxmt/mncl4amEy6OryT+hoE22WSUhoIMvyUDajd4ggFN4BCcRB7Eqyk1WN5H3Qm/0vnt7yRhdzbitQA/Go0dJ52yMJoGWEtDBtxSvhUug+QRwDs9T6j+bX7IltpIA4zaD8qe0sg7LlkAtAa/B19IwLgEJSEACEigJAWfwJRlIuyEBCTSfAKfV8/hhHktcBN2Kvlc5g0JUk0CxCTiDL/b42DoJSKBDBHDuC1L1r9Cd6Mcor4W9hvz5CDUJFJ5A4Rw8P5550aKFJ2cDJSCBshM4hA4ez4z9FHQD2p30BWjXsnfc/pWDQCEcPA49H2E4AuVjCnnL08PEZ6KpaM9yoLYXEpBAlxHIRKP+prg/kZdXz2oSKDyBQjh4KOWjC3kJxNZoYZR2LYv2Qfvj5D9EqElAAhJoJ4G7qewddRUeQHp6XZ5JCRSSQFFustsSOhM4BXZvDaU8UnIVzv3jhIehE2uWGZWABEpGgN/6QnTpK2hdlH/0f8ox4duEwzLKyct/Mil4F1oA5ViyI2Xl7OBw7DusfDvFZbLxe7QTophZvyTUJFB4AkWZwU+F1GYD0MoHGR4YYJnZEpBACQjgRPMu/evRDBzoWwnfgtYmfySX6L7Etmugd6JNUJzztylrWMc72vEo2y2N7kEp50qUMjUJdAWBUezEHW8oP7z1aMRE9BjKG57yw8pjKfmR5ixDvnd8B+GgRjnLs8KrB1gpZwLymc0dBlhutgQk0CEC/Hbj0Pfg97l3tQnkZUZ/Nnn5LG7DxnbXsvIb2O656kbk5Qzgj8j7azXPUALtIMC+9y3q+Rn7XvbLtlohTtHT8SlAiJOfgMah/Nf8CDoV/Z7ljf4Xkuv4+U+7P1uZzJSpSUACxSMwhiZNr2tWvsCWU+zDtfs5ZLzk3Csbp6x81U2TQM8QKISDD21+kE8RXFolj8PfjXium+WxlIaMdXMqLupjlHcsmUv1WWCGBCRQBAKZ3Xw6/+jzO55SaVBusu3v87hDtfdmyjmEcr6WFYnnn/58pOewpDUJ9AqBQjh4foA3A3zxOuijSeeZ+B0Jz+HHumfdcpMSkEBJCPD7fpDf+qfoznWEccwLo7EoN7YN1/L8+t8pZ0XCfEb39WhV6sglQE0CPUOgEA4e2nHep6DTUU7Lx7ZDOWV/MMrpNU0CEigxARzw9Tjl/KP/f+hJdBV5+bDOsIxtZlDOSmyU40cmCseQN5IzAcOq15UlUDQChXDw/Piu4Ae5AXC+i45B+6EH0eMsu4NQk4AEeoAAv/eH6ObFc9tVysk/BlfMbTluL4FuJlAIBx+A/CBz5/weOPqdCSejq9Gw/3tnG00CEpCABCTQ8wQK4+CrI4GjPxMnn+dNT0R5LlaTgARaQIDf2RIUuwfKte5r+O2d34JqLLJCAN7bE10b5V6AE+GdyxCaBFpGYFgvfmhZK+oKZsefjrZFuVlGk4AEmkyg4tzPpNgZKO9XP4m8rzS5GourEIDtSUTfja5CuS/g8coYENUk0BoChZvBt6ablioBCdQROIH0UfwTfWHycTa54/x0wjeS9+fkac0hANO3UdL6cM19RrGLyMtruT+NPpMMTQKtIFDIGXwrOmqZEpDAHATyYpnc6zLbcD653yXvoVjhxRz/NpFA3q6Z93DUWi6HyLqWiPGmE9DBNx2pBUqgKwg8TCurM8pqg99LJHexa80lkDdo5pG9WluHxHA/flO7vXEJDElABz8kIlfoRgKcAt0ZnYV+j45B+ZhJ11v6gSai6ehelOvoI7Gvs9Fktt8RrYvy/okHmcn/cSSFuc2gBH7L0vEwvhitg3It/jCUU/SaBFpGQAffMrQW3CkCHED3p+68KOmj6D3ocXQ0KoPlNa7vQJtWtBr9/RXxYRmO/EY2WBptiPZFk1Bm8FqTCcB6FkVujfKZ2bzj4/VoF7LvI9Qk0DIC3mTXMrQW3EECcewbcACtPoZ0KE7wx6j2PecdbN7Iqqb9G7PlePqVr6zNNvLeQCSz+TXJj9Nu2CoOxpu8GiY28hVhHSd/8shLcEsJDJ+AM/jhM3OL4hO4g+Np1blXW5u7lvO8dzfbq2j8nbUdoJ9Pk84ZiiVq841LQAIS0MG7D5SRwH+Y0eYFLrMts1siud459cWcrv17OS1fhP7kEsRsI/41IsujUnznnP6MQWuh17zYw5H/pYxV0OtQnjsvjdGf3IexOlqlNJ2yIy0h4Cn6lmC10A4T+Dz1380BcFXC+1GuLb+Z2W7uHO9ao/0P0ae96MAFhFsRjkKbok1YNpOwq40+5R+VY1DuOs9NaTkzsRN9e5awYWO7OPTcc7E4yra5kfA1lPMA8a42+rEIHQijfO8+/+y9jjDvLggrTQJzEHAGPwcOE2UgwMEuNy/ldPw1KJ8L3Zm8PxN2vdGP39OJzNxyNiKvcl6LvLyJrqsNR5X7CnL54TT6sy96G/E8RjaSewTOYrtbKWN39EHin0AnUMcYwq412p8nQf6O0rf9UP5xvQCV5QZSuqI1k4AOvpk0LaswBDj4PYnORb9E0wvTsCY0hP7cgj6PDkG3N6HIIhSRO8uPpD95pKxquxLZpJoYRph/7r5bXZ8yf0w8X6Vcu5rXpeFKtHsy/Tm82n7iufS0GM5/2WqeoQSqBHTwVRKGEpBAJwk8T+Uv1DUglxBHMut+iu1y+aLWcmo7dXSzhU9OzddbGfpW3yfTTSCgg28CRIuQgATmmkAup+QlMB9JSYRx0F9CI7n8kJsRcyljtlFUTvNvg657Mac7/zJbv5WWP0Z/8uGa2Ub8Z0TmZZnP1FeYGPyPgDfZ/Y+FMQlIoEMEcFBP4qxyvfxqws0IM3O/AuWNe8O1bLMm5UwizE2Wj6LXUke3z+Dpxux7EibRt3OJ5wbCOP1w0yTQh4AOvg8SMyQggU4QwAE/gONahbqXQU8nPZJ2sF1OY7+PsnIXfe6ov5e8+tP/Iym649vQj7zfYUP6tiThC6Qf7HijbEBhCejgCzs0NkwCvUcAh5VZ9vRm9LzMzo++5cyEJoFBCXgNflA8LpSABCQgAQl0JwFn8N05bra6TQQ4FTqeqj6JcqfynegrzJ7yethhGeUswAYfQ+uinGY9lHKaMlOlLE0CEpBAHwLO4PsgMUMCLxLAKS9KbCK6EOXteLehc8l/BWHDxvr5nV2N8k/Cwegn6Ezy1yLUJCABCbSEgA6+JVgttCQEDqEfmbH/Ft2Ffkj6MrT7MPu3I+tfx/Z5Oc2daDLpOPoDh1mOq0tAAhJomIAOvmFUrtiDBF5Jn2+p63deEZuZ+HAs6/+hboMbSad8TQISkEBLCOjgW4LVQotAgFPgeYXnUqj+rWaNNi9OuH62/lXykj8c+zcrv5NmzFez0bbE765JDytKWa9Ciw1ro35WpoxF0bJoru7HYfsFKuXknfKaBCRQAAJz9aMuQPttggT6EMDR5KMcn0NvqCzMl8kmcGr8sT4rD57xPRb/hW2XIzwfvRP9mXJq35dO1uDG+pdTRl7eMo3wU4SroTeid6NhGdvnZr3D0QpoYdJ5Leu7qSMfZhmWse2H2WB7NAOtR/otlDPsfzrY7k1sn37li3Ybkf4Y5ZxHXJOABDpIQAffQfhW3TICJ1ByXt35LhzNLBzOzsSPJ9yHZMOfHmXdvF0td73vgpZGE8l76RWopBs2tvsKZeUVqrmxLnfPxynnbvrh2iQ2mMK2s/85oMxcx/8a+gxq2NhuX1ZeB21FWc+SnkD8ZMJdScfhN2SsP54Vr0Crst3NpBcnfjrhA6RzY6EmAQl0iIAOvkPgrbalBNag9A/jYGalFoLcsR4HFud6bfIaNbbNi1dOb3T9wdajrEtZHo3I6MNKbJib9PavFkD8WPJ/gZYkPpyXn2xNGS/9w8O2V1FG2rYxGs7sO2cAtmf7mwnD+kHKyT8cefe7Dj5QNAl0iIDX4DsE3mpbSmAmjqb+1aS5Dl97DbylDWhR4Wl/f7Pr5A/3n/Xwea6unfmHaLhfb5ufbeovD6Tc4ZZT1xSTEpDA3BLQwc8tQbcvIoEbmEWeVG0Y8fcSPwBdX83rREg7Xoveg/4fGsmNf9No9zxselC1/cS/Tnw1/qG5p5rXYPhH1jutui7lvJ74N1EuIwzH8nTA0Wy/aDYijGO/Cl2UtCYBCXSOwHD/6+9cS61ZAo0TOIxV/4Cz+TXhf1DuNl8eJ5gb0jpitGUPKs4/GhegPBf/ZfJyU9vTxBsy1n2ebfJWvX8SvpYw9xPkH4UN0XDtRDZ4PeVcSZjLFq9Ga1FH7l1o2Fg/NyHmlPzfCX9BuCr6IPkXN1yIK0pAAi0hMIofYksKLlqhHHyOpU1L0d9di9Y229MaAoz56pScGeWNjHv9aeTWVNpPqbRjPbL/ghajHbPv5K/sj/nK2dH9bDJoFtvmKYH0LafZb6KM3CcwIqOsOOSxKOXMbttICqKcJdguTxvcTznDvhN/JHW6jQS6gQC/jW/Rzp/xuxjW/T/N6Jsz+GZQtIxCEuAH9a+CNGwD2rFfnQM9grwfjKR9lBOH/o+RbFu/DWXdVJ83kjTl5NOukSYBCRSEgNfgCzIQNqPUBJ6gdzkFXmvLkBjxzLu2IOMSkIAE+iOgg++PSo/kcepoJXQUyjPiu/dItzvRzdwL8AYY5wa70SgvujkGHd6JxlinBCTQGwQ8Rd8b49ynlziZFcnMXdS5PjQd/Yi8TTjVujdxrYkEYPoEbHPvR07J74seQV8lv+3X5KhXk4AEeoRA4Rw8B8K0aSwHvxwEtdYROJ6iD4bzlZUq1oH9Keit5F3Sump7s2SY5ga23EWvSUACEmgLgUKcosep5LTlEeguep27nR8mPhNNRXu2hUTvVZKxv66u23nzWO6G1iQgAQlIoMsJFMLBwzCzyTzXuzVaGKVdy6J90P44+Q8Ras0lcDvFbVEtEsZ5njqc76zmGUpAAhKQQPcSKIqD3xKEeYzoBvQ4is1AeSPWx9F23Yu4sC3PG9DOwa9/GL2N+NnojxXmhW10tzYMxkujH6C8gOd69M5u7Ut9u+nLmuin6Pfob2j9+nVMS0AC7SdQFAc/la5vNkD389EKn68dAM5Is3Hkmam/EuUlJ5uin5B3EKHWZAI4vAUpMm/UuxTG+Wdqc7QP+S+dQSHdlUYflqLheSb/+/TtHYS5z+Ab5L+OUJOABDpIoBBvsuNgkDd9TUS5EWkaehQtgtZAuekun7S8g3BQo5y8SSun9vuzT5A5P+Xs2N9C8yTQKgLsl7tR9nj2va9W6yBvVeKHkJdX2Hat0Y98pjZv5Dut2gnyckZuM/I+V80zlECvEuD30NtvsuNAMAUIcfIT0Di0NMqsPe/LvgLFyTdia7HSpgOs+BryvTN/ADhmt5RAXpd7d10NM0jnfpNut5yduLeuE/mdvaIuz6QEJNBmAo06zpY2C+e+PBUcjnZAue7+IZz+LYT5OlVO+SV/56QHM7a5gOVRH6OcY8nM6UStBwgw3tm384rY+dB17Bs5O9Qpu5KK856B82nHfZVGfJRwSqca1MR6L6esQ+nbZPr2FOHLSV+E9m5iHRYlAQmMgEAhHDztPhDlGmUOyLuiyRwoNuWAcRNxTQLDIsC+k1nl91Bmks+i7E85RX4b8bYb9ebrb0dR8W2ERxKuhHIZKvt9Vxt9+yN9WodOpI8nE+ZM3DHkn93VHbPxEigBgaI4+K1guR4HhScJMxu4kfBCwjeXgLFdaD+BzCBz09fpqZr96EKCYwl3J68jM3nq/Q31r0Y71kR5WiGz+lIYfTmGvp1PZ3IPzO9IX1+KjtkJCXQ5gaI4+Dj0zN5zuu9lHCB+zgFjWaI53X5S8jQJDIPAU+xDs517tiGeWebbia6N/pS8Ro3tcl/IQWhR9DTak/Kqp9lJNm5sdxdrR6Uz+vZvOhVpEpBAQQjkelkR7Ps04pccTA+uNoYDxjHEc5ov1841CQyHwCz2pXnqNsh9HnHQDRtljGPlr6OvoLeh7JNnkL84oSYBCUig0AQK4eBx5jmlujL6XS0t8g8jnUduLqzNN14sAji8edFu6PMoz3fXO9d2N/g8KvwW7ZgPxb5JOqfG/zbMhnyS9T/HfpgXMMX+QPqXKDd+ahKQgAQKTaAQDj6EOHjORH+vp0XeJPSj+nzTxSAQ70lLLka5ueoatD66kewFCDtl36Hi3MR2Bco/j7PQBPajhMOxsaz8YN0GeSTMR8DqoJiUgASKR6AwDr54aGxRgwTyopZH8J2fQhej/Un/DO3b4PZNX402xA5FG6It0KfREyOoKNfrD6luxz8tY4hPRKW5Qa7aN0MJSKB8BHTw5RvTdvco17a/XVfpuaRXqMvrxuQpNDqn+S9DBxL/OTqIfxYu78bO2GYJSKC3CBTlLvreol6u3v6X7myELqvp1tbEH65Jd2UUR/4CDX8Pzn1jwsVQzlBM7crO2GgJSKDnCOjge27Im97hzHL/hRNckjDxOPsvojjEUlgzZuzwWQQYn0PrVqB8nnKvLQUgOyEBCRSSgKfoCzks3dMonFSubY9H09FeKN8ReDX5eWmRBgGce16XezPK8/M7oE+ifHFtAqEmAQlIoCUEnMG3BGtvFYozf44ef6u3ej2s3u7M2mfAqfpOhzxlkJfnxNFfNaySXFkCEpBAgwScwTcIytUkMBcEFmLbKXXb30M6j+FpEpCABFpCQAffEqwWKoE5CPyV1C7M2uPoq5ZvxOe0vSYBCUigJQQ8Rd8SrBYqgf8R4NT833DuZ5NzF+HnCVdFS6A9kSYBCUigJQScwbcEq4VKYE4COPmTydkEPYouQe8n71lCTQISkEBLCDiDbwlWC5VAXwI49BvIjTQJSEACLSfgDL7liK1AAhKQgAQk0H4CzuCHYM410zDaES2KbmEWli+KaS0iAO+1KfpNKM/RnwXvx1tUlcVKQAISKDUBZ/CDDC/OJnzy4ZTXo0fQMeR9h1BrAQHY7kSx30Nx6suix8hLqElAAhKQwDAJOIMfHNgnWDyFWeRRWQ1ncybBGYRbk3d+8rTmEIDpSpSUV92Og+1DKZW8aQRfRvsiTQISkIAEhkHAGfzgsF7L4l9VV8Hx5HvimdG/rppn2DQCq1HSUVXnnlKJ/4Jg5cQ1CUhAAhIYHgEd/OC8ZrB49bpVNiKdfK25BPL42BzOnBn8MuQt1txqLE0CEpBAbxDwFP3g4/xdFp+Oo8k14RtRrhHvj5ZCWhMJMFu/Es47oyMp9gSUmxq/gT6LNAlIQAISGCYBHfwgwHA6t+JwdmCVOJ35Ua4J5xrxM4Ra8wkcSJEHo3yU5TF0JKxrvzNPliYBCUhAAo0Q0MEPQQkHcy+r+ErRITg1YzGsc4/D7Bsam1GeZUhAAhLoZQJeg+/l0bfvEpCABCRQWgI6+NIOrR2TgAQkIIFeJqCD7+XRt+8SkIAEJFBaAjr40g6tHZOABCQggV4moIPv5dG37xKQgAQkUFoCOvjSDq0dk4AEJCCBXibgY3I9PPo845/xfztaEF3HU2o39zAOuy4BCUigVAScwZdqOBvvDM49L+75CdocvRr9k7x3EWoSkIAEJFACAs7gSzCII+zC8Wz3J2btJ2Z7nPtEgnwp71/k/Tt5mgQkIAEJdC8BZ/DdO3Zz2/J8nvW0aiE49fuJ58t561fzDCUgAQlIoHsJ6OC7d+zmtuX5elv9R3PyedyZc1uw20tAAhKQQOcJ6OCHGANOWS+BvoSORR9GZWF2Ml0/nv6siOZH+Wrbtuj8IZC4WAISkIAEuoBAWZxVS1Dj9Bag4F+je9Bp6M3oavLnI+xq45T8hXQgX207FcWp5/Osa5L/HKEmAQlIQAJdTsCb7AYfwC+w+Ls4vZ9XVtsV555Px74HnV7J69qAfv2BxkeaBCQgAQmUjIAz+MEHdHEWX1u3yhWkl6zLMykBCUhAAhIoFAEd/ODDMZ3FW9Stsj/pu+ryTEpAAhKQgAQKRaBwp+g5BZ42jeX08SMFIJVnxe+gTcsQXoR2RqNp2y8JNQlIQAISkEBhCRRiBo8DHY2OQJkZP4MeJj4TTUV7dooejnwGded0/IPorehqtBXSJCABCUhAAoUmUJQZfGbKS6Ot0a0oz2IvjNZEx+Hkx+BsZ79xjfSAxnrLsTCz7f4sz3zn9azDMurNPxzHDWsjV5aABCQgAQl0mEBRHPyWcJiAM723hkdmz1fhtD9OeBga0sGzzlpoM9SfxfFP72+BeRKQgAQkIIGyESiKg58K2DjmM/oBvA15D/ST3yeLfxAuIDPqY/yj8B4y86y3JgEJSEACEig9gaI4+EMhPREnfCDhNJTXqC6C1kBpo9e9gaBJQAISkIAEGiVQCAfPzHsKzn09Gj0BjUO5Hp9Ze07LT2b5LEJNAhKQgAQkIIEGCRTCwaet+PCnCC5tsN2uJgEJSEACEpDAIAQK8ZjcIO1zkQQkIAEJSEACIyAwqlfOfnMJYF345KMqU0bAqZmbbERhueTgZYdmUu2/rJyhmgc93f9ic5tM4BWUl0dctdYTGE0VOYY82/qqrAEC4T15hCTGs90W+Nq7R7j9iDfrGQc/YkJN3pB/NM6hyPcz2HkMUGshAVjn5URvgvVXW1iNRVcIwHsSrDcVSOsJwHo/apkB7+qHsFpfaQ/X0K37tqfoe3intesSkIAEJFBeAjr48o6tPZOABCQggR4moIPv4cG36xKQgAQkUF4COvjyjq09k4AEJCCBHiagg+/hwbfrEpCABCRQXgI6+PKOrT2TgAQkIIEeJuBjcm0efB63WIoqH+DxlhfaXHXPVQfrMXQ6nxr+b891vgMdhvcysP5PB6ruuSphPZZOvwBv3zvQhtHv1n1bB9+GncMqJCABCUhAAu0m4Cn6dhO3PglIQAISkEAbCOjg2wDZKiQgAQlIQALtJqCDbzdx65OABCQgAQm0gYAOvg2QrUICEpCABCTQbgI6+HYTtz4JSEACEpBAGwjo4NsA2SokIAEJSEAC7Sagg283ceuTgAQkIAEJtIGADr4NkK2ivQR4KcV87a2xd2uD9bxoVO8SsOcSKC4BHXybxoZj4Abozjq9uk3V90w18N2Fzl5V22HyNkVXoNvQr9GitcuNj4wAHJdnyzvQ+GoJ5LmfV2E0KYTpLugSdD06Ha1RLZr459ANKPv256r5hiMjAMP50DfQXys6knB0SiPsun1bBz+y/WAkW23ARhej1Wt0z0gKcpu+BPjxLYq+y5Jvo5dmlOQtTnoi+jBaFd2GvoW0uSAA173Y/FK0RF0x7ud1QOYmCeel2f5YtAuvpV2H8I/oOBSHsxPB1mhjNAG9h7z/R6iNnMD72XRlFJ7RmmgPFOu6fVsH/+LAtePvulRyNVoSLcGP9Qk0qx0V90gdm9PPJ1B+oLWWH+U/QX0Depb48WiH2hWMD48ATiQzmp3RVqj+Pf/u58PDOdTaOUbvzL57X2XF6wk3qsTfQXg6y2age4mfgbavLDMYGYHw/XSOFRHxG9GbKkV13b6tg6+MXBuC7ByfQheh2zlIHt2GOnumCn6MZ6HP0OEn6zq9AunaD6DkQLkI/OevW89kgwTg/Ax6O7qpn03cz/uBMtIsGN+DJtdsvy/x8yvp+n07Tj4fs9JGSADW16Bp2ZxjxCsIdkXnJY113b6tg39x4Nrx929Ushc7T04T/x/6KDtQ/enNdrSj1+p4FR2eWdPp6j8AC9bkGW0eAffz5rGcoySOF3uTsS3KRCFWv2/nDFackjaXBGCds1Q/R3H4Z1eK67p9Wwc/lztCo5uzk3wEXZ71CacQ/Al5qjhAWmsPUvzCNVWMJf4UY/BITZ7RJhFwP28SyLpicDj7kfU19DYYT68srt+3s5/fU1lmMEICFef+KzafB2UGP9u6cd/WwVdHr4UhO8wY9KWENdVkBvlATdpoawjkYDiupujE76pJG20SAffzJoGsKwauua/kyyjO/Z81i7Nvr1iTHkfcfbsGyHCjsJ6XbTJzj3PfAd7PpIxu3bd18Bm9Fhs7yVNU8Va0V6piZ9mQYD10UdJaSwlcQunjYb45ynX3g1D1lFtLK+61wt3Pmz/i7LMrUeoJ6L3oHtKLRZWaziT8AOll0TjiWefXSBs5gQPYdGX0QbRghfdC3bpv6+BHviMMd8s8o7orO8xNhL9HuR7/+HALcf3hEYDx02yRH+05KDfPLI8OR1prCLifN5frhyku19UnoYeq4jiSM4AXolwX/ge6Cv2M/f2vhNrICXyCTddGudRR5f2LSnFdt2+PYoeotN2gHQTyHyH1/BfuL7SjPut4kQDcc+ptLNy99t6GncL9vA2QK1XAOtfen2bfzj+zWosJdNO+rYNv8c5g8RKQgAQkIIFOEPAUfSeoW6cEJCABCUigxQR08C0GbPESkIAEJCCBThDQwXeCunVKQAISkIAEWkxAB99iwBYvAQlIQAIS6AQBHXwnqFunBCQgAQlIoMUEdPAtBmzxEpCABCQggU4Q0MF3grp1SkACEpCABFpMQAffYsAWLwEJSEACEugEAR18J6hbpwQkIAEJSKDFBHTwLQZs8RKQgAQkIIFOENDBd4K6dUpAAhKQgARaTEAH32LAFi8BCUhAAhLoBAEdfCeoW6cEJCABCUigxQR08C0GbPESkIAEJCCBThDQwXeCunVKQAISkIAEWkxAB99iwBYvAQlIQAIS6AQBHXwnqFunBCQgAQlIoMUEdPAtBmzxEpCABCQggU4Q0MF3grp1SqBHCIwaNWpjdG2PdNduSqBQBHTwhRoOGyMBCUhAAhJoDgEdfHM4WooECk2AWfS86BT0X3QHOjgNJlwLnYF+gu5H56IlK8uIjjoETUd3oy+gUZVlSxD9FUp516O3JD9GfAf0d3Q7ye2Tp0lAAu0noINvP3NrlEAnCBFKjJgAAANPSURBVOxIpa9BK6N3oDjrpMeg96L70OroXvQDFNsdvQ9ti7ZDu6A3oNgpaAbKNsehH6M49/EEJ6Mvomy3OdIkIIEOENDBdwC6VUqgAwSepc7l0UboVrTErFmzbiGMPYe+Rvphwm+grXHU8xO+H8VxT0P/RnHq27JsMcKt0NfRE+hsdA/5axNugf5BWeegvxOf7fgJNQlIoM0E5m1zfVYnAQl0hsA5VLs++hEajU7DIX+aMBaH/FgihLeQH6e9Ino1yjqfQFXLDXPLoVnokmpmJcw/DzlD8Lea/KuJ71GTNioBCbSJgA6+TaCtRgIdJpAZ+bfQl9CW6AT0DxSHvRJOfR6c+/OEK5BeEN2FrkGTyZ99yp5lC5GeBz2Jcnp+LZY9SJhT80tU8vYj3Dp5FVupGjGUgATaS8BT9O3lbW0S6BSBXGc/E2XmfQHKKfeqLUxkO5x0bqD7APoTjjtO/Fy0J9mLVpadTvpAlj1D+Ef0EfJfjpYmfiPK9fjL0BvJWxWNIb4z0iQggQ4Q0MF3ALpVSqADBH5KnTNRrqffiV5AE1EsN9hlZp/8nVBurIv9DuWmu9vRzSiz96NR7Ei0G7oNXYm+ieO/ISJ+SCXvFsKnkCYBCXSAwCh+kB2o1iolIIFOEKjMqkfzu3809ZN+PcFPSa9OfDHCh+vbRf4rksey/IMwh7FsCTIeZNkcBxLy5yN/DNmzr+3PsZEJCUigLQS8Bt8WzFYigWIQwOFmRt3vrJplfZx7Wk1+H8de7Q3LHqjGa0PynyUdaRKQQIcIeIq+Q+CtVgIFIZBT7F8oSFtshgQk0EQCnqJvIkyLkoAEJCABCRSFgDP4ooyE7ZCABCQgAQk0kYAOvokwLUoCEpCABCRQFAI6+KKMhO2QgAQkIAEJNJGADr6JMC1KAhKQgAQkUBQCOviijITtkIAEJCABCTSRgA6+iTAtSgISkIAEJFAUAjr4ooyE7ZCABCQgAQk0kYAOvokwLUoCEpCABCRQFAI6+KKMhO2QgAQkIAEJNJGADr6JMC1KAhKQgAQkUBQCOviijITtkIAEJCABCTSRgA6+iTAtSgISkIAEJFAUAjr4ooyE7ZCABCQgAQk0kYAOvokwLUoCEpCABCRQFAL/HyByKW1tHRa9AAAAAElFTkSuQmCC" alt="plot of chunk unnamed-chunk-2"/> </p>

</body>

</html>

test.htmlWhich would you rather read and write?

Wednesday, 11 September, 13

Page 26: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

knit2html() function from knitr package does R Markdown to HTML conversion, i.e. hides the two separate steps

RStudio’s Knit HTML button does practically the same thing

Admire, publish, etc. as you did before with the “notebook” you compiled from a plain .R script.

Wednesday, 11 September, 13

Page 27: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

What a finished piece of coursework will look like:

a readable report, posted on the web; default is to use Rpubs

+

a .R or .Rmd file, posted on the web; default is to publish as a gist

+

a line in course repository giving your name and links

=

Wednesday, 11 September, 13

Page 28: STAT 545A Class meeting #3 Wednesday, September 11, 2013jenny/STAT545A/2013-lectures/cm...STAT 545A Class meeting #3 Wednesday, September 11, 2013 Dr. Jennifer (Jenny) Bryan Department

Instructions for how to submit coursework.Draft exists. Will be revised today.

Wednesday, 11 September, 13