From ab0fd34fc06fd4baebc7d4a3a905d507a008078f Mon Sep 17 00:00:00 2001 From: changwoo-lee Date: Sat, 13 Jan 2024 13:02:21 -0600 Subject: [PATCH] ver 1.0.1 update 2 --- R/bglm_me.R | 2 +- R/blm_me.R | 2 +- R/bspme-package.R | 2 +- bspme_1.0.1.pdf | Bin 169089 -> 169084 bytes man/NO2_Jan2012.Rd | 2 +- man/bglm_me.Rd | 2 +- man/blm_me.Rd | 2 +- 7 files changed, 6 insertions(+), 6 deletions(-) diff --git a/R/bglm_me.R b/R/bglm_me.R index 48fa0d7..0021a13 100644 --- a/R/bglm_me.R +++ b/R/bglm_me.R @@ -4,7 +4,7 @@ #' One of the most important features of this function is that it allows a sparse matrix input for the prior precision matrix of \eqn{X} for scalable computation. #' As of version 1.0.0, only the Bayesian logistic regression model is supported among GLMs, and #' function \code{bglm_me()} runs a Gibbs sampler to carry out posterior inference using Polya-Gamma augmentation (Polson et al., 2013). -#' See the below "Details" section below for the model description and Lee et al. (2024) for an application example in environmental epidemiology. +#' See the "Details" section below for the model description and Lee et al. (2024) for an application example in environmental epidemiology. #' #' Let \eqn{Y_i} be a binary response, \eqn{X_i} be a \eqn{q\times 1} covariate vector that is subject to spatial exposure measurement error, #' and \eqn{Z_i} be a \eqn{p\times 1} covariate vector without measurement error. diff --git a/R/blm_me.R b/R/blm_me.R index 3e1c309..8d9cc72 100644 --- a/R/blm_me.R +++ b/R/blm_me.R @@ -2,7 +2,7 @@ #' #' This function fits a Bayesian linear regression model in the presence of spatial exposure measurement error for covariate(s) \eqn{X}. #' One of the most important features of this function is that it allows a sparse matrix input for the prior precision matrix of \eqn{X} for scalable computation. -#' Function \code{blm_me()} runs a Gibbs sampler to carry out posterior inference; see the below "Details" section below for the model description, and Lee et al. (2024) for an application example in environmental epidemiology. +#' Function \code{blm_me()} runs a Gibbs sampler to carry out posterior inference; see the "Details" section below for the model description, and Lee et al. (2024) for an application example in environmental epidemiology. #' #' Let \eqn{Y_i} be a continuous response, \eqn{X_i} be a \eqn{q\times 1} covariate vector that is subject to spatial exposure measurement error, #' and \eqn{Z_i} be a \eqn{p\times 1} covariate vector without measurement error. diff --git a/R/bspme-package.R b/R/bspme-package.R index d1fbb0a..0408c32 100644 --- a/R/bspme-package.R +++ b/R/bspme-package.R @@ -20,7 +20,7 @@ NULL #' \item{site_name}{monitoring station name} #' \item{lon}{monitoring station longitude} #' \item{lat}{monitoring station latitude} -#' \item{lnNO2}{natural logarithm of daily average NO2 concentrations, measured in parts per billion by volume (ppbv)} +#' \item{lnNO2}{natural logarithm of daily average NO2 concentrations measured in parts per billion by volume (ppbv)} #' } NULL diff --git a/bspme_1.0.1.pdf b/bspme_1.0.1.pdf index 967c064dedf8395f69adffc1cd21414438fc7c29..2c91986e2cf801221520eb79135e3a36867699cf 100644 GIT binary patch delta 6464 zcmV-G8NcR%sS5n43b5P-0yj66fdVLh?O0oH+_n;a_pcE0q&@Yrq(t3-lLMq_Vl;4i zi;I&sEds90-DRuGM$|fSe|={()UG7U8|MHyPkFF8q=w{NKF(#ixAC}larVvC+1JnG z*vnWJ2eEgx_5|lFl(Cn@oJC>aU0r)`{9y6N)r+s6C*goGOe7OLf@3t&2O+om@=AYl@HwzK@I?uO8z0GT! z+`x@#@WuA`dg9?O$Yxyt!6+Jk1*3$eQ3|mrU_v;zd{QmJousx_*hWxWWv|6y0#eX|I0N*lH67nUeu_c8>_qKITEb=5AG zK?EW+%|b?g4RUm80ri|LBJL?l^+?IB$ zX<^rBtE=42>NR~ncm@|WdX_xn*vKmts8;{ZBW~fJpqc@o zKQwZ}!tj8**TgciVa9lPC@SFbhdlFts5D*HX#=BEHU>jjCz~mMcXiXMUevW#EhuEC z`;1!|I7BAC0ZUsV+pX$PbfYSY8r6I6wyMrZV4|3XAUfG9CV|1sjRo0621r+JV?GdV zs@=XNva=i#of1gw8q=?V^KxT5xu!>e4tEV2@Q`t+s7@FZLPVAwLe}ASo2x>t3MWEd zmmr`*i-668EbyOy*u%BmNExd+$|e|R{G4nN2nk|DW8BnLaf>BgMw%iq2sAJ&>W#6j z$!XSfsN8_h=nhv4X}mRaPDA3+x8NDSRi?0yXvOOO?Z9NE-SeO+1%@m=?(c&Ykk^BF+A(#G!Y`No zwYsGd1SElfP@jMYvhzx}$z1_F)vRi}U19QpMJU{LzUb}iK2!_I{TtKNRcFM)E$_^= zE=>*iyj{c*rJiZ|fUzFh7#<6{^%tu)dNxViO&a4I-Y1Rd4iG!-kTV=>`0M$c-mhVo zs_yxOfo1m7Z?mgwRT_M~Sb;sJVP&eigcO2;K`#1#YxMCEzUB52`e=Nfi;LB2-&CU= z@W*GbuB^Jz4`hXSK0Z~<^0@#z;QRiU%gb3-y!`cu_4$*_SI;QS4PD_B8^nrvco5r@ z28*iA9${urqU0m!`1Mp-m&8!j1${)jfmd2L|29?opP6VT@zA*EJtyS5cd>HuNw}%0 zR3EZ`AUyBt%DRW-+8DTI5k-?4eWD1oKM_kMsuMF2jxhVlI#v5}PNxX!KrRM4g(C3# zFon5@lTKOIL!ygMx=ToWq|RLJl68#&?;Rh+jteSKmL*||TruMuA;M%P@6X8z)T`1P zFm7tmT%jDoIi%vI(YLA%q!ezxopt9A_j~ z{EaGaFad-{2Mg2osjCNw9;LiOR~u{rkbc)!tJZDmE6{mt9Zhae!^h!Zf*uTO4@VD` zQ}?1LX0b?y=z($&-4UdO9vI!!p@+yr#fDQ?vCY5-qCuD(9^6ke*nu7{i9a1XPSW>( zTa#}|h+)Nw^S&xY8j>f-=)jU?Ov*`eSvPfgss(-neoOq92R)C!anP5Y%Kw5?h^2cd z9HMc(M>LSCC*IL8Ky;Lhx)V!F1K1`Gw#h?!mj4r$-L&(4uKx!h|MuB190XY>UJ#PB zFz#sRdkduV8mb>tQJHyUxzpQOT6aBvQXritVX3NHB95*?#@gMbm5WbKz&r72>adTd zlQ}cGs^!79%nN1hX{>ioX71Q&G8J5oLPv^#3WG__i^m~|@B?*4u;AM2jv1iVm=*qV zvOn&ic4oMR3NTzcHnO;m2FKNvCUZN%M&{nz{wIa;ZwH8AEYV7F8fkrj-!Jy8*cBWgGlBLHV zj$Au|++R7*QP+UvzuZdK9L)T{Bv~v6zc9@eu=jO5dZGK zLTL9v-fEVVs8_cRXxF5`kapLBQkjc)G{$ zj%V(jdE9$Cd%T?8KTbj~V_6);-g4~;&RHm9FNrye!oXXedhh&IQGP7-{KNA1`^RxS zl9oJXNgRMKCl^_f87Jj`_m3s_Qjm?Mm0g4)VLTMxB1o(Pyi%>UO*LOcocj-Pfq#C_ zDgq@7!C4k%B-*M={gtBKx!Kn$-t$N}(M2LyV$B#|Sex#Le6#Hgx_PLt-Ik?lu0RL3 zlAF5ibYtpfPl?p+6bP1y6mQ?Xi~b!| zr}obCy2E|5t@!M<@aZV=25Sl`IL8Llli5lE}Sl0n1X+!|Ab&j8TxtAPgarxe&r<*GD&s112-^EzxyRjH8KKj6yIWrYJtqZQhtoXZ;Ae zVs+;li492VDyLn&#&wqx*f@=t5)#97=*3hXTz(dTjBQ5h2FZ0#}Q@yEc ztvYb1J-lash%-uvOnfsJev4#xu6m8Ms-mc^d!Oj4!pl{W>{w!fqo7(U+Ss!*$f zEXeB;>{B)%5b_`s{&624r*OaDNX|r<_A0T{*sSk=lD=oaB8)U1$O76gDA1Vjs#Rxtvfq8#|Mrn_(EoPOa3QmT;eNfmxi=Ne*KkW!*BsdrnLK(s$-Qmn z8H?UB^Nj4Cn0Y2uQfZImp1(RdnH04jf4wJvdi?0*#8%=w%)e5 zhkm}PtCoO&X)D2fY90};uXAQ0f{3Nkx)d{!#$%X$vtQM=oboG*Sj5H1uTTWeuPk$a zeuW(Uf-RxKlFkBmHgad`;mjBeY`ymPJcvmOhD+%0(2hYWW}I7;z%(SUe-IKDfs-R4 z=teS1{lCsNnF8(34F+JJ>WK?HDz`;vu=fD7A0!s^Bqk{n*#WHS6^4Yh+tuITAZ}b+ zhm|G|8ZvCdazh5RGJ=L7e%`ml73E)l;Zo@8%nhDEUl7Z*bjdyZ*Wzmz^7dx<5DpH| zV>{8;qX*yeYta+4SR`Zgz${|Z5u}J7SZPb>L2?)u?3Tub3HU%m2$P}cUQA#oVQCV7 zJ9fy@Pv<5-MNUXY3p*I z>;C}cx8DuJL6AA~f)K5Rb%RSk*jn0ld1We7o$ZNy@x!7;J`L!Bq)Y;9qPh{X1{pm1aw>~FRVZ=dyX-T3c(Yl_j4f?`H4e{`Q;U4_TD8UT<(H z3O|@4hvdv~=9`ZUgLM!FuO>e%CeJTr6eywu6T#vppkYWNE&?e+5=Sgptb)(e+Zm0g z*>btb)8CgBe%#J}Tm1a|Qh$ivR2~W<1%rk4B!wcA(C!4XaZh&cgy&3RrujcJLPn#Pwf4w-E?d7Cw+ z`PKYI)U9^kfTaDZE7b7Lil*htG-C62i-T#Mr`JW+{Nt4QA=phxD1VgE98}JR;(be5 zL>P~pW3{y2FfYyjzkd1f)33|cT7B%mq4Vu}b1FzosaX|W@Oz04JD~`dBv#ZZvaQOz zt?Im7TQ9b0n^#t+v8DV7Fl!klYycE6`^5G%AxyHfwuhD0dAnbo*1<05>u~HmPX|H5 z_&gP${kR5tIbVVD?SC)qlnN0XAMJD~({^8{#cUo431+WnItJU#*7Q_2_%9%8^5Wif zr85@VoUYS#W;$_C-B?!TGAr9!Thf@x+bnJNb+$4R70 zJ(>a5QNone|0<3hrhpPvv3D7Y>A2fnYjDvChYq+{O#VH=cz+86+Cft+6^SAmEVq-- zzlFgH27eAG&73h>2_bV zcF@8o!W;^eh-ih8!Z3OB>8*ie>{#UuxZ)b4X>YuUQw-kYxT{xoP3=LHIC${+CF`Fe zQveFAORq0sC`O8%Bvf@W&tS=fP1{Z0m`!i?{wU2!8t;~M(bQ+=^KO&X z*_hc}IfLGg9K0gop|+v6BTGaik-gAh;=37Ow{$Lqj1xr3D6!Ql3I$9NRCKZA7L9s; zf&$`Ax<%=TW*q#YAZGe2xr6miO+bcs$C_L3n3IA=7hlkZ9Kr?$ZQ5Grnqz}A%Z4+d z#I5rZa?l)v&x)*WhTb`GoD%8d`Mi8*l zW9#9@DcN!D&P;slV+Wn6YPj5`Y2jHs6Co$6|;lB!P}z3`ssM zvU062+Cu4gl8&R`N_Z+(4)n zs(GG74pzGm9R~FsTv}}O+nznFPM+*ugL-E>xX{FYoi9ycJr2`YdJg8+FI-F)@6nqd zFKl<6!5y$PL$9w)->S#70ZxPT~aqJF>+*BncOTJ-ijXJsxINtm6w0 zkxM3sRMB9o!Igmv>5(%yt?^qAo*&$L@KQ6k=)ftVapDKHz7;X9-18P!2d1XKM#~WN zAFdNhbLx>(g!^xNU~XqScWCi|j&kKL4yw;LS5~bDoS_EJphD??3+ae6ytlps;SIeD zRaeJG*_Pt^7Vi;~Pcy}Uu4li|CQZ&VCecLy`$ z4CPVZbWegK1-w@mEN~Cl z$yG#kFhlj_t zzgLxFk5(!ASClD!C_pF5l#1nP&=!Ylv<7YAUIL);07MI2*1gi9l!ahnUJhId(|ao) zwQa4EV!T$#_<%+r^FXZ%p-mqr0`Y0HW$agH*w2^^N(^e1e6VnXZ088+-XZl0rlW&h z=sLizeI*wpiH1e3D^o{R=VCGWAFc%634=&Ohe$&Kw@5<*&N~G&FflYTmm7Zq9tATn zF*Gu_K7Rs!5<@UXH#tQ!K|?V$LPI%0L^LxwHZnmoFfuYlI7Tx>Hbp)lJTOK#IYl!; zLoqc%Lpec2G&4CiGC?yiGBQOtMl(bw3s*$y+Mb>E)bL zgz9(+r&3BPMB$`VN@fg(VVE(~Of!ZBYl~~f-XGx~Ft!%9{0qjy&Vu{QXUFC9wA+X8$9Y-96WS|E;t6AZZ%(Z(nD~VB|UHySQ5Q_ z(DvF?AM|T*@T|wJaa9{=9e2ylLdh-nS|5LOTwU-qoB{pf=b`LYe6Dqxj7FweYRDsS z0Y=?QAGB_a6KJa$7CbLL1#ZDSEWmA8hP$xj zHdfJJ+<_HXg?k)GlcAGACxa#iO$_=M^e<>#(7K>&p5dC(?}7zK@+3dA^2>R{5Le{segj5`>2Fz%T4Y~yFz a`Te~81D$?ug_l^A0wD-63MC~)PeuwrM}D0E delta 6490 zcmV-g8Kvg@s0x9p3b5P-0yr^~u`w!t-H+Qg5`Xt!A>>KEi&;{lzJSXCX__n=Y}#wH z#idP)fY#DBuljIAd6T`rzB3ximb6-Db3pDX4>pJ7kerVnhcn{d*5lsA;@hjmH!tGY z%UBi%v3Irc1m`T2v6sZ0MPcAwU3+i+VEM<@%WqyJ;eatrBojP>Z910qwkY3!l^Pir z-0eQudKHE-3uNN0f&>R?;ZCcf-MzD>T&_YPe70OgoPTYGjm)xBZ)jl-|D2^>gIywj zzF2r~S0W6#pH6tGg};>u7*Cu7e)utX_fk~>7Xh3UgkNp``sft^RvM#P{X37ig+D<#13-Uh zjcrX% zy{1Fu7JNo$;F!ZXlIS>W?O~6O?qU*jRasYCA}=mq*~d#vGJRFK>7$!FyNqIIC{vaR zG09&My(jQa(pR@)q-ArcEL|~0&4pR7-H27?zR-=kpf$CBxoVnQBI|<-b#N71O(?*w zsW!TCcNn?{Yfwid|5X>aYIUK?(w(LbTfCoY7+fL_Y_{*~cJ5NEf}yP-(nH=CqQ>~p zMv?zo->hq zWiN$p-Av(ss-^D+f7tU$^9GC)KzEG!VT}u zwJuE!aJ^l|5xqpyYyv|+R5d&nB*rgRZT0M%arX@wKk5A&65YXNr!I1aW0ik>JZFeE zuu4_;Oktdvo%(HdbfZdxml+eXC)BJ>RhRIhpkVNSH2pc6e1uoKeS|(5@8`l~wLUb} zXa)TE?9r7~H~N9B5Z=eTidjAvZwGuo{Bn6Y%Zi7;8L_^2dinY}MZKXboMMB}F%J)7 zd)i=*wLKur>{*n20v*4e%kq*KlDwb~a5r#E>*gP(YX3PC+$3%q_q>k@`R-nIU9BK8PI`B%~}$!W6k;#yLWS$xPn=L?@8G zN^ii_#k$o0_f9)eAbNJG8px-5;yq23hoUuq6w9(xYPWJpN|G{>O)%~9B7C3@>;Hg) z;|I#+f}=znOaP(L!NRotOV$HKk6yq+S6i$Ikbd7MuhuQ<6ViEXolI`e!^h!Zf*uTO z4@VDWSNEbPX0b?y=z;7I-4UdO9%xlS1{_TrlI0&*% zydWfLq1{u{ml{atHDo`gqO9}8a;MaPvo!8Xra(GR!cbMWM4ViOjI}#UD;J-ffOq24 z)M1}|PUg($vX%$iGB1?1=ds>hnYnGJkE!5t6gpA_R2WQhUfd2rgcDRG!GNo>J7$1d zW7hZ~<#1X-)y;4TC15yos%vo^4UVI0P3Crjjm&+p{TGJurvyYWrs|B!glC6;21=5) z=u+Ux<)KQC9+a*+v&T$fk|W~i6KV!I;^(`T637R87Lp!~m!3XtH8VE_!f_Hr)BNX9 zuKV(>zw6|TrdsvZGElj7)rCFPU)OqrCn*SIPcU8mw36T!mhi-BDTlf5m(()k$w|vB zA}ygE7NpRnLKYRw;xh`aq>TJXt7Qt-qUs5_=XQVf?A0^(;eZ4h<h3z-LE*uJGB5LwU88b+AQC^ zzN6~Yo_St(IB&L9e`I#FoL$TWoaCMu5*%}YAx^x!oW1+Ny;Il(<1EZl?~_w4Jt;s% zh9DQ-%h`(|v4^8$5#)NrA_y>+!i!SLk{~%`(TF06+&dPqEEPSS9$RFLf>6P`VE?Sg7neN<5(r-7aGpfL4JrB758L-`wAv~y6zfkI` z8*&9`M-h6Yn`SN}zi#G22%lY#ZWIShX5w4WbtuNs#9fR+a3LlXpXfGk%%-#Pj-lK^ zx_3i?Ckp|WP=t7Vd2TE=Y3#4JRo)?LFniB=(9w>As5E}eMI-xQr3IoT-%Pown@;tv zuC?kQp!V>bp?{T8LS*8bx$s*QyK~iBq)`<`ZNvLSM-^5!Az~JS>xflMwQHlA7JL(# z^h?z><}=Esf_jIdL(U2=2rF!C&FMv*8_e^e35elKKm#7IB$YP-YRg~JD-0iRb5*ES zK_29F3H~Vy2tXcW!awfuaR&G6jpR&(X>SrMjm`Ql>3@3;EW#*cX%_6W3ToznCk&DQ z)T~yt(yDS(=myujnuw;F=8Dp}?Qn0X*p#tqtsBezpD~xM)_Z)3{FSC$;`*H@x>KfT z@8E8$2TsI?szmt6*_yf~bR|XXpYE#0Jy;1sUsb2Lb>ch%L3cQ+!U3N77C^|RbRfDd z{G!#VXn)*n3X>~xQS|rNfrSgE->7W4zBEl;+1sngF}*RTx-@lBpIy!42y5y=9bDKV z22BIC46j;srYCFPm$h)8B?q-|M-zvhRSfs<&YOEvxW0x{s=8*&UYE(Ox09>4^*>|P zTl$}o-xK}Mq)EEwBfaOZPEICO?Z;p5sUAN%d4G9?g=Y*M>5p#S4Y_;-FDX3YTT6-d z#+=@g%QrXm6-f+j#n7c*!j-C8+LlvsMG=d*7{wKez{Ld@ z-G9ava*PY6#0o1q58U1;oT*1NW6-?y*5C6WCM{?*q1Hob2CbNJZdn4Ik-YwKNLU0x zj+CGp$td;zI@jb1lsnfXfPbnNE^w>d7M;P`1B-oAv7qNMNtwtFXiYCNB&^-8{)PZ? z;z~S>)P+#P0Snl0Kq(_?7+TL&PZG+%!hfOA)tPHTLB60R)6yZ=1YFCn9pvrJ^dTM` zu*X)SuV)Xw>({a;X0b@d?18t;rXxs^JuuRg*n{HGI@l?-4iof&f)FOdrF${KorI-H z{O#N!Pd}ZT{2ZMaw@!B6w?|fn=M8)$vSb;P@}Q`#o4O=O;g6+WbLREppHFArJAcSq z;y)YW9H6J2jmrOl$IJGRp&&-v^`6lns$RIjVZ`Vt*=v(qTpH0fA=)Ni(>(uAT&AhZ zeXjolk>7qd4hKQz+=Fh&c?;tPhrY?RvaRyURHi!HQ~BaYMvLpztOt@Z2}@Pop|DGp zRjlb-T2Xx>2fwIKhY7p&IhleH8-I(L!RuS%p#03 z+wc2fioX|ZfLVM)fC~xj`*NNFTG$eU_LTK!zdid6NkR_okz@dr@NFl4A}UGA(%sKo zYCEm=Pn0?AF0J?@t+*y~;^(kru^jy#p6c3;ma~5Ww&Nz~vtT8t6#+Suu`wur%~(;7 z+cpw@_pi{S4aAH{QKWcn4oK16pv@+|Bu-EuDGIHnwbj~^Pm*`D_i(@chU8EdWv$m6 z9E!pZrpO^VGo1P6Bg0@7gu$!H4~xn3OBn@@Md3rme z@ibp9*G2mKveF;7v)>j!Kfe@zqBoU?f=Iz&VKYf-kx7Uk#nq(HX9bD{r}f5AAb6E*`n3Q4jeik*PBy8VoL3*+6BLt+F>UY;gZCPIz_fs zS+rGMlq>XNo3=%TI*l#mM}S?+DuDs8f!QavrwL(_owYrzv?|(tc3KBq(AVMEd7ch} zgzKkG%yPa0<=bC>*eMkvI6m6xP^Rs^PB*i8#5tI~nwc1E*IV0D-Pr#uEjIVI zna)^e-2`0La7s`O=UGPJ>u5LRo^kcHnnA zI2g&jxs0f!%a8YajD$}!N~s1=Z5*tF@Uh#~GY)|Fwb=H;W>wtg*hE1Mxo}__S>R}F zAIzwWuP1!OK`o|#p5i0p3okmm#UM{;g5mjN-Yn~4XB7Shlzvrcl(f;Oqa+sE%qWhh z+x@1+LBvswb6TL(khZZ>I!xYtdTTLRcVtBauDI6Gw6|W=Q!L=?aaXW#P2)i=aq!^t zOV&R{rT`UKmtSAPP>fVNNvP^%be1G++OCVnZhEsXmj`8kU$K6yuXlGEwQJP0i*6g)8ekx@A*yF|+66ZBYPE;5(Uhtrd4?#FNWw*oJ`r~3*J=Bo zn|)(SZsbaTZZv(hFUm4^-3%(ZHuahPyj$mWK4vyo&Y-s=2d_wYC~jKYktHIM$X;kL z@!br-TRIm)#tEWilwi4vLIG0*6!tUI-$VN3@@|4>*Zh>V<*?(h}fxu|{M~H>(Qkr{P#g_FM))(;;>1rtF%0pH<5$ z!*7OiH}?OsC^xyz%zOKtw{9og+fdPvDy=U;g*}!a0e=+)087AG((3)O(Y1RFfsYZ{ zuW$o@MD7(gt;mOBwM`+B;!szfs?HkxC`p+YtMzqNyB-*{jTiRs?Z>Mt9D03~6}R5t z<;$xN7rhf15Ha7lm$fc=DDQyQ54zT9vBITARrfCQ91 z&^JUE!%yti@5VGFG^4qYU4^r5U|o?>*2~N*>Ua8ZG8rpT^lkBZ9|@QpR>ERLwkUJQLdvYrwb0xZbWfAF-6jTAQ7sQt2lWzJeBL@DRCVf=Cq&wi;X+xR4$>gVP$n_2~0|gIkZj z)XWhbI3+Yr{Gc|sBF2?_-{R`P)Xdjt8G`=9bwU|VJ#vb0|BVmK?QG``EgsNOuH3~z z_4#IoYCYf#HE;$M%3MfCoZ&t89SCpeU8uS`Hp;dXH@A3?kbIgc283;&f6JYBQ#jyp zl(3_YH`Qii9~A-Fp7Lh+?&KqX6$94Y!OS>AdDJ)E8>8I_ek2L%Rp6TMUrmC~M&`{i zslkym9wSj#?=kYYO!?;Nw~elV_v!)y_kf)|t%we^vs?LGY`1TV(+{N)oWAM~z)0J3 zPyaNu#N8k&d{;htDf$LganGZIWFC)BDhru@04DRgaAWhrr`C*wtOGcIB=}mz7xr)Q z{^P|M*W`eUv4^lwts`5e;Q#D0C7>Q29@qX}Rf;`YrRZN#rud-%ohVZ(mZw2m9Inv@ zw1s;KfW`w5Em-$Tr=<{rg?%}2CCu!teALETCB=BHlJNnJK<0s36+)XnP6Xo9X3N;G z&aj^~8MXf7SM^)!yG5H_$Q{DlC zLqmr{Ljku#Ljuk_1u``rk4omARfpgKg=4MmoRsV>4z zmbw-P4_%=bjzW)H{kl$i01mRH4~_s^qMsLPucmsTuff5yKDXvaYM_1GEjtZmxBNT3 zf6#Gt!IN+b^uR=Gp7-YCx`nTm`ubf8rCwClnuM6belvG=WA&{(}5PtJ~~PwQ6(w ztnU)-ZtDZ8@@{|ssOoSpR#bJmSH7w0a$oqMs^GrXud3Vq?Q>N{_b)kBrF6|vRXv_{ zzfsleSv8mLG1roqFC-SLiRJf+$4iN)gNc=~#Iu>ii&|p!Nn-7#8W4*hIl;_>CTzqQ zXxu~~CV)@}lMemG8Fw)5VBEpDW7@M