From cdf6ec95eadc756ef9f4f85f88778b9cc74078e7 Mon Sep 17 00:00:00 2001 From: Maciej Nasinski Date: Tue, 28 Jan 2025 20:52:47 +0000 Subject: [PATCH] improve VIF and actions --- .github/workflows/pkgdown.yaml | 2 +- .github/workflows/test-coverage.yaml | 4 ++-- NAMESPACE | 2 -- NEWS.md | 1 + R/VIF.R | 7 ++++--- Rplots.pdf | Bin 0 -> 9284 bytes man/VIF.Rd | 2 +- src/R_funs.cpp | 2 +- tests/testthat/Rplots.pdf | Bin 3830 -> 0 bytes vignettes/miceFast-intro.Rmd | 28 +++++++++++++-------------- 10 files changed, 23 insertions(+), 25 deletions(-) create mode 100644 Rplots.pdf diff --git a/.github/workflows/pkgdown.yaml b/.github/workflows/pkgdown.yaml index 2ebed96..3b79e3e 100644 --- a/.github/workflows/pkgdown.yaml +++ b/.github/workflows/pkgdown.yaml @@ -1,6 +1,6 @@ on: push: - branches: master + branches: main name: pkgdown diff --git a/.github/workflows/test-coverage.yaml b/.github/workflows/test-coverage.yaml index e3ee635..3853406 100644 --- a/.github/workflows/test-coverage.yaml +++ b/.github/workflows/test-coverage.yaml @@ -1,10 +1,10 @@ on: push: branches: - - master + - main pull_request: branches: - - master + - main name: test-coverage diff --git a/NAMESPACE b/NAMESPACE index 68e6896..86b7618 100644 --- a/NAMESPACE +++ b/NAMESPACE @@ -10,8 +10,6 @@ importFrom(data.table, as.data.table) importFrom(data.table, setDT) useDynLib(miceFast,.registration=TRUE) -export(corrData) -export(miceFast) export(neibo) diff --git a/NEWS.md b/NEWS.md index 4e72035..3c27c1c 100644 --- a/NEWS.md +++ b/NEWS.md @@ -2,6 +2,7 @@ * fixed CRAN Notes. * style the cpp code. +* VIF() should be more stable. # miceFast 0.8.2 diff --git a/R/VIF.R b/R/VIF.R index 58df58b..1783974 100644 --- a/R/VIF.R +++ b/R/VIF.R @@ -88,10 +88,10 @@ VIF <- function(x, posit_y, posit_x, correct = FALSE) { VIF.data.frame <- function(x, posit_y, posit_x, correct = FALSE) { x_small <- x[, c(posit_y, posit_x)] x_small[[1]] <- as.numeric(x_small[[1]]) - xx <- model.matrix.lm(~.-1, x_small, na.action = "na.pass") + xx <- model.matrix.lm(~., x_small, na.action = "na.pass") aa <- attributes(xx)$assign - ll <- 2:ncol(xx) - VIF_(xx, 1, ll, aa[ll], correct) + ll <- 3:ncol(xx) + VIF_(xx, 2, ll, aa[ll], correct) } #' @describeIn VIF @@ -109,6 +109,7 @@ VIF.data.table <- function(x, posit_y, posit_x, correct = FALSE) { VIF.matrix <- function(x, posit_y, posit_x, correct = FALSE) { x_small <- x[, c(posit_y, posit_x)] + stopifnot(is.numeric(x_small)) ncol_x <- ncol(x_small) VIF_(x_small, 1, 2:ncol_x, 2:ncol_x, correct) } diff --git a/Rplots.pdf b/Rplots.pdf new file mode 100644 index 0000000000000000000000000000000000000000..0026ece541b7c66189438f1107b2e8a7353c707f GIT binary patch literal 9284 zcma)?2{_c>_xP1i3JFQJX%UJUGqxd1*(3YTkTJ$InlWO=R<^Q~En6rdDT?6BuS+Z}T?EGg`pZb0J{(j&8Jdb(ieeOB;o_pTcb6@A&D`cRiAr1ygF$)C@2aE(z z0@`4%%wT{70B7&Stf0UQQgubb@o1d08XS)Vh^Wa!Bp}ifUl5oHGH;vD|h zHWM@+L$xvjXu85(uA&ic0DS`9g@C7;8Mxvc2neJrtyKyr36!A9cE`feM54?}O3X-S z2bvd0FWvi569q80FVJ3g>(Z*QWcBBw&gf_J&$-F{)Ag^p{hzLaCH|YQyBhy? zR#F1;-`%ZvsC$o(#2!4qbj~x5)p;phIS+EAdT$J1fZ?XZq3Z6JVh5}a#P2&@o2Q>Y zK&qYU>FEfsXojNpa#P7tB|EOlQkH1C}g^ zt7aS88k*K0^xHVucT@i92@Uj=l5)SrYEm}GP&RJcwnArhWr&=ra%e~X>4MSrfPhgP z;}9-{T&^6*Z8M?3f@VQqu9gorYL!5to7cDKuy3%?|CAaLi;T;_~Yb-OFW|`E-)C^yz2i>8D@foP^u0iCLH8uy>8) zxevto+~lvXci|zQplrCrD~0KL_9nK!mQD-~*(Wxt*eYT_R;!g(Z*hPlaFzE^i6 z$4}=uok|JUJ-AHSSg6!jm~@!I^b9gseEakyA4xgh^~d+%yY2#%F7_)KS27ZmUp}<# z{qoyLc5#!c_PaDO{NR*ptmx|39*=qDLElaMjit4wkd>jSCco5%VrGhS>JDe_6#+^k z$4-w%UTX=(wQFk!yCi&GIlt~;8j+})2|Mcn{4}D^-6Fo-o2?)|)Fq!eUBI}rW)Q%4 z4jD};kvKoDIedwjBbws$c&{0AOHQtz*`ABUI8#zr&K%66|HP(cMpPKA=H-pgVj5X& zj1taHZ80I(Ylc0WMh(XFrXM*Ykx^-qd5zWWS??|w&&C?otEO~$yLonW==Hu*2l>Sd z$2Do!g(jyH_Ep!Bo--FINE>w2rFLrX89U(Vn$28Yhl*Q_8ZVw3&0DmP@LsznV4rW4 z9TUV~<+3QAwO5=i*o?3^eMZ4LMWOy>&e=o4sfN$tDP3i;srFfR=iqE}+-)C7jEm%G;XqF(;pRx>sFf ztV5*1_O+s8dCNTA>t+qYF+x)$@zIYHGZ3a@feNh$O!>=OGSyGq^U|Bvt3+_~IULrG z)|;r!KM48!_`-kGyJjtYZ(dC6IwqYyT*dq;cOaB_@eCgy-D2ZaO#1C`6@kg1tVrxc z+moyhi;p_{XFQ1ym~N_o5!~yGnA@D*&tpoPoG|Gj;VNe*gX*b@-A}UW7GIBUN${07 z(r>#QbKnRW>aKymMZAmdefc4GY>!`}j!Zrg!M#C*62M*41NUn38x%gXO-Jj_Q1ZjP zKZA2>v1SPKC4c2cp(Kmm{4gKuM`elAgjCWvr)AT~XtxFvkMy*Uns?Kp56%^HuF%(* z^!(=eIWH$D@nc%Jy}(3GRx);`?MW0_A(P8k-5kZs@4%@YZ9DP2Fajb*T?zUuWEuH5 z!$$fFi{iovA1ltq#POoJQ4jYA^oLYLTMzIP+zo@(O&@=$nwYsoe~8-S`>bReYLB}i zNWMwV9UJ1CQbusByUsj`+g&MUC~gJ2!^Rhd&E{t8$wy6fYO~>EeFGsfTMHxd2(w{heEwzHj|8{U$PA?2=&xH%O7g4AxiCHJLYBz zDtTmWYWpY8lp6^p_!;{IO1$>tBZUl_khtj{eflnBRsI5ZXZFUaJy!>-@zI?niZ{T{ z6i9e4M+h*~ad9qdGFUc)7ntZA_HYs#BvH}tlvUYlvt`VeTckMLM;Z{O@GrtBRe_%x z0zC$HK3mVR&Ex!hI|lH6UpG}g9OueKmR}hk`v#KbzW*=?dkr~UVU#zxxxRkODBD*E zcJ4u>j*Wp3+X0Gb5EK2*3)OXzXZD`o8=A&@?{RZlPw5(Sw|zM0`)$Z3NBUCE(1!9B zX2_#oUIFltgJUxtm37cQdG@h*pvzSI)Qe`d_2RxA)~e>?bN8kL3ZoGkTPcz;c>-o5 z_p}yHYTXkys`q4nbsivRzI-a^lU&UyN;t}GZFL`b-7W1U+xJ-^J#BbS;QJ`{>C#YJRi+5eJS68;fT+a4xCnGc_S#EU8w+#RI17G7;a6L_vlhh zXCg_-`gnxY$E@dSHI1+gwO_XJ?L2300xbHhs5yJsaoWZ7kY%D-ShL|E1+1N|{$(<4 z`*Wb7)k^T<;j6$M(F?7G*B^YsR?H;OEk-Ho1HWz#?z!x?HP$Lg+2Oh{Fg58%5fi3- zO$~3Hw!QUm#0r;Ntke#_UAncInPMxr&BUUX!|@lFfqx_L{~wsqAk_%)e?TPoH^wr9 zeo-VuPU63VShcyNz!+F82c_e>RtK5e6{l2&Oy-+@Q4$-=vr`ie!CyCzGH1t ze)EVDpY0|-I#Yb}`VG(P8%z42>WKgn3 zY1-QZ{4_|*3khN1f7bx?NeM0rmdb$htBRr%3yb(du)amhA=midO}jC1hrC#>su`)u zPuhhvJ}RVDz4P@edG((lO~>K~q-FYNWgxqcQen(oC;`>;=xn#GX-aTU_N`gd8tkz` zV}@Q*w?ngL*Isw*$8pIYm?kwA=j35xOU%-fa=sgz_-{-b2BI4by;w{Jtfdu350~)1Ts@jfEK*ZpIB_I%gX;8(`rN zt(xcD8>e5@bJ-A2=sipl?=VwqPa5=_Iy;#2JgMc3`S`NWes1!L4YX|DDVC{)&^g%9 zV0`F~>Rq{2{h9=u5TqY}JC)RMx7Ts7W%>8U2@80qTsg0}4SY?bq72EQ6n+x8gK#64 zv@<+z+n(stJ?i^o@{Eh=zE(^8CpiAlI|P{3Vp%75_t@o_s}7PIdk9J@{_5 zaNH34HIDT8>yw`TTi(dd^{@S(9M1^BP2*8jW(iL(*XL2T7l)RIQ*g?@`?+yD zE3jNEPxq}&lFQnZXKuh0(g+-(1}tj7PkX;Qob#yYINTC$jmTwJSfo zE0`$}W@XhIp(4|uQ`D}Ti^)Dz9dXO>s{MeTEB}WCnfpElaveT#f*E4&Xk>L<=5k!G z1bhu{#YbOAvtyC0lLsOibbF(`u@=rxkM8a+1NuE8l(PO?`UoB(akx#y!P3g_)ux3Q zYYq+THuCp#5s7og{ZLYne#5<>>uo`O$R;I4=#Bh%gNflgh z|M_5VlqPlOLA=Q8+=pizYm0o}&vd>Z4U7BsR8dh+( zn7_Q6RP!RN)ON*JEYgMwZWi}2_I{Bz?ezE1>dv%nnQ0Ra?0Nqw(q^W2xV;9JxA7XK z7u}(t7m3A385Ry`KJ5}4WECb8A=9!xVfCb`l+0kN!@r3`oY||<%?fOaClEhg)JWr; z4El*2EsO0gnDpNe^9yNxQ3-Ny0dXgjTMSX$NDC7ND}nH4AoML@EI^OSFf~fzSx#aU zbbp=4jFUif+P0)j8i)+kMqfAWaYYl{tLPcQTzLC0%|LgIVTXg4?V^|dv4%;^HGze@ z)u+p$jm!8O zP9#-~d#2!1GPd()z`QWLg?5Ntrq0Ua>m}725(8bc@x%xQUyXwULszL@%eaaBCaT@U z#Xf6dibmn~qAg+k1R@PV0sj%^2*Pi6!7YAoMYkb}mlg-eJ=#*R<6OD`GMmL9Qo~hP zmR@&KQI$AWCCOe>3;h12TSHVl_pEr=%MUW%Z$&b0*}m-S`4GtybY^V!^nA$2UGU;% z>~^T0m??l;@V$F*jYi=jeq!7h>eOsHWk-IbFZp

