From 2b096f508cd462a8be3de51ba55876cdb4945dcb Mon Sep 17 00:00:00 2001 From: jgabry Date: Wed, 23 Aug 2023 09:24:03 -0600 Subject: [PATCH] add missing link to diagnose method in CmdStanModel doc closes #829 --- R/model.R | 6 ++++ docs/reference/CmdStanModel-1.png | Bin 15795 -> 12728 bytes docs/reference/CmdStanModel-2.png | Bin 16009 -> 12864 bytes docs/reference/CmdStanModel.html | 44 ++++++++++++++++++++---------- docs/reference/Rplot001.png | Bin 1011 -> 5039 bytes man/CmdStanModel.Rd | 7 +++++ 6 files changed, 42 insertions(+), 15 deletions(-) diff --git a/R/model.R b/R/model.R index 5fe9f70fc..37243f0da 100644 --- a/R/model.R +++ b/R/model.R @@ -190,6 +190,12 @@ cmdstan_model <- function(stan_file = NULL, exe_file = NULL, compile = TRUE, ... #' [`$save_hpp_file()`][model-method-compile] | Save the `.hpp` file containing the generated C++ code. | #' [`$expose_functions()`][model-method-expose_functions] | Expose Stan functions for use in R. | #' +#' ## Diagnostics +#' +#' |**Method**|**Description**| +#' |:----------|:---------------| +#' [`$diagnose()`][model-method-diagnose] | Run CmdStan's `"diagnose"` method to test gradients, return [`CmdStanDiagnose`] object. | +#' #' ## Model fitting #' #' |**Method**|**Description**| diff --git a/docs/reference/CmdStanModel-1.png b/docs/reference/CmdStanModel-1.png index 27fd88217456b4a7d0da11a09599ea685749e06d..0aa41c0e91b00f182931ce9cafb34ba5eb0c92d0 100644 GIT binary patch literal 12728 zcmeHt2~<v3Dx@%Znulp|GlAL?a-DmG_ z?{DvY@BMQq%dMMtZ-&8OTajl@+QDG5kuaFd)z3D9ZwAvI&4JU#tJan$!5JK)6E#?H z+7xugF$4xvxd;8OOCFo|hQa;CKdl9(2k`|NH2bmL~@jHkjFpn(9<(!4YPL z#iwndq3dP2jpLJWtni)W3EV+AUmkZXot~Xda+8j{ zw(f(LfLU##O3(Ql+@N!u-FjW#<;&smK~ceP)AZawsw@MR>l%zvIv z2Yal3z(q%WbVf(7T&g^6%cLiAS_ZL)YK^1hz>MeHBiC(vxS-3e;u^DAq|>_0K|LBp zk*iCcU!3gibOx9Fkl7#~^stCOa>P6Ky*~5ptCiJa&igShU&iF8aS|BsIGf2syrA5U zC@rN!bCT1Gq!JPM8+4x?-f6FDvG!ivGogZ^)RB0#E12%FPSb7Q8H@NL#qq=#rJ>ma zRje0jSb9kJO~FpEz-7x3S%bgL9-yNc{d(iRc~#S6Zx+QDaqh~sdCWx8VE;vauz?Xt z2S-GtxMR)H=>I?KwP0kPQ?6mzewp`vj5AO?gov`N85S)q%KPB*_AHBM%{0hVd@G)9 z2=0?@1MbtbbSapf(y7(wY*HaEF)2Np6qxKb-^oP_AnksPYw56PsSzy_)kcf)R! zn?AbZ4mc_`%Tkw8Zi(VfKA^FaA%gJR23*`|h)aA!FcY0C%Q&4Zy51jw=hbF=rEJfG zXe;{^i2G-PLrs|j9S{>dOZvAk`_HV}|1SCc)4-~sF)T(Sb3A%JMScOB$SPx)jdRSR zHh_4l0S5`d$~B+Vg9QiGOj0bN5Ym-6`1Ni;!6^VC(@$r-xt6v*yK6m3Bq5Yk{E)8y z3%SCdQe7Z*J65y6CL*dXb0#KKrq0rsy%S#mSAk8--&VAUzo*-ZFx9}8n?^w)7j_TV zAfNlUtBlOEv{4`2*j?4g;L%Hb2|dpN8hxg=!L3>zrD9WQcx0{)ktof<8e*9p@J~4; z2{0RU{8ri7pHkz$QD*vod;hJowwuOL0GmTvKDoK>Svy6fPSaVhtOW%!LSF(BZ$K{= zfL`KGcxLfA^TmEKT+W30lxtm%3ni1B0O;qKoOF2XjUGDkjtp=103|p?01{s(qdU?h z^B=A~plS*hbjWPJD_U^pDJ4*lR9l3b>D&V{o&boZsjAV0ywMYothDDkeYfuJ^=69+^wm25ktZuGMi8klP7@S`t02$wT-N1l!iuf!f;tS7K?$zHi`yf z$9>scg~ty{!%Frm;~T$+y#K3c5rC?}#K%1Y8=bXKgDmKD?#o)BbZ&?AYDC@OX3%$) zK~@(fZXaZ-)OwS$YaEB15~Zbe!xi1Jgu!b(NpU>FvaTUV}cM0+<1f;HQlv=L`?=KBSjuuwx7z+mfPX zV7uBB^&r3dy>$JF$2zZ>#axoxo?2O$M?FLv4GLq+ zvolloy<@K&_ExRo=CvUKw-lqbup~~qC2t+1uQut4%Z-{klLkp%2$JJD7J&V1$hGcH4nc++Vz;7s&K*Y)T9A(#` zXba4i1<(HJD=0>8tD&21&K`4F*aTD%!Q>cg8(ATRID5;c045zPjtCDAUuXd7#xY$8 z*+$QjwDIUJo@JK`IotPN+m2@m4hl2?%h{)CSCD?O=f*g)^_8cdtm;VoEDc)sPl<9xAtm zpF1DK6>qhON8`ew>pg7u5QWe@*xLtC>)^l8ME)NMmA|LaKLr<7D4PRAj654~^!Cl- zXD7R4V8aAyVO(ROcS`XMMt-OBlAf)IG3A~c8vr`Q)&0<};*N7QyGy-MBomxH2;U2g z71|MXg`77M=tCnacMf(wmM;@!DgWq^;3$O*#8)!I94k0+jqh`#pYqeB;iEefFJcIL7CH z*?!kxw_fxHyg~dVeDt|`{*8=RTfP|R-H%eTy0NMIs#S~S?#&xL=j7hL{b=udLs%#A zBwszns#$)E37e%XORkl)GtI}TLPuXAsk}W0qXUJ54!h0Len(uI-6upmCwvEDCFmIC zfAy*Dp`8!9kKnkjDf@{ytu>|52Ige-sdxX!F7aftK>7f0^vV=wsszQqV>51)laCYe zsQJ-PCY52M$D~rhkyaXND4RlpAd8`QhMfKc#s3aoWk8bkP~H!TALo6PWNldBLjBs_ z#!`u+rQkhJdpVk|FcKP~@E;Z8d<1O2;jg_fs0WX(MGpE`4EkkV>>(8Fg#`smm2~K; zyfvNvINBhiD!5w}sQKbs-Ht~3QozUcP4p7V`@FA>I%eK5q5s9! z41Z9WLC<1zhI=1&)*~sp1zig6+hur}9dn?z$0hw_qxIgVkIj(kMktI-l#Y4*et&yKY>ejB>;wC~^w-`Q5Zz6XGfS@$*e7MB0Fe>XWTigp31 zGyFCb*O+5MFcXIwEJiUOo@_LIce8DW57|0~MbWhJQR0-wgCCbpV3bcGX$$j$M%aS! zsB?i$D(78cHmY&*dB6c?G^$&3qw}ni=_+YkAIs8LdB z-&nJekJ*?0w=_Lhe3rz7df6613N^C9?zh^G1_bK??(p=1r~|1rn!R1^5>l6Xjv?r&7T zOUR>T?IgOxpLnj;yi8&5pg;0^!Ev>W{l3ZCMyE(=jYorP(wOPgr*1$)$D?S=w-9n4 z62~I0m?rQ;r(Yv2hTRdmH2mJ^lF~9s<#MPEBfv!NIa0GdCqqV~seS8c!Q-`i1&90Y zCJ-H4$Ab8+V}~eK{0T{*W;Q^f0oeO(+|WZkmSMaOtr9Kg10mLu19nI&=W*sz8@!pYMi07r6G-0=V@D9mrxlk-Zb)VmOz!>x7+SCPp zVnN;(J3I;*{uobwY*;tWt<8E`aI@~BDY{_tpiHWZopsKKp6L3?{dam$25I@kbX=ez z$@?L$$o*$@Y$3-NnDi$+&aLFhQFm>1Ynx`J$2wyzU9>xzN`t1$y`1-U_V)a908XWk zbiA!ITg+|blpd#@O?4}=qB&#%Cv^*6yg)Or;0AR&Zu~G-Mlic%r5b3`a!ftjD)>?{ z`iIH=GGTi5iXq5J&MXkM2RW(D&0+{5b7guta@HEc{Pfb=?bkG!X`IX+vn~+rZnlZk zDIP>@t_`lo z$b;c0zk?@}J~ZZ0zY{jX>+xRn&vHH63<6nH-)q;nA|VUI3xo*mOV`ME@%{cRs_1TE zZtU`2IR#baL5z&d^WNpfDZ>C36)gn+AlBaOO^&8cn!)hb>1o?aa#H+@lGw=`h*g#(;&i0kOFB`gcP*Z_`(C+GiORsqH4gxkVYyV&~4 z&tdqV`ds`{LK@X81C#j%Cm=zdl$rk>5Y48vc28!|dQJ$KLd=B-08f4!E4j)1BaG2Fel!p#_gz_JIr*FSIHKhDe!cyXgyAP}_oIFYsU zqi;?174xW)r9lOSuMHZ$H?fi1RlsC2b8Y$gZe(Gm|6rx7nZAB`USHH~S65ezWIOY8 z)B23BUgqra@$tcy^ey+o>mPD{ge&$s1-d#pEe|kArh&V{W8(~VjiD##ZcAm8a5>;o zPEJ|X56oexS3hyV@=V`rK0dsxtgHd{)Z#2LJA36x%*xWt`Pub(28T7mpq(S@2%IGf zh0?E9@g%fnxc*GpEr7E0+&8aGT1#R3UOyXe0vnV;TDbbgv@OjH@}dNtmKXpbY8`?9 z4(_ne@l0^%(9lr#%T*^>Q}0jhCu)XgF*tEFuK2-&2k9zCK_i2{wp%+?^Rj%v2pNvH z#Yw;&Fcoz$Bho(ZW|q^`46# z;|zCtFlFWFm`=yP+}|G=NS$xe&L=2p7(d6Ztp+OB)OJsb9m`zy5gGrw-{G1Ffc~M{ zKz1EZxyFBORg%0{yC0^cmOzr=NXDXKQYC1NQ>>7W``JkMKnz>f5iV;4^q`W}o!9-gM$0tT@L z3aRl=b*<}gk6KV0NjMB(6(1kJcFu|t)}CWtTcII)z3@tVc})b*EuHh`={FWe2*g_N z6>nK1rS0eE`kctN0Gj;{1FJM!@R9*t7~J~q<-G`E$-nN%9J3_2RP*8yA@~+nlxe?t zOUJi)CJHm*H3(-bif6uq4UCAqYu;`#*C2yP2#^DWyL*uUO&v!|ISM;qo zFlU^RB|7N3p#@b?tZPYK76`?E;^Y0L_Z4cDu_ zXe}pC>8g<{a@aHBhz!mS36DS(J@>&zGHocuEurgO|OCdhyZ zY-EJy*rO?O{i>qRwQD*vi&0oZ8vEyp*vi2flS0|k0pFC2*ziNK#uY6VQ~1XOKIeuq zbfYxNEojkpvxV0BDFCdx(InH>(uT}yKZVRR-jQRwS&4dqj}#_X={6``y}mG=UN0Ui zPfD=&mBfogqQKU+HWj9r;!wv=MRZPv&Z79M&Hd7--a|kimW&l# z3?IyZlTbddF*6BQ0R(Qn;5GqUv-jtVafUy^7$;v6eIcKJPACfq2$(HDk?>i^&6h-< z+-N?nKB~I~Kq^qRCI+C=9e9QiZXEq_PRmyFP{XDyf>vZQb#OUVrZOGw3~jHt+h}LQ z<6WYmYTyOVzQwt#Z<{Tadmkz?3#cL`|!waa?{DW&B+4;r_2$U2_ zZ{X~3jd0}?ygyKGg@-Tg$R{K&EG%H7e^y#I5U?S{UxW5yhpOdSbw>wZ65)Xb?9$Qi z&`Wv9vg+#UtGAu4(AD+P6Y#wCL;~W2T8>v*&~^YuxMCuUYd;vKhxtWcvWQBOElx_Ez*60*|TSJ`p8&v|CDPL z2vNE!nAQ^48iw2D%OA-|gH$`WuC)g&%|zKllXuP9-CV6i-=?|iMElV?XngB?={pq_ z4V592osV^eje2E-%vb>uz?Szd7Xb0-+mwywW>R9lhI!ruegjO^Zw|bE1?N{(F`Jc} zTTia%VQ^G#$Dk`XCY$#mgHgEhg2F;mUnO~uuC2Z~wewAT+5t$>xJ&_W9lCY&&7OlTKRzsDmaawwHOz`UT zeNYs|ab0paSdGk^CT})38H0!2;vGYlbx0!VLe)`gRuFdD*iajy|!r zphlm6-~%Zm?vjaJp9QtW!g4dK87B_M3DuctWA8(K)8q&Cs`DsX%4-aG?wgPYePKw| z8MO$998AT`d>JXCXQ3`CekOUbsy1@O0m$RaQTK$w67EE;Um@mlZ1GfMnvBF#^(=om zZ;7&$;TaA`JLBA4-x+W7p}exeP2t5e_|~ORob6(|u;`iOfn>~jqp|{6rvUVoJd_2f zp!!;^=vjHoVxG~p!I~9%yTUG>y=5`m%p>NMyMvh`Z@Q#|YbIZ%lQxj;Z^YIFj1m6tsh zeLSrlqiB1JIO>f#k3*d0rz71`#=bRntcwTON+4jQ%l8G#^8Nz066pHU=&ARL9MU42 z6M}q!t{jqu#dUuR+CB|NyI$|03e_rBg(#1;Cq>cqL&upqw*1_cDYCmac)xd-5`>^p zkFNqxE}g*zF7>lhN7*^k{z!w8Ini-uxG>#T11)$AT1*X}L#?dqUG|c(z69=&z+L(1 zYCVIW{ivA3g^$D06Y#RypX$vtXm>p;!*DLc?!r3I+Dk})aJx@3c76&m><-En`ib3Z zFOhXD_cb5nHaIiIr~OWq7iGh0s+se~-0PKjRcisbkd>))xab(}r6if$9!$erA=E;k zYS=(IKo7hcE*QOiP?3z(!+WirA1rBguispytlz|1J#1NHL^IbiYfot=X(@+%!_-kb z#BxRISn=DFMPlp&`GZnuZqm;lTffu^U)$Htk`q?xR}+Rayt)qzvV&Sfh^b>3ZmbFI z5nMq-^b(Fhf6q+%_z14XOBUUSlUvHwXDP^4Yg{v49EA&OM*sSN0!j2MkCCq=*uToG VOuf;aONW+1oKDfQHN+A?_Ee23`BFc8_N~X-M8I^%$xx$Wtq&!(VF}v-MIJmm5mR%TrqH zBa!fLZ0GfZTs=H8Pj46R)9)=h4&OU?izRd|@y(4K*+Lg!J@u?$U7d-cUhnYx-mvN!BbtBUKM; zhncpex%fmTzR9YMma?HP9}LFfDD?Zz$9sEw)zVrBqQ4Z7$@b~3Z_9GU7sVtLd%L>c z7I{ItL43s4zsBRYZ{Hqv^3)v?iDZXF7VB*)C@7elJn84C*b3Pn_Pysdk6Gc_!2J>$ zERv`<`x;H!Sku@m6I^g)(%{TWBrY?Qy>|oGV-(OIw zKc0&x*x)!xl%dZFOTKAC;4+n^!+)E?2YAS`)nW^~wQ^}0DcSg)@o*}#O_LJ6{q zmg`SS`y8hzVDOtc@~r9GH}Z((&B2x|4VpHA+QAVAcfY(d_@xwZ8=VhYFUI1at$+?PG(Tf{q1lM%p|dV4!o?dZ!mQ{f#?iC?)nJtwHH8 zQrO_|ovPj*g84alDI3C8&iVl^U&?xasgp7orO?J?i-<^96-jo}doeGLDPb$0*Uy2* zpZ{XWY4!)`Q$*)qK$Sm0$+8!mS;d;(|8SB4mk>nH1zFQs#58RNW-GV~2DQZ8Dik5L zz(6|Ytu8fDD{C*+;Z`$5QK%l4E*N`Nh0Lt5{)<>!i)d?nWqun+#zy;OZs_WZN!J0( zI0*3aw*nZvBW6SOuY%Xh=YEEU>F-?#VgKI@g?|kC??lo+PG5Uw%q~jA92tK%Sv%*!G0ai>eORRSV z^HeVTp=jei2=_OY!)9&e(gD`sjMTBBU{BRgCG}!p`wv%uKgRtJKK4I&`kK`?m*q~9 z&1zY392>6Ws1Xe}rm78|F!<_vr-%*7&Dr6@y=~lwXfvMVD5*n&3KDMc)DvOYp2bo_ zV6#QtPt?lg^h$ZiDg){__d1`Ne(R|S_0Yd6^6(_`21hS~0_l2z~eh4W*6{)E3>+l2p1y4JmF{FK1&GODcEhDF)Y z@+Y@%-%d>a%#$2eZ7RbPyR}(+)w3l<8F%{ObTY;&QppsBX4=xymGz-&DfST0DC?oo zm>7vmb$~d4N^F*vG#zqU^W@35%VEz^;w9Gu0!Y?2HU!brIRrxgl1w|G zY28z`rg!%}Xw1|Ja?}8pZyQ@!{=myYiVNR!R!^c{SO2VWiWnbRA84tzpCTy5im6#q zw_&iKl)DViwFabcuAD7h8Sqq3+FX9Nq|0nSiV4@S1;z5Or=b7MwZEIa|1X)&iv}Fb zvEas@~!R=$z~g4FOAmo(Yd znHg`IbQS4jhx%KLBf7?K6e?nDTYf%s@0FZT*^wDHTl>?=K7DFPc5n!smw ze@?tBr6Gwcl-oBCqz5?|K34JLVT1IGYE8>bGoW<8LQ0oO0oeIFw)YefQc_ZYAO(0{ zU0t0jf3P*I(%$rz2)CW33qw=8Cv2yu_b=*hT=(w(qiRq<`YQ< zXFHwS%#&OPi-1AdAo`)z)?#%WKHXloAtvPUV#?kJkcbHsLSe(;qFPZ{u7x zrCU*Ww-G5#NtK4BU!sUZ&0q`~P(R5BwaOy5ajALN*T$uD61s13;0bNGHan9>5ZktR zwma6AD#d@xJPZ)c8**}ALfj2%@Qa(p6Wml_g_+cU5j6FH^o0$w%5ptdOR`qyD&nBQ zYjA1+9|Puvs*ZCF<`Cbuvbb?d2>K^I6t?GUui z&O%@+1O&Q&mCdOiONjub7eAOwfXn@-)$6|zZ~xBRgj!i~adCf;-8{qe))yXbC5dnF zWtEx0^GR)zhhy{Ym~Ax^&56Ga!5mDy!yXLpDk0G&tgupcJB z2s5v%+jqP{PtPB2aqQd{YJJc6k`V0KONMu)0V!AZ=gxqEs2!v8visbL301~-qWJ9v zY=Iy502E64P@8Ux2Gxq%3BStdp2&yeP9|89hBe@j%7;DsdC<*AKx@0fpT`YMU><0A zCgNKN9No1huV_M5=$oM@y#q)3NE*CeuTl+^vG>hZ$P)5vdsGG*TI~Xk2OKL z^iQ9lSrb_>+p_p*3cmu8=T~0;j=(JqGDGz(>TZTr%|q3xiWC;+D_A;8A@wtV+OqV8 zhO@)bl~MuIgjId{c7Q|jd!8zjtyP0cY0j>kOULic#F`1i@*SX};3M58L`nLwi3dva z;Iz*_VeUe4`2_p?-q{{yI-f1jx$y;BYLAa)n3 zV^Jn;?N|K=ipaUy%VD5N;-H?UF-in=eda@$US9X-ykE98ec~NwemZv?bUm&$<8YpJ z73&nq<4ATl4}I-s5E_k!`bb4yC@^1w$D5j(PR`Ee)9;r|O|%u^0|8lbo*nrvLIcB| z0*pKb*eR{Bpx~lcZ-+n-PD`&kT?C4ON`j>o_5J&sqVm7w7Zq7E7>wKA{$W>KT!yJP z#>$fAHbdV2Tj$ySi$;G!BS;T}_P^V5jiDETIxbTO&ffh=yWXS&a@jvUL;p%z<;jzbC}sftPUBs+U&w`KjpQSTD-EF`MQjpg zu3L~2IFrKp+0hW02d|QKV5i>+l&4BAx4ssddCaV}VoEK3EBag0s(3IpD=?MATxqra z*4NhG_}`w}Up?jh=Tq_c$5t#@6krIg1^8PCQq?M&DY!k7%i!IF2QfLa2R1)m@`^aP zX@l3PhXq^ZEk{TCOD+V8GLsqB8EO8O*x@Yh@S_C3Dz za$nAO$U87|apm=tT|ewtScoy6VS96UwKqH-JCw+fi+QCO9d~cv(ta$BK}>6@WG;lv zU}$I%HZNYR2eA)c8)UyB>!CwJ0OJA0JVnwVgKfkF{@KWT z<*AU<34I0V~6V81Zzw9`CEQ|hri5xAL##X=Ae+5mzN7IzH#(fE@>c;76`_k z^}3=$w{3j>5yJ8a-;-6zX0zKpNpjZfcB&={nPm25*99^?*jClT1KmXX;2*O3{r&xS zz3B)YpESB+WgBt$RM7US*N5w*Que93yNWpuw91I=_rE|ij=-v<_~sqVwq)a;{C00c zd6EX&+p_m7wQqw{CsC#+ff_vaD*XAv-?_RQ2x!|iRCrc;p;v7k|<*N6AowEYkKmsL%v1lu@ zF`{NCMEoyYO((vT)Ra%??rt|{PbQeu=MH_zrtN^SZ$jRVIQ){RFyTFv-4|o!1G@^t z0onv!-x)y@uHfto!fQJsZvF60MeKSCo%LtqZ}VCs^S_sXa`9hj3r}<))bi0sb!l2Q z3@a67TB8DGF`aD}8}BdNGxrQU08`Lk^P&u*awRr3Z0N~=d;IYO#c2hywaxoI1NXs9 z)^ZdUBcpv3;bHGPb)YHNHB37)CLk#$yjxNKhc!gWPx3ts0_5EbJFmbEFv6#UD9Hb? z4b%_JdME(fFJ{XsO&y}&WbJQKbaP-77psxi!`M&O;zaI!)rG=IoHh}>lN=kI0+w== z^j~iJ7SvB;O8e_?*)}P-BU(293r7U_ocfH^fj#gc0_S=j}aoskeK}=vb22)TC4_(2%WpUcGdrw zM0(41UrY4-yzs&TUGhDn*m<)1^`ja5jP)6LdD2wRfYZ_1HwG?*_Mo=qT%XLA4O#hA z-`mN{eme4SA@U2=9jCgyR-*JTZ`~Tbf}yTcO|pESvzfN=Rf9zbZ;kwCSIg(cDwB|Vlb`1NfEcQ zw#svK+Tt*&6P|&GF#(^Ghx+DHo0=4FBYgTN_m#46Wl4@TcUI_{NxtWbH6-P}q4EU} zl1${H+^UBG@#sF6fOv2W{VE?n3q{vI>_GEyBsQ8qa_X=G!A9jm^E7Ef@= zMCU6useSYbaDBzG&b{H`HbUDU_ayq=yTbRvQ%uh&Ci7uUe$8XrdA&A6kBUV) z_0X{E5cW~q5B#Vr7i1k%xkpzn_DLt_0=PrNxWIczOvFhmb1LuM%<+T=<&`4#VK`B% zSi7h+mj_N~rgCI5_DDXtWXwRGu|y%%2#zkR#+~15ix@baA>>fOJ$+(|d^;A4z3(kd zRn<|nHyU!&t*uPv8;#H%iyU?GQ{;qM;cVojuJr{)*=U@MW=Upu-1%IijF?A{_(|~e%R=aNiG(EctP^|mJ3UF3 z*{qN)S9#8bFY^>l)XtOCRgMsf@6zoN2C0n+DNVUon1PltsgLSY=&u25Nqf{i)Km9M z{lK0yFY+cH8K*okA)zUcwmL(7vXf{XiY+ER-1?-KU!w3qFmU|&=c{vKk$7ojzk+}} zlPyEE8-Lyxi-gM=w(OV;2)`hjzSl8DL2uIEP&T3aE04UXcBJZpxnlI^-jh4hv0gD- zMa`Obhbl-*yLh}?^0U3eukbrGaTPTkqMS{7Vy4x;kX*O}JIxkeJf(d^8b|K610LEc zud5LL=y*4e?4X#r_+@tXKAY53;@&59oBEpgDM~`ZdwSu8>+70FMIB}C61v`o(cLd6 z>%>I~eIl~;r~2zrcELA37-U&z6ZXz}_q@v*XzDNl%Miyt@3U3$ad!R*i&aQ%_+Y$U zZ)tFY&Z}p5hJrtOswIVf_lCRAKZ;0&`816B zdXntaz5BYRe4S6bHzy~6rKWU|T{Yma+FS7^N$R9_@%PH9;I+tRYrzlwlA!AGhD^p# zBzwOyzb!mePQCYF&}94dD}x4@nNqs%Xq=54VvmB|(=Cc6FY0%{|M(zW=51W2{>_V^ zqr99=;bzO!?=X~j5{;d5G0EK|XMi+rb-3|qcXy$iaT%GsP}+jc?H;YVGmbi(h{~%_ zH`8TV#*E;za~r0*+pe?Zze6DHj{&;Ot9WBu?N05cX-dnT$uO^cBU6xwnM?ON3y*oqMdM{1(u$D4M=-j0f;BQaCgEFQrHOAaBUao5mkm!3Jc z6v5uagmR+c{0TF^iH!J}f(n`CaYJ}wteIX2g|Xq#D}%HsHj6cSYl~uCRh53T9&aic zCs-CSSJXd1CFc<=kCaXoHWmg>_vj9U58d56_CU3{iy129WHPjF0APn`*x}8i96Svo5Eve0V?eT+B)>ZpBbqW6O(-6I_SIu_SbJ4+#aGpk=fD zf&^;@gQ@rXZQL+ZfdtM)6G9E^j*#XcA7(3<;zH?{W?B9jr}xDPsMo8lrJp%<(N>yt2k6_-!DnJd>gCgt-$5=jSd2D<^ZjF71uMZ#bz z{cPj*4mltRH21fU!rfZJCu|U(oY`I)6)YYx1yH! zAh1WNeD=d@(gx+Y6R&V9eRRFz&RzH1I|nl^WNmMA6AG zpQwu3aOkbUA}RgUYeWlaIvq@t476>MNrsWbBUj2r<9 zo>ufczEnla!O_~I_Vwu=%$9aqZ`Wgt9NM)?&ptx1ULfH zhEwumi&FSp)LNxEe$7khyzIkCDHO@hI_J_8aqp4?i1`6dr3MK0-r>n=n1y2nuP2Q^ z<*^U6KR_K>0AF(4imkG#GWDI6Y*N^;IPF|a-;LV+WAqJg%7J~zxSs3fBxzI4(WiZ( z3<6Y91W&grAGj+iNo{R?vSG^}nSD6nMMKFD8H6wR;fI->JSjR750u_@^E}cUJ3LA3VvG96WfFm1mwEkrKqBR>C6Y&Kh;w9QK_Oo} zh)Rd*f7paEa&3qrlEHw5VKwJI?)2kMWO}r|`Wds;)sQa6V})UGJ;KkoPf|N)Dm~}2 zn|N29&N-K5(L)qq7;9!Fdpte5`Ju+Z^W8_SU4e<0ms@bDAbv(($C@^PimeK!ZtkIX zDp#1?-|?;6+4kgq@fwHtmp7;)o)rT2?;jZWhDcNe-tY|6gu#CM(v#!}?!QuIl<{Q> z{5;~t_YXij)cFf`Of?O1qtilSClG))cE2$Qhjwf<<6EKsy8md#-mcG7@<_b9!K`#E}w?HPxPe z2l9Xx(%r(OtP1*YtNi5^Iy@~Btjz|i;_`y0?rjsFZw{(dMRnvr2UYTB`mmh-HRY_XMt!D7cI zv-3NnaaKW67j`wiQLruZbFwwm!Q3AWz8UiFP$W;{`o#f24a4ql-<#M_6&x*@s!Ek27 zsY`&3V2<_i$rK*>de6IS*Dn_VN(u3_b1`A8Y-H1PG0iE`Llg#cLI7o5^;&#?liiUx zb|}S@lnAAQ)N-96M53M=WdPkU_9{xTArihaB(rY|CMXExnMZGtm*!M{^*s2g6qGC1 z&}cJI>mZN=V5FLv8tNG{0@&U;l3hPtEU#=5FC%?m5fo?DpqV~!(Rang^tb7c8Ylt= zW$eyQQ!p$5(a-DJfa7@Q!I1PDwOH*IkaKkh*dn-D5k$;>?Lb#>TNgNz=x$=2>yRK* zI#(mE>QeK-MyJR@kK%n%YmFO?e9mwNmb6cC7$6nUy&P;SvM6%9RM}yP-Mu++`lhK&YOvOT9qatWCjK zG(dL9u_ekG$167dy8Ecuymt8PZK(%$$vCzL?itW^fzjy?gVMBs!^GZ0RycncjpOyk z%7ARYP#stt^0i6%E<30ibD`pL53uT~YeB#?3(he;^FUUwQOy(jf;eEVq z+qRAVc1F|$08r)1s1noSAyg+~MjoQWTmThhFO>`pU5)~)gn#pat{~U{GEy?p)%6ZT zpy;EX8dSk>On~FUV@Bp1m83Y~5K#iM_&JDY)#j5fX`VXZNn>QPZStB22L~m4k(FT~ zTpo{u5(iT56=d`CWkHn-DK0dnOZh(mbsCUl{xUsXnzk~g5i6rX3-i-k84{(_;j1A) z145ende>vhQWT{Q1a+(oa=F1q1tcsMKL;h5O5uxMAoJdrgK$H{ynp{*hdZ8bTjqaF z4@fDHR+@+yz@(eQz->?5fdUc1DeJ6IEbjAN<6<%hk{toXlTYzCpCL0Vp$v>du3QNp z)(DpdwG8}i*>kK86xw*-(qtU37}p3OlU-U91aw6yxMl2v^&f$pH3Ndf#>QpP2a8v z>(748c7Ex8+f~QQg~~d<-7X7`SqoNcph6#@=czV zbfDT&SH0zvH(j*LqkW%%pxP!W!ws$X>++EDCuxl#NE`$SYAOQ7=%$rB=?2&Cfd9}) zI{<1#bVydKm)=e<9N|3zo#Q11hWo;(n6|6 z6fKXrGnMKP0>W7laW-lg<^TAmBadvGb96)z^k_T;dZe}mlHbg`@Vlc{Fn~*Ta8|!7` z{JT2u3`sSC5o%Taya~c|pR}X;3l~%ihlf6j+bg+eU-_S}eQv|t#(b_O?7`8N?T9yh zhnb5ifE4>vBZ`SzjQGJRbAGfPiS+q1F|n4oUw16DD3PwM)Oj?ubRLJzT=DZNe;c{E zYej-J;Ok!+=V7pGYCWa(*xE?Lzuow{&PrWl!C6v%7vURpdnkE~?FJmf5Wp1_`tE&c zN-~0wzYoVm8Ze&{xPimYEw(ri7}VGR4>EFUAl=t9hYOHa5iH8oCvP$qWt3Z3bkIJOYli_|5zCeSuszZONM%L!avi$F>=D zu5;lC5#rA2I@l=}M1lQ(x4kZ`IYrBz@IaK%kB2qGni-BUh|kt!+UO4r^%=t4M2R!(f4paWI zQcyX{p;fxrZ->iuu4)?}ojyORZ6TCY`}=JQRb2;JcEBn1+g;j{qd6(U_gF7DDmhCU zv~`g3T%{&VdKX6liq3)a@0vh&#|%nJY=X?4Fy$`1J`E%L#*B2JMUL#Zhkqyeju)|- zx<{USFsRdHp}u1bUX&BrMm(`Uq#7DA*rwe* z%S#}-veYknTwse&68$pUB56)1V+4E~3>eZ0EZWW*|DmPxAtZa+G{J?+LIgA9CS9(wvNSPr7QAT7DPER7(A$McI?O9#?V&^Bc^MnOROAb#<9%VKXu$9p)%7sII*)TGzQZxRJ@tIIY~ujnwB z93}M#nI%pcYF*n67xFOe6ZMXXgbzrjWG8ZQmUXg4bhy@lY9jA}koRCpzhw{SH;*nE z=Ci?!c-#2JrFfWBh=1&N%+kL^Nf8e9K@NFNn;9c!qBVVY5NwlJ?tx5Cb4fGCzS-6W z@TwaHbHcfBtaIq)O1H$9TEof)0bJ62h&hAlhJv$G!@q)hYnt2iYbD|dBNmAA9UEg( znqb_Wolt@E;D46r?$5$6>Jz2CZMV1$)yALm*FiFAQp zNMO~bvwa{-q+Agiwrdk=5fn=En7?^`P%d)R|F}>ZD#Ui_pmULmX+1|^42>YvjAwWs zL0V~nduc9PqB9*ENtw2GIA(mcrj9HgY2{$V;;e6-lB;iY$!X$(O@0EYH=NAp3)^fM zAfb%*rZTs~fco7U{-+Qo+I?qaZKHu?6jyz{tCLgP=^68*m)uk~*};bZ0XaGy3M*FP zZb?BtGB{t&A4N=OnnG!YFjZ++jb?R~RAE40vw*V|=g~jW{CCLH<_NPxgc}8+KWfbh z>Jo^Xe29@udxZL#v9~e!1NSd7gFdYP@$hZWvt4O-%D&J;?b-7ENm5AKjW0-Lc|on8 z^_YJV*<+o1_?+XO^B;bp_1udx{cpt0Wl!H^9aL7_;TGulw(Z%o2G(ywc57qq4L2|X z=O;tBgvzeakhPWs!HV4qa$1xtZV_{60A{T}yTXnFg3|BJmnn3oM+z!1sLoC&o~MIT zmGl0-*iG=mZSG%kus`T0Jx|oPC}^yjLNp`KmrQ*IKQbBos)<)~eFw3^D9rj82Azno zvKx+t&l;g^{`{x%^Cy#%uRc}lv>pD_hPHh9N?mdFGy^kQLvy45z}Z_mj~pgKEP^Ip zqYX~vv(~q^?}g930$MY=3=JHa)fV|d=ggRw$C9u_LT~I=$o;r{dWF6z|7W<)KJMyL zXNXsO=q8=gke|1#kG56>5QiSzd{{fT&D)G=!c!RWYWk#Y|C0@S<(>W)zWl>AY{yu` z+%-AsP&G8)vfN%4Xi(bQUQQSF3D5ReL6?8J07oL8-x0kRauAqv*6pqnJfRW_Up1ux zJwGVtkZt5ZThMfpQ~BgS`ZvB)B^!!Es1Hf{C zkKgCRyI{CZ$5Q-|Svpm4pV$73Z}et$?b|Hh+M zEYnwa3Zce-T;I|-W#D*U=gOvLLwccWSdUj)UAa60n!Ll_f+9MN1xiXU-&s6>?)h6` zx*y)rZt+!>ZWa?6 z2n&nCi&8haViyF4`FY}bIYAL8cW|f@E=L=sODzVx&|_XwD-J)x_1F9uE2;1Xz{Q89i6YnniH`Iu+xtx0xNpjf!CJDU6WM*Da&*dOq z-}~7Mo9v5h6AuJB@wOfNKeLT*+Uuy_1ed(=|V?=AN zy^ZY}((0A^v}wnWQbMkZ%ocpCx|swFYw#NvM$}2+RsM1rB2#sadbvHrXPRV|G?h< z`k1%m`28cR9mUQ@MuF-TW@wDHgNsgAel^=XC#7_1Z}XeGr-_B<#;;AO2&hTr6CWPt zt%U2lm=%_m6Fb70TPpJ~<~vKt@7fiE?eblP^{|YIOABVhufLVV^mrwQmEI+Vhj)@4 zm?=HylkYE3|Mn=xQpd#{Sg|}`@31fj%_!@R&8+GL>fQA5ma>bBu8wUEx*p}+YfhnF z4!YEk2iz8KY{`@XGER6Q%huf?11wigG}K{#7rx01{B+;*!Pm&R5}1RmB5%Dys_(#d zEk3q5FD3Nc_)oF>0loeXOcBDRbmFxff97QFv)j_7hM=JoUE_%C_p+ph2I08tm5&jL zA3Wlpe|nG;e?vSn)PUu&1~N5~%f#nbR8{BBLEFp+8X`NM#7)P@*2ES&71PWf1$X2` zynCdlYCCsuV&bxd^&Pa>Q|6Rk3)*cLyN_7K&v^En5BSRO7WZKtGVzAToF?bdxgIo{ zj$c!_qN*jGgD6UAyRh@Ws>xn}iH)c-=?bBI$^d#$Sn1Q)601yT}Ba;A73}|v~r}TOgk|QE!Or1Un;L|Gq@%KnXmMNP$6h*-3jqj}QVx;Ibn8tzO zn_S2GzKW2L5KYfiZ)iEVyRfs$j~1VinJM$ETzcZ#j*c8?nDuKMiXXU&Pv6KZ>2pXZ zC@5IHkDyG)T2lsFQddg1Y}sNHUb%zE<7rJ66;yiKGu8*xEZfRtL6u96DEuB~AHpZR z5>r5-$b|eM8hYQ)t58AVp`tqa>bS9%ON*jt;K)w@uYE9y!}_v2}6%_uuEw_$2(%Aznig*GG|Mj5>YTlrc)>5vF( z0Bxw~T97G|QcC)g7czF*@U3IGCwXrl2KJNTQtfF}2q~hls0bhP7ql{@Z!w7#fhmg+ z_Sr!LgM)+0n>1}d*7ufg2&ms+C2Jg*MWcj82&D$hu++L)b+PN0s*%ciP7b2bnb%m< ze80t%DFT`xYp&N+PJ8-cVO17hnx%CsH&?=q1Z*1>4PYn~@UO$Ud|?!^fZF4rXFM(N zfLhtxxK=Ap7DJ0iO0$3nS3L5!q)E5GS*7E|d_uvc7;#6t-y}U)nX?1aG*5JG--yc4 zvWyI^y^^?I1749=(3~LP#Fb{@%F4=G@XHMfN}DuL!uN=&SZj|eYfGqN z|4Xb0i4#T4#9GtmOVpYt7e}GA%NxU_00H){4(xtpJUK>f8N;#(8b;^D;zvWDc%k3Kz&48<$5V?C^p^@$Kpga5}LANR?ZMQ*Bx8{(`#AB zPM#n|e87L?*uD5{>%Nx3O{H1-q`Q5~Bb#YM6mFhtdq(|sz$5PtD_UPNmZ<9N;$mN$ zcOLR=KJ-dPB9UwWzPO8vi<#ASyP%6E#S1bD-zc-s+X9GP6_FH2yo%5=EFIi$W)|1a z*4g0fE&7t<+TMmZSMSZpc^het3ai{A-*vq^Rs_Tc6XL-T34jwuN3x^7&!v7K(q5_cTSJ^8 zhNw&Hc&H=tq5uVE^hW~;gpXiF@T5h_AMV>v##&np;#X!TQq$gd_e*2|;&gWhNiR?g zUGjih+VQ10sce;JPk+z#tnULT%$2kQYMGrBd}p}-0@1x9&V@Lfx;BlL=_U-GXWXvK z+{A0h9=s?n;>UDsrup4KRJw=1UU{T)24325VusCT4+G;>V&MaA(_AZ%>FVkdinD5? zl9QL+Qqw>0lH;)I>>Cf%|l1yVU~&yzKyL=6M82?Z^;`<@t~+87K$ zc;f)Do3Gxyc{AO$-xR9o?_ZEzR8qC;jcTNvjKS*#zpl6tUMaI;Is&v69P=*@*Ab5( zB^eXO-VrOny0q~>c!&kVvQ9ZbwsQ-RR#uV^uI-AZrlzuUX6oxE>FZv}DiPCLlL`*x z1o>6Q>8Md4lW!EPWpr)S#Kgo_!)Qmn%viE2<K-4M#R!4Sb0O)I~L|sv_xOzE7Dw_b}s(q`?#3UrpN>->$ zX^J9LaQ#q$jO>bUr1>sxRjKb2@T>5PDkba7y&6VllO2ianVIyZ#f?Cyj>4BOUc6A+ zz74ho(=ZuHk;Hsp3h>yQox%N8e*N{{mD4P1LuhUP99HD&pbPKbx?@(kHE(Kq)gD0G zLvP$AtP(@$t?-YN-X9RnuuXHgx}=mCAXWb9iOX1Tuc-`Lf&d?=H$J>4X$-7)g%lJY zFJcf=YI>%cCR=|}pxCV40yf4L1Dd;t5x^2WS{W#=1KJm|7%m-d#)-{Jv${gwcg48K zer(CUWdFL5lAmuQ8+m%sWuvHKeSp=PtqlDNNXGtteM?+a6kkzERde-SYs+6C)=4^m z08U(G6Fb@c15r7w5(j&%NxrNIVWDM@UtJM2EwE1XPj+17-OcuYnh$`k3aA#i2;Tv~ z-l?mrQ{9sUERJIHmu!Y6?m;%U39$p9$frUb-d~`UKHs*(g6wKX3fet{Gb|EWO3IW4 z9{lwxSMWtL1rzIKq-;2TAX2xiygZ@T8{;Ho%*L;L1SW5&zOUNp)G6OPgYq}b{$r2* z$6C{Pmw$u9i^V z!G#i(Xz62$lTc#z>fx=_jrR8NxYq9EIO3RNgly?Rf}$XX4f4$I!P~}AG0J04?m6qD zsvaSUb>IwTGn}b)Qg<_Dr$B+{#xWJFId}~#4@5L+dwj`UN0iYAe3OxK<6%pQB_VeV zeUr;VWa>)LC~4GI&W%XD;cPA%K{aM>K}m5$3ocuSpBpl8^U}&9C!2d#JYPO^iq$Aj zm3R*b2W9tW?c*Y2AXx)VY|>{6LT^7ZL>W2m??v)u2Rb7&r~whaC@6(V3~N0|R1&U zzJ?QW6BYNEYnQV@H}pL}jvwAs8B^7lIgiux+8aAex>RvtxbRYmIJhDaWW=61upob| z${zO2d57Q1l=m>#rg0ZZ7Ziwh+xZuI(KpwM4fs4^J62 z>LmG4+GNiW;`n~qA>u_;{6y};h-2+{>tzl_dKW$NX5Y-xDM5h=YeIy;C8lI$e#vWZ zlJPIK@NG$}6Uc8M&_xufrjj0&ecgh}1usT&sp)*S-)@>}Rda zgWeZWfoHO1YNMoknW~iF22zVYlM>HZJxlDL$laU}NR;Jq&k~1E)kfv=K1Iq$LS-lx(H?*nzr`DoSQfD8W}oJG*^ literal 16009 zcmeHu2~<<(y7oqjv<}c#1(B&%Dk=ghqRhc#72`k@1!Sg`83bfbfCQzLDgrVn%3M(p z0z~F91d4!w5CN0O7?dGEAYn=Zge3Rdc>1q@{r8?*?>#;D{%5U|71j>fdw=`e@B2Q_ z`@HYBbLXO|$<|G~HbD@y6?x`+GYFChfgtf8Ki>$x8F={X9N299@tn!`U<+O=Cu?zF z^M&sjD}M-5jDr8f;>PElA?O9V09x@qQKFma%k$;w-$$B!$Vzr?d4%x@R@JMEQQoBW_zaiZc)Cj-A5R>dnhN0X7ndoK>DKSaUqq``PlG4*r zh_xRIX-ulWBO#~PJv}4C!kp6Fe9F$wPTT%7K84L=&X2ZomZo~GHIfS%4fe#UjrsXC z_4P%Ii;FGMjv{-UfTM>%5SUCR^|B*jW8?`#%+mAx{D&RSTA|pU(bhQ6;;V0)o0|*k zMN3Od)SEPuqp$&>pZeSML&V;Zc(TM*>Ogjme9Pz$7jns@`g*0ee!9YUndn1(RNC;> zS+pEDdJwfiv%zrw99=ID$0xZKx`w@4iWC{KEQQ=qPIEPR92h^M7TBJ3r#Z(`H9H@% z_?(lFnoh{C7e(fe=g=^401j@3N9kdkUyxe%;r@&Z{t|p?zhCrf7ubHfZ~ji#{*BY0 zzt`}oqqAwgRlO%08s)ZZR9@M^eMTpvJG~F>eCzwp(@@;LT9((uU{ra26Wb9*2*(gJ z_>f5laQG|W@Ze;be&&!+cfpi>V{s)c`M9*??95C7=d;+?A+`?^&O3z8mSdJ*a5lYD zq!}yJwmF$4wms!+ar^ODD#hy<7(GQU<}~L`vysWP2F(~xnWnZ5xsMR*>gsxR&Ky;c z@-)h%t*oH0A4{eB>wMcWUQnU!Z8++W=FV6V>1=Cqgk=$!eP7>WFuR^;&#e6XaX%|d zOUoG9MxkNITvkz0kw>*{MaG6u&WBykr_oPq>gv2ez_P(4E%4@c4i5Bbes1qnVH~ly zp|;Pg)nN3- zPu9PSQiLJxQC`EjL2Qrj98w;$(a?h4cC7&IkQ?>ceJh0Z@g$ze_}EYxte z8$G2l0U!EreR*khE>5KgFpx$v%wg$iX>a|i3pkh|%%&N05*@eF#p!gp9Nqk6{@Zt? zJ}-;Z&Iu~bIC$mxK>EXoala|k&Thb{n&O1b`doJ+%d%bC#n!=rHa@<)m09$M>xD~t&^uli1pkJ2V0#F?F>_y!)zJZc_Fvgy3nH1 zNQ6T9h@`TGuA;Btyah?c|ACnJJHfTpMjARZr%BQ!(6NjcqV_c*KXSUF&@!xsQ;~0J zj^J@hQ>3JVXKI+`8Ss2EJAR;TiV&YyuGcDLULd$gu2oL-o}xIt!8GUE6yPh&IYeHR zw$4cAZr!DeQs4wl1>bx~=U|+xxr3f3@q-H;Tg$mcW}^@VDGGCY-W;cU5#e!-GR3i? z?eHfk^YlMPcKw}}`p2C9{0CTc7JA{QeFO{Ht4-0Cf)t*o^580dF&MX6Bb;}r3wc)g z1q6+43lWRRrd-Ek*@6>>X%^|R4xB=seV;lp)SOexuJ#fKM(om(fJP?MFXRtDsKaW2*oinfuJvic_3aA|1dU)Z6*CkPKix{5R+v2KHGly=X;1{W-(b4duN*fmD`FsH zB+ZG`YudW_UfatM zf|lC~Yw{6#Y#M^!;?CTE9J3hq1_N|t15tHugK<`K&GZu{jUa3+<0Q$0+k|!dAP&g* zpHjBsv8zbtj^q5MTyvb?Fa}fr4T6bJG~gd8|Nr8x{K<94?^nNFYh+BI4RCgke8a&Q zG7u0Usda}`UtDHpmu2h&hM2;hII+1s>RqDrjl-p-rFO#x3T2YgN@W=hYHt@8r4+tN z?re1iHAbQ;3P?k;vtBwiK<_#bnUBZ}?M$)_Il`r7wAfQ7mg;p3-{%#CXWPNl#xfOW zQz+&@_pMX7}uDVa#jn zqjN4ofn*V^k)NS%Lx5x8)!07; z-@WxsK0eImdFl*K?H(D~vN1Alt>28q)>eq%KcO+Ji#}w9FZZ{uv>%q=vat&n8`}zn zur$=mXw0G2Dt%&WCu4D)+5qyY+dPis)ilwi8*cjZjFeJ~raPm55TyDbybQ-LhfkH~ zlSH!Cgpn`L?@dfOzZaIVrF{i1rB<+vbvv3LB`wWQkq8?9y|(5gY=RnY+IZy)N8Bl- zRSxx>cSWB?VQoHXg#x2a8wgqhAZRV|=7DyV9iSyJ1yr_DoM*KyHK_J_jdQ7-g*iX& zE;Qm@VIHVgwCca^x7uWLlBg-4+5BXqhA7sIL??7lqKD{S$I9&q@__wP&+i4FJWKn^ z)(dg8uemubFc_oy?Ngw=DS<_J-dns1VE-U3t#&Bp3c4mlsR0i3tzfar_q9# ztxY5K9{jJr4l4mopy}IPUDuB8ikWVVlr-n8#V9ct3|KCMN&#*Y{hz*2V@E@+J%QO* z>bl+#z^+vRezvo<9h~N)fSwKa5^dN9$JHly$bRR!e<~RMN!`JJ-TR+bEjLj?N{djQ zsxVJx+B4-c-H*JZz8W|MK`C(78;Mx!C1g@?LPR1}%cxh@v}Tn;S_H>IkO>NuAQ9pV zs>uR9IdkG+A5ci5lFiw+7`_J7u=`8+?1p=i<|tEZLX}lz;li*T6wx%|91KVKANGS` zg=!{P2+rlm`{33z-U13UD6-Z-2lPHsSQwU^F(A+U8cwpO6q-kxQ9Vg(|*TIS?{Uj`amqBDhDSDa7cy2=*j0zs>`d+s;&>-g#Jg6R%vI3=*_E|)KC{0`zys-f~X5~*rN zeL475s!gFKI2u1r42N&*&-07j@bc6Y+>vH2eEH@)iHqEJ<|412i4156ZyQ*6C6TdY0wZh29)N)B9Pi*)DN{ zzBZ2%q^}N2poN1n(KU4jgPw1_wepi1Lec1z1!lL~QhWpdJw5W;4L?6Ka0NiWHmIA# zU7!j4MV*qk>H%B~!S!u>G^l#((I#!@rp2PmAFH>WUBOvOp2xt8w&nswl+sEv6-= zUKFrik)D8aTN35rk%dJ0pG&H)EKJmv`8Jt3`Fz#?8{ z{u66zk9Qa-ETK|kHL!TPJYznOYG+PArvhX7sLGt5@34XLTL%Fz4nW~N(Vij;=nURB z^5}~?g|tQK;cEW4qEJBkmH z0?PLW*2#Mb>HbR2!Cb+VO==846c z<69#@UEbT<3&bn$ZR7l~RyeZ$5Vmlb!-fPH`jibVSYf%I$sB3nnmNenAEUD> z_?Yr%!-5hwddwrUe z)Rop8TeTl;4ip|r&)X4w^vS+6M|RxVu{ZKB^WXk7V)N!pMfLBEH#n8Hu`-7IS<2_Q zPn?doO{9f)C59PkJsz6TW$?mK8U7h`|9UhRUDI=UjGaeK%)?&dKK^wqG>)?`M9gi! z1mHsGa=t}E_JnNukACDWsspzZ_PD4nVxmp8s$)nRdSZ~|w~UXi_I)1@hm4W?;dT%N zojv|yP=4xzb@Fj|IL}0eS8`tDe2bw;}%&@w|{I{efR8>Z{~kxi2U00NOFduzCa+D zwf#&o3w`aulRZ$~Uq4VG5{XtXvT`T)y*XxRXlV7q@kIAw#};!uWr@Ri zW`V!JoR2wut-r(Dvr+%pv4_RQ#X%ak#4i%pyWmw6E*` z4VC{D>;8uUtL53WZD?iYw(w{Q9CeOFvjC7)A!wh31pCAU%fw@)Yft&>!X4>Oew@?Y zV3h@evJ|%0N&d?6@EiQ$RPLj}-ji)jHgJH0S2|huZ=vJg-EvJdy!Zu4M!IOByA%COyqmR+ZEAxDSb!SlCa{O>0E z{|9*fYOa9Wzjlv3e(xdV%SyxMa0xzMI=z&;_*-AhlWb_a>AFQ!kki+Y-?h~ZpY6W1 z*IfCW68%yeA?V7EI_T0mGmyK@k+8j*$Aa#ukcKY>Y>Mc5ukgtnF?322_S`$0PFK?K zLT~-xqFeTxgO~HJB#6+jeeDfg3nuq9h<%(&Qp?=rGR@`~=NaD2>tWu#d*A6#KUia1 zO3Nrp=rpV|SVWI5>xGjks(b%lUXyBqi$9qs{?6gkF+;HnfTt57-tD%hPuF_CUQ3dP zU1vaG(D#ngAAP5O?8rk~QMOgJ(4_Qk?2%`*wMvW;Z^HF%tk{g?dKxf})=mZX8hazf z3eFsiHz*EA?QD^-6o8y?6|7cXN4N}o4FY4aEMxR!)8{Q zm70#pPO!G@M(#a+G{s=~`?y;?XdP0J&KXX?cbfWjliL&JcOtP1%qN#;thCzR3)-S5 zWB}3;IE;a+WQLScl_Gr6h>vLck8`>PiNASvsQSyQj)U8zZnk~ShQa7U{1fz*tHdSa4V>n?#5uLr3= ze%5?zXnLeH8J-7pXvfDr|49<|6EwH!mBZ2VKJfY0HJy1adL^@yHD z+iizn?;F~B4`-gMvT@_F2h+4H0@8Xe&Q-kI9o}|nG1~NSb*6!?@HA|Pkr6YdX)K}) z7&z2rUtd}I^u;4;6qPl28j?AOxmO5-*7_E3I} z`f>(yfmNU?uDF_VzJ*TDR;;gQ)>4BdT@|oFnSLfje}}}jWq#<>CK}Iek9*LvQA-i| zOMNEt8I1Fzy(^i!cXd9E@9L@=?E~LOwbDAoWi!kp2NOG;M`Jm}KQa#v*52IfB)Yvid$nG#gh3&8_} z;N@;uQ7GY)-F-&>Hv8&;NrJjKNMhF6y4Wk$q9uDxEf@ZoBe@mu>xG3Uzt`3^NV$SX0=Q^)zg~_nF3nu@ zSsB#L>hsF!-X}bSKT^O>&qkswdxjeIjGQICB5P`Da?->%GDrP`Eu#qR z_r*>RIYjzi{_A(hmR+NyQ*Z8ahA7(g**&qkS=*VRDpQ(`t}U2oYOsDM^PIg-?c{cT z5gzYJvKPgi&eqyIeDmhbEz;7Jn2<>~Aggr8?w?#(>P5K*1e8tc9Lio@T|MI3m8Q}* zGU6D-8np8L_K_8kOHoP~w~Ffn598ydghR*OBJKFb$RkIOCe031%BFqpGWH=+&-^~TbngreQ1uh(KSf=%QLh6YfD=${#47VQNvplo=r$YPdv7& z@Um*W<2dnHT<5A_xfRjnYP9%RrA>Z#?smb|;wK~DA6h&>wtwGg(9k$h+^W>!oNV_{ z1>LM>AJ41j5hYY_J+37A8(Uh^$_Vb%g;H;?zLM{FbQ~_h$4DheJ$qC=+c&K*y4W|e$xBXP;Vg}(t}hT zuUqUom!J^%_EMpeQiN2{OmQn1PBkOs$3}PP5mFP=9lduV!m@thfkdvD2Tj&+iWT>nR&Lf&r3!TUc3a|G{*Li&4EwQ{X?Ptxw#pb1ILE~F0<QGV1LM;~oBE z@g@~ha?l~W*}-5@qcrw@T=>j2tKuI|bs6#oWxm^(Wo~UK3D(RCY$`P>hk|Ud+|tZv z1cD}fi^Hqb$_XA$#2c8~GBR?>{a9*e9_s~qaUyK;q;flp74+K|s$&Zsny6r&J?Evm z;^qjh=Xel57a?%hs9*jt6#GC$SvVI}>eQAXXy0tiEc2&_J34+NYwLRg#F;<0&y}&Z zNT$=Ppt6g2Mg*XQa~jnusmHGO&(6x)hy2i`GkVXR@pI%*l7tUX`1@pnFPcjH#@hichKY~w@JMwa~8e3vmq1tv{I8R;SyvXupHahpBHuR>)GYLx@>$GpxYqy_c6*Al~N zQESR;Ba*%Cn~XinJ?R2IU4-?nnaZ&ZdH=I6=j~-4VjrYV1Jfv~Wt}}Dl$QquB)L{F z=uhrU(e~WiBEcPcD&bSw_Zo#=NK|QKuXV9!xJGciC5-(VrIg~g`=j?bL3}`3rU?o+ zUVe3Xm|%gQJ$D}J&$8pQ$Oen&HEL%{?nHZTc!TagEV)hq+pFH3PR#f!l+T?-CV<8r zC_-nsMzzB;afaM?Nt@l*?*(5u0OA>k7@w^i^gEX@R!#_?t$+Dh*4F-F_ORPg*fUF9 z?nKWy%?r=W+%H^(tn(H98;ZVNme3420Rl5A)w^atcoYleR9Vj{%(KX?W2r(~PapU)f2D~H#HkMcGzdDtRS|eKEou_)9p_jT%{F}BOhtavo{!tEc^RQji zptQ8KtCN#crEEQVZN4>81x7%Q?HRS!4OxQz;vw&2sa8ZxU=fkquhRC|HRJIPHH#VU z*&1kjF1S7!CtSRD@1B!jm7}}(S_?N?aXLF*Wo03mlt#dE@W_NzpZdNPwAiM%(_}HS zU@Nh$sPabp{@k(K^4_mb)7<5~=e2A)yVG7A_L-RFS4S?|=F|AS`?; zQ#)NiHYAVrdlEXoP))nfEUMLDhcWaY)Iy5S75OaSP|yP*&AZYFVycgq<)CH{j>5&f zD}JvJx2mj&ay^%?*G@!&hZxxY_Kt%|BTs9hSE75}Rep3^bv_-hADWr-`QR0HZ;C1_+t+33a22>_^LwvxCL% z(x!zjeI?Gw1Rv~TDjyb72P^lade z3==RN$e{6)Fuw;wE3`?L_o+J!T29+e9eTk&(8&;bzx6g)dUJARE(%dMF!1o8RYJ68 z{qk{PF0s$cG|hoIjVE-)n<(kd7DM|i)Fxm=)d^;OP$dlfyoObY3Qowt4L@0eX5uWc zQt#5El(9O>t8b$G{+gG3?0g}v4L~9cO>h@&b8~ZxEqDw3N-Ak3b&bUqdT-7YQ?T7j zeATa^8QRyeutLTe9Ul(4#;Rsyz=>*H6NffM2$l}o*M(I6JX5sjDhLWTNC(p~1=C^J z-wZ}EECYn6BD%WU$zo8_IdHLSvATO{cENTQ1pD+rCUxABv}!J+VfzTC<&I=c^6n^w z9xmfE$r3}Qk9{m-SHIuXK2g>DaohAiEr+tQX9laa%B|JG8#r1`5A=)z%2=Qe2Q~QM zK_kxu=o~ZO&p<0zi*-43^$KvLJzzH((0O%!1?o4r%yh8@a-V3=p4=1IaUVQ1{vS4U zv?kyL1Wq?vh>uf2W?-iC-g=aS_h2ZSSQ~ymyZaT0({)e1IVs(lrt1$A>e?bg~o^y*Mk$!TZv)OS*DZQrp>h!gVj)&@Kkl!E0$R8KBkbQ_8>L7+$xY8hLbRW3Y8IW6V z*K#s$C3ywtVqsk-J(wO|&93J!0p0|3r)Ng62my3K*K8#{!?UP*FUpW}Wmn8;(Dq#HBV&&tOKyy|ge_$& zhYyctge|nC`Tq7rU)2o|YZH%Pi-IPsoVXBQUfpl8!N%rf#aJ+oP0az;o1wyjn-az0{&u3zD)^rVnkn%z%1|5W*COZowA*dQzXete$d+@H57aDO8gEI!gR)p%=fF(XI!J5h zO!cOos3t8VnBB1fWBqcwR8|cR&^os-aes=sq{1zIR%B+pzhg=tp}DVu&X_W z^cxDKIdfuV3=6F>E}qFP}ZrAl+KBlyx_qbCJF0QQd6r+ zA1b_Zz^W#|9hP^JN+FIUKOoAz1}fC!5W)~8Fk!FZI$gXaAwwr&Cj|AM0uC{5ot~W? zVzkcad2kygmCqMiq%PEh2dR-jGCmBKP9D1{sq{(}vC*{N=i;9&N%C`szw)Faa`AG zZoelF?FF|bVa!YIGf@r?IOaY0XuKoY6$tO^RxLxpbi+$5;Z&y39Tqo3O4rv0F9QNG z)?e=mOD4hD1(^V%v+{b!G7Tr7sUUGUTv7e%G!pgQJOFd0Iv^euc>GOoWB1sSjrcPv zqYO$okOZnk`sAdW6k=@-aK_pINibmK95(rDRDud2q*<+pLMb0;jfsh=h8J(H=Yj}v zY*EB6zel#j`KGfZy7n&S94GEKQMmDa!WfW(=}XvT2P_&!4ycsUZnVhpljf^NW6veP z!V}2#T1DJ=`-48<^Rf&urx|s)Gp5Cg@d+D!$=fWeaUT5ku){R(4N%}fZ+apVVAak> zlS3Eq4uS(8VQf?3)h3BpZY;2^K7w2T>~t6utwXu@;A$yhesc|xX~%8{As4SACnwh? zE^Tlh*pjg=bNTxTa%QH=3PH50$ruH(t`0Nm>#A#UhR5Kd1Qf!VD(**jsVj=c5>Z`> z!GrODp~8SJ?2NWwdU6F%0q=vhRfjAc!UOXu9!9)svNNDfe%d^TfD4BUs_*19e~Hui z@C&F}U~z0tz?gg2&Q=VDjK*w}Mf0anxD=1_6j*wwhi8CQHFv!vkO;)z!b(oGLS_JP znt^(r4cC?MQT|jb*Ks>G#&V)FwbEgfexp@!1@_S(nvo9hY*pff90uhR11=drP1Os) zR~#39FjKgcdl^qD!EmP2{Mz@cV}ii)2fQi)7D&`{-QX?s2`hvDX^Q)^hY%okc!*B! TDKcAs_5^v#^!uW}INtg<4VB|G diff --git a/docs/reference/CmdStanModel.html b/docs/reference/CmdStanModel.html index d8b10625e..9d9cdf30d 100644 --- a/docs/reference/CmdStanModel.html +++ b/docs/reference/CmdStanModel.html @@ -125,6 +125,10 @@

