From c7f9bb5b524d150fd9f6dd3604e76c91833b3d4f Mon Sep 17 00:00:00 2001 From: Eric Bair Date: Thu, 26 Dec 2024 12:06:16 -0500 Subject: [PATCH] Fixed some issues with tests --- DESCRIPTION | 2 +- NEWS.md | 6 ++++++ tests/testthat/sim_multivariate_big_st.RData | Bin 0 -> 24627 bytes tests/testthat/test-Log_Likelihood.R | 3 +-- tests/testthat/test-PrestoGP_Univariate.R | 2 +- 5 files changed, 9 insertions(+), 4 deletions(-) create mode 100644 tests/testthat/sim_multivariate_big_st.RData diff --git a/DESCRIPTION b/DESCRIPTION index 0c53451..edeff54 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -1,7 +1,7 @@ Package: PrestoGP Type: Package Title: Penalized Regression for Spatio-Temporal Outcomes via Gaussian Processes -Version: 0.2.0.9038 +Version: 0.2.0.9039 Authors@R: c( person(given = "Eric", family = "Bair", diff --git a/NEWS.md b/NEWS.md index f96c5f9..f6e86bd 100644 --- a/NEWS.md +++ b/NEWS.md @@ -1,3 +1,9 @@ +# PrestoGP 0.2.0.9039 (2024-12-26) + +## MINOR IMPROVEMENTS + + * Rewrote some tests that were failing on some verions of R + # PrestoGP 0.2.0.9038 (2024-12-23) ## NEW FEATURES diff --git a/tests/testthat/sim_multivariate_big_st.RData b/tests/testthat/sim_multivariate_big_st.RData new file mode 100644 index 0000000000000000000000000000000000000000..88e27d7963d08d97369cc203c181cd287142ffbf GIT binary patch literal 24627 zcmbTcRZtvI)At!XxI=IY!GZ*WI|SF@?ykX|!3jRNGgxrfpc4|@-C=Ma2G;@hd7p2$ z_ICF@SKZa8>zwLe|8KgQChq;a{~jbL+)r;j2>zPT?#c>%=6+V@h2}CIq@}U5(?^`K zrzyUdNhy_2!DF}LEJ_*WR*)eonw5>?b57x_1>fJk1SdWJ27sq-4{x7$#m*mhg{ZR! z1$kh8`BR5})OEX1r{+m5i+p4n8W1fA2z&=BH=j)W3yLK<3?-k@29_cZ}G;WVSqxkiBw?A?Z!Z z6n1m%VoKq+mplY0dK>f0EgPEi{U67?km^m<8jAEJqn)dsq377?xzibf5n`a{irhTH6`Syb0f!vi#dJy zG-2;|SvIH@fLR_S>Jwmx967geyv%GqzIQ{dj^nyry@mBYJ+c?FKk~mKE?ThpMtCcb z7M$NWbZq#_BWC(f@$;Jhow%i6mg8b#M(!I~b;-E5RZZ68)bf!}DF{AtBW^hs=jiiM z8}MT5WNKOeajNqM)L6Z*Jg9203>e%5KGh@vC|`MApVk6q{p$JvvIXM4uUn7dhq(U> zWA{z;b?Ys689YAm^S{-PNPuQVKB%Bk`|Vwt9$D~LKgMn8|JV7w@|?VgW+RX!a`nKu zZeD)U=0}$Q-^TE-t-kV(Z`t{8WAtt+Sc#Wh|Cj6Z+dg9Iy2j8otmw9c?OvVTER4<^VvJe`V$5YO(VmFM4OF9juN z_Gr7j1<|RgmaGLOBccFQ$JmOH4m1x1pllYfR>r(T&sG;WS(7Q{k50YXo}HKuZ1~ANh+DeKFVwHpqy2|o&a1C! zw9XYM`G1A&|AzelND|NcN>h3>b$|GO>$RQBb?MCd|5Sg##d{ZmWB&ipLlsUA4F4~% z^9BZ#MP8PR|4OvpOS~q`knz8&EUNbr`d8x+^50ue;IF1{7hOpI1%}n48TW(t`Okm# zjvJVFhQmU(_1`z&-U5;y0AC*--llfxqp|LMp_p*-*nYMPg|}-5vox|_owZb3cQ3X0 z!)z)jjA2)PMyFiu(pDY+O4}IP4UM<@`tPLVol@?y{D1VQ3;5iqbx$8ibyQswj5hxh z1e3v9ql=Mgby~zmMtnb9>Au)_y#?V;AP!)2S}lf6Wox-l+VlF1+lahg$EsIs;5p=D z6uGV57%Ql@ofiI{Exp{#bYUKKYy)7aGpG5So5$D2gwOuzdkLImqMN3__Hb%P7)D&V zNBO8L-QW`(X6i~_I+fd_%#4r3GyT!c3ze3wxyIHvDTiIQXtItCr1+ zmJ813bMTNdwc&xnz|O<8*85u4Kzs4wd^rn{@ef*_H@73o2a84 zM$tWAD|hJ2{^AQ2>tRPyJg0)kn1ZFvdkVL%Vwu*fLUtAB$JLof67$DcPf%B+n-kBR z-uEe1<&R&m4oCF1#8z@d#QTO_&^`aCz*@`FUR(Gvh@EMK)2WH#H@kTvF4d-P@T010 z5sHs2o|)Uj(L3}Q)X*4AzlCMBfS9?O%Kb%&A6#s2h%`^$a~md;8FHh(V?*14dSVf$ zm3Cn3UUlB6Dt@xPyUE;3B*;oV3f*-w7B~KmxdMeimVNdl4&wDU&YKe=!vvLxzdIma z+J4~&E!9P4*&1HeL|II36)`8*h{{P?)Ov~J)uXQb=}=fh^eTdiLA*tZKg3>+#R=Co zOOk1b2>O3QV)Y6BemC1IdC{~#>bNf#bFVu^qHX#w->rkrU2zj*RXfdiO>?!7zk*fX zQ8k)du-2FMDaBT?mSyZn&as2Nq$G9f7bdk0!FnjouAdur z>-7d$9yS#C^H7<)OciyqGzxcURbzCV!{)D&EJ8{#oS9DRxg_= z{sxqNaRf+gprH5s)sJVZl@gs*|6TY}sl|rLO`b+)IV3Vkw95iHu^_b5&GYX{fKxfp zJsJ!WY)%VVe-+=fe}^b*4s#^$lJIzR>-PVM!Af>bMo~W~=_2*XxKR@E#71r5)z{Aa zZtCct4%>c_o2iJ>ui>B1-@C`JEkDzm{wfh=$C`ZDx&0g0Pc3=*!M15w1c6?A#efsD zRPvh2wCDOgQuAZ(r}ldrDcK>3vqgMH=DlNZV9&{iJtXg3S@!yn!9Hk@fb)IT3Ov`c65vi+E}WDN8e*(hxg!GOb!z^+t@z!p{cnfzR|bTA+2*oKq# zD&W)p4*K__V73fZLw5|%)khb`GBK@$wVW3SrsP%0!hV*jIdczO`@(>MkTt15)VG<1 zXq<4-)$*!}wz?9Aztf4iI{VKj_|O^YHuE)_;|^P0gKvTy6Mq~2W%c#UFo61_fABXk zF(2IGv=m#z-KCHkS1$J>U+%QrTbGZ>;vpu_bu~CqUt4dNk z+&oQr%9$SSIVn04n3-%$<4QV9nG*FzNlg)OqSuuK%xK zJBU7g+HN<=D1S;uk9N7F+CnxgE|^6RisQ|l87m-a5SVy#!omEF*sZKRPQvfe@3xiB zgVsMqI>TW$bogY#3kO6|^y;royrGh8$j75r?&ps=p9znT<6S|F{cFP6!(ZMz2p+dXRK|3;x}; z*JaG#vhc^N4l#L{q z`i}~wxkhkeUz%@c>flF!?xp}(zl5wx(Dj5&$fwBddL?Er+04#CBxpqsigp|I{<;S8 zl6MzmaZXPfC%A#(IB;CWt$tBd-@^O!E`7gA?ghlT?^y-2muDyXceba7md+14#F|^9nKYiv2SuhR zy>`Pq=mVX9z^&F3jm}+bE#_!S$nsdaV1LDy_#|e*eG#jT3D>8azivQOK6I{T<=qft z)c80K9Ts{A^;{WFBkD|yY?MpLngeLyuBKgMgJvO#~gn=)D=YHYkmPg{}8pd8=Fq+4DNQ zD@MF_S+#JLM%bb9u<{rw_3_rcaj8sZ)%f)w(B?T|cnnR7=R~K`$&F(+37>-$AvAnCJ zl6nDWP`VO?_e1LNww0)8$v9mCZSpZ@XLV;ZqKs4lmkquL{f5B1bZDB&hreeC39i4i zbbsP1kPoi*OFB=wkC{gj0F|`?ak0jcN5aWI)y9m4TxJfi`*loZ%;@*n!ARpz8|7Ny zOcaNm6SUnuU!Z}~D(@}s4}h4)h);bo!ei-XBQ7OXgL@Mv%@Q8c?Y}y*AcZFE{KNDHMPJpi z3G3}W;yv^HX?YiIHxPg7Gd3I}re#fDaPQ zM{!D3liLF1PeTkZ7I2NpqwE(=Y$gSDxci{P)r!}o*7Zt~C^r(YqfMQ@z72a7(v#t! z&~q5*K$EEpL$d9Ue0#3Hv8tGO)q2=|FdkM-2eIT8M%a}S2_f-mB(LB)3br6 zO6`waTi7Q&j<^B)Cc_;Z&YRu8uvu4eK>kX6?x;(rfy9P1Shpvs_U5S`beTsN6+45% z!>IN53A1(Q&$3x;nDuG*lR z!!MN&_9ID!K+9s708yQmx678B!uMr=Gtpvs3ky9Q)>ZlAvSmz#mc_R zNen%Ec7T;lTUY9J?Bph2BMiv)a4{nNz+B$7@TPq$s5WfeoxG#hlW<+m|4bcMtIqrt z#BtqciR<%~ML*qppiEu;Fz>J@4kRcD6YULBXZxJTp8vt?w=x#%^{~4xy9)`1EqF~R z&h@A8R&OKkh94F;tfVuuom{H;Z%Morcl2`rlHkubZ5}JByUEHzUXP(6Sa+{ocj6ou zyl~16@Q&$l@Pjq#>_Yj9 z(7Li+Iq+sx0U+s5!_qx-$8j{6X0`3~efs7et&vHW;PAzG#E+e6xYxfs%k>iBy^xdh zK9-$?W;(Hnhz$bo?Qm%&h#IJS=btRPl=A%0@Ddv06A6q{SP*=YByCc1H6gZ*Ck){F zcr*KhF-@jk{CQgp4IVwl@8yl}ag%b?OQlfq{f1}>9gB0`|N$wzYu>bKf@B?~Ni z%vMjb9`@=adk=BzgZ`04SCus)Rjp90FDUN_^N&B|8x$|NF3K-pKY59YMj3Wf3;vL=RuxkBscFNn3 zN8LO@Kq0}wy5LXMfl<7z1YikXMR$*YlZ_s&!0O*?Wmj3IjIkz zWd|`?>ocKUJrRS}2x4i!nF9D(b`a+kbyf{IfZCo_J9p3B<^DN^KmfCmw0h#9K|OB z8_5S6Jntpem1qlktxcj&)Nk9IV(*$Fs3+wEE9~P9E8JLfNJZ0uw5JCO{KxY(91E2& zvpKxDv;i4Yw@Yu zlJl`rgomHwgn?8X5jC?*v!sFD;72$t@$5|CPE-m|><0i$PN1ZfQ0f0k)eMkHt zwl3R$CiXzFB;+*lm|ge+1@!liendqQolO20)yTiqRHQzPo6+;G6oSqan3MInK8oG{ z*aW4W!@raTHmXC#JcjlEkP%&buLz@+VMfYd8#;z%xUZe&i4z%}g{38p$wl_!yWHRP zR9hi#O8v&j9pMPxm@BP+^f_I475n3Sx_z!0VB-gKSAhC`Wx{AdG2n(F7CG z)Dt-E;Mw1Mu%4Vg51xwxH*OvN1Q$>Cdqj7UIWIH)zMq)C(7R^sQN#H@|nnK&;V7wrTRz`&btDt$u2QSOYr!amveo zC_)RI7$0ntCdAB7hqFBcdi-5|ILj5GiH-{^94=#FlWeg)MR}@A4Wv}`f$Htxy>7*- zCiMfYhNbm^xjdm$6!E|(5kdBWaNnRVvL=87Y2rvK=HrzJ^;fk|QpIL__1ket>ytbj zph{6zbtpqn44Hf6F#ax1IUM=?u*K3wnUlb$^Q)s5;E?(PX0q01=S^{sxp7)@*4P}Bn&4-$%opci zPkU`sRxY$BXqp%!fB*+E(}`<#m3q&2JSn}htM9fr+BdF)r6k++iUfogs@8&2mUZ{~ zT(!~H;R!75iBs$mR!f2RtPo!p@P-5Q-^(d&aIu4K%evJ0QQoE`u06JD=3rcZkTZ9l zdGVd5unMgf2kZk62)Qy}$EWk~#!XhYREIym8+E(GwbPlv`v)4+G7Fgj>$+-R>tX)y z$VkM3P4j26{XB^4i0tWTJ5em_2eG|=Lk~9O*5Jm42LVz1B$sS|Z?|YoRZ2cj-2R%E zNNm`_s{81i?Iw#+!hUd(y`%V721(V17LKQXB~o8N44W2RI}-j^QZ?C}8PDddW{0ap zKTuOS;Ones*tck)Yrw-V5)4r>_;-(#iTvz{_p=r%tJl67=Ymel7gT^ms{4Or0t}Fs@B4{F~q24+&9J44MTK(EU|}MgL&=_Q(Kj zEyQnAr>tP^={mzTjTC(+`ra13+FsHnt_K?Sfrmc2x0~F5nPzmdn{<6>eQmm(O8=k1Twd|o>O01^ zgJetU?e*;mn%U+_vd|g(ZGJO^OY#wL4DqHdCb;(g?zSic~o8lkHX=AXkoR z>^4`;<{gSYI@a|lL{{uY$IzWS{`{1$$z41ub!K2nyR?qMoquqgfh4^SCd`qGJ{h8Y zSoY#7phsNEZ8o;oIdu^`G|)zBn%+CclCg6`K4tTavc0w0duz{ zRh}FL?#nsnV^b)|J0^kpjmJt)=s|XMhLoNBbhyn!e9TM^o%ql9Nuq*(%ha82tM9<+ z2ziY|hE4(n`|GBs*F#HS?+Ax9Zb)a~r*9pkh%WxDc-IF$1J^3tz~66OYYu~`SCDgd z#Q-*>lW@$TU;uD`eq}iL^;#9h;u@>uoU&9Z7!o?>m5zirgHfpP)Wedt@5u)fF~V|A z?`{lZM*z;e%&__}lCeTd2p5eG7BvfuiwGsc z%>FTr9QJzIojmt>>VuR8`WyN&*nY}sa=7nRQ0L71+kQbm9H7=d3-XrQkC>Jl>;3`+ zDMpKu;f9>q&t&2jyhW`49rA89tf6yT>SMC{gENZu$r7fmMslTcXD#y-`$H(Uo3V$0 zsmSicp4Q+5L2}!6urtjkPxT8FxM9m$)atWXuFF<0ZQpBRpHwc)9JkE!T{g3W_lg)o zQk`5l(xyr&*p5)FJNsE25{QBBvzsyALObcvXPJ1d7VU1ksov(37uY`{@>?TqfeHir z-^fZPDZ2UHoYD|RM7%)03QPHX`OqI2^v`DhVpzZ2%Ps?~z39{} zB>AJk`NdA^U3g2aF?4|eJ<{-%pngeCqRCq0^f{bm5}0drs-c00+&bA!G{Ls>>h{*d zQQ1iv%QOqG#JN>}GF}l>U6VM=XQ)C~`b%frHsP!Ei=v3 zhLO!Us$c7FaX7({uFK-qA#`mZ4;+`Yue>Y}*7WV;i8<7W9^xRPI|@Q_R#-Vmd13nf z>tE96;KK+K3|#1Yzd5ZjvLBc2TMqa^F*!;dLQf=BU5WFPl=)FW&9UPkYR;2V}voQ^r!lVmSE9a z&CbL=)STh{>?6FY7NP2EY6+isVmw4=0axs`s5>S~oDf>Czj$zPu~4R_*KY|}U;t|=Jtp^{m0jLQZTVcR|5H)pzh4GFejarZio-1x zSJd5udQLXCq%NvGbr>@;h#k@o$I}-<5#6BrjF+kBJeFEWNy>rKWWC5$Obc6OB%&kH zE9J#Vh&xKpsb3M@UMr$?2~{55jxoMfcXhI%AKl=Df+gkyMG7X##OS8>Dt-}Y9+4DvkVMd;^5LNxEI=72>s1Pd)xG3xTjqxH@e+iXL`{Y zaZ2>99w5hrJ@$I>xQg8Ly=UkYmXK9@^Jm-11MAzSniZo9=AN>QPq6FiA$kY?C9vVX zRsC47;fO1M`o`uB*pYjLwLmiTQs>}7doLMmDxbi|Y)_t(dw?M|p)=;}7mb`XxtYl@ zxhJiyLOwPKw z^hY;RG+=fcbtWi~0bO58$J^|_!9uX5VkDsWnyPT4SC(C6gnTvAMQ?psA63)3aY3#_ z@*Cjo>fHJ_fu}aLK!u1$Xq)|JaW=wV*T*JV!RK^I-S%EyM+ux9_UE>dq(J0_cmn&L zT7;V)s<+)ptkX>{bDDmvD1%LVvI8`T--7+SeX)gc<7FS^pT%e@_VUCuH1M&KlNRrV zX$GY|s_i!GAAbRzeOdW8bTNA5SZYmZ4kUXx`L3g<5$s}uYKE2H-c3;ZX|9*m@yGjQzlkk2SMy~Nx760$eCXE zTctN3neNU_CKRQ^X#%WNov%BaIR&RB;|qo8+dm%oEjbl+Io*Bpk3rj8!%gG?4l`+Gdc3e5frc`*8(KHSH9N3P9MxaoKg1;#wTGQdIpg_x!~Q zk^{p~cKkG{WVj!T0xl9|(`uL!Lr>uqn|X;Aku=oxa8`sedA8`D8F(>bbduNa&ztB= z*dB4t4JG_G*(0Pl_;3*_717QxD%{sJWh!9VN>H_6EAz`*Oq$0KaGu{uQz)8yf9Tc3 zV^;&y#6Vb*{BAJM7yTR}z|`DqNxau*)4;)BR}xTTHTCvL6G$(BS>y$OmWsv%az3VTj zs3Il?0>7JFWY2gnsde44P^TYv7XSD$x%{b|+Hem3Jgl{iAeFaEma|9Cq(Ad^C zmL!A{yfKQ|n|@8pcv_1yQVw@Rdz;%Nv@K+UvoH8T44FynjKy-GSV41f6iMu}?^c>Ign2uj0Q@q**H8OOX4(&bNb2lE=F86j8h7OeF~Sw zF-4tSko-%l=Y^W6Du%-Nu%9%Go`V3p>`%DBiP1qD6GBnd)-Qc1mxZ}PA6<3R-|0us ztmGNurY9@LSdI=!>dPs0`77a%Kenp~L5(F+{6Gb&V#pf|V3P!`(a?g;q5(7fj=N41 z2v3!eo!sm{Zjq#LI)S&Ogk6;3AE+6e_=Y&f5ZRZPd5a2VjZeDvETKk0toIkYBwSZU zQ|!Z3Xmd?HfnQ7T#+NJ+qv)Lr>sMZ_)ah`>Z@V-Iaq7dmyPs1`dpFfP_@<7UyGV)V zJ|=#q0N^TqjG66ZdGD!4noPX|OC?yC`c|74Rdla1HwluH)U1ECMqnNG%68QcmhGpk zqRDI{4!`SrHZ0a_j=3heDvE9;+Uo`tKN`kdo0+^=B)saYFWmKmg};qY#74I;*0UGR zdVM@PN(U8Q+xML7v7%tx2J7Z*1~XF`7Sp#BQ_m;Qg(`;k-*R<6D>EM>%cIPuCw0Ya zDl?JhL3-X{fUoADe2o&KgsCmj{hx?zD z@)e%7IM_j$SG(z)Dqyi_lZfMsWpn~90&0$tuSQ5!Kvo%BPZs%dj+Mb%P`Hcs=?;_e zWF&!XGmc2n#wQes%>=e6OeCrvF3Rg#5d`3B``*pudrID0EMRAoM zMjB|Vzr47cXuQM7Ryn=i8VKe^jGPdkW@(jR6)g4j(StvHH{(K#= zZ<`lUs2_**pB?G(G*SikShk;93B1#?h~+QE z_#E6*8PS;w&2?R6n8rVg@vr(WzgOaaozOH@WI?=oDNz%5?J;S4+?q zJ$2F!@Nl;W(oF;?-9slRVtffL(_oZZvd)#!kJgg9Q7y7~_gah2%K7TVOopNYgyGHd z<`#vqFfw$7L(x*5eTobd3JDA}mNMOEIfB6ncy?OTo_!t}x_ATHB)pz$1anCMq2u(p z+Zz9rD<13FRRCfb`*#;EGU4pg@D$@J_vnXL<-KpcGPaU!40}1G&UmYkq|YWBRqt0$ zDV|crzJ{ws`>}s3qvo!nczq=fMxMQDN3jSmr#%&VJZ-h4L;lR5m-cFx$YcD zAQ=7N*IckzXl0?(1h|lt7FFYF^c{l3d*%d3Xeu=75BOn24F44uSih*hMc+)|uwVk$ z-rHCI*$z=*A{6_6P|_>-8DS+MCOkk65H>=2A>x3-j4qnqL#I5R9%NZ3{t{WZ=x9hd zT8F6cWyP-@`qd^g&Zn)bHvDN7=bvH>C+!gD^8KA3`=Qsl*;=t_CGFU6Sr&aJ$=~*) zFLd@iR+$xCis8f(Dk3*I|GcA1r#e8dU--5RnN5NYkMO>dVo|n;7^5skS;aJb6w2~q zt+J9P7{1!RdTKwAv{{ikW+E^Fu8`I*_ZHKXHpk^99u8Hgh^b?hEj1e3mrq^a@*u!8K5yop>B>$6Z%tV>G@V=YL19nKUDlQ0KSv{gizB8d46Ff&+C-|Y)l(|Dg|lFaST zf)|T51RnBvBemIjRSaoM7^&D8r>dE#ZE0QS7j$TtbbU$xH@biwp4QSQ0UN)FnqaSv zcR5ip87`HpcuGGF$j(py9_Kf@QPDreBibtqT4P>5jUxpP%yAZ(H68lS_i%Ax6Vi{; zGXpz#3BYNd*hFr9F}ia5MI2Kxr!~1~`&gh_B@G4%&DbZM^x}85Lp!=J9RFE5tJF*E zgm@o4PMTYNCSt6sk9nAiG?nz93$nPZ8R&nt|J=1)(Vso-*?z306VleQ(Q%1{gO{@~T_iL^1jQ8`BU;iNq&O*tMV?3UqlIREV4jE}l3(+cy? z)^r|8l?3?_zB@OIDe)hRgoB_=yii#-{G?m@|)ly zcUfpl|I!jMr@4&kS@T)l!&v_F)=?cHXp#4U$wmI)J%X5*Z9 zWiBGbvhh~7P$cw1rXWkP6#|@>UQq*kA5w|?vcyC9_P2(tA)^7T0a{Lvz|-yl1X{;$ zFt0HaY@}2z_Q-aAbF0PP+JGjEvl|~w?9I(+)ppk+OcjzdANs&+=OsyTH6!No-+1%% z0Gn-5^rOs3;d0z@l$Tv9pBTSHt`)%Ru}if)FtX%aXw)>s=5-ET2@{?9rDM~Nk@MY6 zTCpX5SZ*IqCAgKZ7++Qz_vOGxg;3B@I_tbiDC<*qgOC*rJI5^8g^Y({yr~H0`inir zq)10pNjH0#&9M0%XGltI0@gpW>t;4az6S}?ehs>gg7XXqf=PP7e;G>oWKe>P>XAo0GT#Xg>+>S=BDkneFKE=tm!t;yC zk4AAWr>>3>W&Tt{BLtI4L3+=7t{>fHyhjgtUP%A6J+r!~z#=DiF>}AA;!4{_MiGW< zVH;~QShthWqb1s`*RhYH4e-Dt-G38*khE07sjUfLZ1f9g!r>~sWCYY}*wJ9MQ)%Kh ziy$f{nKr(Ma%BaVXa28tbZi9}RSn+hJR|YC znP8IM{kwlT`#7rpo#k*H9!NpHot?xn7#KOquC35n{jL#ozRBEt8hUo}x^D8U)zIDU zd69r4F(?1>dpJ*jderLS9^NxZ1rbvNiI3o792UELKPf4^yA_msfg5lYx81F3ic2~R z6{J?2q7~WC2sXRe;`=oRte4uM%Oh%tx(P27vPxrk8?WP%J##k?KZ)jL?=K4b8wPZU z>nf1{>r3MKMLh^8Iwi08UVOZZTOZ>v@eEOatgObi7X$(tuO^)C zk%4E~m`H>Cl9gdyzk4f5M>9T+GybJ096!#@nq9}xDC!;73gUYq1G;sOipV~g(mi#&@q-t1GLdwQ8WP?F?Hs;v!<$XME#(Q!v1VO+KJCZ zimTE{O6+q6rzZIAT%X`af2btp1RVC;YdqB#1G6&07JmlggW7 zd^T3ImD2a^g%!!YObfQv%#Yn6+K1@7iZ=BcstOB^rW9SVxlJ1XOo?URt%#-b zTzLdr>nCDb(SELW3vQ92C< zN}d{0OA|MREt1N;O+;M+ZNXjTTlzRFdglvIL|-b_;|@!S1l^;P^m~Y=Nq4bIvBWy@ zGLbr_a7l%4&qEw~lZ2X+5f_LLO>O1yo^_hzmH9>AL#4B6rl1v}T-Qe!P1OYLUhADM z2TLT(`Pb`vJ|U=g8^gJvQ_Z3i&Cb)>p9E>e*DpAqUEe(#0bE{C`0dtXb52tmF#J>4 znrNUk5sO_zb6u{R`&epnbGF*w+y^JADq|iDQLg%6NyjfJwCu(mzXHdExTsiIe<^Do zEG`wOa^P8-Tk}^rYL_izGkguKG2)O_Hrl9|k`R0#cQ@ zURBsfM{NlM)coHuyUf6}!dfN3%~QdPFOt(aPra(L7~vgp`&=cv*eN&5(p6U+4SUst zCif_kqI>C^r!;61UIrq{-=9%#@00%mVh#JI#ZGR8KM^Sg1ghvJ*-BZx{j~^g*0YPp z6fg`nlTeZ$wBwRA-4_d z)64I2$+tl#3Oihb-o_()^Z^Hks{Pe5s4+1ovBHCI)&u(%jSHPkzl>HIW9|cI>&3SJ zvoGp(nWZ{SHXOub=yMPIIvi7)ftTbt37N{SCy{9l&8z-)53E7b{_)ATZA!4=1=bzx zR+oG9W*AnE+gvE!xROukfWUs82Ge%_hO_k6_`F%=2)uz_ovnD0HMfI)yO8MNq|X?$ zr-xUqM}^2b^a{E46fb6OJ{zF6nU8Hc&wK-%Q;y^L+()0cji!N&YiWD%(nV-Hsc(i3 zK?!+A$b&m+hR?AIK%ZEiI+EXrzf1@h-A?1wd1ILZq+(&-+1AqE+^0uE^zQU{H%hM1R5VLR7-3jOq)QOYOJ?EVkkNWgNXH$VSLMw}bXCQ9 zlXi#E_fRM4>NKg-;B$HwDF12mZR0e4`RPImAyaqNdt72zOdkCEF&2*qD`^g!Wh86= z*K8e&%o$`RDSB!aQKEYIIL0wQqWhQT;L!v5Pxd+k4e}3|S7C(06J91E-ua#PGIu-` zw%z_ta^nG27ZuRdmTz|KUBorD{PKGjVVvw}%j&I76(DKAy9sC6r_RJhGdS_ZQ}-i- z77=Nn;N=={&bo|(u#6az!Y z`n#J0kYnRH+w>V#OXYriiV6!=A275$aIyfT`V8pWxA`)^zD$Uk#1IqorOlL~gyw`W z(NfM*CnVUaNzh1yl6BzzggYmuqt?<_W}ZKqhpyj0K)VJn4DzcI5u?ZKuC0Zd_aq`j zd5|yl15X7-K-#9(^ZTZ9$4cYR$(S^CIeF9dp-)4&m^x#4*W0#Ra-d}vJ2CAgd{f1} z5vmac7XB_eh<^j;>Jv}t_yLn^=M1V$BTi*vuo(RzzllvpX84MF<5|s-dqrNtDplG( zcJ1rLCd%x9I)S~=i;p((k?CkubQ^idXD#LiUL}0;Nsq+)4#~Ds zTv>~5x^??vmia$B_GxZ3b4YxD*#@OzgaciDSqXGSa1x}V^M?Qs+B>&1-<>B9T7^pa zy86w(RH7Fz55hd+M|ZUi83BE08Pm@eb9oHvQzxsdW~~@E(#`A{HByA zTA%$NJINok`GP-Zc3>sMn4M3J0g8?Zf(LR_P|IZg&D@uF3x>l!r7jQlujAYONaJ5a zACyO3@Yyf~dJj$|;Xb2|&dVgH_ZOI})?E}D`iEIrRWV{&;MlC%qKmDfcQi}&o>nnR zh-EO?)$k#nptY0q8tQNbO@KeEej(=;5{9v&r4(wJJ(mVs<9yAfj>ZtS2@QfNl>Obu zVqRLh^5Od^;o7wmF;`}$9<9sefMqvLt*+N|7L3Hquk@LXQ#RyQS(_P}g@ARQZ(OuM zgywM78J@(V#dVLp5pVZ`wRmVAn@}%mmjt~(#1xCxfNXnqS_*tn#J)nD`;~`pYrB8b zjRh=l{CTKf6TXwJT`lIcijWW;%4eEZDk*4R&9(YK;+Zg^gj>*o#vu4wK4_~^=x2U- zL;Zwomqvgm5j}=*CzRxSqPo1kW#tk7`*cC! zP5$K?gfag9Pk+*JypPXm4q|qv#bDHYg|}$qv;Q)6yPS|1gdY|W1ANqT{ngMG6eNO9 z`?89~d-~Vo1^zwRzZ7zB@BZ?q|X z2=0E6>g2?j-wEy5r0P^3x#i5x+{PYdRdQHj6;x*;V68w12X@uMSq_zOJS{%tQGz0A zVZ!IEya1%q+OOmsLFq_ro)_2nwtQG4k=#2$#DxtWl+y9GlC2?b+p?^eL946&$Pu93 z#hXAUJN+hYRqyGMq z>96mLCsX*Nw9|h2(T@4sT)dN2AY!WTnVW8&$hl<74X49<)4lKKnQq^7B0qjW*x9>qx;_fv_<7e_Nn+HPMr#Sw%|!zXFyXY2oK_0@zAi|HZquP{XFdZ2Gu* zUn>c6rFSx?qw0yh@OYisi9j%4N}fD?kJz(xYUSB5um9`c(RG#nG^GoefRxMjls~Qv zTH)T~@3M+Ot7VIt1q*?8>J;;m_bAcVeOh$`&@xB#G&Nu3w+XTH|Ti;rYu)wOXdYLUd z0H(q6li%LkgI@8%`fANM81mn_-qX)Ol6hNq>GCMVs<&RfK;8mkiFWxedIxsQ(*{q2 z{h-v@C6sx5f%>J8#RoZ=P#JV_5hD5=R{2tQx zh&glmeO|@$cVPNUJX;-|02b{^n(6g>V69ka=v(*}GWma~Gp63)%eDV2{kjdz>=Rbv zdlQKF^UU%|;TO;(xwo{vi1T#Sl(e7w47y(ng3``v5}#|MJ*SW83-cFJP7W!6rI1b6 zT)F_%D&K7W4+c;_8RR3DYXgJleUFx9PJ&YL?Z}ly#C-Q5UpsMO3XH(un6JtJ{nxSE zuZH(Ql^Pv4Uw#gBHl5b`oH7Kh*^9_gA7Wr&#*`C`cm~S)&u>3TB!a%nsKz0BKB%6* zSuU4Tpnpnzhoboa*y{W(egcO;35##3Gb#b&uH*hq5tcCM=xOraz6Ua1kM^#fdBokK+yoc_LwVj|V4Vv#ZXPo@t4W_c9-qx&C z;(XPv*d%%zjP=(o)3rK4ttdVJOxzv}SLV$u-9DImJj*xVc^3M+FAII0R|Zz@$s$SG zE*Q?g8C*y7NVBjKe}iXMU_CTHy(g#&I+xZjNnkpGu4e8qvr`!K;rCkM91e`~$HUHO zY$fKK*I55`6R_8C<3c8*K`OiR#%WP7Xq#_e7kW*D<PsR62KsP$` z=l9%bmy_5Q;C1w?6^sqg9awg?tCi??Hm8r>(>V;C5`pDqyKlg-WKnRg&3BlB zSsP(zLcH#(hqE5pK)s?uAj@h1Mx#0BY>t)?Irc?YP)RWHK9A`i+8YYyoyugr_4%OP zot5X7(+L)T)oinCo?w4*S@}h+75Y2V*m-m6VJaaRsV`R!BY^^;7a2#PLvJn_F7AQJ z!bYjhM#Oojp_aCKWPtfrSLlfc4Rm{I`SWmn7!1ECT-p2$M)S00to(@AJ%0VB%c?~% zx<6C5J#r)1-@cxuWmdrOgzgS8h7OoIZGvIxW*{FQ$+-2+7N#^Vd;7|dp!hYE3@G`7 z@$GU_k$ezX`;6au8Bah*rrEu^s}aUaJhqmb*}~wa;l`aYk6`TcY*KRWQ5dC1pU$}{ z2$N^+#o4DCp|kXMiH?el9O{hVs# zzcvMH7?bah&)2oo&isZjWCWy_EQVpFFHBdaO(2 zVkg5NexISAV{zfw9Af|9cs5>rZYgy1WuMl#OM~{N%SHF2G~`U@4Ctw)fhK?QtFzBZ z(C?4TJbDrd+Uqy+<-gw&b8qLs`!zG5U2&BY%|H1l6=hpxy_tYfA;SlcJl={3(mNxWVk-SZ0CB^+1X8hsd^mNqy z@?x?l3=^UA@PSoL?I~!o9=q)+R|ex#rS1A0BLCPe+hOdG2-@u-TK~^r25V^N^Ps$;*1eRx}LVg4CVCcVkG`+}isnJaiig6N$?GUdhzH0Vp0h|UQx zhSt2w-S<6+p8oaoUUM}8X#VN5c)<=Q`qR6kw}~7{cHKj5Ie!?8`GWjw?bpI^U&+=J zwSv%k(!j5Jrw>K}FD@pN9Kd+H=wwcw0!)Ioj5A*y0_%eJ=#PQ@P|dv*;V+j0M$6$Q zpNVuB{mxF%TC*NTD=%-oH|__0@=$sb3dMaWO8227)G16?n(%G|DX~Y%ej7vQSa_@lZtuwH*-WjZu^DInn5&m#% zdYf?g8YH*Bki1CrN*!~(H5UR2KhAYMmUExTn|prv$vMHJ97;!(lL4SM=_vb|JJY?xwz(`yB=xS+}s=#QA$(o87OU1CxI1Cn~Sl zV61R>AOC7TSob0~RIXA1HC_MJQU;OhbKh3HX!}l_bK&@t(K8^k&-k|WRTG~>pzy2x zV;F7twdBT|U>HbL<=^lqgW-sy-mSP^==NOa8V_9k&KWZf)tI z_739jllU&JG=;9pso1B-Bw?~fQ_=5G7|dMorG26cEC$tj!ZnweyIUWuzmf+^zvw)_ zfWM3nM-$iO+7v^t;{Jz4A041ymywxZ9S(!m{E#W#oiM%rU{t6<4cZG2&E4L;0!%f@ zOUguFW7Sv{nMJ&Vq@4BWZu3YO(*rUu9jt>+h~{JKCl18C3`~nLcm@65;DPp|6EJ5j zY8{qsh4I=4zv>?IL#!z(p1ggG*gtDHJx=$btNC_Ch)f13T4VAL#0Ef_`)$d~pGIIm zO26atp#(t{o+vt`4);Aee!>W>jG{;hYY!MjKI1&C zqrh~|xzfO;#C2QNcXi(arG*qdr-hNZd;WWg>iU^4ciqN|%+@ zYAuF_@|VAta$Z5v=8EHD#oN%C2^em=Y6PR)j7slp8!%0t%LkhPu;hdxF4FC}vc==mP!c zE#?=`EufW2Uq1WQ3hLoH9hHX~pfP#l^%_)g>{>g$Y(nOMRP1G$RcSjj* z1^r6KnTt-xpk1iW*?C3c|daE>1Mn zPgz30BJ+Y>(Q%Mom;`;-C-&@+4Zo~MI`mF&dBx{s2fe#-^=oYMK^N;0eq1IFR=Ro< z|5qW9Yv0cOJ!%T}Cq{ipx)s>DeCMZ^SwgSJY5S8?(?su2QvI&(23^Eh&ALkT^i8g` zZ68;Ig{SQ)AKCNyAVHLj(16s(=2l8?6y61m{$N$sEk zm=d96j^t(-)rQN>bKU{=4PoXx(f|w}^Hr(;dIW|fdI#Jt9Z1bQeWqeL? zqzXoRs5Zv_WUxzS5|fNgK%RVaA~wt!v`tQn`QqcCw_x|{BhzM}%eA`aA9)IvuWiXiT3LKN=hq-7v|LHy@7*8J# zn7FKFi87-mwuBL9Me5KNo@RT%Yk&pbO;I?d3Zz9)Qt{ zBYQhCioovKu{DGn3X>~oyXWM|f%N8X?l+}KFpoxE-^XwydYiM{mi1boZJ${9b%zY} z60O%CzaS17&rMR#6?38L+YvKNGlB81yB{p~?1!;<0C($#InaH3#d1G zSwnxrk*QPcPSA;*r)Wd;)1ct><6kyIZ{fk$Q3fkP3Qy@DyfX>&nzswirYJy<AolHv@dwY=t3s!JnWED8L+H+pg!*j`*aE`&5igx#oXQ^BV4_Ck>dl8^ zmQO;*h*UpI^pn?xuzIe zOUpm!CwGG>!Mbzz(jzdGET7b{OQ7A~|NM1TJJ=Eznl3bOp(~t~d(wOl$O*;?r!2XE-9 zfca6>GX5ZOu8JJI^(uluTWRUF_tiqsC3AupxlE{M*Ra(tGeB{j%})E)2PL_qT@Lj7 zgpZWSKNGZJ)OqhpdAAv8k4`Ur{&XD}tJw}4?8=DmvqEF$2+_Bs#}2*^-~%f`g!NUL zA9`muJ{0$-fy9~HrQ%@>>Wcob_)JSEFwUPG-8=x&^4icbxpvSB*7+`4dj<6Uo)5ow z6ZysZ<^5YL+MvH=z0;n%x**kQsctz|3$sR@je=K+-nZt(WTIRvbZ%|n%O5*UoU8EB zq()yDZnIbZI}`!a$f=Cr%^5`gb*7$Kvj%c45pkwluR||$i^w;rlTZ`govl%#0{wvt zW~;-Of%=A4r588=dPtAxo}i|nC8?Z_XC zbTE>=Rt3&K1Y^@Bv5c2OV9s1|x3>QY>ilEJMQC3^5en+tB-{e|H#ZzN+#Lj^-gHyo zTwkycH}u3#4T9FXeYGC<81!cz{8?%G9J*QQAsEYC5#g}_5~#`ymB2hatn`C z3Kc>B1ntGgE6Gs2xO&#&ZzQNVeh4_7xO31f>945g2~ z*;1kfo!4{Ezh>Qn-n}%Pjk3jr-vSgan)HIL8zirHl17}1M3ZwNYe3mPS=T<33?svb zBdcDA!p!@oaQkke?;SFHJ0?^LdfIk>KgF|PeK$9reiu*pqivZWk#`Nh$CN*wa|x`8 z_D$=Eo=HBn^d%#8DHvDjZ$}<>KwI+A$E)@gFu0d+-=Lfh-N>aA@3<_Orp0wE__z`# z0WTkk_Pm6-VZo3nH4*yj-OGQDUV!nuzdy00 zFxW5s*KzPNj1t7YoEqqYeidCg-rNj&FH@Xz^_RlX<wTH#8}X~j3cs7mKyxNeRK#~}cQ(Foex774>zk9B8X zIl;*MyqAvdJm^myEb`;*A#%d?mj*6G{~fmbQ>uR%2Abck25O#yXIw=D7S+lOx z*-K&Un0!&pP6Nh5G8Gm?pRwE#km^>h1JegLY9Fg~K)2+is{493%&m+QZF)@1r#WjE z5ub?g^^pnTZ(NuNtr2h1BJ#=gh{^=1G#GU(d@1Sm4#taQ5BrT@0?BUsx5_tqpkLnj z?n?pDr$deM8WnY**PN>-Yp{ptN0O)VugenG_35g~jS*nEA^OE13$UW*)~-yP2U3Y_ ztDuMru`kC0Z`XV#a=Py4{?dhDoT_D9J989bXMXC>iD7UcN-W zo!nPlG9>`kQOADiWe;KAQ7-E{Y2gdDy(kep4@P0d#IU7(TWp=Bixd8|Aa7 zzz*%s=nPi_YnijpPiNvB+I`)earYuvhhDvYeR?_c%+k`Vx6X&A?!cuJ3njsxS!-?d z!W;Ve{*QL;*Z?|r|I)Li%fNI?TZn&A z44MR6W_~d_q#U~EJ+<~9SqtN>iC5k~^@NV1HUF;Vt)O*@zVkaC2xZg!F4-%$VZ8Uv zhngBF*&vzervGp|5AxxgtJdt@45j*|gTrj1N40f@ExB|9OradUXHOEr ztPT5UK;$!e;}`A#!x{{G$wj1uL}*hl?)f2g4YYkP1HF#NaewjmWLKyIQn9Lz-dA=& zbJ}9@>!T@9xKdOedcX#X4dT@{-ZU_O{hfU4R1X2hpM-TCHz66@@==d92P$sGUCM5; zpdT9UHs8erb@Q9f4Mz?@);CR*;{Ow>-h!KB^>#q{9@T5zseEWX)x5dnayR4)nS4zJ z(_mol;)e=3&|0?L_OdgPPm{kKd+AldRwUHEsmwd0PlE2zbEnE3szG}m zeKTg=F3z94Z)=T&h`g&mQTfwuDQ9Gz_)ha1OCe#i?dkn^8WcjiWVQy>L3h*V%WM3c zq2+gM?!{;Zan5QCZl%OP*CISajL#S1$zIlr8YGDxu)=Bb)^A9W7MQA@v%%~Vxh>xp zHbcSQsbAq*J2WM^?v;&=Q1|@vw82UihGw<*8Ru7m(cW1a`uh^JE?t@v$yWrufTd31 zQ42Bq)F*8@B|aF5=YKIaCGyPZOp{cAB4kd=CD*|j8l<#OX2lv1|(;AE*54`P-ocJ0I=7{|`SMIy^xHt7d|AOWkGvd4A}(_eTHlSE!diiC|TWrtCk@yER`$;-AOUtoM5U zuis0s=l%0|sX0^kANKXvsDWU~D-XGUe&6cBm-5f+mTj5N`serAv*vvI=kZ?umDc~f z@4lra2maCT8S(l0&+pr>oVoeW<2f5-@Xzx~WE|G~!xpcJ7~=i@&s*#BN&KH1-~V$f P{ptS?1VrX*<9Pu9fuYS} literal 0 HcmV?d00001 diff --git a/tests/testthat/test-Log_Likelihood.R b/tests/testthat/test-Log_Likelihood.R index ff6bb26..061485d 100644 --- a/tests/testthat/test-Log_Likelihood.R +++ b/tests/testthat/test-Log_Likelihood.R @@ -200,9 +200,8 @@ test_that("mvnegloglik", { }) test_that("mvnegloglik_ST", { - source("sim_multivariate_big_st.R") + load("sim_multivariate_big_st.RData") P <- 3 - Y <- cbind(runif(10), runif(10), runif(10)) cor.matrix <- cor(Y) cov_mat <- c(cor.matrix[upper.tri(cor.matrix)]) logparams <- create.initial.values.flex( diff --git a/tests/testthat/test-PrestoGP_Univariate.R b/tests/testthat/test-PrestoGP_Univariate.R index 344da24..ef5ccf3 100644 --- a/tests/testthat/test-PrestoGP_Univariate.R +++ b/tests/testthat/test-PrestoGP_Univariate.R @@ -267,7 +267,7 @@ test_that("Simulated dataset spatial", { expect_equal(as.numeric(beta.out[[1]]), beta.out3[-1], tolerance = 0.07) expect_equal(params.out[1], params.out3[1], tolerance = 1.1) expect_equal(params.out[2] - params.out3[2], 0, tolerance = 0.3) - expect_equal(params.out[3], params.out3[3], tolerance = 0.3) + expect_equal(params.out[3] - params.out3[3], 0, tolerance = 0.3) expect_equal(params.out[4], params.out3[4], tolerance = 0.4) # Missing data with lod