From d66fa06536d0d7a6727c74e20ed2ef1b097768e9 Mon Sep 17 00:00:00 2001 From: Alex Gebben Work Date: Tue, 4 Nov 2025 17:12:36 -0700 Subject: [PATCH] Remove old data --- old/BVAR_Pop.r | 74 ------------------------------- old/Bayes_Pop.r | 35 --------------- old/County_Data.png | Bin 43686 -> 0 bytes old/Popultaion.R | 106 -------------------------------------------- 4 files changed, 215 deletions(-) delete mode 100644 old/BVAR_Pop.r delete mode 100644 old/Bayes_Pop.r delete mode 100644 old/County_Data.png delete mode 100644 old/Popultaion.R diff --git a/old/BVAR_Pop.r b/old/BVAR_Pop.r deleted file mode 100644 index a49ede8..0000000 --- a/old/BVAR_Pop.r +++ /dev/null @@ -1,74 +0,0 @@ -library(BVAR) -library(tidyverse) - -source("Scripts/Functions.r") -#source("Scripts/Load_Wyoming_Web_Data.r") -DATA_TO_GATHER <- list() -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c("WYPOP","WY_POP",TRUE) - -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c("WYNQGSP","WY_GDP",TRUE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c("MEHOINUSWYA646N","WY_MED_INCOME",TRUE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c("BUSAPPWNSAWY","WY_BUISNESS_APPLICATIONS",FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('ACTLISCOUWY','WY_HOUSES_FOR_SALE',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYRVAC','WY_RENTAL_VACANCY_RATE',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYBPPRIVSA','WY_PRIVATE_HOUSING',FALSE) -#New Private Housing Units Authorized by Building Permits for Wyoming -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('B03002006E056023','LN_FIVE_YEAR_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('GDPALL56023','LN_GDP',TRUE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYLINC3POP','LN_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('LAUCN560230000000005','LN_EMPLOYMENT',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('BPPRIV056023','LN_PRIVE_HOUSING',FALSE) - -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('ENU5602320510','LN_NUM_ESTABLISHMENTS',FALSE) - -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('GDPALL56041','UINTA_GDP',TRUE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYUINT1POP','UINTA_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYSUBL5POP','SUBLETTE_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYSWEE7POP','SWEETWATER_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYTETO9POP','TETON_POP',FALSE) -##Idaho Counties -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('IDBEAR7POP','BEAR_LAKE_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('IDCARI9POP','CARIBOU_POP',FALSE) - -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('IDBONN0POP','BONNEVILLE_POP',FALSE) -###US Population -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('POPTOTUSA647NWDB','US_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('CE16OV','US_EMP',FALSE) - - - - -for(x in 1:length(DATA_TO_GATHER)){ - CURRENT <- DATA_TO_GATHER[[x]] - - if(CURRENT[3]){C_DATA <- CPI_ADJUST(FRED_GET(CURRENT[1],CURRENT[2]))}else{C_DATA <- FRED_GET(CURRENT[1],CURRENT[2])} - if(x==1){RES <- C_DATA}else{RES <- RES %>% full_join(C_DATA)} - rm(CURRENT,C_DATA) -} - -DATA <- RES %>% mutate(US_POP=US_POP-WY_POP,WY_POP=WY_POP-LN_POP-UINTA_POP-SUBLETTE_POP-SWEETWATER_POP-TETON_POP) -colnames(DATA) -TS_DATA_ORIG <- DATA %>% select(YEAR,LN_POP,LN_EMPLOYMENT,US_EMP,US_POP,WY_POP,UINTA_POP,SUBLETTE_POP,SWEETWATER_POP,TETON_POP,BEAR_LAKE_POP,CARIBOU_POP,BONNEVILLE_POP) %>% - filter(!is.na(LN_POP),!is.na(LN_EMPLOYMENT)) %>% - arrange(YEAR) %>% select(-YEAR) -TS_DATA <- log(ts(TS_DATA_ORIG,start=c(1970),end=c(2024),frequency=1)) - -MOD <- bvar(TS_DATA,lags=2, n_draw=15000) -opt_irf <- bv_irf(horizon = 25, identification = TRUE) - -plot(irf(MOD,opt_irf,conf_bands = c(0.05, 0.1,0.15)),area=TRUE,vars_impulse = c("LN_EMP"),vars_response = c("WY_POP","LN_POP","UINTA_POP","SUBLETTE_POP","SWEETWATER_POP")) -plot(irf(MOD,opt_irf,conf_bands = c(0.05, 0.1,0.15)),area=TRUE,vars_impulse = c("LN_EMP"),vars_response = c("TETON_POP","BEAR_LAKE_POP","CARIBOU_POP","BONNEVILLE_POP")) - -plot(irf(MOD,opt_irf,conf_bands = c(0.05, 0.1,0.15)),area=TRUE,vars_impulse = c("LN_POP"),vars_response = c("WY_POP","LN_POP","UINTA_POP","SUBLETTE_POP","SWEETWATER_POP")) -plot(irf(MOD,opt_irf,conf_bands = c(0.05, 0.1,0.15)),area=TRUE,vars_impulse = c("LN_POP"),vars_response = c("TETON_POP","BEAR_LAKE_POP","CARIBOU_POP","BONNEVILLE_POP")) - -DATA2 <- RES %>% mutate(US_POP=US_POP-WY_POP,WY_POP=WY_POP-LN_POP-UINTA_POP-SUBLETTE_POP-SWEETWATER_POP-TETON_POP) -TS_DATA2 <- DATA2 %>% select(YEAR,LN_POP,US_POP,WY_POP,UINTA_POP,SUBLETTE_POP,SWEETWATER_POP,TETON_POP,BEAR_LAKE_POP,CARIBOU_POP,BONNEVILLE_POP) %>% - dplyr::filter(!is.na(LN_POP)) %>% - arrange(YEAR) %>% select(-YEAR) %>% ts %>% log -MOD2 <- bvar(TS_DATA2,lags=5, n_draw=15000) -plot(irf(MOD2,opt_irf,conf_bands = c(0.05, 0.1,0.15)),area=TRUE,vars_response = c("LN_POP")) -?plot.bvar_irf -plot(predict(MOD,horizon=25,conf_bands = c(0.05, 0.1,0.15)),area=TRUE,vars=c("LN_POP")) -exp(3.5)-exp(3) -acf(resid(MOD)) diff --git a/old/Bayes_Pop.r b/old/Bayes_Pop.r deleted file mode 100644 index 21e7a70..0000000 --- a/old/Bayes_Pop.r +++ /dev/null @@ -1,35 +0,0 @@ -library(tidyverse) -library(bayesPop) -library(bayesMig) - - -DAT <- read_tsv(us.mig.file) -DAT[,2] %>% t -sum(DAT$`2001`) - -# Toy simulation for US states -dir.create("Output") -sim.dir <- "./Output" -us.mig.file <- file.path(find.package("bayesMig"), "extdata", "USmigrates.txt") -m <- run.mig.mcmc(nr.chains = 2, iter = 100000, thin = 1, my.mig.file = us.mig.file, - output.dir = sim.dir, present.year = 2017, annual = TRUE) -pred <- mig.predict(sim.dir = sim.dir, burnin = 5, end.year = 2050) -mig.trajectories.plot(pred, "Wyoming", pi = 80, ylim = c(-0.03, 0.03)) - -load("./Output/bayesMig.mcmc.meta.rda") -############# -#subnat example -data.dir <- "./extdata" -dir.create("Pop_Output") -sim.dir <- "./Pop_Output" -dir.create("Output") -example(pop.predict.subnat) -read_tsv(file.path(data.dir, "CANlocations.txt")) -pred <- pop.predict.subnat(output.dir = sim.dir, locations = file.path(data.dir, "CANlocations.txt"), - inputs = list(popM = file.path(data.dir, "CANpopM.txt"),popF = file.path(data.dir, "CANpopF.txt")), - verbose = TRUE) -pop.pyramid(pred, "Ontario") -?pop.pyramid -summary(pred) -write.pop.trajectories(pred) -load("Pop_Output/predictions/prediction.rda") \ No newline at end of file diff --git a/old/County_Data.png b/old/County_Data.png deleted file mode 100644 index 15834d29ed7485afda773356b25386eb79a42218..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 43686 zcmdqJby!vHwl>ad0fH_R6+~cxln4@1Qj1hX2^9sU5$Tq$g$lToT7ZB^ij<0Uw~BPP zbV*6~qQ5cm?CjmNDjeo2HJm@!VgmGn~!XXh)&xf|NS*+h;||(Vjz;2mAdZ~H9O#{sks)t zzqIh~JOeW`GX<#(!!K#+pMl3snal4;-^-m=?V&bPqr{rEw0(H}$s$74Q7z@%363rc z`$Bf7O_Mj7cpd)-@jEmM`h5&-z5-WerTv(n-Q@r4`>*4Cwpfp-Cq9}uy@@DczVu0* z3Sry`pONHbbo0i=+AL1EH6qEg$5P3VKTyctWI%p)$Q=3$e!QnYO@jP&El?8q@ey|f zc_2w~W(wr5zZst+KklGNkHC*-F8_c3Ds}79SasrUdZ*#%chlUqrjtYKla!)At2^u8 zi~H<2ys2D%*bKP`tZ1Qh`IE zu`7*2R)ynXQbM=7J4*b+Ne0=!S_MfG>Gnq{O$wDQ`Ah6hD>G?jI4%sDr4>#ndu*xL z+2urd`zwWw4-`9c;P*H7_l7?2_Ur6#bqn=ddqzmKCg&9SQ*$FH;yq6GjLL1nbtt=U z&Y`$8IazK-y<>4_rMY;tI)q(iZ>D@NB`+^8F>$<6XsDr~;eLR}YO5S_$=sB?hm|p# z3Cm_JLqi^`Q?0~kw-$TT?e#gDb4e#XwmF>EdyM9oUitZT6p0W?`5!s-m5f>bs`dx_ zD;USE)#>epavf(e+lhwW-d@**QpZZ4v$m5>yAwj?o7+SCd%siCRhuv?4Sa%;ip96e zf5Boblx_@u<(`?DQI5sTRia0>=y={Kyng)}CBE7svA-fwjEx$w+!$?ctGkrnI<$A( zkEy+gU*O{+8S1M(*t<5*ENwpyjeRIwU_Bz8+`jpPQknBn7PFX`+kAha(14waib{LS z#* zx3O)(;Jv-Q$yg68Tn<(^n%|@?O4#-)A78A?Z}x7s?jXx}YJ1kA$aKfi@~E+KDLv^| zYC7%$;^N}M!aerio*yP&URg;mS+0Ba>QzEQf`(<`_<-XO-PTzD=3(^8l{+U zu?s^FGjye#oCq(DlBE^T7roVXsXc2a=bTJrtbFX*sRE1<8RLBU{@&fRFCRbN;qra; z%9_q?Dh?OoNp_C&VL5j9{XW6V+Z)_WclgchK(P=Pm*td4Qt+EMc(lZP(QE}uVy8~6 z?CNS(U|=9(Y}0lw`JdyuQZ*)v2;g=y+!=jyXr0H|JyjYGSZvu-)6UNi4#(zw{(MW= z^6EP=Ct>j9m6q7OLDz-Oj*fz|{zBX2nRnY;Tbi*692^{SMHkR+i_d5JZr{3fEcyy3 zr(N_eID&ZDS)JYAYQ)%x`)=bLI>L5Sfkqn^KZpERCY!l+%G``$y%KM~y=*mpGv6j* zlD2y$n!r7|x4W&?UB398EM8|;N>5&1zI>}Qomu|0bVIbLnlKpn?oeD@oJ{N_o5W`A zS2mkTRDX27WZ3v~l?LL%Q6G*fH(0iSBZqWl2X8SPmk(|*%&uLOC z-Me?wP?BnzHrHKTie~E4J$7g642+s%q9r`aW(cssVOKQHojW%nvfOXm?40>ouW)Z# z2sX4)R5&e5{5HdUN!0VP_uAM^gYIuXe?A-^etF`8CRO=LlW1E07@VJ}G^draX02N% zJt{)bqF2_^dpSrgIHAyX@}!%6x3+8QpGDZ7FK%9I@TKBzYHErRFbkWhfdJHg*K+J; zvv69s=#oy;bJ8cwlg+Uk1sRXB4GYhy-8o&efZY)7&j0c-b9BAWN+-|!-Aqlh$DZx@ z`Tm{5jw@H3mTTFxh^elq+qEkSUC}Ai++;Z={cC*Oe^2SJ7TWj>qK|q5MfR~Cd(Lul zaz;kUGly?!X+_Vxz|p&WOi8&-!V5uVW5{FQ%453^qRfK&cl+iTakTqt5XEB%PqA*x zuPMIk@`!=sRR_>*uf*>6x3xX(e4xgzw)cxo;yizd@x=!xZhulC2?`INnm{dE_Ao8Y zj#K+rf$d*+*3{CP+iF>dpt`R@ujx2opVk~O59ey8cm8#%yX#U~A#cpU@ke3h%c1Ex z2OqbkN-8A!hK^a4@2-E%+@oCWHA`nlA$bUo*_w)W9Gsq=wNo_MgnOC&X6v%k<+^UU zRg^JSjqfwh17AwX%6ba7eiAo^lh}!S`g?DCAw4zKeZFvVZ)Z}XnDtiSUfRIH;#u!L z-Rmpg?TeGaX*-XH78MmOFFze0o=Gp~s?W-|JVT_z@SJF|u%)GiO2^G+qT!wNtEs7} z1&{shx~0Lv!8dHmm+?9tyY}PX9Lu){Iy&x;FI~y+UtUZ15NAg|^6%~}MGo!U5VoI5 z%^$R>rmNG>?%4r9;n+XqGWvtGTVhbT1gs6yU%;0xIZ~M=H@Z}RO4|HqhLm&5z`z_$ z*}_)-oQ73tXCuvg*|ueOjH8W>&0)t(H^K-ml$d_V2SSzhYq$1fWe-@kii(PZRm+PU zM-m#i{o?pr9W2CEEyrfV>}X-z8_sTUSyW6{cDHDdV|%hUDI-|*cve={&d$!-Ot+AL z0G*iA9TGhrYV2w}x-f0CRvG&t$5JU4Q)<*qKQf3e*tYrx7fll*>@Ild+`Ze{)y3*n zFqhXi=sXEq<80KOrgaim!Qn{y3=c4hpy}r3CM)YTbGthn&FHnvJrCw>Wnu9k;EQJA z)8UQ1dB<`U!Mx(NE)9q+m_362I6+ZpfPXOyR?{Hh3x)J%b;oNpfzxU-UmvVYm1wNr z6vGNjze-6jbNTk-Xq~>nz87pfM@7o_@856GSTC^32hSHTU<_^*Rp(!1WE}C~NV9eq zHtXW-rrnsBn0Or`l+lDlWsRXNRk$^0=go4++ewLuS`Ix19y?=bajX(vXF+yQ@okXm z*aQu~=9r|3L`Fm)UbugS0Garc7-6+td7)?~ zP1tsVnT>6KN=`z=X7opvL48(sHt|WV!^1cN<RRRi(Q+q9@x3C2|3cPIsezP!*<7@;M03#3Au9F)=Yc-D--qAIuJNg}n1lZo0cg zTY|hcrzx8rNXoOZSwdpQINc`Pl(0oFfAcLLJYHE@xqK$ASTo*O3`>k_ z(sgj(26!jVVf7J(AutFDscp)+%LhYlu^#gmT3=sRRa0|a_o6Q!BrS{dxqEJKtsq0n zBj#|o&;@xlwYQNT!{6ufhrzQPuhE`8dp3OnoaY~3nBuW3j9DY=c;s{%Lri`=Z&-a|aigXcVCAa8_ z-uc=nVZAZ*0vufqLLXsi%!5#Y20JQy&hdC-eZ8n$P>3iM%-PuLcgLqRk6#XR0%?S| z9-D6=nQ>_qt%j;;lb$$Xf8w0WTBjPZw!`A^H-MuY__v&oM%ztSR#rm)+LC?+=eBDh zE&zm4P4BAUOnd57P)rPFFiUISQhX}A*C^I~qVesL0VD{ARKY4~kAC=jnLehYugX3f{-ZSLCi{bJlxe zj}v)7NiKjL!b&s8aY# zsD}x?T2C5F9?gXdTDK3+4HR#@W{X)IO|fpAx3|faig6kxHK?;vQfdN_T|$uYPF{G= zS)V_IEnc=-fk(_AaEQ#3&~kDrm?33z8gg-U)xv8n8E(uEq>Y!n6j&8JIf^RJEaTv0 zVKEpiah4_lINKh(*I~coO{ERq_uAjz-`7{=hH2pn@%)5FOKWS8%WE^C3lOSH7u=h3 zBQoxFO6r%X5ZzClXe%dLHJEC78#g*OHlFD)c=OgR8bZ!tD(>I2v-eMeX&0=QlaZ1Z zZv9TtN_ladv21tFisLOrUPFM*X=4a@%Y5u@_K10DymWbuldP<)gYFvxNJ^f$uc6WRewsMFZ0q;b zh7`S-L||N%J;$**;>SGH=eRVlS&ueC0Wj0!eqpP)wWmjTaU5|0Fjl&2Q{t7) z%`Ws6XfOQm5FQu`YT@v%&DyEsNG~X3xpHL=kVg6M&ncRY%f&;xX;0sYdJNB(A5DvM zyb*6RUKe^n(7a%8YetJQ?D6Bru-pf*^#FgFvQj(>Lo#CbRMK+*RTEocjP=Vy@KAch zI4>`+evkd#`Gwb*pd%OSmm#b^efso9(UqC5R5LTP3r?lm5Guv&Hb1eeB=y6S){VVp z4)oaV_Mju=c;>y*&NU64>Buo@2fGuRaHqC9biO;)nkDcohbq*d_h&K9I$lMnrJ$3 z4e-%lz0@{=@^4sop+vCWZE5uw`ODb*?n9!T>q<&WP-UR7i=JP;d^yi&7~VLM?mlN0 z8X5}aR?Li!vNr|z+d}z17J~2zp?)66`F?mx!zWMNCZVtei&ua6aH{9C6n!|rM1UD( zP>7b4lx(g{D#wU91HgcfW4}!`dY@1YRa0ekHL3)wX|O-EzwP1cR?@jZlnO2u44gydcj(qnX+Um;RgWr}B zMtu*zc@3sxS#}g+73S%$8g7rx`cMwFrdSCHI2axmxd5`9OmYF_nv0H_3ulNrZ(iMg zd5Fl1LfX#WzA;KzWpn>~q*d8;ue}t<(C)lr%ID8s6e+5yz4=zLd6h&&CEoNTM8g*& z8Ix%G=A!SmQJi!yO<6h#rGa zLgkjmO?Lb^Kj$$bil4ygJoDh8@xOlkx|+6gVwX1w;}3}KpMG6j+n7B}L^#RU3MYMM z6eDc^JEX^4yYT5aR1Y~^#6%P||G@`Us6Rw#qW@1)O0215i+8N7(i;u=sL8Vp8?RE= z+(o!7%+zdlai{4(flYjmd4QjvA0ChQ^(8xT!XJ+x%FvUsut@FdL>Mw^Pv-ipE`b8G z`mieiNwPR28z%tm$jZr)F#>w-?L$7^6q}oy;1lX^N#QHHvPI{;AE>IT-n*A&*cbt+ zW_5K{$_UaA85tRD=DmCO@+=2F%X`DIL{$lq^fwI+C%_<6Ra0vD-AscGghfUBrtRj3 z%5}zjkz)$vg;pfLsLm{p*UR~!rAg$20lt1wO^JALaEDX_%z+&p9a~#2xC-7n0egG< zPoF-04n;US-4eSul22`I0k%Qlp`oF98Dn8-+0@uLkZ(2An`_pYegckq7HP+ak|nX} z)^dDs|2_#LfRxjel<^-vFv|n4UEhyIE~s(y=1m0!a@P-s;rA|a zk-L(|c&W*cc}nFH?1sPoRXmW^)QnJ&tW8g@YC-%~SL3C9qy{%`2~-;tqH4Pkkk@B> zX1ddd%00^Q*eRX8WvUZH49I2pz}a-%R!N;Nu67*m$uhWYV)BWz-Y{F3iHS)>M4hCy zv$Jz>&@%HO=egThEY`pvvD?E^M< zA>m15gaEk6^z^i6lc1m=l+LOT9wa6tm=6>^cyk=t20wWrzE(ayJ_u_jx$1Jg=d{Ac zXJ^fm+Y3GRu?u@sz#j_n@x1}pB_djto_^lZ)s^-zQ+h`ISbO8;d#Wk+o6F-&l_*9p zW_h6aR8t-_eug87{Cq7>-!wHfb$2hV7@;1I0Goq4q4O59NIcp_@(hKAg^G%b!otEQ zPMnZw78DYioS2yJ%eMl;@OSqyIO?ln(Yx1WChbWtU%vG5@%jDxH>9N>KYo1tcu79E zsJIxArMkMh7lpI4b6=iCd(l1QvP^A?hfYa{ynRb`=FAx?s%{`ek|HnoSA~a%L(v4k z1O#iQ!3ak_Ydx@owzf8h`F?hUzFxZNhCvtOh)zygMI9+lV8t+_#`mnaP zdym2z4Q+4SxB>hn#fOxXl+x0H$FP@SP{`HQ-Ab)T44&sAI0JY`aB#4IS=X7fXA_f? zQ6ha{Iq-S#AT&5Q%YH^x@$$j3tn%C{?7(#WKU+4f&RS0gIIFU<(xKQr#4Lbbd<#<9 z)*mw&z$%UWE?_5`JXy8E|)As(n0DJa0te;+<-YHZMR zU}e2={knW`$H2fqS65eC+eNRM`g+*$$gWVtn5cs+Prsu00J-`7`}YL}1xZO&0t5ZV z{8c}Ge9b*;|n8)vl&}h zMDnpo&0OY>?{1z)cAjU-fL#y%MY6TGH_4uqrv1^OO!d?YLYDlIjX;zxEiFNz4^9Yw z0JH^sc_;9~U;z9-P;u#ZA8O>8vtPLqYS3dYXaw~8e}>;H=X@(4*ETn&fBDi;q=0;b z`UlDzoCdr|UcNR$zzmEYU<7frg}FH=C+8O+{N&|(itJNZ&mzxVanY4LA@c4+QHR-^ zN=k#Jt{W``Mz4oiPhP`m>*^p#zTL`2jF^cLOa3A9$m3k~x7<3GSN(tH1TnL7aMYXT zetvWo{+;LeX!RX%NC>`bYilQQ`XOfEYp}x)IFPJBeRN6x;lJXqdpBCO9C0Bs7U0GPhGxzD>z|(ejW;=u&{HCvfA2%c@}*+X=mUPQvS2NA_L46jI!X% zQc_Z|pa6h?YYg?g5zouYYPh?#cJ%1co}QjwRR}}R802=S-^9SKU%Ya~0;pfTVs6(#zM8CH;Z7X_q_HQ|xF3 zk8z0#3>7{cc!Ni(@vQLAM})F)Myo-RmPOM^K(tIE#mlQ#1i)qz7QH#ojTWF>h16|7 z(`E929*%!;VyQxZ>m~yLAosn^$@KEQtrss2`BpBE`f_8}QgrO2OBc$v!4?3!YMjEu0but1HlJ?x3%D&KjUS`VuwNpunG z6dgeMYBU1r}@mp^D z@WgXCh@l`*Jm>ns@lsxdIAP#Bh=@qesaK-1swMw?1*4AM^MvvNesBE)RW5yS@N>T) zwY`@IcX8<4QU$8@F!D^gZOqs}@hVxQ8!7odZf&HQ+WRP2&ok+rMcY+LI0eNuYUlC& zse?;PQb^c6KRAUQK1fYO8u`|Q@K+`t@s|*b(KAmVU)@I}8#V zAP*>sKIvS5yKz%Y=^|dLn1Q762MH=Bor)9;W;QS|c$BItkOJHP?6{rKvJZjx_<-N{ zTZ)RQ>FGdK0bEM>k2uu%1?gX764PfHjok~Rgv0|OZzLrpMMXscG>9{bjEuCkwcVKO z%XgSlUo!)Q_KbUq16we(b~B(I(NwInlj3E*N~eLa*O zKEA$Yrl#u)L+P~qH2#1T0M2n=x>VED)RiX<_vqvuFiPB`U$Pnk-T?KbsN?(@X#mkb z&GVqS6CO^zE(_p*gj%GI{sm#=R5S|(6_R17pFHX+LORz@ijq}ejtn9*^bb?J{4WSW zpAPsL0)YUD4Ja@8!-)%mfWS<(=-_a%Hh5eG0pls1d4!=kS~w2lRkblp30j7p2geSyV`y@`)^u+2M-{T` zxay($436}c#!Fx5L>++rT!1WOyh9-k?Z+RFGoKh30G7LbklBaNm;OL;&#&3+z1ZD3 zT5M!+mh)xq+g)i| zDOd|Q+8UjTXGo5uuSrsj1blw^^5xGdDU|=Ha7W9AS z<`7Vl`D=1gO;y!Y3&~|(BwvEmaxf(~|66R3$-i#(Kj9L|Q#X5LcxEd4z~LQce+ETG zmAG!q>*Zh9IY1`UQ&T^4&AOc^Z_3G41DPUMI1C)=SPB3UhF3@L?JVA1zvjyI@og>O zcz{`uOD+gvFEPv4)zw|Ta)qwDonH7UGdp{2G{?KOldvnS`*?J39_EFZ>!p7TOQCbC zuA#BLy)9Q*^*9q@wt|`cXWqfCzVPQ_jV`u5pNfiw3)Vt?S>kL9^rb#tJkM^rO$v!E zXe?m_Q0iGKD*VPfFq9kRy1<8wxho-YS$alBYNmn?@KL0og{iSI@ERY}(m+lD)uNw| zPcZYR;^N}e)YQbp>m)&M-{xEN38r#@(UNfGV2#ou5F% z8mKbd@;Jg!O?vDZRr#HVrnRk|-72)!Jgvdj+iS6cP_?F zeD-%k&-<8|w8$e@@969612XPZ)7kCPu(~cM$3YzHI$yW|j2+aC9CrYFcCSCHKQ29QI1)`t++Vt`60 zE!867wbS_y?SlEAnKHfzNIIm-^8_U|BXhlI_Bp3!!5tEt*`JwG)uN)J06tTsq~MUh zE>n-<^3~Bs&1;YdjhcFKRrq$J7cP&BVou+*wOwJti!DQckyTiDGpnA<^;eMGv6$U@ z-~c8ft++Ju6jGLv>-dQ*RaiQ@P)L*%40qm4cH^D>z4~sC0gslJi_d@%1(+1=H2hb0 z5Yj7j9oioJ-I?M|0UC_`yAz5kuT@dAew z2Thw*MsVI$>q2$4C_@`WX4qi%H5sK6#$9|eY!r^2RFg#MVPbl5QTOUxDO<}3sxnW@$n!- z0vBDq@r*$Af;V2-+>s+jj?HQ2GAu0?1uB&(%XG4i3G0;KMEFb*fRN`I!FxDgB?J4#btH zzUe=kMfwYx$m(Too0Q68nGi?J;7*{ug2;lP^Ln0&Dp}JS?lw(5*A#N*pT*|B`~$y) z_1%UD@6z?g{TNPk%6UfLc!xT1csMTd2=ElaVPO+n%ggoLRMN~WEGJ4>Sy&Pa3m?tZ z!5XiAf;cZRRvq3A`GK0nmE?gsy221O3cEh>JiO7+V7>6-<;%YjAy|%a>sK?({>}KD6rgoaLRZq_4L45k(P$M>k1l#Ie-Sk z>?xft&EgxH(7=b@@;($1UmhT2Vv!_t2Q|#QJ@wQ|leb@?xgeZ^3`N3v_$#myVDxx8 zQNhLgz*qp|Enw1?fT#|DSZ@0$3E%dZoOz9=|KdE}lwgZDa`$trDimc(7kfBodT0Hb zebN4X>t@o{#!H-xve0Ox2}C`0$TY9$;*kkHA;h!MCP)zR@WYn@+Ej)qZq$2x$03|g!R zh%0Cw8kN({$GmLXpgto}WWi*Z%pipzev*eSw9bh`vy4@y?t$WGkecJ!SJ{~gkC$2} z>gd@#MSTJ#;roXI=K;NIpr7EEE@XP0FDoQ*aI^HOxm zdzP`Jd&@YtR>qh7^1&&IM%0!n(CTRHHN@Qxs;cMT9G95FHPB<~CU0j)SdmicS=mvj zDOTWS+P4^7{h%iV4zIsQDX}U#>%)$|%!FKeDs}r)qaUOY? z@a2)_gOhKK#S=!FQ&Ic04%gQgi4e`Dq_VL)`L&u14y}VHn7X0+17)v{eu65$jx0nA zu?MwtqVr4zM@3pyZW~um#mvE(_5!-96FtZ19Rz~d(F#Hvrw(pZhiteDr=8>$C|=@= z+~shxu{+5SdGX^qT;89rz3Rd(ChMs^L2(~kjzFc3mR>{tE`an7h<09SkKz{8(Y+^; z?;QX0or^>_v_B1PQFmu>B)Jg|A6#ku(EGJ-85{-j>7)ty`Q3)u=gD6}1UOHA?4Xcc zjFR|)>QqN#B5azt&>|nQ5l#h>K7tt1&cWd_LL!@+&(?)<2w4qAy?gh`!WphUmgpO# zhV9(at-}oh<7nQKm6(%rKg=MzCn_lD6dm0auPI1HR#y4I?)tJwNR($~-J7`ztJHrs zaH+!AcTVZa&$B9+ov?qI8Y)nxVoE<0zb3W5CObOqZN4^KfnXqxbs3s-}D`yD~v+U}0)nS6gdhY)tKvS6C=4B;+_$=5E%V_Mlu443RfC zFK_cA_UD|E(3qS~x$PWowrYP5QqZ};#RZcbtUD?v|mK!OLm8Z&=^ z9c$M&%T$S1ja1lKr2MxlJY60MEB7N`VfQ0vUj26l(Qtonw@xa)m70<=8@mf+%nM+| zEfC^`x=t(hmk_C_q#z4d@r#7pcEV#HlVL`6fHfcg>?bbuyqq>xopLqkWN7+90HL_X%n zHL8+(Mu<_oz{Q-df#xmHwFZv_BOA}3Ki}HgD*Y-gEp2_>0fF)lLF5P`z+5g$2qdCF zi5&Nx+nPkYIbiC)mVKAR#ob$as;a7>h4A_F=a7&y2gNGQvda!hZx?Y7(;e7!z&%XD z!p_jJudwtjds2Jax3aCZHPft{3vr+ppP?$m&o7e~Uvxb}^rymu{-@p;PC6LrM*taq zQ1Ra(Sgcdf!NftB9s?Dv3Dj-)ccn$d+7?;5Tvq$b*E@;Myk& z5f+-)l#1*#DpbwY4?_Gy~FrN`>rj>WAQ0eXApwRSe7EJlv$ z@zWp2vMHsjx$DiKTLv7@7+#_b6q`*=O5?9ujs5)nwV?Mp(HIFmrUQI_L`=dD{Tp)V z1&6cbqU5*URQQomDz1!rEsZkJP1b75srbgPB;bkx#cBzdVN zukzjBcOrZuqB4>6P!py2I$k|0mUtmTjj2`f5K!M7?CcQM6+vt2;4g6$_h^GA3c#`j zsK3g}%I>Ct4oBIGqQQljh~e5%jTujN^gNFOo`);OC`w#~gvGTcK)F1_!|p!DUVW@? zRZ@tQSwJAStZZ<^8f2Y8X`eths;CH6bA$pAJ0C%C<@dA^RIuI@X~sa^zQ7@j17Fj- z<+#8P7iI39K9Oyx?L(gxziTVke%EN>ygx@f3IHIuSYxEnp%4oWr`hy@y!t& zdNAKg2lNvKF@jgGA`qNrHzGL}yM<_F33X8AZUI5na(jl!IN6QbP)R7Uf$v8{0U@`2 z+E~=mO5}1ZCfP)5D;}2f-*YjNlCu*0RtVDD{DIRt`gr{ne_{hVQ2dG9w$1jpbDQKn z6_o#mmmS-Gji$-Nlw2E^N$-(4{qpBKM&HR>ciXw>ZLZB;p$nK$-Rw8kC-91hvp|bN zUA-4nE#4F$M-Ys>a?qcuvNjpvgG}~Xg%4}|b8FH3=Ytl7%|P;OZZi3I9M^AL(y8Bd zzy*tQ`Q^Ph)emg@C}@Akj~$*-fi%Ccd~0Vhzl&G}Qzdh?p^WoMt(eL0Gvj6~u0K(9 zXWc9M-JB$iNPn2-B2^nyNFb@b9{{`i^XE?xXjl%G2=no&^g;RfEbY@diJk1VQ?n`< z@?c)Zu8^opiqT(a1G!g4F(&zFOUm6RD62M=L;q5T`G6wSTE-E0gSv3XE7A^Ka}7V; zSQX-Y&nf(=WLe9kv1KlL^ZvD$9&$XVV6lPEg_6#!>kB({dX$yXc6lI-f>P}B=gTZC z-r^v}_*y_%so^+*)ADz$!_jBs&OE-_6ltLu&Y*bwK|qaBb03P%mUfOpX}6TNv_f-e8hgwp+k47Me(L>Q4sMJ2fv$V8km)x1`nqnlT->UE{Fjb=z)zMO- zv$uRmgt%_;HwoMhPN21k5i^jm>)7m(w;&@k{MX%xAY5+YPA>92Bjq(8geaC~fgd`o zUd&z;S=J6QQ{e56Yo=Q28SN_2pHE0koa--CA!%)I2bs0Bv^1#G$d3*9u>c=K(D z9^4}NFnfy5L_NwlTtG$B90FY)_bg8Rmx@u&74qXij=lwDBrsZ9&iw-ebmz~@p9EK8 zBK_V~?F13w=6d*f6ZH^dp5~H+XtUoxXm>v^x-*B%?KL1K;2lBEG6*~%41OTQ2L#3) zARvmJS0(}i0)Xd&Mpe?6FPvOlN(dBx=0OSx5dTtZTEO+@cZv&iSx#wZ&N?r>s+hjk zt7l~y+8ooJbCp~MDAMF)Ha@-;pc3ldf;g)K1mTFF{H3>ds&?s?FL#*}RGl!WQ40Ea zV_?A!*sdMi@udnY4KzKfxJ_l9BA9pNs>-q^`&H>p^*qt2Shrq}0>!BE_Z9}wTRmM9 zkH@%xXf;(Me`0)`P6E3FLn-QzVS&7dwrxmRHOyp{0BA)h4G@VHb|FBsoaY-MuC@kX zUorVH{qLMAPe1V6FIO_{I&4>41mOF|%l7V;Q_{369&(1;4dd3?gDD~s#xO3eEyzG1 z;9c9_jNPB6aX_911kmlk>2~Ho*Ls-hP}v%9ke*sm?6O{??uBb@ zmJx^nS(A#G!1Ar_`o6}n!WqLD7EV@-mhDRAkxxb?brjV}r?!4nwIvu|9 zi#_$WXZIBQ^8#ve9U8GbT^JD{J7@d=c*K?#*->0wkBn1M^=9}IlVnzP!*tazS$*@> z#KT$HyDNK1`mosQVMoy|1_+x#QRIWhn7SXSva_vdxT#X-Cez^Ot0Nj6u--hUSpB5$ zyyDDqsIUr4>g^V4gq)bdE6Zt5UvzzyndCOb=z8;*Z3w$Tj9yiaFLB7MO~!SMty24U zH}Bavm?$A4ibj}1B&KwPpgg~Y^VElk$wi1<7iQ?A3y6>HffPNtePdrW8q^pE&0v>N zwxx=62o>zE)F-M6ZIZ6rCe*$OIlSzRt**9MXPW!qv3TVKQN=O*)TeKIq+&mHfFf-< zHB^a*$X{}GOl&wC_vGq%ss=$Wg~0F9p-jt#WPm`nTTkktv3K>gE<}iv&~x=ewe=oF zUyX{+!yGcK!=0gdtRAkbNx?16g;B-G>B!5!fF)~_nW&kkj(R+?%)Qk!IXF1 zlKyBaDho#;9dX7`GhA1lM}yiQoUFrSdqmL4ESc6?G@`g?m1KD~ZW!qfcud9<;*@48 z&Bzwb)kiZ&!u*edL>0P)5STq(_YQq~5=U8YS z8yz2iQ^*+=&9)U8P^*>9pVrVwzbIXLFI)-Xhn{&X zabl-Z*5bBjA%fX9B>j?Qzqm6Rge~5!Nru??i|nu5Plb+=^iuR@imx*+OO$zT>m_5>DW&6lBP3CZjh!Yg3Mb}G7!i>3Uc7U(S~d$uZ8XEAo8)%$ zjkN!YQDpGxz06OW)3yOjbL+~zhfP#gq9M#0_?v2KdpME|37jj2t8L_H7h zc+W~LS`#Bjf4p)>16>qt$#uTrZjwiYRXApC^!eA`hb?1?_N8sB3{K?Zs)sfkt*|Q) z7Nf4&JpP#teK74Sq)-4~g#s9_jz*FTKQ-l;g74`B_QOKw%pGYh-)wKm?JBu1_9e;5 zT2}{WQr2WQVO-t*<@g7`1ua(2+94yXhIxt zu_eIBIXbS$>p@9A_TpMCzpgEQLVsO9_UP9MdKy9S(fJEZ&ug;W2#X`2e%W z3(|8ihg|D${4D+9A-WM%Y<^Ti8-dEvnP^M=5@+dlXlNx-eBI-H$kiC z;Za`u%aO9n_3vaMBW>q|q0W?W=EZ(&*&-=nN8fwm#qONOyLWY@w|`6)Ki0xqz4XRq zNbL?@5zqc+cRoJ1uF7c_8f+YDX>Kt30uyK7pwrUIt@^yv)E&mRH)ea(t;>_m{!F(rh{|Bo1qOr&B-?RYWs4jKY8IeeZVjI)3J3+5508dOQwC?W&Z5Jz; zA_uYuW)%049fhvrdopM-yQ}o+{#P~JC5in1M-FN410j8$+G;1friKHhjvkC&9+B%A z8=atEA#ZKxq6?0Wjm|-v+pPc0;WYnl>~iE;rd2h}&JDt}Hne47V9ekrS;zTn5C5r- zSy~)45DJ4n=rT(vCF^iH3p~a=)RczjsTz3|BA)`f8@aaiuC(w0##Ot$j@7gw35OTF z%I4%W@N;7v+IFgO0lahn$%iT3dh~Z4ZYsp`J%)=3yjuyj+~6k=sNS^l3pEw+NJY5W zY#-^dbflfjL5r|9>9pkKyk30VE)&f@p??ftJMs-9<34+2lx-{m*c3#T`HqtN@5$I_ zdj1?(g7L4=?Nn4LsyJ9$#*b{&pKLTS<2Kr&d1Fl|Uj*+khe!L<}o7npYg|O%K= zIlo^#^w|PoHn+w!cMvqvFs&cRVZOSywmWMi;R-!NL^c4Dr0s$l-+v|^ zQcXS$QihOiXrOh;F_0Grd+Dj>unV{Hth9BJWX;YJzGJ7tvo+GmW?fRIdDDg&#Y|nC zIGhwt8a?M^^e+sXhrtJS`8K307MBd7WZIL}8Od&`f5!t7T$EQ4W(N9yIX&j5PMyN+ zN&y%MnO8@P)YD5}4Ubn~Um@q?vX~zUp_{9oQ z=G6&DMcKdC54&x7it((Ai`z}z0i8`oavC&+kqP-emcmL9Hjy8*&$|c?9e;-rt5CAW zTbIF6-b8<_Hk`stT};h~MD1kFu0wUElH3aKW)v6IPdZx@EggPwep-v`)%p+hhe8kl zE9PYVQ-v5s9|bc)`L}8O|8cys1!-plH0hqbu0}27VOXz?^?PJhw;XjC2YnQ-56B$p zSXjiw>2JmwVdq;-B*@cb=z%QYjUG_ode6eh8`}hh*%7*Rw~u zYP(4|t{WOho!9U9xo1OSQw2>-cI++)i|b5tFB3zv7KR__Auj+M=f5{M_`j2!gx(`# zxQSF$c|AQ}a^$&2-XBe8-ePh@@s2%_cwvc!u8hmgA0;ZIhT}vn2XY}OfV!OxCgR|| z4q%8TnwtIK!7|eKwSNCJ%}s;O`%gk#AR@R}9G-R5`o4 zZqN2)6E|VdVIsFjK{bIeZftCmnbq_fMRwEgv!TBQZ*W)yi5aK^A6z;{^lUf*yA!TE zQ-M0*xElQVt^ktfwp``hnnN%fRNG|k>JP@xhepyL!9U~O87P7IFK>!Btcs_l`$5n9 zS(X3AzeVRzgtMNOV7)f&b`dbT5ugk9jsk{nKuDPdacs^W-3}5ox{dpia@66lZcb>o zS_hrgdia@h9cF7M_eXWy*B*3_fz%vO6VP>WZ;@BKd_TJ8P{Gd44cf*_TwH0NKJ~!N z%F~wrCaTYU44ZUW+Yha6EyP;p-U}1$C?0LiU$-!}O>y~09XfT?gcg&Wk2o({ zpWo(ch!S1_HT2H{o3R!ZWwE6n;Cjw5=Z-!jSOVhWXi>-GjDI(8<2TagV-Glo>i#SA zps71SV?FR9lY;LA^=ZSxm${;@<^a#>Qlq6qZDzGc*y2-fkkWs@pon^bGrdtlhbWbs zqlDwdup0x8I4mu4}(Sp#IW zmAj7S^*;|z6pLa3Xv*`>dg@((I}d|$c?AUS3kkO3CyKA0SF|02=F35W`m}t&Ii~_!pomExEa`a-uv~U= zG5XX3@>-T;v4VAYwFBtVkr_nOQ2vttbyVeh%oGukLX_!J1&Bn43{mWy9R<31CQ>g* zN^@#DH^~q$V2V91&iQmalzoCyr>`z~rJj@?dp)A&z7Gh)R&u+wFPJia6_G06;0ZJy#(t60@^d0UEZuT!@$ugeDe(8t`VLOK<-oqTne3U${E0$z73_ zz&=t19W$})f*k#E_pqS5^EkavX=YrIc>PBbOxk~ zpzv2y>;CGgl#H0|HHX=75rm3I!c_uHTl_6v*=Z`ff~pQigTIeh`?<4&B>h`9VPvhV zda5FgznT8!`F4d`TXLJXuW$QmLNzQt0Ofs^-F@9vnn8V{FJ;gk8;9=sEU$jeM zvdQ&lHuQYg9po#Zd-ET1=&@H;+91(1HT_)>;{aqyX=!Qvg6Z;iGH#wMG9n6Q{?};B zrfVBe^9A!X`3Q1F=ogq*RWfB+(@W?{bo>@;jn_UOLnH1wCGJELYXC0d=YC zf#KD&$A;8Cf$OHN_=a4f6Zei(-qc}v+i>=@ch=O=Z_H0UZRGnC+kSWpJT8uGy@F+Xd;Gw-F+^Z zVKY_*5_1Sn0WptmAAHs+MO;9BNyKWJjkR+bkxRy5)T11(D{ts#{;WD7RQP*$`IYR& zbAwc5UcpA8b-u0ab)59FgY{x}#EN3`_SK;k2`^IFefEqpS%nS@LJtPX_-0bYH>&IbO5SffQ^FUjh>&$Uz)WR4c1_J{tw@84v zZrQc_fW2fqCDWOtb(dDKxK=Ukcjaht;%NWHX6Gb|~j7D^u9Y!>RGny?fj4 zJA}-4?;cg8<4#;~04d?iPj_HU4yMToo0GAj?UE#n(9H)khpG+UgJhy<*<_dv0~+FD zVmrILTAG?6Z=lF{A9?rA?$l_DoNs$tNUu1*z|V74w7%bD=W-)8$Io^T-Cc{+{Pc8* z;c784<-$kb(bl^b`&T00y$cNs6WnpB31kFV40;vwhnBp$U!Q>uUJ@MMe`;+F5>^H; z0QIntP>#F1yF>q_0T43LfM4(oZZ-~e{@{9BtB*CTg=PFmw#N-#=49Q9YJL`}Y5&pn zy%WFW(w(fl{UT&UMFw{=8tOn8bkEMq>+0zEwd4f9KEa5fco7j1A5_l^psZm+Px)dH zsE@llJ7Ef{wm=afpmLIY*W#o3(m^Cq6Yh|>!-D1%r^^c^?K37~X>y*_$6 zw5&F??U=Z=KF1k+n%VY(I5YDv*iWR+h1?7fAF4x28X6aP-wb3sh><`Vvoo?VJ^^i; z%uDijFhNsc1w(Bh(LF~?s|iC`B&{tii|gwr$;jeUP63GJ3A*@OE)s5p+`UC)&gb!k z>C9AxM=t)|r|;R+>Sw;cvrPDO?~vUMxrG$pQK^gunWH(oM|lJU8XFq?UcVlBTIdz{ z>J>@c9VbvxKy`w=jZ2S5mK16dc=yhQ3pFsUu?KI;SzS1&WFZZ9r|ay)Bps0CX_|!} zV4y1MCai~FSL+sI^*gmW+O0eHHgxZrr~TDe=hd0N#Hx1BX32eNPP#UGvW}#M{b{Vl zb&)plJ7OY=S|{p>d-It_fHX&5)%H>cCS&1!4pz`XIqO>qM(5Me^%>?VfDy=tmvO94 zx5MR65RZ+HPBh1sLvR2uX$w9M)+TCVv^zD`g7Z5&$!~`eI9c5vZXq1r=C(FgU$nI) zS~*(c+GhI5)w}hkfssa>Q_aiKrJ}HM<}Al4CNE=CQ$dX4rv6D7jkxz|6_Bm*7qnL( zoF25#&^?TMj%>I5I?*e?EynaXl*>X6jKXH&p7t2tI5Ro^AZhh_?lrYbWf~V#sp_#; zv_i6f&{iI*YV==Ma1=wDfXHK|wf0o12a_AU6WRi3JBcr%UVIg*q>T{>(Gy z)0vLys|zvsqI{Q|dn!?xW5$Z=mJj;UesGfbRB8U%(#rVwRL}QLJhJsV2{15yNmvd? zPA3vVaqr7euttCi(RjIe(PTya@`dd!4J38NXD!%dTWZgXdVY_RpNt~X4rO(@3 zXU4~N4V7`QyMBq-C~zrtz2ux1`eGg&@UIivI2U+_ z8sJzl`t+q%pie+&&2Hf;3)f*a_1F#EZXvc9@T`1%VaE200(*3o_rz(Nr3x3RxK)$H z$3sOb#!(iPUq)M}yQJ=jO;v@9^A~7ZLTetz059NYhsuYVn-@TO%IF0dM&#nf4^~5E zFs=Ch%0(!aOy9Wf-dT!6`Bdng;2k9&zgMcA_9(Qg^Y(mL*Cc5o_vf=4Wzk)GQm6+H zC!UW!A8j?U@F4A@eQ>AgnLKDHKsW)f6p>`80VLemc$Vade=!BUxQpG_XMaJL$fCf2 z2cChywQ}p*mpIfdm1ql0xPC2*DL41o1cBU9(T(3LUI}^H>t>^h?My6) z5*oJ{T_-|)6pBAw@-g_*OmF)Aclz#f;Qh;m8}GFlLP8`|euz3mxEX1%=QsuC6%qRo z+g+a|&I5dnz^29jQ`?)zW7V$h<98`k#t50Ikdz@L<4r`Rk`$34l_8aRmPm%oijbKU z%1kmhky6}9gF+?TC^98Pzhi0d_kH#|?*03GUw`y?igmAbt?N36<2=se{QkY9L}t_7 zx~A6Fi}T+q&?OIrNK>Sr_^8>RKT()c78E46*`ve$!SQdo0^ia+I~OyUPZOuT`9<`} zXW3>Ku4=@JOh~(-!sd=>1v)$5`2!Hx0Mf7M-!PLl1B4_K_h-XCQxo z2*WGjfd}b5e}2zEZw?dzAb4>KY6X|^RHX2ZhPk%s3}&tU)$Im@404|sr02#oCWpJ` zepxZTT>H`qn&^C1SZg3&3u-keH+S>%<4PBmwr_Wy`0|W#I)Ft6 z#9 zvel-6OxAhdA~5J~-3nqUIr~m#Hu%i5EbS=b#m%uIJw==G$B!RDgMr+^7E#nKNr%f6 zNbQyjHmxu2xG`N2EATFwnlBB$QA9Vi;GQv?6^c_j{MDvs@!e6US6^`s@hQiYRBNC?09KVl0i8YrHEFIzw1}09o);U{Gn9)&|a>K$tP5n^Hi$ zlkA(ehLY~GJ2G76YLqiA*SxHDI&4}gIjQF|{=h+7JHw+a=i9-?D{P{%(OLT;V+Z%) z@9*F1JAbb4Tcy{`C?6Z!3kwTiYIb&a?i;;C zfqENkx{QnyDxc_)kEn%#I0t9T=v<8L5w{-rzTf6kuc9JbYV;e+w(Mt&tUNTy_p-sW zA+Ly1(@~nQr)=yggwIHJyKr=yd0mZVLnyXrt*^tcpW|p%<=Ji(Zh1*AJ2%%8e}RvK z(U}4@r%${jBeG;}n4Z;y{-S#;D zuo#3lEk5~l?#DOp*@@H8$>={7w_iLCItv4XgFB8D^|*q5YjocsLucbK-thZY?3eU| z-ox61R%UW(Zkt+1)mdMsFXP>^&CcR?iOEQ?z%?nG@;jN| ztx=2JF_qEyDS^8rUpL20lJmyCk>Q5G(|-cs*&40eIXdIiDF?bk5ue*sf|PIsEX&3{ z^YZi8tqz0=$Jy~i+qV;$kGLSwm+qo_ZCUh>H&GncW}9880iD&H%~zYHw{J|2=ge&% z=5|!ybD<#mRj2#Huatn#ZJJShqvU?QOY*8FyiC?`rBM6ID=0v$>9df5z#W-QB%jjJ zm~mV={TyGOz_&!t={ZWUa<#VFxgXH&RLXJ8sGPDreW5^4qUPYnvTrs$#xoY0<9adO zG>K-)6vf!P5zOnu8oCgsp??MB6Mj}!X%{Lbs1IH1pjCTLeAxqKtIxYob|QPbX8M*X zObNBx8nQRYaOYSxWV^VU@yDzkUwDf3cor~GkXln(bzxlZQEI=F1V7g(^-$O11E@Zr z<4knZrUeMLA?<)Cm>KJpl#qzk*$8SbIDrwZuR1$be0~-QcMzX*ID$s--1)h-o4#Tv zr~#jg5GkGAcQVhtFg|$WTE|Zp9L-|ByqhdYuIwcTbu5hbLpO%r*;(orppXm zQhqG@$Bre&$7>TJBbzRxlYMgbSEKdt$jBCP5iC{Ht3_B#!fywhJO8BkA=&L>KSEux3*Z1Cbz&EV5`m~SqFY%vYFq2lGm^Mk>wDchrI#;hvOL2Zk0{bzrrV-H6Z ze;o6HBrW6(7#In8jdkqAChp3e{dGrvL^}+Acl#8Mkf2KBx_CLgit&UmhqP+bb@b;F z#4ebqvq#aY6nWs8>)bAgX$iMMKU08>4H$2i!;=ROEZQ=)n&TzMx~nU`zQJnxX}A{n zHaSxRURow`IRrUv{uXbzQa{%~|Ce5}$;QnWU&pB#i@I1Q2J?n#4<5aq<=V6|lrErb zVOf4&9%O1hPtM|f;gs)?m92qz3dp5MY1r7@VQK=g{9>d=`zr=ex>=pt1B0 z02(N@>v&N(AQ``wBDJ^7wfVGpA`uolxAlh}+fv+P{B=CDO102PUr#vSu&3ZscZr8o zjw{QD7poNmukL;S)7={S#_U=P3p2eh1)N-5*uxsOs=i7=dU6u-P{J`eNDo}B5kCul z=_KjEXiX*>njW8i*gn>}{`qM}ZnBx_uYJ$ct~iHnMDJ5gPBf<$J9A)9`|8Jhrf(KP zgRiVI%ov2AXD5ZCvwy#}V0mE-xCFNXp@EGo;mUU00uI+p!DB}j>)NRI4U6vLVM z4ffAg)u~ps-w8=c%}T_&G_}Jmry{UOG-MBcEWiDfNupk@S^lqBohJ>D`?R4|Y$BSd7HMh{`ljsDy%(25#{!xNEL>bUP0&UseOFUF zq^8)b;x+Tx$nl*>(B;YD(0f()RLzYp4)W(7`Dw{)wT`*Vgid{&G+;A7S&RGzJTo*D z;shhr;^yM|P=E0^G+B0>|G*`43W;{eT7BX$Z&*K_H46o-PUhDnzQ|4Sr)8sS=00U1 z!dIxjUzZaevZ=)Uy_U*7Y-n<8+YjqnM8*KS(Rmsg5z#+6%SZc8@+ziHQif}y_sxyE@6H4oym!iSRMZBZwA-xdCtGZ<*<4sC-ZKJ@6_WX zf((-IdX$3JT7UnCz}7Yh3$I?Y28RUtCFr-ibxXvKb|aQSOyjZgo!y0RqP8(!$Z=E{ z*nQS{TlVys;$S0va|4lxVV$s;el|7_nihpse(IHphy?E!eV&hM@YHw~AR21)wJqawM#*D78@{a?{WiDx&34|fnNE3~~U5S8|dMT1|pz^m5xOPoP7VzohsAO<+X(SYg}RPCk^ zn-6^XA|_4fG?d--T!>bg3?1$WJNK|ZVZq%x$}mhNRL3>nr&Lao|68(kcG;)M^V^@T z&fmr6T5$L8TmU%)nz{YF`{j^z0ine-Z|(- z5rgY*>tgLxz2h4P8~DFT^7EN?Kb(^5=Nh33=_z*=cG?Ru1#jr9Ci}9t z%i{0fu*xbLkB6Nn@_b`ITojV|=&7@y1?aB?ZK_-$>{Bub`zRREe_nQCbFc2r{@v6GVnl z0`gSmk9s;eJ>6aDEZQxN+h*^do^#wVV|y*pP?FP|L3-m?ENq@#!^)W_vkmPx?(Lwp zy>0ng#k+pU=MQC>Lj1`b82o4>-r3(R5KDjdn z*R-9ToGt|fpe-5_L1+WIb(v@pdj7^~A)VeXeVeyfyo%0Pr!li|ZVvk($Ca=#cgJGw zB=76Zj`su9XIf9&jF9_nM9Y}8?55}QiO%($92L5_4XUa=U0ri%6WCgyQ(Ru&Kkk7z z9wJFLf_R~(M^p=|b;Lt_Tt1|Q)9{N>_8z_XSd!_ba#4pN6DblHD{uV$@rEbD|?h1c_jT%b8jd_ z2@RYgR zq8rPn`BXB#?m7Ccw>qoi+Fgp`7yA5N@&-(u4-?|3W@cvCC~ofVy^vk1t-Wwfd*$+n zjS)Q8u3dvLqsP-Ymj#tdN$aiPbw^bP7GrY~-51erlu7BN@^<(4#<6Q$b$UAut7< zgW$kG2$>T7!Jh<^UP3lw23h~ijD1Oz5*DK(BkP%EX~_Lsb#)(86@I4zRSOrbP6kf4PQ66Fi|;o#7a&yd$Q3i!GrYP!sBgH< zCuq6nam*WS`H&YSZRs6hiV4-R$xT_V=}tKkWlRH>Ml`mN12K{~8hoo3&y>SmY~(XD zYU)pRm>k%TKj(j(+pr-E+9n)n%ITbTH+BhVG<7^bTETf%e>=A;b@a#z8XUg7)%jW8 z7DhCmCTF2DguF`|5bQx!^r%Fn=QCo-X#Oa!v3?12Awr_0_Tj(B=wcUeXw=LkQ|AoVSA>{*(|zbQ6e77+;o_?JesYb_bMcwwmyfV$^5j0>rGW9~c%|ej(F|YtLu6 zy4Z)?3L+wk3J*_~Ilb4$hMIFI>xS};YuDJY_&hvzCr<_t{S#zYb`B1REB*$NLtH|_ z2O{3W8#Y+KNNB5eAG>G%Zmo%3xAwuj)cT)k7u3Yxx9u^o-tN?Xa^Pr&y$ij#cGM&6 zJplnIdS3uc!cn1B@gI&hxf-}Q6uRN*XAnk}6&H^~4*ACqc_}ICTp}Iy2Y32DGlS0h zoZ$x+I?P{rR>bWp&v2SLg}e7ry$F?Ugm1L$T0lj#wX`;GM+!$~)6aE*Djuq3-+Jro zZ{Hpy!~`v$&{x)n9hHncD%ld2_v8GfZCSRmsd_)yy58-e7pE&*9!JA1earvtUk~%m zzZ~X|uBb;4?cty^1O0(g<)*&h$DG>qcodDax)#kZ&NtI$xC+s!_v@XY@(9wLjTVWiq{0Dz%xFOmd9v?@H9ayxBp2J?f5p8Iz4@qinRF zpWoN7Um+@Z^XA5F@z4u^hRnW02A-akdSVh{VuF{kn^&w{*_MKk-zyWHeh_W;dufFT z3cn*pHD#c)J@O`m$Jf+!vez*xPkemj;D}ARO=idbzRu>h9bPPDOw<*PvMd@p5Dl!a zzi4Fi7UB#j1aLEfiGpGbGc)tWix)XLIlb7HlgO>;ZRv;NrM0#D&+jMD+#6VD_@;2I z+DK|~#$Y1ZDyS~wvz|QJEa4SEOeQ8KC_G$_4h>x$qk}HqTZM(aE#5hZlEfYX%0T3S z?!ASMgKD|C<;5ROBDo6}@7O(Yl2}*LyDze6Rr>I*e$vVle?QdM4aBnvTw%%vtRf^+ zzfoe}arSNc<%+|GCp%y7xxst3m|yn5>oW&fT{HKWO)*}2-lm?zo6gH%hiwCK*r-Qc_ZALK2UIybp>?!$U)?8Ym;{A3y%0Z`?yU zwXkpuFwsyX0I2R>*mo3$oRMOpBG}#8dF0tx&o4z2H=Y`GD-x~E4t%-Ydbu|f%v4&^O&5$pG z#FRP2zS7kLVg=IotP;8``y={1&Hk(r$_{cI&IuAk@O8?K)VTH7A>0X#_s)!t>({Th zq^eXyGO*E^U)4v4{$ZXwiG&W(Bt_R@A$&e2V#zku5;g5SX4)MuSUfPEJT>4pRzb_6&L477-PN`*Nk085y2=~B{*s}AC-Hmvb^`~V03u4*W8MJow`<=fc&G^rjbX`QzMRFuQCFNUR z-v;*8(EY&^dZI5JVUC{7yDTV=~= ze0i^WzhA1Khek$9jY1bb zF>=$B$TJQLbFK(LP`T0Cy_?uW4TAuM#@>~RDl6CAy!qtmQ&F^-R#icezXNQjy$9Y4 zi3y1_ou-LB7{^0@n9l4NL5#Dsxd2bU2~YPP+XiqU1T260_}T%S&%`GV4!R?l#epWE z2c*P!kK$%(^g<#c$B*~rMFvSe9r@jsLQgK2vf|UBQBA|_pG6Frf|ZlwJyH4&C<~j2 zu2GH0IA9QJ1fno;F=(L!0YC^v)SL!2m^i|qqlzXMtIdi@K61y0ER{8%=}D0kL;ofF zApQx_9`qH_L-PRSB&!z*F|o0l&9s-EzwuBH7ta*T;U$?Os4Ia$6l$kXL%$pp6v-`@ z_bdy_whtdZ1UVH};4K>;DP44nn3w|6M|3oT&$koIaBKz4&3W_a(WTMjvJav|j;Gs` zn)ps~>-hN8l$TGV5!mDd2SiV=Kkb7$8-!3Z{b6Xy04KCDWC4^*ut9+K!d!nF8nQio zTFZUg+zA~`d3FX;KcYcIN|^9qSL`dBF7_PC^00}CxyXghmivFiuc3AZg_;t45HY_C z_MTBdHF@-SUe0W`+m6wV)ODo1aB0RTEiK#7U@E~;`tIxN)N84!a5A@F7{Ix&QBCXB zQ%0{kk`F8lF|**~N3#32Alq>n(>S|TUMtBn4-#lU_xG23j9a?9&j6x_mdg%L13{Lv zAMQ(Ifg~ituOnVdUxDv~nfMevzFR!zj>1b@KuH2Mms?7P7430%ejTXbG+5)uANVgG z-`RO=2@H0az~RI!5a=Wl%2O)Nj+LRiVn}{JEdJMs8ADt$b_gL3Ptg-Qe(&rtOUq1n z57cV*>*`Vu0A4YDbdy$t-toVrOXMn5(~$fa5Z)uYGp}~bM*}L}7|^vL;#Zv>6yaA1 zON_7~Fi;DlqHZVykjnyE4dlLKE$S9usJL|o1_lC7bZ~I6un_cpyMjbHCp9E61$z>~ zyVHr1u8t6^ko;Q5KmkaxmxhK{SK-hI&$WB0{0IQW;PX=}Z8wNPw$#f(bWW$BZG z;Dljrj~i5VAV+$qpul(Z!%m!hgg7n79+zIS>^3^*Cpcffoy%!}}j%M6DzR6`vGP6QR2VE%2b=U=w3w z;NnmddILqqX&O>G4|cle6qPl8Ud?@dMM()bti&uAz-`pMFyVYqHY|7Q-;1t)l!S1G zVQ#EWpI(Ph`fc&!?rzil`y)pEum_D0$0Co^-M7!q$pHo6|m>gW%qw<$FsSJ5yGOh5ow4G0Rh@kR5=3pC+eU&aLXBKrFA{ z0X!3S6DpwQBfA_qOPSVWC-f zsr5n348pJ)@ST@#YHBILoEmu}DJwe#@EK=zKmWv!A9bFSR_f{xF%QS~8Sb4Cf**3o z$B4-#+}v1dLUI?!4pM`o$MN*yT%`g!Nq?TJP4+044%dp}Zj03tGQ@ zcU9mV!l;rH3@b?KgU?^QaB*@%wg!e1Z+m!du5HphjTvqYT)NgYwS&64J`1}8kvrzA zBlx0L)?sW0A~^?W*R3pslqE52NsRq3{wgUd)-6Uu_Gc;Y*T}&2ugAuM z5j6;fH9#&XFEX|cb>4Y{7l%n-^C0h2SOw5Xv>J<&9z2Ig8@Q9n2Z899RuyjU57C1O z`EY28pt!R?Rhi(*eNtaWO7~uRqAPg6#94s1h^*5yHRY>NrGsX`gW4iQ*}RCd(WeTe z`k0XsB?afki#MEbHVWW;jzIhY3S7J0(EvrgB1OBLL@5nv_zbp|q;XkDbAS8N&kM`4 zko;0lL)i!Z3GG8Sm$Cvq0yI+5L4g~DbF>pH+E4_mN^$a;1BvSdJZWrH^_(~bNdgx- zyztaD#G~B?LX7b*lTf)M6fE4q0;5(^$hfIcKvW1qA3Pve1>ho4LBvHi7#uxHbXXrN z-ECr$G>Sbz$s(vkQ}l?$X@Vb*&&=Fjx*3d$ zMqItB9t|UWd?^^65k89BCUFrvVH9VC3&J{(WQVXLaW5odWnApQS(+Sp3B#M3n);_o zvVyqzYA>YMP*w*@3|Sta+z|3#vYBN1r2y-&i3tx(s;w<7|3jil1hiT&Omdz2a_%D1bn3^3Kn3Z~kEkuZf&IxU-fpjQ4iRIM;J=g8(&}O5VBAT1 z>biV1p8o)!QnEVFOhMO~gQ!~IGAHKlS=*m~;TKSCUoZ>v^ZLb8N$ z3a2rKh6s0viUZ$+gkwg4Fi`jUHwa^L*5~xC(M09&H~PvD2bG;}=J*i8iyas$Jwjf-d^w`5la`h^-S>}! zO_N(tuuba_R6{i~iqg~lSTNy;jh$W7U7tjntwyqpZm*rZ2aX=SUS2NO?MU1~)E+gB zu*k?duj#ztU~SmRy?dX-rQ%dyzj4DD%CDb3IU)sys$^iMSk=kpq^5TWQ$ZH_b|`B% zFx5<7T%|+Wis7bzKmZy=mJ8kdOZspUa-%MWKxMFD-y%)MFj*Xc@4JX|OUw0|48V`j zHjCRk#I&50E(ZGmf`a(w&6N(HE&=6Qy2UR6C$i{#RC_U%ghJ7SmJT{bnwkj1fp0oG zw#&)Ijk+Lc0>@_$9bGU2fZvEZa$pWADSsflNG@qhm`NZ|ii%Yb-KeaLSie0eQ;a`q zhwrbS=q{(Eu#(d6Jh2>yQ^R<|4bWQKZr#GbEgLebiN~9Gf^G1><8XeJFbJ3y;R)&m zA}eC3>eB9uR%;CRM|xyp0(mOW63+#}%Wu27COnj_tp}ex(d@UTAw`Fw4v1tnR5KF> zG97D5c5g_aeKa_U)UWwdnWD|EL$Rf$+tXSN2}Hco1IGhk-!OI!K5zW9*5CXhNo@5= zs}ui@b|-arIsZFAmOgSIS)thB(<6Wow{P770@BWhV@fXO+D-)H$R)QT7zcbcfQ);& zq+}t;)d;%8)p$tcL8Mri?(5rWZmdU`C~dzO*9MB?_s@>PR9wJkmgLlPQ!?3WHXc@HxGf21{| zX~$Nf=^u&B)cy1%3O~SbOgq7L2L!^rdi7Q5N}y@!^C$Fyb1Nt)KzfHo1MR;2$znMpJ01DQKnbR&V9s@AY+G>wLmY6>h3|p_aj{Wd;{cxTgqd5OQs`cgP+OqX? z0YGR9+t%zq>J7&Zn}j-YlSd06M`AY0+ABn|UrfQ$(g!{Q5{C;ys=G7a#=i3vY2V8!hKB-6DJzyg56tS!8@#)!#|rJ9bS3mE*Rs5riBg-Xy51^)B)mpYKDcF}jJJsXyx5qi$y-Y*iQX&Q(} zN%>DZ|34G@iGAOF1Rxg%D-nx~jS(Ks5SodY5z(2P^MEc#m~_B1JXxOK9>&WopFe(t zt~{AOHzx=F13kTJC*Y~TcJ*j_X#*w@&Z1@x zw*EyF<#9aWo@yIJ!!$E7!wMEccL0y2$Bk3^rvJ+q&E31}Ff0{BvDY0Pb_!gaoI6mJ zLM(0n;eJey2W9(}hR-NUBVI@U4O$dnwPQv#RayX(pt4cyJ#$_o1dL$VCip9n`At~# zy;qQ6k$yKc(LxItHufSkuN#p!6U_@CKj`bP<)#eaOwdPeA9WnVqk zTQA8Gwn_-@09k7?H4%C^Ea{ zqd$MP)CeIepnaG)3Xs!f<0bx#y1K=$hnyaQYa>Y7qp5E2t`z@*hFaro`mn@Z^yi|$* zf27_MOhzvm^@<-D#)6r95VuLMsriLSASx_uzW(>T+IqFq%uGx$HH*}&DD*-*xL@9l zeETKPx#!P^-xQ|R)-HswZcuc$%f-sWRebHKRg#pn$H_pm4SdF9r1GHjBBr2UaS5mE z#N(i2CMLuvn2j4RKkbuKy%84&HghhKrX{dcCBnR**n{1K^+xOkrq&6#1Ng>sb%)HvM#~XG7b7z8CO1!3r zkuh)GnvXG_&lf^Z{S{)48B0`8)D&(@u^R2EMbxsT?+8xt43sG+CS>H~zEPEc?_#^m z&CZ^KXziMZiMLUfM#B@5PN2+iVK`#;^)f(rf~q0$YzTqCQ(!Rq8J0K^A zmT1&(QOI2i>MED}r;k6)#x~=VvUq?vaGd6ljU}QDygm$ZgkusC(!|V7ybgR}p8-30 z@#1dzeA!wJAZdtu6c=6-eg7C*t03TLzKZ=;%9v;AyXF#^U zE!_y3JvJDj)Ce5LwRy`|*bpE{K#`6fKJ1LM22){O`flmCkvap_X{f0;^@oro-pbEk zKoHWN5(q2iRQUmrEJ6ZGidHJJZ(xd|khBNBE&=6wK^8nk29Y25owX69s9+!sN;nXf zH`=!^Y*a$?@Zm|A7sPU4pSxJcBw^{_b#;N-E+Z&NVBv^_&l*`mb(18>!%gq93;3fs z+z~w2ha%T-aV{iymX@D!JP_1STiTD{FrQ240fJQqI0#c!SKoOIHxZ0+ zr0Bn%@e5kjnlTOnB7^`4@7;3~b+hv0*O1`A`X(#^ksT@7D{&Y=(kDbZ3(RWISXhwV zlOtKf8qmD)uIR*LLT~|#X!(KI3yQ((8K1Pb%LofUdinAjuI}vUxYP0Tita09ExeB( zI4^+ot`|4c)z+cngGwqs2|50C^{{)*&!6MSV^D|NSPzi?J?_L^T1biH(%Kc0uohNL=K`{ zsLy~AMm?}BPe+U$(}xP7OAiZ&7>*{sxx3ce9R)W;i$g8#L3sVBei5^So13?6+SEYQ zE{Npmfyax4z)l;WKQ7v&Bba%sgx#T}RE<9Nwv=pBoW!GdjSyY`!s~XZJYLmTMwZ}6 z`x7gHS^|*wDmE~F>Nz(cdlrVR-sScf1}1spxPf{qTu{)l2!K|$oG#BXGh2m}G_ zKk9et^&V3}^$;mS*3us21@hBpdf!1vM&|0%KIxHBeO{*Er+w?#*#V?n9hJDuLTnLW zbcsnxpz|$a-LYz0wr<_8qeJQJ%vOTp4asE-LsH;M%4Gv_#?`BDD7hEhxwEDFnBogq zGa;c9z;|G^fSxLDT@N(?0>{BRgCf-2*5(MIN+h@3${tc%PvA=FqGXz7R{L~xeqt0m ziZ&YV(FblQIDW>*4k8o)7Xm07|U&r-RXOi>Pp-=!~~Ii6YlEv~t2TfK3U0&dS5)< z7MBy}l(Ld>MZz!;ef$%m8JO5xO8X2Mxuk$6(ICQu^F?7Rrg`UKeSMUG=YfycAp?Zc zKqs08PC|4gemij_Tvm<-uGC$}%X`+sqKLqU0-A!+8w8jH3sH}|L50d~P8l7qh`k&T z1v)y$=xhXN0^$z}rJ4j7Lt5ugyF5HLDT#-V4^oRqnQqUA_btr#{#HfB4lYJw>sG)O z2+Dgh;?=N5D0P5r3MCPun*)t27tpp1!mpUvR$PlbnC?&70HsE+5BD7b%?Wj@fB6c% zMS_5UL<}DRof8VYEW*8alDYZ$waBso%~c*BMbLy{&PE0X?CaO(mMPg`sYxyX{G<+| z;)8`mReFdilLf5_h=SqAU<{FS0M=0gmES=8!op58MB+uE6$`yINd*O?H}_mabO?qI z`izg}8waYaz%|Knczus$&#F^DeHsT=)SjZ{59aLALbDV|kzlI;6#DYTW&CpktmXat z3YaBNLqh|+1Y6=pQqtJ_iZQ%5q-ed3Liz)=*$pvQvd%1&8b3LYO z5RqU5r#Xgs{Xp*^)*f6;FJz+?XWva=C;4qRRGxiT`bxCICr6Cg*x5$VH#%|>2C^(Tu$bQmr z_Gn11Iq{hGZW$5pT5OYLu?UZjWw&zy_Et zE#TuHlas{IJ_O%?AX^f3ZMeT*o;Owo!k8|oK7aXgiU6WrIHW6F>|ela_IXK?_<77g2C2Ho6NA3KLo}mcU4o;KBPm zk~!DW(sS^J?4JGmzajJ^U@nkVm=E(b0rSLt7$T;R(E)bt1a3k~N(!C_LElb9W<+TO z$NU$dj{gn3P>=u8-(PwEKG1EV-2~8U_MZn6TXBtB$HL+x;$(!fz@z|MU{fNG!l4lh z{XFgg`)6Kcv8;Q2#x*eGrI4?no&k$k1@HW~1<*S2V(-*n(Fzv?_WSu?K6|!`i3zW~ zytw$-y|e4LuFB{@APz7XR0g(n>l!{u4QyYscPS4Nr2m~EBVQghS9i@q$Z-OeO=?KxGgqTYNKuOF+`uVtF++$Hmw= z*x6Y$e0+Ss8lxUqeeNg9tcdp25M@Bs1Mt(|K3*R=EQ+lrC;%1E`nZX}AW<=TdQ+OHGhqI2t^r1EfryiH`?fYc(RbzO zC=DoW+yjgnNH^S65{U*Rh}tVL7qPC+7bOpKpK&ds+(Q6q;EVW6GoJb@pZZ+gRxiL- zsIi|tdp9EDC7Qf!_Ruo|Sw}#%uh@uaK19_7L>2H-h$|S-1A;{c`rpF60knWW5+Hhj zLKa+u2SS|w=J=%`1;ofdq)BW--UD|6E^PqJT$GYPK>FqmKU6k;f+{P;a>oRT1pCNf4O* ztPNuqmzFQ8bYZS6KOYl1&VXbF4tdejoy2c8HG=-UY}qm-I)Ge&(H`H1T0fXT@tHXj z6HY)C^z_K|NcUH*UY)6AdvF}bYTcHjPZ}Dwn%5EaO$5|SwIvdX_MfC#oM;8N=F7;) z`$F@2_%sxcD1Xc)SHORVbxClX9JuPrjN3pY4wQ?LiK&rk$vsenz%Ps-$_@=z_h!4$ z4LQQQ!NR~N^xicpzC=_B;VzUJcvMithujG!Ng~%phYpfj1iAPISpgiS%k7ax--T&I zEG|6IQq)YM@`^Tuj zhR0&h`Tw1JKTNR=3k^luc^uQ0`*xz~YuDgcE1-9n7gq;vEwa3&dN7c}sK0E#zp#Pc zWOU$V7hrbe65S9}+v_1Mky`nF2J{knaM_Tgok3i&R3ioV#BZq@W(QA!2kzXuB@YRfcddK=VCqU76(&Eiv2X8Bx(6BxY8Zces;X5$+x}@b zZdnpS?Z_oJvG((XK=aNsvoH$~D68{^Pd|HF$ z3%*(d^9vUi5F(;x05OfO&cNhv5Pb^O`-RYLzn5klGo z*#nqxKwi`0;+Al05@|JvNo%>d6fB>F0OSHS6y_4~;lO~r)eElBOb_LrpW}U;5=LbA z-M27w^#TAxAWvWzV>I61VxMjfn7DJjbt*`!pT}-SwEwa~tMxzm@5C+Jww*_X0Y?z< zzRV9FA`l}{6+<`g+YI*H%+K+YV2X8ZDlrguU5el%DxO`9#|cMN4N$?Qe*Y?fV))7W z*#TnM;tA}GzXwsoIPJeT4=w>hv;I5MTGQmu-3U*k9R3^MHT#9Z-cn%BaJE1eh_APA zr*K?G27N$00-Avu=NzFvfU?*K(G!PKF&kLn|Fyr5?E|!!vfVlosR$x1u)8a(tF5iAcU%PN z&BN1klc*@s1cjJ4If>2%B!L1dUa`r^H5iBr8VTfSBihC|4m2|NdP8 zNk0_dg=gl}vk>V021-kmZV`it^x1-plMD+33Nhgh~;%}HAk%P$A66@=K zLnVfa+ReXh>@tgXuBMTHsQYEV<=1NFYC=$ovb(*FO+sw!t3^5f_TTnnqXC)Xy=sx> z)gMC#E^SaHsjgF9cT3^SNx@9H%bF(F71yh7zi~BFZs$P4^@NCH#>U5tc4r)#l8uZD z$_%Hy-j|bQb0*#D%gNm#Az~#(y=pB!Gf(Y@c%F~;1P6@0x#*pilr-XEzE+^+T`AL3 zFkA>u1NPjWYa*x~0JaFz!BlQ|+9WwUVgMkzbG%C(uzJ)eJHga9@XH{g)2&?TspLyV z!UhI>xlNlTyl2G}hx4^(Q-H8X-&aV&mSZykf0V5G{s@k?udfgJ`ZR?8S%`o)UGlN= znvowTD&9Xu*%qelXD!VI*^b}tOXj*f@queNebE&Q&xJfPServO zQTGV}B+2qO{_J)^MO?CWs(t@~IuPUx!^Tqyh|s+57vxUBccWAj+?fYk!ehwT2a&*D zz!wK~$2cD|0;`8jV+KH7O+7N*lV?SzH?lAruhqhO}=J*bA_-lQde_;GgC|R z$H_G7`?g8kMt!Z85c45X9B8Y-Cc+(3bOedNcqQad{6{aC7&=EH^AXzGB$Dg@tABiI zIJTVOfn{-s4p>Cf+X%CCe`VKLaYO_P0#wt0WzfNT|IMxT&i(wCIPj1+1}Ndhxx5cp zYy?QC*j4zg#!}QafvI8o@`x#Y_#kGx<7g0&RZrjp4+2R9r}3_Nea}L+&xH$Uh3qze z+IOGU=`i=;K=)CcXioGQ-k@VF0;3KEAAIKk1BX+M@AM1|Tu@e6 zMng+A4dEPD1s9A&8jdCt|H;`!=EdIa$qXh#&6fdTlu4K3}pdQes zY$JuwaZ{3o1Tq$YMz0?C0A~HzSfRU+R4Qs&rSCpex!Mu8dc4l(-T4YSTH0Ye^lg-2 zaC(rx63c^HDB4K?nzUE5{Fw)cS6E#{+|mEP`}|jabB>5F;CiNsE#?aIi|I4&125Ls zu36@?d%3@A5AsyVg)XX0JjgODi6ur?#Xw*xIDaR zVPITm{$Yr6^)xp94P{R!uqgrz{ro*{dcu)laLZs)G4d%IXjZw8nHd`oC!?WpqcCJg zfNv5+9IcHzcOGlI^3V=6xL}6$!or{tgIWOJDVL~=iNfIe`Z*t{H5pr6!4)(AvM=UG z0IR3l2ndIE-wX>8E_d?Z_Pnt%O172L+~-8G5pk=?@d25_Z6OcDD^D`dXcP7)=OU7T z?b$cadVGif+q3elv%_}^%i@t!Lrae(1~bHLHo%=`hSY$p(<0fn?W%-vnU1&T#|O(; zuY`r6dhDdY3p})|eij89@A+3Bj{uA1jKRB5)Ru=fo_rEcPI>(F`Pz9l3%$`DQS)0K zq$}T_CR$KXG6a%$d@cP;NNn^fgWQB96hN8NVxlXeUSS!`D@{6{A-w_jDVpj58~^UC z{GHC>t7L4CmsXPq?aweDH88+5Z(t~Pk^-JzSXnDQ#xtECfl{mW1*OGs%Aq?E4Kr6y{iCQ0?z;$-g>?;;E58p)+yA@T2>LvEyDuc{rfx3xT`#M1SwBZpBTAHMl+nK zFeFaHn_uAEBblo19(7Pa*N%iRVv@>9DgaBZQ^?c1ngJEM4ZnJJIuH#}a3QT}G!RH! zgGb#KfJ{Yr0%dc>^OhssX`}$Bde0B$%V_&^^N}f5LRY!D83mKXeakaEZ-~z%iEi9@ z6M@n1iLDIX${SgzS?8A%JTW)v;vebV=k4F$%MVFL4e^r$2ADIUmIiYOMgW4Jy*fJI zr>7&XT!DtLvjT6INjnG#Ungp58M!wKck#-u=ae)9;d%rlYrN~FeJ7c^Dz-X%V7z^` zR1WRg1q`E>;2^|Z_?Uu;+d03b;*~NAJghMwDe{8L099|q!xxas!x>-r_0#d(t+u3K z>k}Ep@It%rFQuiW?!dx$b@%Mqylopt3{quQ5T_EPPF#BYsR=$9gkxY#&raRIfP&65 z>sY^~*;!bmH$U&vfecHDvwsDM>BvCM>U}SO6N{a+wC8G^xG9hG<)d~S^;4I#7%4QZ zmPZ5!Ma2 z!Xnp(1UPpj4gI+F5ta@MkE(hvj>^$OgQt;}$UJxnE`XT|$TXE#V2Lh?o(XQ@eyyBp zx{S8-V}2R{8YG6-p&hhf{$k}zXC)h>NrNLtt~8;djK*s8oOZOg1Ny#qzS|222nVeD z#p<4&Jgsf>!VKpRyN*3-XJ=;*N0SY#34S0fDoV|T8vUcG8(~HL@5T3!x=Te!;Mf=B z>e2>2wOI=$tMg*Dd*{}P<^`pi`OcnZ9YXId z`XP;GZ2s8Acm!fG>tJ%uE}$T`*njW=F@3sr%LB>>Ta;qYfWnAPTct};C7lB4X1ajV z6hlsrx~MCNAPVROPIgNS^1r5wYG&Hr-$#My9FEOQ;QaXoobjuRYTryB5xjVPeL{L|*mPNA z|6*$dN(~6xk$`nB=4I>(zVUJZ9wzD7lGPLVoj?WySW*~uwhl8h;CbTW&!0XW1PQK7 z2ZBkU4Z6AAo7{(2{=#Y!4)S&_IPkYa!UfdKypI7oMZ(I8UdYMO*#}9$slaLIL|Hao z;7NIB_4MdLKUSDsep@xJm>u`(Mms!gQ$?O_cv*^K5NqO2=UbE4^N@b%S7(ff7}~t> z`Q1M!nEn=Hpi6+2VL~Pp|J+7ipF_q11e4H-DBqD4)`Ze6nvp)jAOq9{Gz2krWl_8skTUxkpa)8aC|s|qJ*9v1uF!TpE=rgP;H?M z)TTcUL=z4~yNh}EO{NtqD$#O6*B_z(mU4zV%qmps!)&SsRbe4l~ zhd&WERdODr*$c+!uV&_AKn*q>!`)spq7NKC*z+GO7hrhGV`4tT-IQNI+V$OH5=a@Q zhC>}5^m?SO@yM%G&wLn=-pLzzan^n*ljeU#P2)hP5-EMv|82bde-wuQtCw{Lh}V8f zW*TWSm^q_R>Xf+l%muypaN#+Hj3mliA|@Wa0eU@H3`_#5HK+#1?{cfY0EAw+|_DN5YU*lt_f@zJR zIvw-|Y+ckpzjl@1xqJ78RjF+Nw9dG9I-P=}M-GZk{JwE}oRT95p2trwBbhar8j(n? z#OS%M3j5z+gMgQ(#{=WPdGlt^Pl22W;0hvu?Mj+f$%3DNffmsGM$D)PHtaVoRDz=* zBZw`464|nC`v*h9504&PRRgc?%Qi6D@orI0Ac)5QbI{` z0lXCefkUryJ5YW!DYFq0hW3#MH6s{kA&R?i!NtMYVH_s6q4rJ4b@3^}PDrHH18Cu^ zAyEw?7R)DjNDb(^1)l_MSHf2l67GX5gT4a2Ynn-q;y9c-#Q6S~M3Yv%&nK+JJ#<(_c$lLv;|I zLq7~1Fa*0qTvd|fciRE&B_HJYpp$9}$Qxp13uT#g+UW@b{ z6x1I>9iobg#+u?th}NvR2ulS+bKBvK7=0KON&5`x#kz^{BCm{rTe z%uHd{cW_kzz(#!cfK3wugQh(HwUlT#fyqQJ2|g$g8ng>hWh%kx>5&j1eKb{GvxfBX zM3;&|xqbQ7-;F^;1iiRNQ}c32IPXpw^yp%+trF@xV5A^JOCICOtyU#HNWYK7TTcjW z`%+R4$Q-uy0b{)noEKW)lIC^E*RE}cCMG)i6r--k#pP~gBYkxJ1^r0Kyy+9n$uo4} zxFBrt3`9MI^-WFBRlEnHmZA+hsu|`Ru)-4_GT%yiuy5&uPJy#Qv^HSk96(?x8%jEg zst{Jo#o`cr(GVm&R^T^dB_hJZeLnirLDmPF1}yyi+l1jG;ZJ?O#Ar{cqPvrhcq!?a znvP=A7_dBECT+Yn6#Q`z@1cB!gyjsIDoD6s{v)L#-W1X_l(z;DTtmC#kbwbk(A^ip w=tv%`)9C&n4F&(Fz{9_4@%|qoAt}E}-^4?mo_FpuAhPN`S_d=>G%POte;##1761SM diff --git a/old/Popultaion.R b/old/Popultaion.R deleted file mode 100644 index ab08874..0000000 --- a/old/Popultaion.R +++ /dev/null @@ -1,106 +0,0 @@ -library(tidyverse) -library(fixest) -source("Scripts/Functions.r") -#source("Scripts/Load_Wyoming_Web_Data.r") -DATA_TO_GATHER <- list() -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c("WYPOP","WY_POP",TRUE) - -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c("WYNQGSP","WY_GDP",TRUE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c("MEHOINUSWYA646N","WY_MED_INCOME",TRUE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c("BUSAPPWNSAWY","WY_BUISNESS_APPLICATIONS",FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('ACTLISCOUWY','WY_HOUSES_FOR_SALE',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYRVAC','WY_RENTAL_VACANCY_RATE',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYBPPRIVSA','WY_PRIVATE_HOUSING',FALSE) -#New Private Housing Units Authorized by Building Permits for Wyoming -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('B03002006E056023','LN_FIVE_YEAR_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('GDPALL56023','LN_GDP',TRUE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYLINC3POP','LN_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('LAUCN560230000000005','LN_EMPLOYMENT',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('BPPRIV056023','LN_PRIVE_HOUSING',FALSE) - -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('ENU5602320510','LN_NUM_ESTABLISHMENTS',FALSE) - -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('GDPALL56041','UINTA_GDP',TRUE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYUINT1POP','UINTA_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYSUBL5POP','SUBLETTE_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYSWEE7POP','SWEETWATER_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('WYTETO9POP','TETON_POP',FALSE) -##Idaho Counties -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('IDBEAR7POP','BEAR_LAKE_POP',FALSE) -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('IDCARI9POP','CARIBOU_POP',FALSE) - -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('IDBONN0POP','BONNEVILLE_POP',FALSE) -###US Population -DATA_TO_GATHER[[length(DATA_TO_GATHER)+1]] <- c('POPTOTUSA647NWDB','US_POP',FALSE) - - - - -for(x in 1:length(DATA_TO_GATHER)){ - CURRENT <- DATA_TO_GATHER[[x]] - - if(CURRENT[3]){C_DATA <- CPI_ADJUST(FRED_GET(CURRENT[1],CURRENT[2]))}else{C_DATA <- FRED_GET(CURRENT[1],CURRENT[2])} - if(x==1){RES <- C_DATA}else{RES <- RES %>% full_join(C_DATA)} - rm(CURRENT,C_DATA) -} -DATA <- RES %>% mutate(US_POP=US_POP-WY_POP,WY_POP=WY_POP-LN_POP-UINTA_POP-SUBLETTE_POP-SWEETWATER_POP-TETON_POP) -colnames(DATA) -feols(log(LN_POP) ~log(US_POP)+log(UINTA_POP)+log(WY_POP)+log(SUBLETTE_POP)+log(SWEETWATER_POP)+log(TETON_POP)+log(BEAR_LAKE_POP)+log(CARIBOU_POP)+log(BONNEVILLE_POP)+YEAR,data=DATA) -TS_DATA_ORIG <- DATA %>% select(YEAR,LN_POP,US_POP,WY_POP,UINTA_POP,SUBLETTE_POP,SWEETWATER_POP,TETON_POP,BEAR_LAKE_POP,CARIBOU_POP,BONNEVILLE_POP) %>% filter(!is.na(LN_POP)) %>% arrange(YEAR) %>% select(-YEAR) - -TS_DATA <- diff(log(ts(TS_DATA_ORIG,start=c(1970),end=c(2024),frequency=1))) -library("forecast") -library("vars") -VARselect(TS_DATA,lag.max=4,type="const") -VAR1 <- VAR(TS_DATA,p=3,type="const",season=NULL,exog=NULL) -plot(irf(VAR1,response="LN_POP")) -plot(forecast(VAR1,)) -RES <- (predict(VAR1, n.ahead = 20, ci = 0.95)) -names(RES ) -names(RES$fcst) -RES$fcst$LN_POP %>% as_tibble -CURRENT_POP <- max(DATA$LN_POP,na.rm=TRUE) -0.0157*CURRENT_POP*1000 -0.083*CURRENT_POP*1000 --0.0489*CURRENT_POP*1000 - - - - -# View the forecasted values and confidence intervals -print(forecast_results) - -# You can also plot the forecasts -plot(forecast_results) - - -install.packages("sparsegl") -plot(VAR1) - - - -#Check a VAR it looks like lags on changes to Private industry could affect other variables - #Idea check a SVAR placing limits on which shocks are first -feols((LINC_POP)~(WY_POP)+log(LINC_PRIV_IND)+Year,data=DF) -feols(log(1000*LINC_POP)~log(LINC_GDP)+log(LINC_PRIV_IND)+log(LINC_LABOR_FORCE)+log(LINC_PRIV_IND)+Year,DF) -RES -ggplot(data=RES) +geom_point(aes(x=YEAR,y=WY_POP),color="red")+geom_point(aes(x=YEAR,y=30*LN_POP),color="blue")+geom_point(aes(x=YEAR,y=LN_GDP/1700),color="black")+geom_point(aes(x=YEAR,y=LN_LABOR_FORCE/20),color="orange") -itial information and data that will be required from Kemmerer-Diamondville Water & -Wastewater Joint Powers Board includes names of key stakeholders that can be interviewed -regarding future developments, new businesses, and business closures. -version -install.packages("pbkrtest") -install.packages("bayesPop") -update.packages(ask=FALSE ) -library("bayesMig") -help("bayesPop") -library(bayesMig) -?mig.predict -example(mig.predict) -mig.predict(51) -library("bayesPop") -example(bayesMig) -?bayesMig - - -