Stan codeCompilation

MethodDescription
$compile()Compile Stan program.
$exe_file()Return the file path to the compiled executable.
$hpp_file()Return the file path to the .hpp file containing the generated C++ code.
$save_hpp_file()Save the .hpp file containing the generated C++ code.
$expose_functions()Expose Stan functions for use in R.
+
+

Diagnostics

+
MethodDescription
$diagnose()Run CmdStan's "diagnose" method to test gradients, return CmdStanDiagnose object.
+

Model fitting

MethodDescription
$sample()Run CmdStan's "sample" method, return CmdStanMCMC object.
$sample_mpi()Run CmdStan's "sample" method with MPI, return CmdStanMCMC object.
$optimize()Run CmdStan's "optimize" method, return CmdStanMLE object.
$variational()Run CmdStan's "variational" method, return CmdStanVB object.
$generate_quantities()Run CmdStan's "generate quantities" method, return CmdStanGQ object.
@@ -145,6 +149,15 @@

Examples

# \dontrun{
 library(cmdstanr)
 library(posterior)
+#> This is posterior version 1.4.1.9000
+#> 
+#> Attaching package: ‘posterior’
+#> The following objects are masked from ‘package:stats’:
+#> 
+#>     mad, sd, var
+#> The following objects are masked from ‘package:base’:
+#> 
+#>     %in%, match
 library(bayesplot)
 #> This is bayesplot version 1.10.0
 #> - Online documentation and vignettes at mc-stan.org/bayesplot
