updated install instructions

INSTALL.md

@@ -54,9 +54,12 @@ version](./period-luminosity-relation/TutorialPLZ_R.ipynb) of that tutorial.
 
 #### RStan
 
-Follow the instructions
-[here](https://github.com/stan-dev/rstan/wiki/Installing-RStan-on-Mac-or-Linux) for the installation
-of Rstan. Where the last step is to do from within R (run as root for a system wide installation):
+The recommended installation method for RStan (working with RStudio) can be found
+[here](https://github.com/stan-dev/rstan/wiki/RStan-Getting-Started).
+
+To install from source follow the instructions
+[here](https://github.com/stan-dev/rstan/wiki/Installing-RStan-on-Linux) for the installation of
+Rstan. Where the last step is to do from within R (run as root for a system wide installation):
 
 ```
 > Sys.setenv(MAKEFLAGS = "-j4") 
@@ -65,7 +68,7 @@ of Rstan. Where the last step is to do from within R (run as root for a system w
 
 Note that a `Makevars` file in the ~/.R folder with the following contents is useful for
 optimization of the compiled Stan models and supressing of irrelevant warnings (see [RStan install
-instructions](https://github.com/stan-dev/rstan/wiki/Installing-RStan-on-Mac-or-Linux)):
+instructions](https://github.com/stan-dev/rstan/wiki/Installing-RStan-on-Linux)):
 
 ```
 CXXFLAGS=-O3 -mtune=native -march=native -Wno-unused-variable -Wno-unused-function -Wno-builtin-macro-redefined
@@ -129,9 +132,14 @@ version](./period-luminosity-relation/TutorialPLZ_R.ipynb) of that tutorial.
 
 #### RStan
 
-Now follow the instructions
-[here](https://github.com/stan-dev/rstan/wiki/Installing-RStan-on-Windows) for the installation of
-Rstan. Subsequently make sure that your can run R notebooks by following the
+The recommended installation method for RStan (working with RStudio) can be found
+[here](https://github.com/stan-dev/rstan/wiki/RStan-Getting-Started).
+
+To install from source follow the instructions
+[here](https://github.com/stan-dev/rstan/wiki/Installing-RStan-from-source-on-Windows) for the
+installation of Rstan. 
+
+Subsequently make sure that your can run R notebooks by following the
 [instructions](https://irkernel.github.io/installation/) on
 [irkernel.github.io](https://irkernel.github.io). Step 1/2 should be done while running as R as
 "administrator".

period-luminosity-relation/TutorialPLZ_R.ipynb

@@ -57,20 +57,9 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtAAAAIcCAMAAADoo2lhAAADAFBMVEUAAAABAQECAgIDAwME\nBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUW\nFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJyco\nKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6\nOjo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tM\nTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1e\nXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29w\ncHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGC\ngoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OU\nlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWm\npqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4\nuLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnK\nysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc\n3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u\n7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////i\nsF19AAAACXBIWXMAABJ0AAASdAHeZh94AAAgAElEQVR4nO2deYAU1bWHT/f0zMDMCMOAyK4O\nKGJigCAoLyioiAuKshlxIuIGPjRGYlTkYSCgEVEJLonwnsElKoosohh98IILmwqCSMRowA0B\nBZkRBlmEoV/Xreqe7p6umq7bt2s5/fv+8PStvnXqzO2PsmvpWxQGgBHkdgEAqARCA1ZAaMAK\nCA1YAaEBKyA0YAWEBqyA0IAVEBqwAkIDVkBowAoIDVgBoQErIDRgBYQGrIDQgBUQGrACQgNW\nQGjACggNWAGhASsgNGAFhAasgNCAFRAasAJCA1ZAaMAKCA1YAaEBKyA0YAWEBqyA0IAVEBqw\nAkIDVkBowAoIDVgBoQErIDRgBYQGrIDQgBUQGrACQgNWQGjACggNWAGhASsgNGAFhAasgNCA\nFRAasAJCA1ZAaMAKCA1YAaEBKyA0YAWEBqyA0IAVEBqwAkIDVkBowAoIDVgBoQErIDRgBYQG\nrIDQgBUQGrACQgNWQGjACggNWAGhASsgNGAFhAasgNCAFRAasAJCA1ZAaMAKCA1YAaEBKyA0\nYAWEBqyA0IAVEBqwAkIDVkBowAoIDVgBoQErIDRgBYQGrIDQgBUQGrACQgNWQGjACggNWAGh\nASsgNGAFhAasgNCAFRAasAJCA1ZAaMAKCA1YAaEBKyA0YAWEBqyA0IAVEBqwAkIDVkBowAoI\nDVgBob3B8zN1VrtdiN+B0J6gitqUa5Se63YlfgdCe4JdtF7E8X1dLsT3QGhPAKFVAaE9AYRW\nBYT2BBBaFRDaE0BoVUBoTwChVQGhPQGEVgWE9gQQWhUQ2hNAaFVAaE8AoVWRodD/nDNndY2a\nSnIaCK0KKaHfGtTprAcOhsMbu1OE1s+pLiqH2LtGsBRCK0JG6FkBzeP/OPBtO9J5RHlZOcN/\nGkNIr4gmhM4UCaG/bEh0VAOiCbdQs9sef+yGxlT4mfrCcoQRwyo13qcXRRNCZ4qE0P9FLd44\nUvP3svK2rbdp7a/LaazqsnKGESNE2AShFSEh9Gk0SwvTiWboC16g01SWlFNAaMVICN2Edmhh\nM9EX+oJdVKaypJwCQitGQui8kAj7KXg4cQmwD4RWjITQRfSDFrYS7dcX7KOjVJaUU0BoxUgI\nfSKt0MICojf0BW9QJ4UV5RYQWjESQo+g8w6Ew5Wdi4rOEt85avrS1arLyhkgtGIkhH6L6ORb\nbmhLV1xM57z1/XeLzyZ6S31hOQKEVozMlcLh4tJW6b/ezTOucl2mvCxGVG3W2ak3q43mdr0J\noRUjI/TBWwqJTl0dDj/bUPg8eJ/yshjR1fhXf7TeHGA0Q9WiCaEVI3e33Z4P9B3MV3dfNvS2\nZSrr4UfHu8Ue+fFivdn3RtF8lXaJJoRWDO6HzjYd9eupi6JCjxdhPYTODhA620BoR4HQ2QZC\nO4oCobeMGjUq8yxsgdCOokDodZFj9qRF1b+/I8aYQZlvws9AaEfJjtDfXNg3Rnc6kPk2fAyE\ndhQFQlcvWbLE4u0VdDDzbfgYCO0o2T8ohNAiQGhngNDZBkI7CoTONhDaUWSFrn7nyWmTJk1/\ndl29R3wQWgQI7QxyQi88vyA6n0TJsHoe3AShRYDQziAj9O5+FE9g9GGr3hBaBAjtDBJCHzmT\nqPGwhxauXLt2xbz7Lorsq6+y6p5zQo82TsD/Xm/aErqiqb5u/y8cLJgVEkL/jejmvbXN7YOI\nXrfonnNCFw8RV0jP66g3bQndq6F+eTX/VQcLZoWE0BdQ4sXsmt6WP1nJPaEXiTBDSmhjhpMi\nCC2JhNAtaXHigvnUzqI7hBYBQjuDhND59Enigg1UaNEdQosAoZ1BQugyWpm44G1qZtEdQosA\noZ1BQujuyZONjqGeFt0htAiLisS8uZW9IXRWkRB6MhUkfImeH6L7LbpDaBGmRk/b62eaIXSW\nkBC6sjkFh87dekR7feTL2f2Jyn+w6A6hRZhA+nwcjc8QTQidJWSuFC4v1fY0ha07tG8lroC3\n3WjVG0KLMMEY6bK0hM5rUS64womCWSF1L8fnFaHaK99Fv/nOsjOEFsGe0HTxTI2Kjk4UzArJ\nu+2qnrttQK8uXc8YOG5BfdMmQWgRbAo9MSEHSBvcD60cCO0mEFo5ENpNILRyILSbQGjlQGg3\ngdDKgdBuAqGVA6HdBEIrB0K7CYRWDoR2EwitHAjtJhBaORDaTSC0ciC0m0Bo5UBoN4HQyoHQ\nbgKhlQOh3QRCK8cQ+tG2SwRtHxXNOkLPF+/2GyaaEFoVEFo5htDDo7+AGC6aSUK/FH33RNGE\n0KqA0MoxhK4o0JsFFSIkCf0ibRKxQwcRILQqILRyILSbQGjlQGg3gdDKgdBuAqGVA6HdBEIr\nB0K7CYRWDoR2EwitnOwL/d5MndlqC2cBhFZO9oU+t1TME9aGqtRWzgEIrZzsC500xzSIA0Ir\nB0K7CYRWDoR2EwitHAjtJhBaORDaTSC0crIg9E3ddIbqTQhtDoRWThaE7jhwisaI1E/SAnFA\n6IypuWOkzgt6OxtCWz4aDsQBoTNmF/UbqtFJf74VhHYVCJ0xu2i9iOMhtAeA0BkDob0EhM4Y\nCO0lIHTGQGgvAaEzBkJ7CQidMVkUenR09o7Rogmh6wdCZ0wWha4I6bMvhfQcELp+IHTGZFPo\nxBwQun4gdMZAaC8Boe3zw/g7dJaKJoT2EhDaPuvpEnGt+zh9KlwI7SWyI/Sh+XNiTGYotC7S\n4D7iD5xF60TTcaHX0Sx9iDcr+9P8T3aE/uKE8hit6EBWtuEeUaFPyBd/YDuaL5qOCz2f2okC\njrpG2Z/mf4xhPuuZfdnaAsevHLrQSTI6LnQ0h/EUAKBhDDNR6eg12dkChBYBQjuDMcwNtMtR\nnR/OxkEGhBYBQjuDMcx7nr4gP6J04S8X16jeQs4I3T+/iSDwkGhCaFeoPSj8bkbvYMTpYyd+\noXYLOSN0rxL9nENqGSG0MySc5dg6rXtE6eC5z6s8L5E7QlvKCKGdIfm03abJHSNOl/1W3alN\nCC0ChHaGuueh3+2jHSDmXfm1oi1AaBEgtDMkCf3pXe0jNgd6NY3spd9QswUILQKEdoZ4offO\n6qXtnJvfvim8d2oxlXysZAsQWgQI7Qy1Qi+/pkTbOfed86NobjiaLleyBQgtAoR2BmOYt03R\njgXpmDtqDwafouZKtgChRYDQzmAMc562cz73xR/j3tlOISVbgNAiQGhniN3L0WJs0pm6b+lo\nJVuA0CJAaGcwhrnf3EO1yw5sER/YocNKtgChRYDQzmAMc/vT45Z9SD0VbgFCiwChnSH6leOY\nuGVbqInCLUBoESC0M6QQuuaPVKBwCxBaBAjtDJFhvq9p06YUaBojn+gkhVuA0CJAaGeIDPO0\nICXzuMItQGgRILQzaMP8/f9OoILzYlw8+k2VW4DQIkBoZ0h1UKgWCC0ChHYGY5hv/UPWtgCh\nRcii0MOvqBTsUfC3+R7MnGQfrwndyTjyCWLCGQgtg9eE7nDsZo33jSnJchuqqJgXCRVJKNwC\nhBYhm0LrOXZB6LB2VxJN0EMCCrcAoUWA0M4Aoe0DoT0MvkPbB0J7GAhtHwjtYSC0fSC0h4kJ\nveFPIjzRrVXPvyrdAoQWQYXQLzfQH4r189tS5YDQGlGhHw021sJN4pBwjMotQGgRVAg9OXrM\n3jtVDgitYQzz+jxqtC8cXkVUdnE50XsKtwChRVAhtHUOCK1hDNENdEplJFRQs6/DB8+jUQq3\nAKFFgNDOYAzRT+jlyH8PNaa7I2EZnaJwCxBaBAjtDMYQFdP3YfGNY2MkfE9HKdwCI6Fv1x+D\n1Ia2iSaE9iDGEAWD2n//SK20cIjyFG6BkdB9L5ypMUq9jBBaFcYQldDuyH/Ppl9pjW+pWOEW\nOAmd8vwvhPYSse/Qi8PhbSF6VmsspRMVbgFCiwChnSF2lqPzh19fQA21/XR4CA1TuAUILQKE\ndgZjiD7Qf/l9a+TltqtInPJQBYQWAUI7Q+xKYSjic5/9kVdLiPofUbgFCC0ChHaG2L0cn06b\nOF/Mzrizwa/3q9wChBYBQjtDirvtqtVuAUKLAKGdIdPbRzt3ru/H8xBaBAjtDJkKTVRVTw8I\nLQKEdoaY0MvGDu53Toz014fQENpLGEN0eFj6P5I9Jh6i5lqw6A6hRYDQzmAM0cOaxcW1U+pa\nrpICi+4QWgQI7QzGEPWg1vPSPbuRR3m/vDEK0fVasOgOoUWA0M4QuznplbRXWdaOyp6MrY/v\n0BDaSxhDVEiV6a9TOZho2G5jfQgNob2EMUTH0l47a80sovJ39PUhNIT2EsYQXU2rbK228WcU\nule74wNCQ2hPYQzRR8UX2rshaf9NRH23QegwhPYW0RNur5Vdusnemi83paNfhdAQ2lsYQzR+\nzCUUOOHs2IOD0ll161kUuBNCQ2hPYQyR1HS6NfdoN1HngND/6tFN0GioaEJoDxM9bde4rGk8\naa696vicEHpRwRRBUdZkhNCqyPBuu+9fWvBjPV38LPRfRgrON7zKnowQWhWZ3j5aP34WumO3\noRonZ11GCK0KCG1FxxkiZF9GCK2KOKFrPvr7bJlfx24ZNcpqckcI7VAOCK0RE3rbjY2J6FA4\nfHCgzXFZl+KsyIY1MWZBaGdyQGiNqIrvNROn6yJCv06hxbZSpBB6UyD+HOABJZW6AYT2HcYQ\n7WlBwYumCqGXtKRG39hJUb1kyZI6yypjvI49tDM5ILSGMURTqeka7fJKROjw18fSFIVbwHdo\nh3JAaA1jiM6kx8JRocNPUE+FW4DQDuWA0BrGEDWlHeGY0NupjcItQGiHckBoDWOIQmJGaEPo\nI8GCeterfufJaZMmTX92Xb1HfBDaoRwQWsMYotJ8TWVD6G+pvns5Fp5fED2FUTJstXVfPwr9\n79H6Ne/SyaLpD6F30hC96rlyfzQPjCH6OS0Nx4SeTqdZrrO7X8KdeYHRh616+1HoGUeJS95D\ng07JqCZHT1F0p75yfzQPjCEaR92qo0IvLaEJVqscOZOo8bCHFq5cu3bFvPsuiuyrr7Lq7kuh\nO+rRQRnV5RgPoSPHgcV00tObib79dM7QAJXssFrlb0Q3x/2kdvsgotctukNoh3NAaI15ebXf\nIfIWWq5yAQ1KaNf0psssukNoh3NAaMGyn0Z97rzSepWWlHRpfD61s+gOoR3OAaENVtwzfMCl\nV099t75V8umTxAUbqNCiO4R2OAeEtkkZJe3C36ZmFt0htMM5ILRNutPYxAVjLC+VQ2iHc0Bo\nnf2fr37/q3RWmUwFCV+i54fofovuENrhHBA6woGHzszXDggbDXq53h+tVDan4NC5W0W/I1/O\n7k9U/oNFdwjtcA4IHQ6/0bb2pF2Pj+pbZ3mp1q+wdYf2rcQV8LYbrXr7SuhH7xCc11JvQmjf\nIYbopcKIlh0GXX957ybaXvrN+lb6vCJU+w+g6DffWXb2ldDFXfpqtHBLRgidKdoQfdaIqEI/\nFXfwpW5EpZvrXa3qudsG9OrS9YyB4xbsq6erv4ReJIJrMqrIMeaUOYJ5/v3pWwZoQzSAaGZs\nwY+3EJ2vcgsQ2uEcZwSaCAJWNySwJTJEnwXp5vhFlxKtVbgFCO1SjqJX0/2rOUHa7wmLEp5H\nsaOAble4BQjtUo6cFXogXZ64bGju/qYQQvudyBB1oIcTlz1meSnbLhDapRw5K3RZ8iPdFlFI\n4RYgtEs5clboPFqWuGxVmhOepweEdilHzgpNyU/AgtCuywihZYHQCXASumCIfh1/drp/PAsg\ndAKchKZycRm/U8d0/3gWCKEfXpLAwxDadRnV5ZiRg0LXReEWILS7OSA0hPaMjBDaNpEhuiMF\nCrcAod3NkXtCZxkI7W4OCK0YCO1uDgitGAjtbg4IrRgI7W4OCK0YCO1uDgitGAjtbo6o0Aum\n6Dxd7yD4GgidAGOhT+okroR3Kal3EHwNhE6AodCPtNVvaGj7iGguKq5/FPwMhE6AodBXRq/+\nXimaEDpTILSnciwqMp7wW++Ub/4EQifAX+ip0R32A9Yj4VcgdAL8hZ5AawTd77IeCb8CoRPI\nBaH15jkQWhII7ckcEFoWCO3JHBBaFgjtyRwQWhYI7ckcEFoWCO3JHBBaFgjtyRwQWhYI7ckc\nEFoWbwp9xLgAHH1+V5XeLILQPidXhR5jXAAOfSGaDyZeEPa4jBDanFwVesRF4gLwUlovmuN7\n6BeE/SEjhDYnZ4UeIcKuqNDGo9D8ISOENgdCiwihuQChRYTQXIDQIkJoLkBoESE0FyC0iBCa\nCxBaRAjNBQgtIoTmAoQWEUJzITtCbw7FPw4AQnswB4S2w5EVtY8gmg6hvZgDQsvi6a8c2+nk\nbhoty/XF/pARQpuT40Jvol+JKTmP85WMENqcnBfajzJCaHMgtIj+khFCmwOhRfSXjAqFPvTU\nTJ2NYR5AaBH9JaNCoTfQceUaR10d5gGEFtFfMioUej3tih+P8Kf6L3fW7Az7FAgtor9kzJ7Q\ne4LGxbD+YZ8CoUX0l4zZE3oXLRO/fr/VuNDkPyC0iP6SMZtCJ1459R8QWkR/yQihzYHQIvpL\nRghtDoQW0V8yQmhzck3oH8aMFJzYTzT9KSOENifXhF5PlwzVKPGzjBDanNwTWv8AfS0jhDYH\nQovoLxkhtDkQWkR/yQihzYHQIvpLRghtDoQW0V8yQmhzILSI/pIRQpsDoUX0l4wQ2hwILaK/\nZITQ5kBoEf0lI4Q2h73QO++8Q2edaEJoI0fn/xCjcg19J5oQOm1cFnpRSFzqHto8cY/kaxmV\n5GgmhqWnkQNCp43bQhfr8RwIbZUDQqcNhPZDDgidNh4Rus9NmzVehdApcyQJPe4XYrQ2fxH2\nGzkjdLPo5L5fiaZXRPJKjiShz4qO1pthn5EzQpedKnY5j3pMJK/kSBK6V6m+h274athn5I7Q\n3hTJKzmShTZyFEHoOkBoP+SA0GkDof2Q48xWfTV607KEHBC6LhDaFzmaiAuHo5JyQOi6QGgf\n54DQdYHQPs4BoeuSdaF3vDhHsCDldiA0hFZL1oWekN9EEHg91bsQGkKrJetCR+87SD34EBpC\nqwVC+zgHhK4LhPZxDghdF8eEDpbo36UHJ7wLoSG0WhwTmoaKkx0jOya8C6EhtFqcE1r/AGdA\naHU5ckfo6neenDZp0vRn1x2oryeE9nGOXBF64fkF0TvAS4attu7rtND3Hi0eRz/lyjtFGFGo\nv+svkbySI0no+frQTnkq7FlkhN7dj+IJjD5s1dtpoSuC4r6xvsGWIpQH9Hf9JZJXciQJ3bGT\nGNMuJWHPIiH0kTOJGg97aOHKtWtXzLvvosi++iqr7o4L7bYEnHIkCz1DhOjXOC8iIfTfiG7e\nW9vcPogo5UVnAwjt4xw5IfQFNCihXdObLrPoDqF9nKPoBfFg2Urj0J+n0C1pceKC+dTOojuE\n9nGOgHGc1ENv8hQ6nz5JXLCBCi26Q2g/5xi5RmOccSqUp9BltDJxwdvUzKI7hPZ/jhmshe5O\nYxMXjKGeFt0htP9z8BZ6MhUkfImeH6L7LbpnTejD1+nzinbqpLc9JgGnHCZC/7f+AQy9zvJC\nhLNICF3ZnIJD5249or0+8uXs/kTlP1h0z5rQu2iIeMpxC29KwCmHidB9u4oPYIgxX6AnkLlS\nuLxUO/AtbN2hfStxBbztRqve6oV+X7/++ntaI5oelYBTjinG7QRHTxHNFwr0ZofRohmd0dX4\nXKZMd/NX0TJChz+vCNVe+S76zXeWndULPaJpN42f0lzR9KgErHIExIh3C+g5JpLezNNzxB5r\noX8unenDsHtICR0OVz1324BeXbqeMXDcgn31dE0t9OolOvpcoOENRvPflqmW6536XS5aXpcg\nZ3LUeU7LzITP1pqvjE++nlvc0kZSaBukFHp39Iz9+Xq71Gj+1CrTxuj/FE4UTV9LwClHktBf\nUOJna8350VvcdqfTu35cEjp5pvjiRSIk3csc3qNfeK3SW5h935s51tNm8TFdMSIhx/je+qe3\nJ5zAj5UJV9P7jjdyKDqw9LTQK6L/1qeJJoT2Zo7Xoh/TTxJynBtdnHgd7nRjaRu96UGht4wa\nNcri7QyEXtRQXHhd0/0u0YTQXs2xUHxMbZNyNNI/vQZJd+yNE0unR0/9eU/odZF/bkmLvjq5\nPEbrVEJXURvxZnFDvVOwmQiNC/Rmo1YitCC9GWoiQhs6VsT8fBHa0TEiNgjqvUjvVRLQm4ES\nEZoYOYINRDiG2iGHwzlaiNCqkd4saCxCM+PdhsUitKGqzE3UyI7QB5+YGePPd6Za6Xn9zXvH\n6fHmR0R48Fa9eeuDIjx8s94cd68Ij43WmxMn6nH0Y8jBJcfM5zMXUaBA6OolS5ZkngUAFWT/\noBAAB4HQgBUQGrBCVuj0J5oBwEHkhLYz0QwADiIjtL2JZgBwEAmhbU40A4CDSAhtc6IZABxE\nQmibE80A4CASQtucaEaGUgL+53jFVqSFhNA2J5qRob1+R5YKpjdUlmpNj2uVpZpNS5Xluvhi\nZamW0mxlua49U7EVaSEhtM2JZmQwfl2sApVzSBi3OqpA2e2S4dgvRVQQva1XBdFbg51FQmib\nE83IAKHtAKHjkBDa5kQzMkBoO0DoOLI/0YwMENoOEDqO7E80IwOEtgOEjiP7E83IAKHtAKHj\nyP5EMzJAaDtA6Di8eT80hLYDhI4DQtsAQtsBQtcCoe0AoePwptCnzFKWanETZanCF/5BWaqP\nA3vq75QmI0cqS7Un8LGyXH+4UFkqG3hT6C3qZuCt+VxZqvA3e+vvky6b1aWqrFSXS2FZe79R\nlyt9vCk0AJJAaMAKCA1YAaEBKyA0YAWEBqyA0IAVEBqwAkIDVkBowAqvCL2+K1GdWR9/GzfJ\nw6o0E20Z161xqKznhO3Jb1TdfVqzglYXPnMk/aJMcsmUtfyak4pDTXve9bWCskxyyZSlsyyP\njlFQV+pU8mVJ4Q2hf5yYTymEvsb+WDxdbPQvTZqebOUxxht90r4ryCyX/bIODI/9vue5TMsy\nzSUxWjqV7aiO0BLDZZJKuiw5PCH0us5UcGoKoQdS61VRqtPK9HKACm95be2CS4gafh7/xpdl\nVHD7in+9PojogjSrMs1lv6yriTo+8tb7C68IUmhthmWZ5rJflsFgqrNblajLLJV0WXJ4Qehv\n8+mUD+5PIXQf6movU81xVLBcvPod0Zj4d66koL6bvZnolQxz2S5rY4BO038YH/krr8+sLPNc\ntssymEllv0y20H5dpqlky5LEC0JvCf7uQDiV0J3pbHuZVhFdq7/aE6DT496oyqfB+qvqRmnu\nc0xz2S9rzonNjYdPVgepV2ZlmeayX5bOxiJ68sYkCyXqMkslW5YsXhC68q1wOKXQ7aLDmi47\nXvnrBuNlS+NBvTqzieYaL4dRQVr/5zPNZb+sWg43oPMyK8s0l2RZB35Gl4eTLZSrK2WqTEZL\nBi8ILUgl9FF0XWSYtlbaPNLWOFhAl8Q1byPaarx8iOidjHJlUtbjRI+qKispl2RZN9Hxu+tY\nKFdXylSZjJYMXhb6cOS761+65RE1uWyN3XSziOKfTTqACqMDupDoyYxyyZZVs/3ta/No0CEV\nZaXIJVfWyxSK+JpsoVRdqVNl8iHK4GWhvyPKN8735D1oL9tnZdQzfpfQi1pHX64kmpZRLsmy\ntPNg+ecknmmTLStVLqmytjaje8N1LZSpyyRVBh+iFF4W+t+RQbj0/74/+PH4AqI5dpJ9chy1\n3BK/oCu1j75cRzQpo1ySZWkSNhvwTI2KslLlkimr5iw6R0uSbKFEXWap5D9EObws9I4//WmR\n/ur1ILWrqbOKKa+VUtvEXy93pQ7Rl2uJJtuoq24uybLeXjz/3u5Ev6hSUFaqXDJl3UNHb9Ni\nXaFt12WWSvpDlMTLQsdxuZ1DpqlBOnVr4qIzqFX05Qqi6emXlSKXbFkaf6aEB9RIl5Uil0RZ\nq0IB/RRgsoX26zJNJVFWRvhE6L8SPZ1mor2XEQ3fn7RwIBVEvwbPI3om3aJS5pIrK1YIfZZx\nWSlz2S9r9/F0p/4q2ULbdZmnsl9WZvhE6BeIHk8vT1UPCj1SZ+k4oi+Nl/cRra3zvp1cUmVF\niexW52daVupc9st6iRLonEFd5qnsl5UZ3hY6NpH6VKKX00pTeSo1/kfdxQuIZhsvL6KSNJ9P\nbpLLfllj+3SJbvJhogUZlWWay35ZC0wttF2XeSqZDzETvCz0jY2DnxovT6PAznSyHOxBZe+n\nWP5DEfXXX31TQJenV5FZLvtlXUNkXK7WvkfGzR4nUZZpLvtl/VgV5XpqXlUVd1ud7brMU0l8\niBnhSaGr163TTpTNIjpX3zs8YHr8k8QYCr2RsMBIpd1kI/ZmNUMo8G56FZnlsl/Wm0Tl+qzw\nb4TouCMZlWWaS2K0YsS++GYwXCapMilLBi8I/bF2Z+GviZ6I3mG4hGhUJBzqRnTK/7zz0cJB\nRM1SH/4k8UmILojdq7hqX22q8M42FPrtsn8u6GO0M8hlvyxtX9rs7tfee/HqoHEuVr4s01wS\nZcWIWZhBXSapMilLBi8IfV781y/tHvDosG4/Pbq0wwdpZbon4ZvchrhU4X8eZyytOGSdo/5c\ntssKH6iIrlH4Z7FAvizzXPbLilFHaIm6zFJlUJYMnhY6XPPCoHaFha36P57mYdxkc6HDe6ee\nVlbQZshr6VZlkctuWRGWXdupJNSkx1hjHyVflkUuibIM6gotU5dJKvmyZPCC0AAoA0IDVkBo\nwAoIDVgBoQErIDRgBYQGrIDQgBUQGrACQgNWQGjACggNWAGhASsgNGAFhAasgNCAFRAasAJC\nA1ZAaMAKCA1YAaEBKyA0YAWEBqyA0IAVEBqwAkIDVkBowAoIDVgBoQErIDRgBYQGrIDQgBUQ\n2i2+bESd9dnx97aj5jtcroYNENo1niKaIl6MIVroci18gNDuMZgaboqENXl0ndul8AFCu8fO\nFnROOHz451Re7XYpfIDQLmtMsRgAAAEmSURBVPIq0ZPhByi43O1CGAGh3WQkNX2vmMa6XQYn\nILSb7G1PhdTloNtlcAJCu8pbROTE4yhzBwjtKuMjQk92uwhWQGg3eS9Ep1L++26XwQkI7SL7\nOtIvDpxCJ+93uxBGQGgX+TU1+CS8Oo9ucbsQRkBo9/hHgO6LhN9R4B9ul8IHCO0au9vRqYcj\ncd8J1PZ7t4thA4R2jeGU/6F48WaArnS5Fj5AaLdYQDTReHkD0VxXa2EEhAasgNCAFRAasAJC\nA1ZAaMAKCA1YAaEBKyA0YAWEBqyA0IAVEBqwAkIDVkBowAoIDVgBoQErIDRgBYQGrIDQgBUQ\nGrACQgNWQGjACggNWAGhASsgNGAFhAasgNCAFRAasAJCA1ZAaMAKCA1YAaEBK/4fp4SRlvRW\nEKEAAAAASUVORK5CYII=",
-      "text/plain": [
-       "Plot with title “”"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    }
-   ],
+   "outputs": [],
    "source": [
     "# Define a seed for repeatability\n",
     "set.seed(12354)\n",