Probability of a girl birth given placenta previa (BDA3 p. 37).
Simulate samples from Beta(438,544), draw a histogram with quantiles, and do the same for a transformed variable.
ggplot2 is used for plotting, tidyr for manipulating data frames
library(ggplot2)
theme_set(theme_minimal())
library(tidyr)
Sample from posterior Beta(438,544). Obtain all draws at once and store them in vector ‘theta’
a <- 438
b <- 544
theta <- rbeta(10000, a, b)
Compute odds ratio for all draws
phi <- (1 - theta) / theta
Compute 2.5% and 97.5% quantile approximation using samples
quantiles <- c(0.025, 0.975)
thetaq <- quantile(theta, quantiles)
phiq <- quantile(phi, quantiles)
Histogram plots with 30 bins for theta and phi
# merge the data into one data frame for plotting
df1 <- data.frame(phi,theta) %>% pivot_longer(everything())
# merge quantiles into one data frame for plotting
df2 <- data.frame(phi = phiq, theta = thetaq) %>% pivot_longer(everything())
ggplot() +
geom_histogram(data = df1, aes(value), bins = 30) +
geom_vline(data = df2, aes(xintercept = value), linetype = 'dotted') +
facet_wrap(~name, ncol = 1, scales = 'free_x') +
labs(x = '', y = '') +
scale_y_continuous(breaks = NULL)
![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABUAAAAPACAMAAADDuCPrAAACtVBMVEUAAAABAQEMDAwODg4RERESEhIYGBgaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj5AQEBBQUFERERFRUVGRkZHR0dISEhJSUlKSkpLS0tNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////cW+DcAAAACXBIWXMAAB2HAAAdhwGP5fFlAAAgAElEQVR4nO3d/4MddXno8dQbeqiRLvZar7eltyve20IL4XKLky+ENU2IFIyu1YLYiiYoCv1CVyuXUrFNFG1VhC5YMbVKuyiWglxxVcqX0CKrGBSNGwIRyAYSMn/HzdlPzubsZHZmzjPPM59ndt+vH9ycLyQzn88zbzfZzGZZCgAQWRb7AACgrQgoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqAAIERA0S6f67zv2IO/6Xww3pEABBQtQ0DhCAFFu8wL6Heu/3/xjgQgoGiHgweP/mBeQIG4CCjcW9/56VWv6gy/5c7ugyMB3XfF8EmnXfnMkQcf47fwiIqAwr31nQs6nVNe3vmFj6TdgL7tt1++8qxOZ+1LBBSxEVC4t77TOefx9Lm/6HTu7wa0c8bjafr1kzrfIaCIjYDCvfWdV/+s+/HyzgXdgJ74UPfBaGecgCI2Agr31nf+ZPbjE52TXjoS0HNmH1zV+SwBRWwEFO6t73wu/OA1nSePBPTi2R//CQFFfAQU7q3v/Fv4wcrON+f+GhMBhQMEFO6t7/xD+MEpnccJKDwhoHBvfefq2Y8/6aw4SEDhCQGFe+s7r3mu+/EDnfNSAgpPCCjcW9/prP5Ruv8vT+zcR0DhCgGFe+s7b+l0hld0TvxwSkDhCgGFe+s7993xpv/xaxd9pfuAgMIRAgr31nd/6w44REDhHgGFVwQU7hFQeEVA4R4BhVcEFO594n2Pxz4EIBcBBQAhAgoAQgQUAIQIKAAIEVAAECKgACBEQAFAiIACgBABBQChugE9ND29T+VA7O3ZH/sIqjk4Pf2z2MdQ0fRM7COo5sXp6WdjH0M1h6cPxD6Eal6Ynn4u9jFU89L0C4Y/OwF1h4CqI6DqCGhAQN0hoOoIqDoCGhBQRyYm7kkJqIF2BfSBiQn/R0tAAwLqyM8te1VKQA20K6Drly27P/ZxlCKgAQF1hIAaIaDqCGhAQB156+j7UwJqoF0B3T46+kTs4yhFQAMC6g4BVdeugLYCAQ0IqDsEVB0BVUdAAwLqDgFVR0DVEdCAgLpDQNURUHUENCCg7hBQdQRUHQENCKg7BFQdAVVHQAMC6g4BVUdA1RHQgIC6Q0DVEVB1BDQgoI6cPLQyJaAG2hXQS4aGdsY+jlIENCCgjizaWzk3VGJ3nO0KKLdyqiKgSvwH9ITlp6YE1EC7Avqm5csfjH0cpQhoQEDdIaDq2hXQViCgAQF1h4CqI6DqCGhAQN0hoOoIqDoCGhBQdwioOgKqjoAGBNQdAqqOgKojoAEBdYeAqiOg6ghoQEDdIaDqCKg6AhoQUHcIqDoCqo6ABgTUka1bPpwSUAPtCuiNW7Y8Gfs4ShHQgIA6wq2cRtoVUG7lVEVAlRBQdQRUGwFVR0CV+A/o+PiXUwJqoF0BvXd8/OnYx1GKgAYE1B0Cqq5dAW0FAhoQUHcIqDoCqo6ABgTUHQKqjoCqI6ABAXWHgKojoOoIaEBA3SGg6gioOgIaEFB3CKg6AqqOgAYE1B0Cqo6AqiOgAQF1ZGpqV0pADbQroLunpiwveR0ENCCgjnAnkpF2BZQ7kVQRUCUEVB0B1UZA1RFQJf4DetbKC1MCaqBdAf2zlSunYh9HKQIaEFB3CKi6dgW0FQhoQEDdIaDqCKg6AhoQUHcIqDoCqo6ABgTUnSUb0IoEx0lA1RHQgIC606aAqpaRgBJQAwRUCQFVR0C1EVB1BFQJAVVHQLURUHUEVAkBVUdAtRFQdQRUif+Atu8v0hNQZfxFenUEVIn/gLbvVk4CqoxbOdURUCUEVB0B1UZA1RFQJf4D2r5vZ0dAlfHt7NQRUCX+AxoQ0CUe0FYgoAEBdYeAElD/CGhAQN0hoATUPwIaEFB3CCgB9Y+ABgTUHQJKQP0joAEBdcdFQGOUsSLB2RBQdQQ0IKDuENBigrMhoOoIaEBAHRkf/3JKQJX0n027Anrv+PjTsY+jFAENCKgjju5Eip2/+vrPpl0B5U4kVQRUCQEdQOz81dd/NgRUHQENCKgjY2PbUwKqpP9s2hXQ28bGdsc+jlIENCCg7hBQFf1n066AtgIBDQioOwRURf/ZEFB1BDQgoO4QUBX9Z0NA1RHQgIC6Q0BV9J8NAVVHQAMC6g4BVdF/NgRUHQENCKg7BFRF/9kQUHUENCCg7hBQFf1nQ0DVEdCAgLpDQFX0nw0BVUdAAwLqCHciKeo/m3YFlDuRVBFQJQR0ALHzV1//2RBQdQQ0IKCOnDy0MiWgSvrPpl0BvWRoaGfs4yhFQAMC6g4BVdF/Nu0KaCsQ0ICAukNAVfSfDQFVR0ADAuoOAVXRfzYEVB0BDQioOwRURf/ZEFB1BDQgoO4QUBX9Z0NA1RHQgIC6Q0BV9J8NAVVHQAMC6g4BVdF/NgRUHQENCKg7BFRF/9kQUHUENCCgjrx19P0pAVXSfzbtCuj20dEnYh9HKQIaEFBHuJVTUf/ZtCug3MqpioAqIaADiJ2/+vrPhoCqI6ABAXVkYuKelIAq6T+bdgX0gYkJ/0dLQAMC6g4BVdF/Nu0KaCsQ0ICAukNAVfSfDQFVR0ADAuoOAVXRfzYEVB0BDQioOwRURf/ZEFB1BDQgoO4QUBX9Z0NA1RHQgIC6Q0BV9J8NAVVHQAMC6g4BVdF/NgRUHQENCKgjk5MPpQRUSf/ZtCugU5OT7keVgB5FQB3hTiRF/WfTroByJ5IqAqqEgA4gdv7q6z8bAqqOgAYE1BG+G5Oi/rNpV0D5bkyqCKgS/wENCKiK/rNpV0BbgYAGBNQdAqqi/2wIqDoCGhBQdwioiv6zIaDqCGhAQN0hoCr6z4aAqiOgAQF1h4Cq6D8bAqqOgAYE1B0CqqL/bAioOgIaEFB3CKiK/rMhoOoIaEBA3SGgTYq90BkEVB0BVeI/oCcPrUwJaLNiL3RGN6CXDA3tjH0cpQhoQEAd4VbO5sVe6Axu5VRHQJUQ0AHE7lpjYi90BgFVR0CV+A9oQECbFHuhM/gzUHUEVAkBHUDsrjUm9kJnEFB1BFQJAR1A7K41JvZCZxBQdQRUCQEdQOyuNSb2QmcQUHUEVAkBHUDsrjUm9kJnEFB1BFQJAR1A7K41JvZCZxBQdQRUCQEdQOyuNSb2QmcQUHUEVAkBHUDsrjUm9kJnEFB1BFSJ/4COjW1PCWizYi90Rjegt42N7Y59HKUIaEBAHeFOpObFXugM7kRSR0CVENABxO5aY2IvdAYBVUdAlfgP6LZtf58S0GbFXuiMbkDv2LZtT+zjKEVAAwLqDgFtUuyFzuCLSOoIqBICOoDYXWtM7IXOIKDqCKgSAjqA2F1rTOyFziCg6gioEgI6gNhda0zshc4goOoIqBICOoDYXWtM7IXOIKDqCKgSAjqA2F1rTOyFziCg6gioEgI6gNhda0zshc4goOoIqBL/Ad27t7uWBLRJsRc6oxvQ5/buPRT7OEoR0ICAOtLEnUixg+WN3UqLcCeSOgKqhIB2xQ6WN3YrLUJA1RFQJf4D+trhc1MC2iy7lRbpBnTr8PCjsY+jFAENCKg7BLRJdistwheR1BFQJQS0K3awvLFbaRECqo6AKiGgXbGD5Y3dSosQUHUEVAkB7YodLG/sVlqEgKojoEoIaFfsYHljt9IiBFQdAVVCQLtiB8sbu5UWIaDqCKgSAtoVO1je2K20CAFVR0CVENCu2MHyxm6lRQioOgKqxH9A3zByaUpAm2W30iLdgF47MvJ47OMoRUADAuoIt3I2z26lRbiVUx0BVUJAu2IHyxu7lRYhoOoIqBL/AZ2a2pUS0GbZrbRIN6C7p6YsL3kdBDQgoO4Q0CbZrbQIX0RSR0CVENCu2MHyxm6lRQioOgKqhIB2xQ6WN3YrLUJA1RFQJQS0K3awvLFbaRECqo6AKiGgXbGD5Y3dSosQUHUEVAkB7YodLG/sVlqEgKojoEoIaFfsYHljt9IiBFQdAVXiP6Dj419OCWiz7FZapBvQe8fHn459HKUIaEBAHeFOpObZrbQIdyKpI6BKCGhX7GB5Y7fSIgRUHQFV4j+gW7d8OCWgzbJbaZFuQG/csuXJ2MdRioAGBNQdAtoku5UW4YtI6gioEgLaFTtY3tittAgBVUdAlRDQrtjB8sZupUUIqDoCqoSAdsUOljd2Ky1CQNURUCUEtCt2sLyxW2kRAqqOgCohoF2xg+WN3UqLEFB1BFQJAe2KHax2stuPDAKqjoAqIaBdsVPUTnb7kUFA1RFQJf4DesLyU1MC6pHdfmR0A/qm5csfbOwXlCKgAQF1hFs5vbLbjwxu5VRHQJX4DyifgXpltx8ZfAaqjoAq8R/QgID6Y7cfGfwZqDoCqoSAdsVOUTvZ7UcGAVVHQJUQ0K7YKWonu/3IIKDqCKgSAtoVO0XtZLcfGQRUHQFVQkC7Yqeonez2I4OAqiOgSghoV+wUtZPdfmQQUHUEVAkB7Yqdonay248MAqqOgCohoF2xU9ROdvuRQUDVEVAl/gPKv4nkld1+ZPBvIqkjoEr8B5RbOb2y248MbuVUR0CVENCu2ClqJ7v9yCCg6gioEv8BnZi4JyWgHtntR0Y3oA9MTDzb2C8oRUADAuoOAfXHbj8y+CKSOgKqhIB2xU5RO9ntRwYBVUdAlRDQrtgpaie7/cggoOoIqBIC2hU7Re1ktx8ZBFQdAVWy2AMaOzKLmfYeL4iAqiOgSggopLT3eEEEVB0BVUJAIaW9xwsioOoIqBL/AZ2cfCgloB4pb/TCugGdmpx0P6oE9CgC6ki9O5FiR2YxU97ohXEnkjoCqoSAQkp5oxdGQNURUCX+A/qGkUtTAuqR8kYvrBvQa0dGHm/sF5QioAEBdYeA+qO9xwvii0jqCKgSAgop7T1eEAFVR0CVEFBIae/xggioOgKqhIBCSnuPF0RA1RFQJQQUUtp7vCACqo6AKiGgkNLe4wURUHUEVAkBhZT2Hi+IgKojoEoIKKS093hBBFQdAVXiP6CvHT43JaAeKW/0wroB3To8/Ghjv6AUAQ0IqCPcyumV8kYvjFs51RFQJQQUUsobvTACqo6AKvEf0L17u2tJQP1R3uiFdQP63N69hxr7BaUIaEBA3SGg/mjv8YL4IpI6AqqEgEJKe48XREDVEVAlBBRS2nu8IAKqjoAqIaCQ0t7jBRFQdQRUCQGFlPYeL4iAqiOgSggopLT3eEEEVB0BVUJAIaW9xwsioOoIqBL/Ad227e9TAuqR8kYvrBvQO7Zt29PYLyhFQAMC6gh3InmlvNEL404kdQRUCQGFlPJGL4yAqiOgSvwHdGxse0pAPVLe6IV1A3rb2Njuxn5BKQIaEFB3CKg/2nu8IL6IpI6AKiGgkNLe4wURUHUEVAkBhZT2Hi+IgKojoEoIKKS093hBBFQdAVVCQCGlvccLIqDqCKgSAgpj9TefgKojoEoIKIzV33wCqo6AKiGgMFZ/8wmoOgKqxH9AuROp5eqPAHciqSOgSggojNUfAQKqjoAq8R/Q1w6fmxLQFqs/At2Abh0efrT+z2SMgAYE1B0C2lr1N58/A1VHQJUQUBirv/kEVB0BVUJAYaz+5hNQdQRUCQGFsfqbT0DVEVAlBBTG6m8+AVVHQJUQUBirv/kEVB0BVUJAYaz+5hNQdQRUCQGFsfqbT0DVEVAl/gP6hpFLUwLaYvVHoBvQa0dGHq//MxkjoAEBdYRbOVuu/ghwK6c6AqqEgMJY/REgoOoIqBL/AZ2cfCgloC1WfwS6AZ2anHQ/qgT0KALqDgFtrfqbzxeR1BFQJQQUxupvPgFVR0CVEFAYq7/5BFQdAVVCQGGs/uYTUHUEVAkBhbH6m09A1RFQJa0NaOwsoKr6m09A1RFQJQQUxupvPgFVR0CV+A/oxMQ9KQFtsfoj0A3oAxMTz9b/mYwR0ICAOrLAnUixs4Cq6o8AdyKpI6BKCCiM1R8BAqqOgCrxH9C3jr4/JaAtVn8EugHdPjr6RP2fyRgBDQioOwS0tepvPl9EUkdAlRBQGKu/+QRUHQFVQkBhrP7mE1B1BFQJAYWx+ptPQNURUCUEFMbqbz4BVUdAlRBQGKu/+QRUHQFVQkBhrP7mE1B1BFQJAYWx+ptPQNURUCX+A3ry0MqUgLZY/RHoBvSSoaGd9X8mYwQ0IKCOcCtny9UfAW7lVEdAlfgP6AnLT00JaIvVH4FuQN+0fPmD9X8mYwQ0IKDuENDWqr/5/BmoOgKqhIDCWP3NJ6DqCKgSAgpj9TefgKojoEoIKHwo2HwCqo6AKiGg8KFg8wmoOgKqhIDCh4LNJ6DqCKgSAgofCjafgKojoEoIKHwo2HwCqo6AKvEf0K1bPpwS0MWvYAS6Ab1xy5YnTedMAwENCKgj3Mq5NBSMALdyqiOgSggofCgYAQKqjoAq8R/Q8fEvpwR08SsYgW5A7x0ff9p0zjQQ0ICAukNAF7uCzeeLSOoIqBICCh8KNp+AqiOgSggofCjYfAKqjoAqIaDwoWDzCag6AqqEgMKHgs0noOoIqBICCh8KNp+AqiOgSggofCjYfAKqjoAq8R/QqaldKQFd/ApGoBvQ3VNTlpe8DgIaEFBHuBNpaSgYAe5EUkdAlRBQ+FAwAgRUHQFV4j+gZ628MCWgi1/BCHQD+mcrV06ZzpkGAhoQUHcI6GJXsPl8EUkdAVVCQOFDweYTUHUEVAkBhQ8Fm09A1RFQJf4CGvtCRhwFI0FA1RFQJQQUPhSMBAFVR0CVEFD4UDASBFQdAVVCQOFDwUgQUHUEVAkBhQ8FI0FA1RFQJf4D+vM/PxTjekbDCkaCv0ivjoAq8R/QZcv+S4zrGQ0rGAlu5VRHQJUQULTCimXLfqf7MfZ4FiKgAQGNKHvlJMnq5i9XeLM2SdZ3P8Yez0IENCCgEcW+UuFa7PEsREADAhpR7CsUrsUez0IENCCgEcW+QuFa7PEsREADAhpR7CsUrsUez0IENCCgEcW+QuFa7PEsREADAhpR7CsUrsUez0IENCCgEcW+QuFa7PEsREADAhpR9oo57bQzYlyo8OWs005b1/0YezwLEdCAgEaUvXK4EwkbuBNJGQFVQkDRCgRUFQFV4j+gp5zyuuYvV3jzW6ecsrb7MfZ4FiKgAQGNKPaVCtdij2chAhoQ0IhiX6FwLfZ4FiKgAQGNKPYVCtdij2chAhoQ0IhiX6FwLfZ4FiKgAQGNKPYVCtdij2chAhoQ0IhiX6FwLfZ4FiKgAQGNKPYVCtdij2chAhoQ0IhiX6FwLfZ4FiKgAQGNKHvFcCcSNnAnkjICqoSAohUIqCoCqqTBgAqvnJe97OcVr0O01H992cte3/3Y2LhKENCAgFqIfQli6WhsqOcjoAEBtRD7osLS0dhQz0dAAwJqIfZFhaWjsaGej4AGBNRC7IsKS0djQz0fAQ0IqIXYFxWWjsaGej4CGhBQC7EvKiwdjQ31fAQ0IKAWYl9UWDoaG+r5CGhAQC3EvqiwdDQ21PMR0ICAWhBeC6985X9XvbbQSq975StXV393Y0M9HwENCKgF4ZXDrZzYcOxWzmoaG+r5CGhAQC0IrxwCig0EVBkBVeI/oGeeeZbwv8Qi8jtnnjlS/d2NDfV8BDQgoBbsLi5gvsaGej4CGhBQC7EvKiwdjQ31fAQ0IKAWYl9UWDoaG+r5CGhAQC3EvqiADO0RJ6ABAbUQ+2oBMrRHnIAGBNRC7KsFyNAecQIaENDBxL4OABGF2Z+HgAYEdDCmQ3722eeY/vxoheTss98Q69eufCEQ0ICADsZ0erkTCRsGvRNJV+ULgYAGBHQwptNLQLGBgCojoEr8B5TvxoQNg343Jl2VLwQCGhDQwcQabKAJlS8EAhoQ0MHEHnDAAQLaQ0AHE3tyAQcIaA8BPSr2SALtQUB7Fn9AY88asETFvvQDAlpP7CkCoED1elfLy+IP6JwG74Wv5eD09M9iH0NF0zOxj6CaF6enn419DNUcnj4Q+xCq4bfwAQF1h4CqI6DqCGhAQB05eWhlSkANtCuglwwN7Yx9HKUIaEBAHfm5Za9KCaiBdgV0/bJl98c+jlIENCCgjhBQIwRUHQENCKg7BFRduwLaCgQ0IKDuEFB1BFQdAQ0IqDsEVB0BVUdAAwLqDgFVR0DVEdCAgLpDQNURUHUENCCg7hBQdQRUHQENCKg7BFQdAVVHQAMC6g4BVUdA1RHQgIA6Mja2PSWgBtoV0NvGxnbHPo5SBDQgoI5wJ5KRdgWUO5FUEVAlBFQdAdVGQNX5DujhmRnLw9M0czD2EZTZtu3vj/zvSzMzL8Y+kopmDsU+gmras6SHjyzpHdu27Yl9HKWOLKn76yk4PPOS4c9eN6AAsGQRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqAAIERAAUCIgAKAkCigMx/bPLL15vnfiefxa96+bvQTLfkOVy59ZeTHlZ5DZcct37euevPIJdf7/15HfmWX9ObkqB9GOqDIJAGdfkeSbEqSy/u/a+WX1ibrNifJZmZTbGtyfCzznkNl2eX7ZJKsfcua5LxHIh3PIpBd0r8koAN7b7J1Ov3R25Nrjz31/TVrvvhS+tSVyZ+qHdkS8/wNyXGxzHsOlR23fP+5as0dL6bPfTAZbcn3MXXn+Incmjw6M+twpEOKTBDQR5IN3c89f7Jm1U/nnvvj5O+6H1548yr//5qLRxPvPi/Jjmbec6gsZ/muT27ufphZnzwe55haLm8iL0xa8q8SGBEE9IbkmtmPlyc75p7bdPQz+JtCSDGgj51//vmrMqOZ9xwqy1m+9yTh9+7vSu6Ockhtl7OkL6x6Y6SDcUIQ0Pcmd85+vOXY7+GfS5Lwr7HemVyhclhL0fk5scx7DpVllu8bdz3f/XB4Y/JQnONZBDJLuit5d6wj8UEQ0AuTB2c/3plcNvfc+mTX7Mebk7eoHNZSREDV5S7fzEeTi5foH9gpyCzpfcnVt1488rYPTcU6ntgEAR05+idI30zePvfce5OPdz/MvCl5o8phLUUEVF3O8o29dU3ygekYB7M4ZJZ0R5Ika38vSdZ8IdYBRSYI6JokfKHo4WTz3HM7V63+9E+e//d3Jsn5Wke25BBQdTnLd8X5SXKJ/3812K3Mkm5PNtx7MH3u48nq78Y6orgEAX3j3Geg7zz25G1ru38XbO32vqhiMARUXe7y7f/cqjU7mz+WRSKzpN/7Vvhs/i+Sq6McTnSCgF569I/g70z+uO/Z73/s8vf93a67kj/UOa4liICqW2D5/iYZa/pIFo0FlvSBpfqZkyCgVyZ3zX78h+T64177XPLhuke0ZBFQdfOX72f3fDv84OvJO6IczmKwwETuTtYuzS/Mif4e6F/Nfrwy+ae553743WdnP77naFwxOAKqbv7y7U5GwkXOX7aTm7+kz+/4YljSh5NL4xxPbIKAPphc0P1Ln/vOXXfsW4d8Itne/TCVXMBNclIEVN385Tt8fvKd2R98KLkhzvEsApkl3ZRMzv7g+uSv4xxPbJJ74S9LPnQoPXB5cl33wRd3dL/8NrV6zV3p4f+4ILlV9/CWkrnRDEs6/zlIZJb0U8nmh9L0hRsTvseVWGZJb0o23Z+mM59ZtfHpqIcVjSSgT25K1r97XfLW2e/GtC4085YkOX9jknx4af5BiIq50Vx37P+GCGgtmSV9cWuSrH/z6uQ8/pxJLLOkhz6RJBs3r0o2TkY9qnhE3w/0qY9cuHbzDbO3xc1d7d9678bfvZzBrIGAqssu6aEd79ow8ofX8/fo5Y6b0vs+sHHkXdueiXhIUfEd6QFAiIACgBABBQAhAgoAQgQUAIQIKAAIEVAAECKgACBEQAFAiIACgBABBQAhAgoAQgQUAIQIKAAIEVAAECKgACBEQAFAiIACgBABBQAhAgoAQgQUAIQIKAAIEVAAECKgACBEQAFAiIACgBABBQAhAgoAQgQUAIQIKAAIEVAAECKgACBEQAFAiIDCqcPjnz/yv5/rvK/wdSAmAgqnDnV+OS0KaHgdiImAwikCCv8IKJwioPCPgMKnt3SO+KVuQPddMXzSaVc+M/vs9Ptf/8rfeMcjx15P0zsv/K1fPPX8f4l6sFiqCCh8+tQfdV5x2XuPBPRtv/3ylWd1OmtfOvLkN36ts+I3XtF5xW1zr6eXdzq/+tuv6nQ+GfuAsRQRUDjV+y1854zH0/TrJ3W+k6Yzv965eiZ98doTT/pe7/XvdIa+mqYHb+gMxz5eLEUEFE71AnriQ91Ho53xNL2uc8nsS5d13tN7/cZf2jr71K90no11oFjCCCic6gX0nNlHV3U+m6arO5OzD77ROT3zRaTdJ3V+FuEYsdQRUDjVC+jFs4/+pBvQX+384ZauSzqvOBbQw4/s+Mg7f71DQBEBAYVT8/8aUzegL3aOOdgL6PZTjjx69UWvIqCIgIDCqeMDmr7mxBezr3+086o/v/upNP2fBBQREFA4lRPQ13cemX2w7+s7e68Pd/519qlhAooICCicygnon3feOfvg/Z1req+f1Jn9G/aPn0hAEQEBhVOHOiv2ZwK67791rngmPfSJE0/+Qe/1szofPfLqXf+r09kV93CxJBFQePXqzm9tmh/Q9I5f7vzC6a/uvPyLc6/f3umcPnJqZ+Sczmm3xjxYLE0EFF594XWvOD0T0PSJy/73yaf9waPHXk//dd1rfm3TZw/d8zu/eH28Q8VSRUABQIiAAoAQAQUAIQIKAEIEFACECCgACBFQABAioAAgREABQIiAAoAQAQUAIQIKAEIEFACECCgACDUZ0P0vlr/HiYP79x+MfQyVHdwf+wiq27+/TQfbphlo0cG+OBP7CKrbv7/4YJsM6FPPN/iL1TMzPd2eC33/nthHUN2e6adiH8W88SMAABvOSURBVEJ1bZqB6en2VOn5Ns3A9N7C1wloLgJqhIDaIKBGCKgIATVCQG0QUCMEVMRTQD+1cuVXC99AQI34mYGez69cuSP3BQJqhICKeAroVcuWfbbwDQTUiJ8Z6Llu2bKP5r5AQI0QUBECaoSA1kJAm0ZARTwFdO/U1LOFbyCgRvzMQM8zU1P7cl8goEYIqIingJYioEbaNAME1AYBFSGgRgioDQJqhICKEFAjBNQGATVCQEUIqBECaoOAGiGgIgTUCAG1QUCNEFARAmqEgNogoEYIqAgBNUJAbRBQIwRUxFNAHx4f/0HhGwioET8z0PPo+PhjuS8QUCMEVMRTQLkTKRY/M9DDnUhNI6AiBNQIAa2FgDaNgIp4CujdY2M7C99AQI34mYGeybGxb+e+QECNEFARTwEtRUCNtGkGCKgNAipCQI0QUBsE1AgBFSGgRgioDQJqhICKEFAjBNQGATVCQEUIqBECaoOAGiGgIgTUCAG1QUCNEFARAmqEgNogoEYIqAgBNUJAbRBQIwRUxFNAuRMpFj8z0MOdSE0joCIE1EhRQDdU09zB+pmBHgLaNAIq4img1w0N7Sh8AwE14mcGev52aOjTuS8QUCMEVMRTQEsRUCNtmgECaoOAihBQIwTUBgE1QkBFCKgRAmqDgBohoCIE1AgBtUFAjRBQEQJqhIDaIKBGCKgIATVCQG0QUCMEVISAGiGgNgioEQIqQkCNEFAbBNQIARXxFNB/HB29r/ANBNSInxno+cro6NdyXyCgRgioiKeAcisnAe3hVs6mEVARAmpEIaDNddbPDPQQ0KYRUBFPAZ2amPhx4RsI6JIJ6K6JiSdyXyCgRgioiKeAliKgSyagCyKgRgioCAFVpZtGAno8AmqEgIoQUFUE1BoBNUJARQioKgJqjYAaIaAiBFQVAbVGQI0QUBECqoqAWiOgRgioCAFVRUCtEVAjBFTEU0CfnJwsHjgCumQCunty8qe5LxBQIwRUxFNAF8GdSARUCXciNY2AihBQVQRUCQFtGgEV8RTQW0ZG7i58AwFdMgH90sjIP+e+QECNEFARTwEtRUCXTEAXRECNEFARAqqKgFojoEYIqAgBVUVArRFQIwRUhICqihHQiooOu00zQEBtEFARAqoqdiULFB12m2aAgNogoCIEVFXsShYoOuw2zQABtUFARQioqtiVLFB02G2aAQJqg4CKEFBVsStZoOiw2zQDBNQGARXxFNDtw8O3F76BgNZQdNh+ZqDnpuHhW3JfIKBGCKiIp4ByK6elosP2MwM93MrZNAIqQkBVxa5kgaLD9jMDPQS0aQRUxFNASxHQGooOu00zQEBtEFARAqoqdiULFB12m2aAgNogoCIEVFXsShYoOuw2zQABtUFARQhoNbHzV1/R2bVpBgioDQIqQkCriZ2/+orOrk0zQEBtEFARAlpN7PzVV3R2bZoBAmqDgIoQ0Gpi56++orNr0wwQUBsEVISAVhM7f/UVnV2bZoCA2iCgIp4CevfY2M7CNxDQGorOzs8M9EyOjX079wUCaoSAingKqOc7kWLnr76is/MzAz3cidQ0AipCQKuJnb/6is7Ozwz0ENCmEVARTwGd3LbtscI3ENAais7Ozwz0PLht28O5LxBQIwRUxFNASxHQGorOrk0zQEBtEFARAlpN7PzVV3R2bZoBAmqDgIoQ0Gpi56++orNr0wwQUBsEVISAVhM7f/UVnV2bZoCA2iCgIgS0mtj5q6/o7No0AwTUBgEVIaDVxM5ffUVn16YZIKA2CKgIAa0mdv7qKzq7Ns0AAbVBQEU8BXRm794XCt9AQGsoOjs/M9BzYO/eA7kvEFAjBFTEU0C5E8lS0dn5mYEe7kRqGgEVIaDVxM5ffUVn52cGegho0wioiKeAbh8evr3wDQS0hqKz8zMDPTcND9+S+wIBNUJARTwFtBQBraHo7No0AwTUBgEVIaDVxM5ffUVn16YZIKA2CKgIAa0mdv7qKzq7Ns0AAbVBQEUIaDWx81df0dm1aQYIqA0CKkJAq4mdv/qKzq5NM0BAbRBQEQJaTez81Vd0dm2aAQJqg4CKENBqYuevvqKza9MMEFAbBFSEgFYTO3/1FZ1dm2aAgNogoCKeAnrLyMjdhW8goDUUnZ2fGej50sjIP+e+QECNEFARTwHlVk5LRWfnZwZ6uJWzaQRUhIBWEzt/9RWdnZ8Z6CGgTSOgIp4Cundq6tnCNxDQGorOzs8M9DwzNbUv9wUCaoSAingKaCkCWkPR2bVpBgioDQIqQkCriZ2/+orOrk0zQEBtEFARAlpN7PzVV3R2bZoBAmqDgIoQ0Gpi56++orNr0wwQUBsEVISAVhM7f/UVnV2bZoCA2iCgIgS0mtj5q6/o7No0AwTUBgEVIaDVxM5ffUVn16YZIKA2CKiIp4A+PD7+g8I3ENAais7Ozwz0PDo+/ljuCwTUCAEV8RRQ7kSKL9oCZ3AnUtMIqAgBrSZ22JoSbYEzCGjTCKiIp4DesWXL/YVvIKDmoi1wxj1bttyb+wIBNUJARTwFtBQBNRdtgasioEYIqAgBrSZ22JoSbYGrIqBGCKgIAa0mdtiaEm2BqyKgRgioCAGtJnbYmhJtgasioEYIqAgBrSZ22JoSbYGrIqBGCKgIAa0mdtiaEm2BqyKgRgioCAGtJnbYmhJtgasioEYIqAgBrSZ22JoSbYGrIqBGCKiIp4B+cPnyWwvfQEDNRVvgjG3Ll3889wUCaoSAingKKLdyxhdtgTO4lbNpBFTEU0D5DDS+aAucwWegTSOgIp4CWoqAmou2wFURUCMEVISAVhM7bE2JtsBVEVAjBFSEgFYTO2xNibbAVRFQIwRUZMkHNHaxnNFfYGUE1AgBFSGg6Ke/wMoIqBECKkJA0U9/gZURUCMEVISAop/+AisjoEYIqAgBRT/9BVZGQI0QUBFPAY3ybyLFLpYz+gssw7+J1DQCKuIpoFFu5YxdLGf0F1iGWzmbRkBFCCj66S+wDAFtGgEV8RTQh8fHf1D4BgJqTn+BZR4dH38s9wUCaoSAingKaCkCak5/gZURUCMEVISAop/+AisjoEYIqAgBRT/9BVZGQI0QUBECin76C6yMgBohoCIEFP30F1gZATVCQEUIKPrpL7AyAmqEgIoQUPTTX2BlBNQIARXxFNC9U1PPFr6BgJrTX2CZZ6am9uW+QECNEFARTwHlTqT49BdYhjuRmkZARQgo+ukvsAwBbRoBFfEU0FtGRu4ufAMBNae/wDJfGhn559wXCKgRAiriKaClCKg5/QVWRkCNEFARAop++gusjIAaIaAiBBT99BdYGQE1QkBFCCj66S+wMgJqhICKEFD0019gZQTUCAEVIaDop7/AygioEQIqQkDRT3+BlRFQIwRUhICin/4CKyOgRgioiKeAbh8evr3wDQTUnP4Cy9w0PHxL7gsE1AgBFfEUUG7ljE9/gWW4lbNpBFSEgKKf/gLLENCmEVARTwGd2bv3hcI3EFBz+gssc2Dv3gO5LxBQIwRUxFNASxFQc/oLrIyAGiGgIgQU/fQXWBkBNUJARQgo+ukvsDICaoSAihBQ9NNfYGUE1AgBFSGg6Ke/wMoIqBECKkJA0U9/gZURUCMEVISAop/+AisjoEYIqIingE5u2/ZY4RsIqDn9BZZ5cNu2h3NfIKBGCKiIp4ByJ1J8+gssw51ITSOgIgQU/fQXWIaANo2AingK6N1jYzsL30BAzekvsMzk2Ni3c18goEYIqIingJYioOb0F1gZATVCQEUIKPrpL7AyAmqEgIoQUPTTX2BlBNQIARUhoOinv8DKCKgRAipCQNFPf4GVEVAjBFSEgKKf/gIrI6BGCKgIAUU//QVWRkCNEFARAop++gusjIAaIaAingLKnUjx6S+wDHciNY2AihBQCOjvQwYBbRoBFfEU0O3Dw7cXvoGAeqG/Dxk3DQ/fkvsCATVCQEU8BbQUAfVCfx+qIqBGCKgIAYWA/j5URUCNEFARAgoB/X2oioAaIaAiBBQC+vtQFQE1QkBFCCgE9PehKgJqhICKEFAI6O9DVQTUCAEVIaAQ0N+HqgioEQIqQkAhoL8PVRFQIwRUxFNA/3F09L7CNxBQL/T3IeMro6Nfy32BgBohoCKeAqp7K2fsxixq9Ta6Am7lbBoBFSGgEKi30RUQ0KYRUBFPAX1ycrJ44AioF/U2uoLdk5M/zX2BgBohoCKeAlqKgHphtsWlCKgRAipCQCFgtsWlCKgRAipCQCFgtsWlCKgRAipCQCFgtsWlCKgRAipCQCFgtsWlCKgRAipCQCFgtsWlCKgRAipCQCFgtsWlCKgRAiriKaBTExM/LnwDAfWi3kZXsGti4oncFwioEQIq4img3InUGvU2ugLuRGoaARUhoBCot9EVENCmEVARTwHV/W5MsRuzqNXb6Ar4bkxNI6AingJaioB6YbbFpQioEQIqQkAhYLbFpQioEQIqQkAhYLbFpQioEQIqQkAhYLbFpQioEQIqQkBhx2IGCKgNAipCQGHHYgYIqA0CKkJAYcdiBgioDQIqQkBhx2IGCKgNAiriKaDXDQ3tKHwDAW0Z+Sz87dDQp/NngIDaIKAingLKrZyLjXwWuJWzaQRUhIDCjnwWCGjTCKiIp4CWIqAtYzEDBNQGARUhoLBjMQME1AYBFSGgsGMxAwTUBgEVIaCwYzEDBNQGARUhoLBjMQME1AYBFSGgsGMxAwTUBgEVIaCwYzEDBNQGARUhoLBjMQME1AYBFfEU0LvHxnYWvoGAtox8FibHxr6dPwME1AYBFfEUUO5EWmzks8CdSE0joCIEFHbks0BAm0ZARTwF9OHx8R8UvoGAtox8Fh4dH38sfwYIqA0CKuIpoKUIaMtYzAABtUFARQgo7FjMAAG1QUBFCCjsWMwAAbVBQEUIKOxYzAABtUFARQgo7FjMAAG1QUBFCCjsWMwAAbVBQEUIKOxYzAABtUFARTwFdO/U1LOFbyCgLSOfhWempvblzwABtUFARTwFlDuRFhv5LHAnUtMIqAgBhR35LBDQphFQEU8B/dTKlV8tfAMBbRn5LHx+5cod+TNAQG0QUBFPAS1FQFvGYgYIqA0CKkJAYcdiBgioDQIqQkBhx2IGCKgNAipCQGHHYgYIqA0CKkJAYcdiBgioDQIqQkBhx2IGCKgNAipCQGHHYgYIqA0CKkJAYcdiBgioDQIq4imgVf8ifewsoCr5LPAX6ZtGQEU8BbTqrZyxs4Cq5LPArZxNI6AiBBR25LNAQJtGQEU8BbTqt7OLnQVUJZ8Fvp1d0wioiKeAliKgLWMxAwTUBgEVIaCwYzEDBNQGARUhoLBjMQME1AYBFSGgsGMxAwTUBgEVIaCwYzEDBNQGARUhoLBjMQME1AYBFSGgsGMxAwTUBgEV8RTQh8fHf1D4BgLaMvJZeHR8/LH8GSCgNgioiKeAcifSYiOfBe5EahoBFSGgsCOfBQLaNAIq4imgd2zZcn/hGwhoy8hn4Z4tW+7NnwECaoOAingKaCkC2jIWM0BAbRBQEQIKOxYzQEBtEFARAgo7FjNAQG0QUBECCjsWM0BAbRBQEQKK6AaaAQJqg4CKEFBEN9AMEFAbBFSEgCK6gWaAgNogoCIEFNENNAME1AYBFfEU0A8uX35r4RsI6OKUs9Xbli//eP4MEFAbBFTEU0C5lXOJytlqbuVsGgEV8RTQ64aGdhS+gYAuTjlb/bdDQ5/OnwECaoOAingKaCkCujgNNAME1AYBFSGgiG6gGSCgNgioCAFFdAPNAAG1QUBFCCiiG2gGCKgNAipCQBHdQDNAQG0QUBECiugGmgECaoOAihBQRDfQDBBQGwRUhIAiuoFmgIDaIKAingL6j6Oj9xW+gYAuTjlb/ZXR0a/lzwABtUFARTwF9Kply34z9rWMCHJmgVs5m0ZARQgoosuZBQLaNAIq4imgU2eeuTb2tYwIcmZh18TEE7lDQkCNEFARTwHlDzeXqEFmhIAaIaAiBBTRDTIjBNQIARUhoIhukBkhoEYIqAgBRXSDzAgBNUJARQgoohtkRgioEQIqQkAR3SAzQkCNEFARAoroBpkRAmqEgIp4CuiTZ5+9Lva1jAhyZmH35ORPc4eEgBohoCKeAsqdSEtUzixwJ1LTCKgIAYVHw8uWvTbzVBgSAmqEgIp4CugtK1asjHK1wp0zVqw4M/NUGBICaoSAingKKF9EQoEwIwTUCAEVIaBoiTAjBNQIARUhoGiJMCME1AgBFSGgaIkwIwTUCAEVIaBoiTAjBNQIARUhoGiJMCME1AgBFSGgaIkwIwTUCAEVIaBoiTAjBNQIARXxFNDtJ5xwRuyrFD785gknnJ55KgwJATVCQEU8BZRbOdHDrZxNI6AiBBQeEdCmEVARTwHlz0BRIMwIATVCQEUIKFoizAgBNUJARQgoWiLMCAE1QkBFCChaIswIATVCQEUIKFoizAgBNUJARZoJaOxrD4tAGCUCaoSAihBQtEQYJQJqhICKEFC0RBglAmqEgIp4CujKU055ve01iLY4+5RT/k/mqTBKBNQIARXxFNBf4U4kHMWdSE0joCIEFB4R0KYRUBFPAT371FMT9SsRrXTOqaeek3kqjBIBNUJARTwFFCgQRomAGiGgIgQULRFGiYAaIaAiBBQtEUaJgBohoCIEFIuL/TRXR0BlCGhG7IsKS4f9NFdHQGUIaEbsiwpLh/00V0dAZQhoRuyLCkuH/TRXR0BlCGhGtdF/w7p1620vLrTFevEs2E9zdQRUhoBmVBt97kRCT86dSBXZT3N1BFSGgGZUG30Cih4C2jQCKuIpoK874YQzZBcNFpvfPOGE02X/pf00V0dAZQhohu7FBSzMfpqrI6AySyegsa8WIEPt0lBAQGUIKBCJ2qWhgIDKEFDAN7UrqAgBlSGggG9qV1ARAipDQAHf1K6gIgRUhoACvqldQUUIqMxiCGjsAQcsNXJ1EVAZAipy+ooVKyP8snDojBUrzjT9BRq5ugioDAEV4VZO9Mhv5ayokauLgMoQUBECih4C2jQCKuIpoGvOPntdhF8WDq09++xzTX+BRq4uAipDQAHfGrm6CKgMAQUWhXpXFwGVIaDAolDv6iKgMq4DGnskgfaod60RUBntgMaeIgD16VZBma+Axt4qAC1VsTKqP5m3gDb4i9XTzHekV7J/T+wjqG7PdJt++9aiGZienol9DJXxW3gZAioyNTHx48I3EFAjfmagZ9fExBO5LxBQIwRUxFNAr1q27LOFbyCgRvzMQM91y5Z9NPcFAmqEgIoQUCMEtBYC2jQCKuIpoHds2XJ/4RsIqBE/M9Bzz5Yt9+a+QECNEFARTwEtRUCNtGkGCKgNAipCQI0QUBsE1AgBFSGgRgioDQJqhICKEFAjBNQGATVCQEUIqBECaoOAGiGgIgTUCAG1QUCNEFARAmqEgNogoEYIqAgBNUJAbRBQIwRUxFNAP7h8+a2FbyCgRvzMQM+25cs/nvsCATVCQEU8BZRbOWPxMwM93MrZNE8B3etuHhd0YM8eN/NY+hnoTIvmce+e4nl05Sk3M9Cz8Gege/YcaPhY5J5v0Qzs2fN04etNBhQAFhUCCgBCBBQAhAgoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqBHzHxs88jWmw8d9/z3z72i+A0oVLquxz1AJaUL+8Nr3n7exR+ZbviwWi9nXW9Ojvph99H3r33TuW/50519rxPQNJ1+R5JsSpLLf5Z5/oWLkysK34BCpet63ANUUrqw96xJ1l60JjnvoeaPrc3y1vUv+wP6b2uSZPO6ZNWNx14noGn63mTrdPqjtyfXZp7fnhydx4XegEKl63rcA1RStrB71yeffTGd+WhyUXu+R6gHeeu6NXl0ZtbhNH1mXfKhp9JDn1+16t/nXieg6SPJhu7/5fxkzaqfznv+W6s2hXlc6A0oVLquxz1AJaULe2fy7u6HlzYnfAo6gNx1vTA59l21b0neM/vb+xuSq+eeI6BHluOa2Y+XJzv6n953wQVfCvO4wBtQrHRdsw9QTenCfib5q9lnrkxub/rY2ixvXV9Y9cZjb/jz5EuzHx9JLpp7joAe+cT9ztmPt8z/LdHVyTfvCfO4wBtQrHRdsw9QTenCTiQXH/kNZ/rCBcn9jR9ci+Wt667wyXzwpxc9MPvxseSCuecI6JFP0h+c/Xhnclnfs7cnf50encf8N6BE6bpmHqCi0oU98I7kmh++8L0rkysOxzi+tspb1/uSq2+9eORtH5rqf+ONyR/P/ZiApiPJ47Mfv5m8/diTPxr5/QO9ecx9A8qUruv8B6iqfGGfvXL2C8f/l68hDSJvXXccWca1v5cka75w7H33n5s8MPeAgKZrkt2zHx9ONs89d+iP1vxn2pvHvDegVOm6znuAykoXNv2XNyQX/NHGZNM3ohxfW+Wt6/Zkw70H0+c+nqz+7tFnDvzd6uTmY/8RAU3fOPd/PO+ce+7TSfefYj86j3lvQKnSdZ33AJWVLuzXko3fPPLhq+vWPJD33yNf3rp+71vhboS/6H3h/e6LkpH+L80R0PTSo3/Z485jf7Lx2Op3df++wtF5zHkDypWt6/xFRmWlA7s5+drss7cmWyIcXmsVXecPhM9K940lyQd3979AQNMrk7tmP/5Dcn3vqbuSY76b9waUK1vX+YuMysoW9pkkeX722V3JOr6KVF3Rdb47WXtkKXdfmPzBI/NfIKDpDXN/a+6fek99c3TWBcm60dHH896AcmXrOn+RUVnZwr6wOgk3I34veWOUA2ypnHV9fscXw/8FPZxceuTR7yfXvZj5jwho+mByQferlfvOXfdc5pWjvyNa+A0oULquOQ9QQenCXpx8efbhZ/gzp0HkrOvhTcnk7A+uT/46Tb+Qs54ENE0vSz50KD1weXJd98EXdxz77WTv2p73BlRVuq7HP0AVZQv7b8nIV15KD35hzeqHIx1hO+Ws603JpvvTdOYzqzY+naZ/kPz7S0fN/TcENE2f3JSsf/e65K2zv+9Zl9w690Lv2p73BlRVuq7HP0AVpQv7yVXJujevSdbcFuf42ipnXQ99Ikk2bl6VbDzyieih1XN/zjw6998Q0COe+siFazffEP7gPfdC738DKitd1+MeoJLShf2Psbedd/E1u2IcW5vlret9H9g48q5tzxz50Y8TAgoAeggoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBCBBQAhAgoAAgRUAAQIqAAIERAAUCIgAKAEAEFACECCgBC/x+Z4sx/T0fSgwAAAABJRU5ErkJggg==)
LS0tCnRpdGxlOiAiQmF5ZXNpYW4gZGF0YSBhbmFseXNpcyBkZW1vIDIuMyIKYXV0aG9yOiAiQWtpIFZlaHRhcmksIE1hcmt1cyBQYWFzaW5pZW1pIgpkYXRlOiAiYHIgZm9ybWF0KFN5cy5EYXRlKCkpYCIKb3V0cHV0OgogIGh0bWxfZG9jdW1lbnQ6CiAgICB0aGVtZTogcmVhZGFibGUKICAgIGNvZGVfZG93bmxvYWQ6IHRydWUKLS0tCiMjIFByb2JhYmlsaXR5IG9mIGEgZ2lybCBiaXJ0aCBnaXZlbiBwbGFjZW50YSBwcmV2aWEgKEJEQTMgcC4gMzcpLgoKU2ltdWxhdGUgc2FtcGxlcyBmcm9tIEJldGEoNDM4LDU0NCksIGRyYXcgYSBoaXN0b2dyYW0gd2l0aApxdWFudGlsZXMsIGFuZCBkbyB0aGUgc2FtZSBmb3IgYSB0cmFuc2Zvcm1lZCB2YXJpYWJsZS4KCmdncGxvdDIgaXMgdXNlZCBmb3IgcGxvdHRpbmcsIHRpZHlyIGZvciBtYW5pcHVsYXRpbmcgZGF0YSBmcmFtZXMKCmBgYHtyIHNldHVwLCBtZXNzYWdlPUZBTFNFLCBlcnJvcj1GQUxTRSwgd2FybmluZz1GQUxTRX0KbGlicmFyeShnZ3Bsb3QyKQp0aGVtZV9zZXQodGhlbWVfbWluaW1hbCgpKQpsaWJyYXJ5KHRpZHlyKQpgYGAKClNhbXBsZSBmcm9tIHBvc3RlcmlvciBCZXRhKDQzOCw1NDQpLgpPYnRhaW4gYWxsIGRyYXdzIGF0IG9uY2UgYW5kIHN0b3JlIHRoZW0gaW4gdmVjdG9yICd0aGV0YScKCmBgYHtyIH0KYSA8LSA0MzgKYiA8LSA1NDQKdGhldGEgPC0gcmJldGEoMTAwMDAsIGEsIGIpCmBgYAoKQ29tcHV0ZSBvZGRzIHJhdGlvIGZvciBhbGwgZHJhd3MKCmBgYHtyIH0KcGhpIDwtICgxIC0gdGhldGEpIC8gdGhldGEKYGBgCgpDb21wdXRlIDIuNSUgYW5kIDk3LjUlIHF1YW50aWxlIGFwcHJveGltYXRpb24gdXNpbmcgc2FtcGxlcwoKYGBge3IgfQpxdWFudGlsZXMgPC0gYygwLjAyNSwgMC45NzUpCnRoZXRhcSA8LSBxdWFudGlsZSh0aGV0YSwgcXVhbnRpbGVzKQpwaGlxIDwtIHF1YW50aWxlKHBoaSwgcXVhbnRpbGVzKQpgYGAKCkhpc3RvZ3JhbSBwbG90cyB3aXRoIDMwIGJpbnMgZm9yIHRoZXRhIGFuZCBwaGkKCmBgYHtyIH0KIyBtZXJnZSB0aGUgZGF0YSBpbnRvIG9uZSBkYXRhIGZyYW1lIGZvciBwbG90dGluZwpkZjEgPC0gZGF0YS5mcmFtZShwaGksdGhldGEpICU+JSBwaXZvdF9sb25nZXIoZXZlcnl0aGluZygpKQojIG1lcmdlIHF1YW50aWxlcyBpbnRvIG9uZSBkYXRhIGZyYW1lIGZvciBwbG90dGluZwpkZjIgPC0gZGF0YS5mcmFtZShwaGkgPSBwaGlxLCB0aGV0YSA9IHRoZXRhcSkgJT4lIHBpdm90X2xvbmdlcihldmVyeXRoaW5nKCkpCmdncGxvdCgpICsKICBnZW9tX2hpc3RvZ3JhbShkYXRhID0gZGYxLCBhZXModmFsdWUpLCBiaW5zID0gMzApICsKICBnZW9tX3ZsaW5lKGRhdGEgPSBkZjIsIGFlcyh4aW50ZXJjZXB0ID0gdmFsdWUpLCBsaW5ldHlwZSA9ICdkb3R0ZWQnKSArCiAgZmFjZXRfd3JhcCh+bmFtZSwgbmNvbCA9IDEsIHNjYWxlcyA9ICdmcmVlX3gnKSAgKwogIGxhYnMoeCA9ICcnLCB5ID0gJycpICsKICBzY2FsZV95X2NvbnRpbnVvdXMoYnJlYWtzID0gTlVMTCkKYGBgCgo=