@@ -170,6 +183,7 @@ 

Examples

# here using the example model that comes with CmdStan file <- file.path(cmdstan_path(), "examples/bernoulli/bernoulli.stan") mod <- cmdstan_model(file) +#> Model executable is up to date! mod$print() #> data { #> int<lower=0> N; @@ -251,7 +265,7 @@

Examples

fit_mcmc$summary() #> # A tibble: 2 × 10 #> variable mean median sd mad q5 q95 rhat ess_bulk ess_tail -#> <chr> <num> <num> <num> <num> <num> <num> <num> <num> <num> +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> #> 1 lp__ -7.30 -7.03 0.721 0.380 -8.82 -6.75 1.00 902. 1006. #> 2 theta 0.247 0.233 0.122 0.129 0.0786 0.470 1.00 762. 712. @@ -305,7 +319,7 @@

Examples

# Call CmdStan's diagnose and stansummary utilities fit_mcmc$cmdstan_diagnose() -#> Processing csv files: /var/folders/s0/zfzm55px2nd2v__zlw5xfj2h0000gn/T/RtmpFBtN6X/bernoulli-202307251434-1-21173a.csv, /var/folders/s0/zfzm55px2nd2v__zlw5xfj2h0000gn/T/RtmpFBtN6X/bernoulli-202307251434-2-21173a.csv +#> Processing csv files: /var/folders/s0/zfzm55px2nd2v__zlw5xfj2h0000gn/T/RtmpGdLwVH/bernoulli-202308230921-1-6956c8.csv, /var/folders/s0/zfzm55px2nd2v__zlw5xfj2h0000gn/T/RtmpGdLwVH/bernoulli-202308230921-2-6956c8.csv #> #> Checking sampler transitions treedepth. #> Treedepth satisfactory for all transitions. @@ -325,20 +339,20 @@