3fUhsQ}&1E2soo5Fpe(bH@y~ zEr8t)#plak;y&ndyLyz><({7}ymvTQHhRQ;XJh93UVr1>O0{dU-xLdcaHr{JZzKX+ zkCQ66M~DgVK6kz{k%x8V*%`b8iOwZ?yQwnFzgJYL&Bxc!W!I3%^gXDc><M!^5W7O&6nH>5D2zt>XAwg+;ll3kUO*uK0PYQ^S*4$vijAd>pWJnzlI2TvOgtrh$3I_ zU(HI7VPYl{yB@@{-TBK+)?nTmbw9lcnrP~?;JIfF5HJ9kS#VHpSK z9rC2JLFESJo1^uV9$*r^-hb8d+n1{A%c9N|90_ynUKIkz`!}@Q+C1Ezi@lf+ZCCIG zC^X%D1v4ydxwXPqteTipZ3=dPTZzkS`Y4k{a*p2bSKXRqPA{4}^kvc)8nt}dCd1uY zaHV`tj{_xEM`6dDz76)Kbh&ZKPcn+MZXJ+fsh6^@%=4hDVxKg*cmQ^-zkO(r+tylA z3Py07&5E&h-(N6BtMC0NY5mD!ZUB%D+QAKA4WQNgi~zJ!nh{_Z+>8J=)FL5{uv_c{ z=^`D_a21>8=o894<2OO*16=ke|wT-iF%HdUXJE?J{mM6h<28rS0i>o~_ z9_~7*qWL2sl73wxQA)R2Lc95-O=tlOfP_08`OkZf{GkU#~V2Dj#6I?4!5>RfhTA53RY>NQc&#Qb_HtDowXn-rtU z{`oI(&!t#s3?8g%E9>m&ouT%=wt3iP?a)|i$f_XIHbLMh`{wD&JD@mqwT6t#fXng4 z*0mgW7#^quUsz3jxXM<&{j{ev{_dw$_Lg{yz0Xb@zu+rCCvuh@aNvf+xt4jKgC8{M z*P6mUK1fX7pxoN_F)z8iU3>JjUnD|P4N`K%uzBX;rqy040G;{;`l$=N6Z9N!_YvQ+ zm(a0&rr&!}fkCy0xr{>a0M%N)f*^}1-cfo6Rdxfuxcx$^Dq4HFE=)f* zdd`;>JeB4+VLc&&KVC-1qDp#fc=gCSmp9WXi31r`RzO1vZIbX-HX2PDN5j} zU}>mUbY`tJC!@ML*x+*9!FoZipy4O+X=Z6M=~;$M2FcHr+D+RuOwk()6); zf(|bx&yySm)G&XX)WSbeIKjRQUFTxForVYpM;`OL%y)*XaBm?`A#YoxToBKPaCk<( z#Y_24zUfew8rN6$7%-9f#^D=BZZMpW^tpX%81F)aUV$m--8uK(^s|+Vu#341M+wjC zBLj!$!`W(<78<@HS2b?3`R-SaOg$mO5q=;yP)p@3Ts$qNMFXl+qT_k-^pS6O^xnHX z@h=IsxTK(_pwA~%DP*YScv?m@;PeUXiCN*~ll@cxjIGmfU~^TlygZ=dd?xk8pzs^R z$7e0h?h{$yG&|QWfg^&?uF%kA?IP_I?TGs8T%!DGydUW&CYK@W@O8rac~fyyebdURGaWJ= z=d$uT#-4jU|JtG5{;5Nwec**;yH9&tJL!eTi@~-F=tOkjwNa;;kI?L|cvfs=1K-@E z7`c$~zRp1DCaF_W=2z$1_*x!ywstPI)wgD~%BO!bNy?o~OU|lL-3a{tl+~21tNFN4 zb^K{FBAB^4HoEz_!!kGnFHP7Y2$g*+yIm&f{lVMNJIH(1du92+y4cp~E&lb2bz??( z#uJPTj5&-XmaoBL)tuEf!CS$35#bT`;WEbvJc68)SX{Dm|5YDyusz3B#S>=e^`|+~voX zJA0}uXKr?{adT(rJk$wmNcZ}%w7PVBJfd`})D<>u^UV1WNsUzDY~*3%5qabG20fb( zTZ(2ZIWY4DDr@{A1!dc5doqtTuM30vuKisU_r8}>{{BWw$S}tj<5{;sbBb&asl%h5 zGJp#j6!k0ei`tRd;oW1kr_X=P|HTFJ+xr1N7oR?6xY%)#?N-w*^(w^^E?RcA-=oTx zSjG>Idj#ET@o$+;TX`2=9UtT)7S#h%A20K-V?}(k@S(I;Y{vzvA$9lg({^q=+OGf z5p}OHpjfBHqsrCQBQ-*ggfx$QybHc7es^9j?Ga=1a;yKdVx7zq%p7ozBZUdago~N4 z2aW5WEZfAgKT4&w?{uaT80!94H1^GT!|T7sU-$kG+m z@f70L$TAAtvKh|%{O_!S#M@c4-WN7^pV>TkkOF+?$@+lR(dJ@MbTun?~w0gz#V#@rpTry z@rUBb9GC2;ghjh|RSkoler^39Fxj|FQi8zVkkcUoeo3pd->NT$+_?8Dvd9OuYJ+7e zPJhZV1sjm-xZf%?8#yhYblkIuM1Gbza-(5>MX*iL%nkO%v**COyH7Y$^~O>kDcvVb zO-_!sAu4pd=O5iwbjty~#w}o4dX~m-Qj#ZEW|X}KzHR0oZ03%%b+DN!8+Cc%_SX0D z2-Ak3JR?8j`4GO4GSGROIi3{Ro5P_S_uI@X@>0r=y3DrLZblDH8+PFav`S7?)RnJI zn62lGliZvseU@(AKBUo<){4-Lhbz3kQd{{pQvgM!DMk4OyQ%E+B%-`5i- z*5EaP&1y@hdNs*se1t*DwyU0-mpWq~WY#F3+ktH~j`n!G@!W(+aGjUlhAnTlalBP(tgd-5tbAM)#1DZON8=9I&0G)Ao2c#na(UEAM0YU=nuWfBltx&GfWxxZzS(7$| -#include "miceFast.h" #include +#include "miceFast.h" arma::uvec get_index_full_R(arma::mat &x, int posit_y, arma::uvec posit_x) { arma::colvec y_col = x.col(posit_y); diff --git a/tests/testthat/Rplots.pdf b/tests/testthat/Rplots.pdf index 7f8ad6742d9eddf6ef05edc2f9afa83dde26e0bb..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391 100644 GIT binary patch literal 0 HcmV?d00001 literal 3830 zcmZ`+c|4SD+ZG|yVyAghZbGAQo55J3EMv)%H9MKcjd>a~)66|tGWIRo7nM@h6fLh; zLbfQ`vy3gJ6sad7OHqXCyT@C-&#&kE?%&KG*Kr-^c^>C^{c|3Qb_6q3tePfDF={;O zbrdIRz>|T(0_p(m*qlYN6eP+4FlP|y0b~*ruw}t?7L0`1F=)Ol62#!%YO3MXun4wZi^heNQ3np7AgV7{ z1@4y`0CWnWLjY((gozZIKLFYh{UIixfe;*k{|0FPBY?}}vKbH+{y_&#XjB+EG69|S zaL@r_(pU@ecZVFA#LOo0*opcw@j4FMw2d?BuCCX7rh5EY!Q z;WY1{?$qvpd4t+WRO68e=W|0k`f>fEH?H*mr2g64e_8^q$i7U#1K_se0C*uS(^!af zq^b?%OC}oAf&mY8myo0 z(_q-#5#~OvzCJexQinogIMRuvbp&KOm>UCH>mra_b2|b zGDMQj^&Ss!;I-IVZ;`IH7b_%U`QQwXJf7f0C)lSKZNg#3A`4B@%jA>dHrSUyji?qPg5va;>*(YwcCI%t9u6yrZn<=&pTJ+CXeN@yO=WZ zae3GiE#0AG;F7IaW!rJ_uBJ_ox@C{NTQtu3Ec2s^N7&~qNLY6P7psr%WcG|h^KVA> zyiU~>Ye-4=K}%0;TCTx0D3M~|L-B^Kv^Qnwu)}yBOvZ@XTc)EgMU|9+!9}aS4&fQI z=r}sx54v`7#5;vsdb-@z#1x7X(Mq0<@u_FUlV5Q-HWSG$UFzGkQ+JPJpJt)mWL`0N z{5WfF5+Q-;JM#OV>gejI2dV}(UD%Eb(wL@>pjuR? z&&H7r4(i)-BbDb-U4|mDaqS|K4hb+|(k^o+&h^ZWH+=jiqINRpHYu7ITkuL6PUkyZ zm$?}~Rp2-2F^PfysO8yeQk8EXAhseID!5x+sHoBXwkUcM$BUzo^jtZmbUF;Q34!Kh-wQ&X&+^|{TDWF_OquVfZD7vKwT+6&s{Tt6^)WDxTw>XWFA+3_pgeLa&r zu!nGx;2GZROb3G zy&;+Dq^&KC+s7!_yQotlr^HV2?@wdL>>h{dcXf_>>K{(r*L-BmoxX?eLKm+_-xGT& zGM6IUvM}HI0a`W-6As^WFfC6GBc3Ak?^p}ty+qZ5j9xPxt7@xYBL%TfiHDo%S0bz9 zT@UFK^lfDn8x-v={1ouYQ3`T_ax;5!BHH^jU zjtraWI_a8(lrxe)^8aNnq0*8#Zt}`L3A2jv$B>+2Z`wkRHG`}{vnjVJovkJsMf|-E za}%hX3FB$am-?C9g0e7A^OA+E4i_8(C_Vtb@vO6f8F@xGg- zkKbMox&Gm?<>0f&W`i$pXbiFk2L`Kd1l<@NFeG0j#~y$4=ledLTOZ(UfoYvGv)LKC z2@@kjv0B}lyER<`W(Q_{r(=zh??^Zq}Qu#(uuK?MR!Ad49Rc z#MK^BJgVtzde3#=MQjnQ#rn!pto>9QQ>zi$8fqUJ7djKVv?#Qq@>SvM&XxKVM*%$n zIRSow5`n6%AL3P-B$~Y9zs8rQrlcNA!Ea}w4RFpH+thossag%kjqxY6-TeA}7~K-2 zM>rde7lEO}Tkz%DR@htE9k_($($Ed!-Ib87B)Q@30mdIGfc(fzsUjRbx~) ztJbOHxuIQKJOaz7s#XF5X*tvv0qnN;W8wjh0RyCkTP0blPwuvUe&+enbF}8x#7PdD z>Q8N}dQJ)&$o*6H*!{+aC*NXX3S-&}x{3XL53qe@Lmvc^TT!GyznmcB{UN#~C-vSg z4h=WDy$^f(MOwPZ>XKD*XJJU|!t%l&6R9;*H4M*buWQsTRfMW~szZ=hP}-@ORX$<1 zaIX2;w%Fnu{x>I#IR4&4-twi}N}o{tKU;oQrZqq3)HR>#O&Ayd%W;M|>cY_;u6i8w zi1U&bH>w;_5pm`_{`-y%?i)rT-$vdrY2tCxYUldn0>i3YH$GG-bxIWz&R&Hy>Nb>ZvF&`tEmjRW;ZM^i?t|bE1HY-oyMM4Jw2yekS)-&*dKYV(yF+cGOIQ#o-0Tbq-D%K z0w-)0YLw3JF|a&o^D{iD+ZDW}c3#Vf@_*BTCQwt1fSS=m6T-hl5y`o$AY z7NSK$tGG&y1zbXTS9wzT)z0h-GoE-wl$aJ|GMs#;u0?a0GZ@)SxZ{nx^g7e=?1GZE z-jYQwBUhc~Q#1V0rgP6n_QHMdDECIMK978l+jo87YZH%Fj`ZHVb}Q6nx>U{GprMu4 z9aR}yNj@HdsZdh^!)nz0ylm5*`z~E6`xHLykt~^H+MEuaux6ApX>}MGj(Cw}Q3dgz2`PI<$^O35Z7#B4TDEoZZoqDDq z$-E=&!u|{I@7WPu?5mP*Tivm8M|%!?L3WjuQ*4&yCyP`Un^PrzGeysy)P>ZQzW?xa zz`6K{i~AC_)Mu_c^!`eZKw*3m`yyL+br;@D`*TG1N^Q-d!W$Lun!;a&zl%!bV|S-@ zcdK4fg-YnRuChM*JZ$V74UQOi(Mq{RE3Ue*i#I_bVOPZE<(W@SMhU0>=XP2J+ke?B zP_VM_s`!-WOWnuk`xR%>rgt6qBe5FiXIDw0TPJ@WMnGC+FSL`-34bKdqH*?vYOO_Vs#S`-Xml z`5-+1wcx6(p1@9l{RuJ&wcviQS#+*;*m%J`b>ZOTKcPL`zN%Q7S(d~ z9d%zOomWaGs+d&Hh#OOyUG*lnzdmX8($bD_&97x%Q-HyNDFeNEpT?N|q>-{c2chqa z;rAxxz7Shtdk71=pPRQSvG;%ny_bX64h@~nDsDNr@4M$}*PG#>`@w5Cb;YzCg@8;JPwmd20(uX@g%b25=kuXs!Jj>NUT6V z3KR^0Fqz_uK!HROgGL3wV+@E(L?VeK65>Hnpf4FI#UwKUP>o81eIY*pWUPxMA?rSc zLi`&Zz@qvS8LU7Gkp%-F&7Vevjw98G>zSF{&E4y7i+V5v@8uMfa$X=!Ky ze!wpbu`}zp0QDD!#UtK?`~3~m(D^5(iAPrb-*{Sh#JBv0;gL+gV>m40?SA9obpL^A z;B*nc@skHLh-3=HK+GusI+8;nq&EOOMD?Q~?ucuaKnEHP2Dqki-4illA_HEZf`*nB M4yCAQYHNo2UvJTC-~a#s diff --git a/vignettes/miceFast-intro.Rmd b/vignettes/miceFast-intro.Rmd index 81888d1..db23a6d 100644 --- a/vignettes/miceFast-intro.Rmd +++ b/vignettes/miceFast-intro.Rmd @@ -86,23 +86,21 @@ upset_NA(air_miss, 6) ### Intro: dplyr -```{r,echo=TRUE} +```{r, echo=TRUE} # VIF - values bigger than 10 (around) suggest that there might be a collinearity problem. # VIF is high for Solar.R and x_character which is obvious - x_character is a factor version of numeric Solar.R -air_miss %>% - do(vifs = VIF(., - posit_y = "Ozone", - posit_x = c( - "Solar.R", - "Wind", - "Temp", - "x_character", - "Day", - "weights", - "groups" - ) - )) %>% - unlist() +VIF(air_miss, + posit_y = "Ozone", + posit_x = c( + "Solar.R", + "Wind", + "Temp", + "x_character", + "Day", + "weights", + "groups" + ) +) # IMPUTATIONS air_miss <- air_miss %>%