Examples

#> Inference for Stan model: bernoulli_model #> 2 chains: each with iter=(1000,1000); warmup=(0,0); thin=(1,1); 2000 iterations saved. #> -#> Warmup took (0.0030, 0.0040) seconds, 0.0070 seconds total -#> Sampling took (0.011, 0.011) seconds, 0.022 seconds total +#> Warmup took (0.0040, 0.0040) seconds, 0.0080 seconds total +#> Sampling took (0.011, 0.012) seconds, 0.023 seconds total #> #> Mean MCSE StdDev 5% 50% 95% N_Eff N_Eff/s R_hat #> -#> lp__ -7.3 2.6e-02 0.72 -8.8 -7.0 -6.8 781 35502 1.0 -#> accept_stat__ 0.92 8.3e-03 0.13 0.64 0.97 1.0 235 10662 1.0e+00 -#> stepsize__ 0.95 7.9e-02 0.079 0.87 1.0 1.0 1.0 46 2.0e+13 -#> treedepth__ 1.4 1.1e-02 0.48 1.0 1.0 2.0 1874 85179 1.0e+00 -#> n_leapfrog__ 2.5 1.4e-01 1.3 1.0 3.0 3.0 89 4050 1.0e+00 +#> lp__ -7.3 2.6e-02 0.72 -8.8 -7.0 -6.8 781 33958 1.0 +#> accept_stat__ 0.92 8.3e-03 0.13 0.64 0.97 1.0 235 10198 1.0e+00 +#> stepsize__ 0.95 7.9e-02 0.079 0.87 1.0 1.0 1.0 44 2.0e+13 +#> treedepth__ 1.4 1.1e-02 0.48 1.0 1.0 2.0 1874 81476 1.0e+00 +#> n_leapfrog__ 2.5 1.4e-01 1.3 1.0 3.0 3.0 89 3874 1.0e+00 #> divergent__ 0.00 nan 0.00 0.00 0.00 0.00 nan nan nan -#> energy__ 7.8 3.6e-02 1.00 6.8 7.5 9.6 775 35215 1.0e+00 +#> energy__ 7.8 3.6e-02 1.00 6.8 7.5 9.6 775 33684 1.0e+00 #> -#> theta 0.25 4.3e-03 0.12 0.079 0.23 0.47 796 36197 1.0 +#> theta 0.25 4.3e-03 0.12 0.079 0.23 0.47 796 34623 1.0 #> #> Samples were drawn using hmc with nuts. #> For each parameter, N_Eff is a crude measure of effective sample size, @@ -365,7 +379,7 @@

Examples

fit_optim$summary() #> # A tibble: 2 × 2 #> variable estimate -#> <chr> <num> +#> <chr> <dbl> #> 1 lp__ -5.00 #> 2 theta 0.2 @@ -377,8 +391,8 @@

Examples

#> This procedure has not been thoroughly tested and may be unstable #> or buggy. The interface is subject to change. #> ------------------------------------------------------------ -#> Gradient evaluation took 7e-06 seconds -#> 1000 transitions using 10 leapfrog steps per transition would take 0.07 seconds. +#> Gradient evaluation took 5e-06 seconds +#> 1000 transitions using 10 leapfrog steps per transition would take 0.05 seconds. #> Adjust your expectations accordingly! #> Begin eta adaptation. #> Iteration: 1 / 250 [ 0%] (Adaptation) @@ -399,7 +413,7 @@

Examples

fit_vb$summary() #> # A tibble: 3 × 7 #> variable mean median sd mad q5 q95 -#> <chr> <num> <num> <num> <num> <num> <num> +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> #> 1 lp__ -7.18 -6.94 0.588 0.259 -8.36 -6.75 #> 2 lp_approx__ -0.515 -0.221 0.692 0.303 -2.06 -0.00257 #> 3 theta 0.263 0.246 0.115 0.113 0.106 0.481 diff --git a/docs/reference/Rplot001.png b/docs/reference/Rplot001.png index 17a358060aed2a86950757bbd25c6f92c08c458f..6823093a50447db6ee651a92b3785b9e33069224 100644 GIT binary patch literal 5039 zcmc&&3sjP6xBj%!O3g0+q-DmOlb88_ODzm16{)5)O0>yFGw-J{ z0KN4~msVFFP6YrPJm3fR9Z$}m?n}FL=X?tH-9VakC3vK-hdeD!<@O|G<>Z{~{eI)c z1D;lq+FjwttbML5u512$|LXM@O>cS~Jb3-$$vfYS@FJI=IkN8Tq31uEF|HqcZ)l?Z z{oc=t1|h}o&UU*;C5=eKO2HUbiXwzQ`E7k((yj<9iNbb0IUUJza?F)cD589Meu{_o zFTrZwnvY9u7PWzHX7w4VC9#({it_AV0njh)lWN=FOH%D;r1iEV`~^3ISq<|+n8?kJ zT4ix7XRqoD#3st|cvo(%H(PMcOzs|R4nWpw6%a@`7bQI+kXsA3N{QR71j@?i?im4@ z4-aK(6=nj<_G)SAWNK^c11pvr1FPqE-rg}|noEZ{7z^4Q4Iw@O1*4CRKv67@;2oJJ z1pqp(C2sMBzrA1^f)hP+jk8D5RjW!T7!_9-HygNT!b|fKrqHPkRh~XNLD_2u4Z@vR zi;f^+suGv+^p_`R zA#YASLZO|1#?C0X;yF86ZT>s6mKKJ%pLE&gc*u|4ha|DvISx-Z(`v(@$@Bi&8F{z| zrLzY)rQLijMPF0s|B;(t!TNh~{HYt73sbdtKPozktTA}N@IO>knC3t;^!3G%^; z-Y=IM_qTSIMkEGD;NA{i-ZJH&%pqWgE*x}ca?D`KdCt{UsZDDL}@yNry^#k1T|Hh|Xnxti!7m@-lfxScpU zCr*%Q9dgpq2h5jv_ve$!+g5+k*O?zfO#`>=zVv=vC$oM1;QbRPSoncjM?HT5%Q}gw z6(n3cc>jiJc!KismbVV_4%8lKP6qI3=a-1p|Hbz|@Aoga!EeC)Ycc4&)&mA8EmWlX z7+Lj_l@%S3ym)$?TU@q2)?m3Yq;{VT)v3w~Um@(xMvvY#)7C$@WDVH$rd#C9ovWx= zAA9VVS?+%_zn5gzk-9TIUSDrBh*X$w`e(!d`IcxJH$< zirfEw>Lca(7#o%G`Uf=@3ri&uIdO|gPaw@Y=FMd+i|GEO!aQ^W#dmX9uttteh(L=2 z85ABL`yaXt7F2#yFS6Ndc<0B0_VKp=yJq=x2TlJz3Xo;gEi{IkV>CJLm_sjbmR5!h zS5#8fvF)^}83*~Bo=YX=hJY{s`l{IuTcbp@5N9}Fjrpa6cyE?E_m?NeU#s4~CY^s> z47NW706*>A9HFM<(gP+-q5rD*#$(Z!z9eEMEibZR`)DA?Wfn zACgj0RlmEgK z?mQ5^zHVJmP>abNF}Mw|Z>Y*DR`$mBPA>#M{1x?PjcuFFy;_+aZa5jum9KO(piktw zDCRKKQ(LKBk-b67jR)5kL29wiqBl?~gVC`ghk!14k*ArRxKIC3y#61cqSui@F>?0y zFdjB0%+;Nua_l4_r=HHOvfX3%Zleb7#G*tvv30l%SA zqNola!rNYwASANpksoV%ZFzDR`=R~7Rai!QdLowOyVcy1JwBzGCNqW4{xrs?2}U8n zAj2M6DQ~u^cCBMc+Z3=SIj0ki2q~$-y!D)h76q(HqM+XIk!VySucVygqM3ckuL`N4 zM z(3HTYkv;vy2+ZWOuVogs)>?^6l4{4Db*8$BPZIET0id~Vc9n7?ArLeNp4`QyZ+dVE zBkB0uCTZbDdr_}5_<0S3_3%$zKF zaX2MTEuBa_NE!QDwno<{k^uG}ZyniQ#SHqHeOD7Te7uzy1zHnd(eB4DZt5DpUacb2 zE>+@r>?8hy>y|zCnjOnQ5d$4qm7gzZls|(dT=y?97wHIz|=>w!eVByEY{8Mi!#26jctpgl2Rki9McS6 znYs#ibj7b7I@ok$R6wZj8M7%Mkv7Sc3yTVg0R0w~5_3>bPtUPGUDD|;{|9e=Tp$p5 zd0`#E&sw|f5S1B98QGoc+&A%HgK!cG_v(M&idq#*On_ly*-&|X$$$uyK6$sqc>sGo zXSz>>lf!UgrX|d)L02g3PZKThj~qD?nDTz2MC~oLen5q$$vQv|32#BpiczAWm#E?_ zy<@5PbSmQDIPf^#iHCUH^dFLVvJ*8eP$ z+dPfyi!GLIfBQISHl9ysy9-4nukd#RGNoL*PJ5T4m#EiiWsS0_9n`8fQ}p0}N54GWkHb;Xp7VkzEE25D0b)x$xOY2 z$%W0HHr(=)PqwJbSIzF{i{zq#CcIj1rT{OY8&Li?{HcZjg9-b~TlR<>h|Ez z4L?1Ls?W)pzSlOh6j83JqpagAL+C5VxZwk$rx!NReyrN#C*y5};Z!qbK7M|xTg6D0 zspgw67S;7KEEmP62U6TDlJGMasL6Tg1%{ivkff-i7ZU-*3_bW)sJWs4vT>8Wzk0m8Uu_J$ATPXGuhNM0q08&Cx zw4^og3H4V;j}pj)e0HRt`Zz%zsFJ`AiLcaj_}_zzF>qYaJ>O-`7B$8^xJ6{+Z?{Lz8NFzdVy{C4JDGniAuwbRz!!jnHYzem|fWUxn0qkaTFAzeV= zC)S6>xDss46}~0fQ2PF4@wz1w{Gj^*+Ww&L{t3L(Ss?%b literal 1011 zcmeAS@N?(olHy`uVBq!ia0y~yV0-|=9Be?5+AI5}0x7m6Z+90U4Fo@(ch>_c&H|6f zVg?3oArNM~bhqvg0|WD9PZ!6KiaBo&GBN^{G%5UFpXcEKVvd5*5Eu=C0SJK)8A6*F U7`aXvEC5;V>FVdQ&MBb@00SN#Z2$lO diff --git a/man/CmdStanModel.Rd b/man/CmdStanModel.Rd index 94ec10d61..2e8176aa7 100644 --- a/man/CmdStanModel.Rd +++ b/man/CmdStanModel.Rd @@ -34,6 +34,13 @@ methods, many of which have their own (linked) documentation pages: } +\subsection{Diagnostics}{\tabular{ll}{ + \strong{Method} \tab \strong{Description} \cr + \code{\link[=model-method-diagnose]{$diagnose()}} \tab Run CmdStan's \code{"diagnose"} method to test gradients, return \code{\link{CmdStanDiagnose}} object. \cr +} + +} + \subsection{Model fitting}{\tabular{ll}{ \strong{Method} \tab \strong{Description} \cr \code{\link[=model-method-sample]{$sample()}} \tab Run CmdStan's \code{"sample"} method, return \code{\link{CmdStanMCMC}} object. \cr