From 19638e5a290b10c4c591ab6587f35daa28d62436 Mon Sep 17 00:00:00 2001
From: Shu Fai Cheung <shufai.cheung@gmail.com>
Date: Thu, 5 Oct 2023 22:37:18 +0800
Subject: [PATCH 1/6] Fix influence plots for reruns with missing stats

Tests passed.
---
 R/influence_plot.R                        |  15 ++++++++
 tests/testthat/Rplots.pdf                 | Bin 0 -> 14975 bytes
 tests/testthat/test-nonconvergence_plot.R |  44 ++++++++++++++++++++++
 3 files changed, 59 insertions(+)
 create mode 100644 tests/testthat/Rplots.pdf
 create mode 100644 tests/testthat/test-nonconvergence_plot.R

diff --git a/R/influence_plot.R b/R/influence_plot.R
index 1f2b3dd..60fd415 100644
--- a/R/influence_plot.R
+++ b/R/influence_plot.R
@@ -261,6 +261,10 @@ gcd_plot <- function(influence_out,
     } else {
       gcd_label <- "Generalized Cook's Distance"
     }
+  dat <- dat[!is.na(dat$gcd), ]
+  if (nrow(dat) == 0) {
+      stop("All cases have gCD missing.")
+    }
   p <- ggplot2::ggplot(dat, ggplot2::aes(.data$row_id, .data$gcd))
   p <- p + do.call(ggplot2::geom_point, point_aes)
   p <- p + ggplot2::labs(title = gcd_label)
@@ -432,6 +436,11 @@ gcd_gof_plot <- function(influence_out,
       gcd_label <- "Generalized Cook's Distance"
       change_label <- "Change in Fit Measure"
     }
+  dat <- dat[!is.na(dat$gcd) &
+             !is.na(dat$fm), ]
+  if (nrow(dat) == 0) {
+      stop("No cases have non-missing values.")
+    }
 
   point_aes <- utils::modifyList(list(),
                                  point_aes)
@@ -564,6 +573,12 @@ gcd_gof_md_plot <- function(influence_out,
       gcd_label_short <- "gCD"
       change_label <- "Change in Fit Measure"
     }
+  dat <- dat[!is.na(dat$gcd) &
+             !is.na(dat$fm) &
+             !is.na(dat$md), ]
+  if (nrow(dat) == 0) {
+      stop("No cases have valid values.")
+    }
 
   point_aes <- utils::modifyList(list(shape = 21,
                                       alpha = .50,
diff --git a/tests/testthat/Rplots.pdf b/tests/testthat/Rplots.pdf
new file mode 100644
index 0000000000000000000000000000000000000000..4ee28b0d21b3659c618ae7b5d66354e3418be5a2
GIT binary patch
literal 14975
zcma)j1yo$g)-Dp<g1ZLy#v!<c;O_3OO>lP!?oMzB!QDb|cMI-r!TmLvnS1BndGCLJ
zuf;muyZ5dwRfqm|ov$e5Ma3BajBM}}ej9$9eq(-9+K%u55@r%RLrZvGUU()EM-u~}
zg`KUa0nmhmT9k){nS~X=%*@Wh!ph0aK?BbuV`u#T{7%^dXbrMbB>CWIU~g_=<U}Ip
z477I!g3RO{?Tnp`OdMaoWn*Mz1b}3}ec8Q=XyEzz;Z1CfU!!38S2QF{%5L^1BupX(
zKm%(#GZH3w12YpR5>}9cBFVoEIR9brs`skxXkrWeQ=Lh~&K3yz=tRQx+nh<!#L3Rt
z(Fimi+kcJss_<7iq)lwifaWADtSlTPOyU;SKu`uu;?|(zh?*GL8JoPu?F0nX%m&^y
zV@5+Jnz9bHy`$qt!!Xh0ZH<kd3#0@c6x^qLg10a%i2ftqu)(l29sDX+%1vcvAxDVR
zFvt%lyx%H{mQ_qzGJg!k{5U+D+8&~Me%LxG8_Nz!Z;|?G7}@!<vfQe49-)R3s@Hbw
z|H0$->em$~Y49l{O>4eS1rtr2@a==xb9^>d*~97XqLiyGqA-mXVQ#)nq5MsK51u9A
z;@5oaXVeu!j7<(F<mY6XXQl_K+h?Qbw+*UMA3|@4z(;r6KNGF!9uQ~^!LMrTmT-Np
zA93i%B!+KlVCWGMI?t8AZr%n`Q}bh3(rFEZ_#3L?c%>_4zm$LB+2BYEoyyR%mHIYL
zgW*1lmFzEEB=7X|eqAUR#bFTQt?gYm_PC6Oj246Kpzm@5<`;$K1n?FbDW}EsbQXP1
z69;QA7bR%=K&ypo&tbS6k1tEK9fnHsdsH^74a!j;Aur=+A5?n*7QV8ilF;{hvO6dr
zV)OVihv?5DPu{=CX6*zim<8>Ie?)mMFzZ_7{5n_fm1e!F$%+5S@GOeL1<<R@O*`C?
zd&c7H_D3Yp-@IMc&Wbm!;Kj~p6v3hUre)&?|3xlJNpeq7_cs0A__P9`%-Z^0I?^+i
z&lec!RIHKUINEf9JFUyN2rTx#pEqH8YwW_Z41N+op2kW=#)b2YWX*HVC=cFqevkmz
zS{w+5D-u^Kf+frMfH_o{&V@dHlRf9OBFh;Rl#b;6I4-LUpO|j2WQn<9MgFmByV1lF
z&0~eOPPlY;5(%y-W5Ir)V1`~xYy<x%4k`3juXUJK<`C!GT2!s=K`9%#)+>!-z!A%?
zLtoyW;D)~xyVEu5M2deBz_+_J@nSarGy2aXux>hK6JPZagHLaYKbEP16-!FT!{gUD
zX(vAGvX!{M@ggkNqJp%%knmuUs{KTI>ar~=0~N+J`7M>{la5jc9Tj0t8r_iW&ssN3
zu$j_*rpYgNOrN64CW)u#c*gX_q2>kG{Xf>ucVqP%Y>3%Qpu)xo2B@Z)PkxMN8J6{P
z#Z)28ahXaY-2SMi3-z7OuoT03dei~xNtz)$JVR3>Rd^NpUiH1-fc4&TA{GmSM*-cM
zAUB0A3+2d?*EGmJh!z@e&Tm0w)Zjw1T%DdWe-QdpR@r0`yIrzrY$8AO-ujTNFMRcR
zmzkqV^;SC=obncXTCJ@O{d@2;cojs%6~|&iJ{rDQcPQV|ZZ?{|lFZH0%ZbPRFzDtO
zB;fN%f|-bEwzsn%RE!`SO$`H06ZH_XV5=}D<n+u#^*a|QP}VBl$~5!EuXb+nl(QQ9
zhhVzc<+KSQ#)3_07EK>D{N8-2qwb@~%vd&sk)Mm0^!w31!gxWSHrJ|~(EFFEg-)?L
zkwP<vJa;EcIb}x3=OJo2iw-b?;&C43bYu+8T&|X$QQy~eSw3*Ikh>M?Lcw}{E>WyR
zZ9!KXPd+9%0Kj?Sr&O1JQOVCmIC_KDzy!7~Y9${nrOCG3^Ub_8IlmNVPm9sCk79n9
z4osP%P%;teh*zyZPO5X*w#Eoj`L;VODTcI0&wSF{4_5>f^tM?&g-zVDB}`>KFGmb_
z^}>Mc&sKO}+CS$=<$kH|y?_;MP6cG<RLvQUR5<Ie39$q)@^MB<%EXspWI@Cno?4dv
z%!Q967qnodp(T^0iPW!4qo@FEfSuh9m~+OWV-6>TLc~^3`#m975Mtt2!1IWrW=hP7
zb&+6ecquZuH;gIv9bt+Z_*r-gEfjv8leSQ{9BoU&wdigYxVtWYpY(u;B@if`2W83x
z&m~{5UDQW061QBbM53D;?niHiUnk}0saydes$wKuYOTN9MC+vqU#KkGd$cH+oJq_;
z`flJl?5o9J37<%5HW8*1{;qTJ#f*ts)7FSpU9dht1)jrlw~^NCpg@@{A+6l69a%an
zUXm}xvQUXBCzDpz84NA-o^z(;Bb8d^0(huev0=`T2Ii-POAodRF!{t3%HOf4dhd{N
zd5w%;;EjeEkn(>HIXn0g^6rx6JuRUjU*Oss;}3Bnt)s|X7W3Y#R(pj64>xt`pWvH?
zdXU|Y)iV8U!K+~E2y>MxWw1WES&eOD0PzF7rhm?GK`m437gb*)$qFSK9ziV)wqgLQ
zx-45}+eJE)(7&2lIXl4rth|!J0mA=8n}}Bn0bPOPS3kgIItTH*OjMeb#(xNHf8%k4
zPfIIC8d5GgD@Jkm(Jbimrhh80;ChzYe2OlhLE)!9&()G`E&+zD!bV+0Ud0f#Y#=`=
zL&bNHbGh_jHl}%oJh#(j#f;IrlIgBDLl(9nZ>nM@le+@DhKq>qQ49{$0nb>n{!-58
zu%}(;>wNY=PdFtN;xUm%_M6r~50?ueR}N0SgdbMVuKZ24WgBnY;_fDHHqSo_48ed?
zo&HKqa8Q&!Z0#Nj7t<4<ENZ<#%{sSjbng6RkN`MN$1N)4GUCQjsT!8stos1iVM-NP
zf%<x<+PV^ba9TYkus)yhqJwvZnjL+)eEu%-QTaJR&Z~L%UF1#Pm|P<pr$9!<#0c@O
z6Y#?7?gjGOTYlZzv5F_8A0os$0(bNw7~eWjXv*;s?!aZ;26xwW_`9_`kkqa;8oj;l
zF7F2#7;auNn-cqP@cK@}7I~G>Qwp%kFrY8Odn2B_mtv&+chVRZjTol%hoqT8!V?*$
z@&bO;k}sTw@1y6<&Nl{CTSTF>LCpEq%oBGIuc&(2&!9Kx&b0+~^u(Y%lj+Q16=_9^
zwA(VsCgX&7>S$)1FeBeJ*Kp|=R6g#zxOX_=0Yk*Ov+o@^{qG&s-cbrHCp4y`#%D81
zF7<>C$lZ`X3%%%DoRxiiRQW8RDL~W3I)g5-R_*B-ocyFa-e=}|tWju^%Ac`LY)96S
z*Llt<z%%E}&!V_jQPJ-HLO)LCbY*qBv1<3?_Ha>S_v32lK{@N~Y|njJUW{?YP(wlO
zP{Q}@)%HAdzx&6T20QOZ6!j%<Y~5C$mx>PFyiKRZW(1#y^OIz|)#WAr{)1>~*(OeI
z|2njNVD`>S#8z_h3j`t&+XV>J{fUqMgm>KkLRWvmyZ;fb0hs@PaP7CoKkyni;D3YH
zq@s1$(N<RV)=jFv-+Fo0gG1AxNyH?X`l$zIk%T%11Vf+S?3N|VW#;F$7_dhyjqbWH
zny4*Mg(Q=PG^}?d+hv`dPPN!|c%s^#GIAFt_`kHZcvWnT=`M?SUoFhBH3tbWu6hwG
z-i}d(B;Vd#$D!l=P{V%i65wOM*`4D|KKsE7Ti44WGU70n`<X}yp?su~^DM|G<R&D7
z=A_`wixlH?;v>c5tvoEvwR>O;Kdex~uBPF^y1T_gWE+mRi(3Jr$rxKZ4hAt|dDr@r
zzPvP|Q{GG1i<+2CAy<#!%8mN{7djv!X8Uz54hB!Ks^bze5`0_O6m57{TlWD3+8)HH
z{KkfYC-XEz9}6;ckM*}1u0hmT!Ckq1hEW>;^(AREsuQ0u)}F{pz4|U{wP0V2>(dic
z!DL@v9H|yRm0Qjo4y0RUUvuYoJad|QEcefC;KpfeuF+>0OK4WnMtakz5A|DbLRC#?
z{BVbC9i1Z>(pABSX*m*0Pr@>8uC9MU3~S|u4V26O(pv^gA`2;>(%hLfT<qd7|AiSt
zI=sDt<=m^0lmRh2_?0Kd2ww_$pF3F<dbpG$sB|SX<3?p^wISsLwO23g$>Mj}&Ufqf
z&<U~BuGY!gh%vbtio-3Rx4H)^;t;9~*JootSI-$&(MwRNF!-j-R4UF;VDWO(1Ga`K
z7uOOy?os*QI~u?&Qdh0F{wVA@pfMbiB<&QdNRW|IMLoo_#of<%Ly6#18C<9pV;^Bt
zP^L%g@k`67Q<Ajj;5!;sSdX~<TRDwZ=-NJR|6?B^Ym|Y&Y28a%jSgeXD%<s0)Gc{F
zV~?eW-1J!cvftK+PqXUc1)g*kZ}$&A)He#5vJFnhQ5JXxYm!%!)x*s6Ux=-mBhYJz
zebs7o&R1Gt(nygAu@(POUZ42QyX5PQmQgAf4nS#Mn_Y^Y!=IpM@XODfC>84w|HRCA
zG)>gSbZQQZ373p8rvW^}<zbt53`w8*$<BL00qm(RN~(1z&pt3zjR$fJV@P_W5)9!j
zd=&Y38(-OnW<g!-@60NCT`XMZJ&s@_7SG<@9HjN(t($JO%~{R?R`-&n*8Wa%f7auy
zx@9=zd7?}+G^hVpKQfGOQB5@O=QM%%2~+;_M3(eu`V;sR-{Vw%+e5<`RJVpD{YM5u
zZVPj?+RLz}%?2FuDU{_>g>}=K-$IX_S5hKI@tp?X(UyOdzvE4(j7YV$*vBGYGF;p{
zC>eC`V8<dS`1Nf|64(3cUGb;efawpTNs+Ch$Oc;S;BULH<<DX+BHE|C)vLqR8|)7y
z*862Rrka|8X74+b|5bc>C+bt?dZ}3LcYCp<0hTH~b7AjHcDr;?7s`sA%NHfY2)Cxq
z7J_~ckc;j)V9L{(Y%5L>v7p*~-&1YYL6H=m7SKx9f40P(#O^N=Jjtmrs_1SiX_Ij<
z^VTE8y`dmL*FxaW2{g`|w5y&6r+)f#`K0TSy)GZ25uufBLRK0NpU;Tb$XD>(Qm{@+
z3S&ELpd3z%V{qG?NLl|}0deDT<wu$OrsS{XKS=&vgz0yK@2w^I;^$!TPDaI&9kIMb
zON17IV%8`#(T*#7@izejw0NH=D`H23dW`YeIiy@_7(HI&FQ{Iz`wYbn<Cc7pZ(d>h
z1|egO@IXZ-Qr?vA05^>nG{6-yw;ga=_y{Img}ga+QJg~)A|C2TW)~0vZTYWRphyxQ
z30cCg&M>kvF8WlPt9v!N{7D=?)SaW0w0drjBN^F!9a_W$NoEn;>dgmL-D>CYv@~0x
zVbBt`^$t540ELi_1tZiMXdzH9t`EDCt1H7Z+8Dz=UBiu3#W&sJMF5-*{nSq@IUBtI
zi<N&Wg@xcxm-FX3*hn6DT?f;7k*R4t2hVsV$}RYD><X;YTO&D2vUG+X36*y_?+C9&
zR28b4JH#8^G^)c6?k6+9o-qBn{L8Cyc5+0OqU(n+T|ZktLW|fKng$wmF>Zkt?*XB$
zOGLZUfcBp&D9^Ii7+)o?;BU8Yp*s9X_dZlV<D6s5TuVK_U*Oq+VS`!||Dnn3#IY&$
zjZzym&wcaA!19=_epc5Bz%j4d<MVFiqTiCzQT-bC*WIL(Rd@ToDjtpflJc4>f>~n8
zBFm>QyOrOh?gb8i@P=HmOsZ(uZE}AViP7!hf$pn1_yz=14dDw>VJ#_BUFHp$fh|#0
zW#tXYY*mCjBL75i$QvI&Gs-z;Koygr_Q}!}dC|iWeqQ&IAZenyfxe5dFONv}#V&O=
z{`d*S(dX$PBUbED8*tl`XvJJl%BMFTK{v7j8NISgwJv>Q^^19EQ(RoYXU!#YKk^bs
z24i!$(FxOL=?K2pD`)+(Uwriqb}mY4tng{0f5do0ZpyKWTL7i>6*P{A_}M&;oyNKN
z^V!|=ZpyXI6FP6Lp|lxFN_OK2vA*(X*Xn4?VQA4sO(!bTS-&m6Z^QBkcDdfCi-nMp
z?BR~Dr!95-L}N#{ehT!mmtTQd0`4-YkwpgE$#;;W!l(Xs2^;RsG%N6Z755HaIFWtI
zUoVteHfS%BrHT1uD?s>g_QRc@)K$t@gNV;2%vu+}{lKWeR`gcU{pu~bzy{>EsUEqz
zrWO$H8yO+qc5-@!hw+|^Zm!84jpx~)PCTMtF7=9ktnYlt=s#3Kdx0GdOKzO-(1+)o
zJNwFq!fD5(x5T@TKJWE1+CC#noDsDeZ}+gk9vn0WL4o>&eU~`QH2F|T(uQjQ_pw?)
z6^%V2o1039OM<27fwGIRJZL+5H9C4Wx;>cUAnkZ6=0(0+&m^U}9Wj$1gYfdoVDB<B
zbWiTZ`i`?}+!X%irInC4r!;HsgH+z-)GqbS-N0T^g~G{U{Dd5Fo@dSr^T6I!+I^|$
z?peHJP#XIQp$v-vkI>=mcXRLBnfh}$r_dx1*30$dTh1tW_p9#x;bk^ypPQpay(}1D
z+}Ja$-i!Cs`bfoC2MQ;X0B{@5G_zbT1EEQ>Uwe^er$YDT@_za7rsMJKLg)5zw`Eng
zqr>YWMZ<LvaT>+cx5_)Dcbo0>r9a+|@ddIOPciWyfDiER?Zy8c;(OfyR3!NyHs}EV
z-Xer&`fG=d9l-YgM|{T`Qg%8Gs364Gsk1GH(>}I7Ym5F}5WW^`Tr4IihdNAK52{WZ
zW3S@W^Mb_9cp)k6%ib3R!sY0THG3`s<bgfimc5ncj?CxF$c~EBm&ao%7cZZ1{=l`_
z)2FM<IK4Kvz{i~#;)&)apN#gGPN<trn&!F7n`;b6K;x<GW3&K!!uIMzdc21PR%6Cf
zl1^0n2Q|%eEf)&x)lVBvJ`OHZ-D@6_wogC0#|E;O<n-`_&)COE17SqoZUe}()8A)H
zu6oG?uV${`bhSUEU_*DtliJ{3BSH7Xlg=T>!vt&$$V*-o$U^t%hmzhZPD7am6tL^H
zFkf!n@Z!E<M6{^LTCzZd#}00_X_tWxLJ(Q#H*#R6*0$ja437-!RZjvSXtzp=NDswe
zGU*na2y3>|x}IBFhlcE{bVjIIP#&_>gm=nE$W4{_&z>jjEN%e(%H*7G40v?b50uCf
zdx{{`vI~!HKv(SG?+UHnRV4hjiCn824kbKl=;dU|{g*h!0SlFPjZKovP!KIBGZVKr
zBj4H8%TAr;@gWAcFe*VnF*tkx(6nKMtBK28c1@#PC%&bsRUyX|#O;!Nz1KjwV+9$+
zc_^FjQRKFx>jcQ*h^_hDd`A%c8zcUF>v>4o;4B?viTVBOWRc<wUaXc4=7pAT1IvK*
zhL=^b_o7aOy>;QUrW3<m3Sv%Pk{5h{q=qxGBfL6Yjityn14`3Y189#!U-ek2*z$9K
zd42euNtpy|75Uf5i=V<q1zgK%W9qeR$q8{ZMlp~WI;=M*Lmd0dyF4bpabD7%QF?V}
zk=Fs}clu|;<P2k?&M}b^`1O45E-${KAm<jAqoy<}R?L)nIOd-6cjf{`uq{$5!P*05
z`!3{dT5L}BJWZ%Q<s7+!Qa7@OpFCI)lO|L<wg`Us)I~P1gII9DBWw)v+H&ggAVE}%
zX{3+CuSL^=UHUHdeXF#>#M~qY!S4a1RP&u<S^0t2Xj4iSqn|MnE<*NWzR(OLZ$Cab
z<h7C*3!7zp@mp@FuZbGzJ&4-`iFx@;Ef$zLPT~BCrL&<+mSyhQzW|@v)jF*3CWCi8
zc8>g_IAr0_>LKRk2$3hD(-z{SN5QPxl47=k{Z(uod99~-G7P=oRmxGDW>ncbq=e#v
z@4S?$(cBXnI3%axgfuuH=XeyH<!kcAW0hip?|#0m0#DceK5D9oyC=~kXUlA*exiSP
z=3O)%0e98oX3^goWW5$yb%s^l#RByc6_e(tp#CYPnQt8}9p{@USWvt;A*#>^UMHYg
zLK-5pR=J>|7*|*7bQEAh!c2{r#VlA(z(d3SsPscoe<}moIRrXh>APmC3O>uAV!Tvw
zp)?)j@_dtAZYld~l^gWo2RZ0qZDJ&SxsBrS{1}vSQeoiaOYn=t6P1{Me<8S?FD_qR
z<3G+H^hbfx<ieurMKtl#+`i^kffdYWf;&nH`YY$)rpY3U%e0Hev&xI_Oau?(if+e=
zqglx$xqbaBO-(Q{M~NdTaE$eDg*_7F+~?>Sl=@hL;@fC_TNrGMCFQx6f;9woCWVQe
z^uCe5-5sho!?umCa5Ev{?t&d8mc~ybL}@bD;L8i-PR2niQ_XB4f|fz^k{~{i$nnT!
z3a{6+_+Qf@V6YQ=Doi}%v}R2`!@9jF+7qlgHsY%?<|dC{-7!F;U#`mcVbo2N&ulbU
zlaCWplg}x8qI8!TbmM(E1JdqOv(_o_IYF9ePRRIuEGgBffS;BFer2ZJZt?1;fmJQB
z6xo}%Q(-*fiY!(I$K%~oDlz;#SY98s)>$jMD*h8j)>70E>2o|knjul4OB7Wv25w1l
z<g=+49{YE{7PHjJ1;$yb4hD9yyf7rPn)xCV*$?)KKZ<fMCMylH7jJuW$KArOxzmJ7
z7cg_M)5O9J@=W<Zmi}?S%PVB3YIl;Sy!M1%jsY57jxbsV&*!&MR&krCOMG1Z?{9W2
z8s5K{-Dex)HJC=PQhGjbR2(JD$y0EC4QC~uyY)?yr*gC&7nc)|T8V`m4HeZ~faM=_
z9zF|-7}%`S$#wEfC(MH1-d9EHFw$;#ck<ERFr(@Ma1wqv+KBH3vl1iVhcJt29~E9R
zUV`7bM4asJzrRLc@zlF)Z>Ak@h{tu-?ghQ6Sm_NRqPGn>bcB)!)IB2Qy)|Gmprzb#
z-J0g^TfCLDybrKPi0FVVd+E+*NPFpi33R!R9vuruc??`=4?QIA)IK|RhO^)J)#AW|
z`Sm*V1(qt;tSgrp`RJ67EY_g4G#e8<`5R6jPXhc%<*kccaMy-ikCAsmcD<jxYU}2S
z*#T9?;iQ1=*rXwGPegaoPrch)fxZxg`u50<GjE=!08bt47t*Z54ZVmLy{+e`A&a%`
zZTF|=d7I6-)p6eDq3k_%)^87Y?{G>;emH^CZT~uA97J%ryS@w=Vrt{^Jiv3(xj{ik
zsn`8^G8t^npIv<NK^65DIVPiCI9w<=#fb--+(1@?<aonZqon2gpfVPO@e2pPeRA!l
zc7IGOQc8f}FCKj<_@E|Ezy81yDOxvzB$BXQBUqI;QdxwkOZ67wVF9?}Xu}T#5S^+!
zU~nfk$Q|Tx=FZ4M=ZTHTYUUdH=(KuchxGJs=Vp7QE2*7Oaa0nGu@WkmV#!sPkiM5P
zx?$X*%TGakLJ=H+=ok0>h=nvu&^l4`T@hCj{v8eurAp!h@T+h-l!@%UAizMIt!Oz8
z_eQF+)&jFry=WwCXan%2c7T4*HNM>A9NMs2PLE;);^1S%?#s%T!qR#X3&BexnIip&
zdv3ev=qEy(A#f(I3y;cjVi_YGmv;NDE<D5s+EfGiN0W~U)YB*ty+z@Lh{Q$R5vX6*
z-Qu`SMG?uf7tBrw@p<JsjaVzmr@A*T9vm=lCKij9H1_eA?+yDZ)T9aJfFal|ABmZO
zAyYr{qhXr6JeKpLZ<@NC2=~T(ww^hcNtk@p4;Niw5_;-*A@T0cXd7E0@Se`78+l(Z
zl#Q+6pbL3p^E&6Z^udrDj2MY*%+5#8c<h1uOV}N&IU3yTB@5YK@SI}VOIEqg$Uevu
zMMd9MGdRu(CaV|iF-3HIe9n(YmPa12qYm;S%3K;=^7UqPmLB(XpY9L!d~SXy1i3k!
zdi!|1tRI}*@AROY{JIi(pKUN@*sm7hmEd8Iw&zI-5z>0UeQ>gi@8YrB-t@HA-qO_k
z^6-#ebi2p>Jb?1pr9OANQq|@)iQ#f}F>}296w+t<SmWN{pLrv&bsPOG`vP&+{gD1o
zu=)CM_>axN{~g%;<KgfB1U6qa{%d=Ym6?_Ge*rc%{|z>6a=6zYs_Le2m}4Ny5sgue
zb7{jF)Z!F|#QdJ`N8>Z{v^A=adUb&7R$duYAzStmS!=q@Yb(t@1kb;=U${Hm#sqrI
zE><6$Zt`0AJv}bwu*a4bMo0zoT&{0+qzA_2qGIH-UU0S2ALiDJvNK&e2>d<ZURH=8
zQXkgs8}<V{I0&UJtOg1yKZE(~mUASZcRb@`KXDiWud#I=Qyx2?-1V>@RVNSpJ7%C%
zi9DASoz+5=AClTOhPuL7F582Qad4P|=h5ytao;cnD?>Mf;3QLU#!%9b;oC^Wyf1*L
z?}<6uhpN!AA~ss)@r>il-O~*?lmUWWF^}_xsfZBeVavSwA-RwXz|dSc+*@K9pEv0-
ze%!<|{W$3{LCtjr!oBl>h`}W(2I!H9Gj`5bDY3^W@&^Z%g7_@*Dt+*{B17ZDwXKTg
zMrjwbv&FD%h^y^fVFgWdE{r4CRNy69nnz4>hgjEiQ%#!nZyZLLg?Ll6k!&uO{kCJ^
zF4yg&$f3;o;n8_yEkkBFx8KC58WaRq2QhE|N_w1@^i}VnRO-P`s%g(6%p!H_!s#Qb
zS6%S3kZ4rxe0#}t9T}e=WrG1_UL3@LXk+dWp}UO{qiUE4ErF=8{c~PYz8o>X$Xuy-
z18Lyc)8xisS;$rr3OcPzKnt%va0j3aM`=Xu8z*E=U)!8_Q;wq{l%L*IK#>zcn{KgQ
zYAZf`%%JmRxRANJl#D@ghp1ljwoKY<A_=bne3-`25&msztPsYLvzR%P!F@zL|KsdK
z0%mt}n3ksDe0X+sWpYA6I;um2O>Tr9S<6NG1DZi8IIMxSg}h!Zg+CdN?w${%2ViqL
z`fW)PW-g!RPv?oUM38-|mgYgy&jh4_15?o@L#@37nxT`&ByMO~jJ*f*`4PIsNd=M6
zpq+x7Mr7f9DEiPHp+M0EEH8cCjB0-gPA4e@2E|yUcXoNDpzR!)dMceqG>yZCNQ7^E
z<H?Ytv2=?NO{5oC@GPmU(pK}-p@wtZX~H2h$PvmWm{6IqKOsEyXp}fGr3XvXW23l#
zs`5NgRKM?VKK2NQKD(i{($nNIZIOo@*%!s7oF2Z1RX8G1*s0vhOez_KCtDt7U<eET
zU;=s1;|y~=RLcN7mohG2WnpAKfTBw)_DW@d3J88+`#K8?N4@qf(;*n)%REItmD?n%
ziY9FtwWVGyl#^X{bxPdi^hl_4&7Xzy_$<}Oj#1^?RL$Tv6sCv<4W&lV+qR~R2BgF#
zC{6oFjA4<8DcYQxu9f2%*BV|;l;oP-j0@mg1jSxZJ!`iw<xhY-;ah-V;xmDT!<Plt
zSZ`qPg%n3$`ysNQ2zw-6m*#|)E)>YxP%F5;nQT38sMib!Xqd-pzNG+}A4`3Qz%Q~L
zZ<j!z7rUf#{Qe#b*}Oto>?d~GurO1ksOmVqx~Wv5jkpQTuIw@^mE58*c8BQfJ%g8@
zDA=PbTP!b!mU&WehIyWU-7@Kk{CRhN=#daaULMkvffDo#C~qsGbcV{V8484mV#eh_
z^?!NZN07DI^C(fr(pg9ncOtq`3)R`54ux`47QGhE+p4&XLU8)C=v$o$+sAs&=!Ujc
zf7ZN=NU9Qm__L(KTB;Y1K#eYcj?26ODQ0F<zMEa?3wJq?r@i<q+2qL5<c0&Ra=jCQ
zsEeKh+M=*phx(^o&KruUSY#K9v92#jxECUY-%0ntc%VhWe)^}O3_06U6Xnr5!b{Lq
zTZ(ie<~oA6DU`n%od>TTP(Bad{#v29K!s@FFXECosEbFI)SM`s22s~J=s!)i7Ks9`
zz6)Qs2`B%0n={Pif-6#UkvPWa7<|AiM+{B0flq8h<vuJMCa2aS!KnOIYE3p7!<=t`
z7^z7k@i~h_3%b@wC>7p}tRTsTJXYx~79G3bQi<rhrlm+#Xf~gDO3qMD`i2V?1P0n?
zR^yPY{P7K}VG3Swk~w2_$y&P_4smTBq4pGpDr}zp@UJAZpzP5A7sd4CoY{r8;+p2U
zHU$wy^eq?o5B?u^TB5rjq~Kfk7xdT~wXl&!GkK5v4_>5S;-z`;cc9ip3VV*<55M6?
z;iMWz#D&<A2vc_J+`%!48W)q%OFC(lUe;sZ2#)G`n}vbD3??+Qf}w|oY!oe&Wskj-
z+JW%{5`ntA4uLVUWwPp97~Ff$k9;CGU`6COa+AJR7#U;+633CO#50C!8L)GgsE9o(
z$Swvnqv2H{LEFsF`OUWzdf!LVpTUS`F0pCz#$|ZHiYzx0e})tHqE~iDZ)TetHD+uY
zISd37Ge@BqPNMWxF<w%q$R_B^Ro!--2M9n<wHeMURtoST^o|f$?%AnRkVJ{uuhlGs
zhQC{%yh+)la*Z1Mv3%-Ke-9p_f7In^zI=-vKP-1nY>!ju!~7x+c!F6f4MQ27)87cW
z;7++cc=YMAgTASFFsq1S7+7WWZWetaazE02g2&wO{JLf#&Sf&@B7iwJTt;_8F5q0T
zLiB7bv*gD25OY*B2G>v{_jOA=>T@<W+3F@-J2~RR26E9T7lQ3DF{D9?b*lZKKl<hu
zU@IKVQrz0;B>>i?e$EYfinJZ|n=W%ZV)rO_0P{Vo{lU`b^O)6>A3^rWy+p`Tfak58
z$0IxB)}1;5Uz#EdAMhLN4RENSXi8mQvq$DES3R6JQs8Jb{4D%feWtf%MgEjZiCvA%
znT{A6)YmjbCd@uDx9fCo9~Lo+gc@BtbJ0zoe7*8-#}Cnl<sTsd+i$uvH*sq=tTz~3
zqoU6sN1dYsW>7Yl)2gD?{1$Gpx4EUCLhh!^_KDZ<WM1;Nq(7sKu9|LO?;}=ilwG`7
zZH218`^3?Oi)U>Z6{;KGz01G0vTge77NXhP#Wt!f`OP{xl80Juiy81dy_oS!*aHP*
z#alEa6wIDX!|@=r$B-}2Q{?CJFj~Zt-d{%E>}Hb(I9St1k?%@a&Yqs5cYR~KOu^V?
zEpofu;#A82rY8Ha*R*#SJ4Gah<<;i!lWI&%PcGxS2Yo-AWu#q#>lT%LlJw5<=DVe?
za$>88K08XO<Z4?*+RgOGhOx38Zk7F(ITidUrur{ywFp|$`E&efW?F7%ap6!kYjF&V
z0cz0vDT0`7#*1CYy=muuHUk^ue9Bw~Iz(S9_s@+k=HE>+=o$*Qp+JREae8xkIk~&J
zKb_cNbDp$Bq@xtgneFy1xhEj>o^bXOwp#yaVBn$2MIT$u(Cy%N!Ja#Lh(I|BxtL1#
zX%{%?=~G<`7>y}G%5c%g8<RCSOC|)zM({cLrHANsfNQ6D#*I}fc64*NySh9xlqddy
zX!KbzT9l56D#s`mzA`*>$;I3Y#bK3G>X$)xTjI?z6S*~q6XPZ|io0Efr(FfZXXI*b
z==+2G2|^p$%K~4LX7G`W&Y&ze%&_18G?pMjv}wlWIi^ek{``sFNbWy*?6s=uLxx*f
z`9Zr$S_0*gRKZK$zM-V>E43q0lff^5V2^@S%$<Cy60kuswOBKWAP^D?8ct?u$$#!d
zS;VY;g6NV)K8-k78P~SH9C$5d;^OTy8WV_mmfs^X=}>lmAhH2~GhCkxM}%uIz6YZS
ziEKiP*h(;+$2pGy{!|Ur-AC0lfvU3INWzj2=UgW{-0#$JX5An)`2B<QY4>Iq5zd=%
z!bJ4;<}gYEFDR?kF?DrELaHD=8*!l3Nayq%kwyW%dnyY#Lb9&{X7CluxQRAGt%eOH
zTm^(HO^k$o2()pE06k&z=r|7boK^kwHEG~>U+qhD5sRJijVg=Y5{I;8D88CON5eVn
zu0#2T_f6K60t@f2oX18>%H>j1H0)ZS;Q??6cI6S+8tzmC+B^Ja2S&`mCgLFXe0ynj
zbFh6(cAuw;PeEFJzwvng$VdI_WMF^0r$b}=j_KsC#l`2PO9{s6Y=C7?)@~_;bN9rU
zMRt)GX?b1130|=CA3ewaz2$HsVUn^ib|TRtdF7`nlDzUB6iI$>hJ!AzHi{%VAa1Lj
zGl-)K+Pao8F}5%ewsR%XdSyJYGqb(^>HLl6`u9Er#HKY;G69lkF@gBJphm~U6{rKx
zBo3mulJN2}iL;OZIDR_-u~V7ET!9~yfUk63US1NW4@v+Mz#kET-h&kXq{D*VvvZOF
z{)qDR8+I-dmOr=`uY-u#S=%`(*&7)BF^ZGo2VoNSzg+yAiu)>h?Ul5cK8lD48#tL5
zll%@u=Pz>a>!-ikGoe>Dt?g@IOrj>AER0NE2W6545qAH*sn7j?r`^7m<zKYh%p@0O
zr3I{zu{O>&`M9^?3FB7giVzY=Z4$AJv;erblHX&%$T&om?L`%03t?HP*L?FtVyno)
zgTNK44&UI4PCM9R;3s`FF2`MMbvqqNUnRcrw#sfBTQ)oOnr#M)c?^L1A*~MQi<qm9
zvA)|A`Z>Z(9|j!p3=Gm5Y|Xx=Mi}c>93mCMb3W2K75x^wtr?l*?tyknAhHweB|D)>
zgA4#BpAxzI>vlmK7khw<Pc4I@L2e+Wo=s+yS#p$2!=FVp-su;;mgjY<2@3ZGOCT@q
zkkjaf$=TKt)#hhzwD!nYeO&Ba*xM$Sc1oiF;KyLWUc0R-Tu)gJFzUD^#vLamV`~PZ
zIBFJuhD|rHHHiVevq>iP;W6|fl%COHO-3~GqEK8)ZKdGP@#qojV`DPA5k14q$ef?w
zZ2;y|aWx3m9l=b;9yKuBY;dc|-p_M!b1~a6*tAnm`{+ET$K8BN(hc-E!4)OG?YdU;
zPDKK%EHtTGxDT<hzy#Ad1eb4l3K`L70t>HhcSiPdAZqs1RqWoNH=`?CTGEyR&$gmr
z+B?0^M)ncW^BhxbUq0P7aJ`7<iUR^fbUAHJ-3vv#`$sqI?+|ygLvD#-o}G!XP#?+L
zLzxo9M28BcNu-nNwEEFQp;Co|1#h#{Zc)0P^B0<u!hYVOjwh)@6p-QxdJ}<B)1s2V
zc^K1;pL!wmet^6iiI_-D$$S{|fA&&qlz#5VBKL_h`XI{Eh^8=ll=i3z!A1fmCJ1>T
zh`$Gk{{2n;chp8Olr=~QA?naK4*GyKe~cJ#Q4(wn-?|)}1pjh<s$yv898|Ud4tpG+
zpV1l{7Yw7{kp7z`a7AY1EMJ;4_+i1ffkFLmF%?6ABqIF;)j{ep_*;-rBB=5N39u9*
z!V(aef`>VZ#RR3n2f3zuT6@$$yk;;Yk%k-vbF_O*cbIp~aD|<km8jNGHTrT~uu*~7
zodahqT*Qq24c`;(AwA!o1j=<Kuc4~L3-}eG`o2Zw`4Yj5iuFTWIM$UM4wbYl1}z`L
zLhxq{Q9f}~xI}DGzZM3xm>57_Issvj7&B<&OH!_CE=OLe0*riSG5@s6G_@uC6{?K5
z)t8Zp(LFFAJr8bX1f#fGPWm|EUd<lrCD%PBbU?0ABp?dMN1A{Vvj(CDw+4SIiaQ9m
zH`1W6O8qO(EWu$oQjcS$p*5f$-s7zYng`VTD6fEb8$kPdE<J5t*-*M}l{HO!Dtk41
z^hVrPw3WA~ktjVEXG7;Ex8j~C-mn5u*`(Cyk#OGvC4^}W7;@vr#kr&!rCf!`(XK*e
zyY0XDHU_JI<Q3(WBcN!fP>?Vs=b-T;C$%9xrphE+CTk;!9}qW2apcvIUJ<sX%qCr>
z>QKm`Ri}MJeTJb*H%$$s)mCsV3{_)OzEsx!no^;j>niFXabAL^J0KBO-mRcnpqvL(
z-VdIvWl>oW=T_zxajS8R_yx870fW9LbVFoaA)NYw+Kk#rIj~gDM5$@od0PBH<A4iE
zgtZWAT+yPS*3LYjmIgCCJX5k(GE4IFpaUiiVJ`jz<lg=z%e}$9^Zk1j1{FD#_5;cp
zjv2brikY2ax8m~|Nf6Ofe0rsXb=qrsYPzAsrDSzV&?3bm&}z%_Xo9Qk9EfZaHAHZd
z5yu^}yEGffKEn2nP0jpdieNlrc5?RD)Zk?4Bv0Oza{9O9+|1H8k%z$Rd}NhAnGZQN
zBD?vcM#1o1@v)=D#+QIXAiMLEGez@Nb3ikzd#}5KdyxCF`^_cXJ^d5;6XAW^y%IDJ
zG$}L`bUAbb(s?j_7e<#(@KbQb=g7~7ksLVAxO^<CtjNrxoVM)kR>B+(oEoN+`i>(Q
zMuRLetjjj;3y2&woKk=?06t5|2;nf)PRmFL!#I{KmT&}g1im!dRD5cewMmrbBtrlL
z0z(UZwg#?xkCsi%LBqYdja{bgvbk4Zupzp+lKGU;MOk?&!(4su^-t|z+N(`vyPw9q
zY|U)@8WxRQrm`%F4Lds9=bi%s@&fvEM+_z>z5^zzX3wD`df|<1Ofy}C-@9>_f8yD>
zoL%VDIP#po!^SR@N|TBh%5&?zxV^yJ{oHiW<fwh9Q(%kOAllGotLUQR66F!_0Eyy-
zlJz0JFR-Y@taMjs%uIJym#hN0V$Ry^TJoC4u6uE;rQ2gXWCML$>DXyiZH#lFVa8=}
zY{f2UmBy#mC+3CY1s_}!e93pmw?we-d!nD0P<{@S(2Nj@|A@a>Cm*T3gns{ZOv?q*
zF2b%$kpH;v_;K#dkI1g1ATQ!5I242mI3FAXbcMGQ9lp)_1O7=+6=ZE>Lu99-siGgh
zxaVQ!p?^uw)9Z8XOX!m@DmPT?)eWBtMela*Guq7COy4NNzNcSmvghsO&*fUbC)*Tr
zi(sUmRnKU@y+!Mx$e{RuHW3C0V+cFt&dq=xy`1zdsFNycv_4@xLC=D*gRzS{9c0>-
zBX6Qiq~eqOB(svmo4S#u82d~4JC-s+R?&UYY(cv$g3My}H4PtS8y{qHY>I=j1ZpAN
z14|w21xrZHa7}nk{!m8pu{H+<JQ}-+$U;PQOAp(^*tBoAXtgd&+Gdhc`~@W^&y7Tu
zV-~ZWe$&D)nIWoQUKbs@ewv*+6IwZ1mG%0-g5AM;mGRPoGIzDZ3Pw%7_FlUYzq-IW
z3o9S$T1I*%&n8Aw9l2Q5iL@_OSKfzO5t!j(J)9hk*2RmVKdI-r8ZAx=OJ#_qTq&@n
z%BAQeSM7Fe+iw#`)Drc)e#X_BR1W7_kEsr;uFoDWE;SHRt1*r-5ib%|w>iNWeHe&J
zexH2QDsP~s&=656RV}57tG%FOB41Z`;N>iNCGlBJ#-rTyxbWv(i(5;@(fRz8YLSYX
z=8bKI{^^K&$NeaDUU0ZqiWm0-5r;VEwa>_XbJNGXlG^<)?{)8ezfef8k*JXohBO9~
za{ID;=U@6iI)_$WeWsRst;_6+8j^`1LdZjie9~`^ueyXnJfbV3YQ4;Eb!=ej^77FS
zv{$%i5+^B+qYjDq@my;g`U;9RJ%&zih^P2eowT=I7vO$`eL>Hb-O0}2Pjo)m-`|=t
zYLjw5%?RUjDrai7JF^~NxY+d^%iO;?5^!6&daObi#g5W7);Vh4vM+J^?wzv<^FYi4
zO$hxygdn7u>AlVgZWgEKM)*VGlv-OwRtuK>@nrvF?E0a?oZX5<BWc?}%iW&pefe&K
zlkM1&h7-0|!&cU0TlhoT4Zb(qQ<csE2_OFfAJ3V7XTW=-r79``lb1{H);-cYgPy=q
z(TjJBANnZ0sF(zFZ(Sch&c>$}^$5_tXg>^ZEx2^JKC&=lzGr{dzI>djNY^Xwn17i(
zF*ya!Up}efR+#E0w<`EJtuY+hytJgaptCleR-Mubq<HV$;ar}d6h+Eqdqds?UOfB*
zvIG7Cs;@`h{y1bNW(zvcWnl|CW%i0&|NTkL!dA%E$>J{oJd>z}si_I*WZCP{F)b1{
zb|xDOTW2Q{oxdGn`;ET;g1LVT{~Mb#DLESge;;Xk6@nPbvIaIF^b5l4zn%V}BVl6w
z$pmO&WB`Khf5r4W-oKxD`wij$dMfc}0)OFr5Ly=^VP<4!X8Bw9@Ehm<PX+)+W@hHU
z8~mH7{jU=J+k=t8A0Gae^zUxr|LFl_@Y}=R3|t*eOyQYHSOD-K&wsv1IM~@)*-1=E
z{v~4p(T84NB)0#O0XRVI(|^bS0A_9wDd@js0B+DU|CfxF9hBF9%l=`@%*^siGWria
zknEp&EF9dRCh0%*SU5TV5i>grD69X{V`qJx#s9SB<OVf<|B!Kk+NS@|1F$i3{kJWE
z4Z!u^F$36Gxc)mn02?dIKWf0n$_Z)&|H}s(%inwe9Stn3O&mdoL`j&GEZj{%^G3q-
s(bm+C<dxs`%Ga{91CqQlvwp8rC!m2N@b?S@*w{EZ;3+7?<iz3s4{gE0i2wiq

literal 0
HcmV?d00001

diff --git a/tests/testthat/test-nonconvergence_plot.R b/tests/testthat/test-nonconvergence_plot.R
new file mode 100644
index 0000000..ea6327e
--- /dev/null
+++ b/tests/testthat/test-nonconvergence_plot.R
@@ -0,0 +1,44 @@
+skip_on_cran()
+# Slow tests
+
+library(testthat)
+library(lavaan)
+
+mod <-
+'
+f1 =~ x1 + x2 + x3
+f2 =~ x4 + x5 + x6
+'
+set.seed(154151)
+dat <- HolzingerSwineford1939[sample.int(301, 60), ]
+
+fit <- suppressWarnings(lavaan::cfa(mod, dat))
+
+fit_rerun <- lavaan_rerun(fit, to_rerun = 1:20, allow_inadmissible = TRUE, parallel = FALSE)
+
+fit_rerun
+out <- influence_stat(fit_rerun)
+
+p <- gcd_plot(out)
+test_that("NA in rerun", {
+    expect_equal(nrow(p$data), 19)
+    expect_no_warning(print(p))
+  })
+
+p <- md_plot(out)
+test_that("NA in rerun", {
+    expect_equal(nrow(p$data), 20)
+    expect_no_warning(print(p))
+  })
+
+p <- gcd_gof_plot(out, "chisq")
+test_that("NA in rerun", {
+    expect_equal(nrow(p$data), 19)
+    expect_no_warning(print(p))
+  })
+
+p <- gcd_gof_md_plot(out, "chisq", circle_size = 35)
+test_that("NA in rerun", {
+    expect_equal(nrow(p$data), 19)
+    expect_no_warning(print(p))
+  })

From fdbf92a5621f8f9856ba6e234c7f617b960795f0 Mon Sep 17 00:00:00 2001
From: Shu Fai Cheung <shufai.cheung@gmail.com>
Date: Thu, 5 Oct 2023 22:50:59 +0800
Subject: [PATCH 2/6] Make index_plot() to work with NAs

Tests passed.
---
 R/generic_index_plot.R | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/R/generic_index_plot.R b/R/generic_index_plot.R
index efb0c23..9d14e05 100644
--- a/R/generic_index_plot.R
+++ b/R/generic_index_plot.R
@@ -188,6 +188,11 @@ index_plot <- function(object,
       object <- object[, column, drop = TRUE]
     }
 
+  object <- object[!is.na(object)]
+  if (length(object) == 0) {
+      stop("No cases have valid values.")
+    }
+
   if (absolute) {
       object <- abs(object)
     }

From 6db620d6a7aa57aafcc648cf38fd0c82980de1a8 Mon Sep 17 00:00:00 2001
From: Shu Fai Cheung <shufai.cheung@gmail.com>
Date: Thu, 5 Oct 2023 22:59:37 +0800
Subject: [PATCH 3/6] Make est_change_*_plot() to work with NAs

Tests passed.
---
 R/est_change_plot.R | 8 ++++++++
 1 file changed, 8 insertions(+)

diff --git a/R/est_change_plot.R b/R/est_change_plot.R
index cd0839e..bed3e35 100644
--- a/R/est_change_plot.R
+++ b/R/est_change_plot.R
@@ -256,6 +256,10 @@ est_change_plot <- function(change,
       }
     x0$gcd <- NULL
     x0$row_id <- as.integer(seq_len(nrow(change)))
+    x0 <- stats::na.omit(x0)
+    if (nrow(x0) == 0) {
+        stop("No cases have valid values.")
+      }
     p <- ggplot2::ggplot(data = x0, ggplot2::aes(x = .data[["row_id"]],
                                                  y = .data[["change"]]))
     p <- p + do.call(ggplot2::geom_point, point_aes)
@@ -379,6 +383,10 @@ est_change_gcd_plot <- function(change,
     if (!missing(parameters)) {
         x0 <- x0[x0$param %in% parameters, ]
       }
+    x0 <- stats::na.omit(x0)
+    if (nrow(x0) == 0) {
+        stop("No cases have valid values.")
+      }
     p <- ggplot2::ggplot(data = x0,
                 ggplot2::aes(x = .data[[gcd_name]],
                              y = .data[["change"]]))

From 2355274a039c6199a4392379ed23c83a729fe10a Mon Sep 17 00:00:00 2001
From: Shu Fai Cheung <shufai.cheung@gmail.com>
Date: Thu, 5 Oct 2023 22:59:52 +0800
Subject: [PATCH 4/6] Tests for plots with NAs

Tests passed.
---
 tests/testthat/test-nonconvergence_plot.R | 28 +++++++++++++++++++----
 1 file changed, 24 insertions(+), 4 deletions(-)

diff --git a/tests/testthat/test-nonconvergence_plot.R b/tests/testthat/test-nonconvergence_plot.R
index ea6327e..569e85c 100644
--- a/tests/testthat/test-nonconvergence_plot.R
+++ b/tests/testthat/test-nonconvergence_plot.R
@@ -22,23 +22,43 @@ out <- influence_stat(fit_rerun)
 p <- gcd_plot(out)
 test_that("NA in rerun", {
     expect_equal(nrow(p$data), 19)
-    expect_no_warning(print(p))
+    # expect_no_warning(print(p))
   })
 
 p <- md_plot(out)
 test_that("NA in rerun", {
     expect_equal(nrow(p$data), 20)
-    expect_no_warning(print(p))
+    # expect_no_warning(print(p))
   })
 
 p <- gcd_gof_plot(out, "chisq")
 test_that("NA in rerun", {
     expect_equal(nrow(p$data), 19)
-    expect_no_warning(print(p))
+    # expect_no_warning(print(p))
   })
 
 p <- gcd_gof_md_plot(out, "chisq", circle_size = 35)
 test_that("NA in rerun", {
     expect_equal(nrow(p$data), 19)
-    expect_no_warning(print(p))
+    # expect_no_warning(print(p))
+  })
+
+out <- est_change(fit_rerun)
+
+p <- est_change_plot(out)
+test_that("NA in rerun", {
+    expect_equal(nrow(p$data), 247)
+    # expect_no_warning(print(p))
+  })
+
+p <- est_change_gcd_plot(out)
+test_that("NA in rerun", {
+    expect_equal(nrow(p$data), 247)
+    # expect_no_warning(print(p))
+  })
+
+p <- index_plot(out, "f1=~x3")
+test_that("NA in rerun", {
+    expect_equal(nrow(p$data), 19)
+    # expect_no_warning(print(p))
   })

From 84763218f42247bb00e369cfaba4c291bf1bf8b0 Mon Sep 17 00:00:00 2001
From: Shu Fai Cheung <shufai.cheung@gmail.com>
Date: Thu, 5 Oct 2023 23:07:41 +0800
Subject: [PATCH 5/6] Delete Rplots.pdf

---
 tests/testthat/Rplots.pdf | Bin 14975 -> 0 bytes
 1 file changed, 0 insertions(+), 0 deletions(-)
 delete mode 100644 tests/testthat/Rplots.pdf

diff --git a/tests/testthat/Rplots.pdf b/tests/testthat/Rplots.pdf
deleted file mode 100644
index 4ee28b0d21b3659c618ae7b5d66354e3418be5a2..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 14975
zcma)j1yo$g)-Dp<g1ZLy#v!<c;O_3OO>lP!?oMzB!QDb|cMI-r!TmLvnS1BndGCLJ
zuf;muyZ5dwRfqm|ov$e5Ma3BajBM}}ej9$9eq(-9+K%u55@r%RLrZvGUU()EM-u~}
zg`KUa0nmhmT9k){nS~X=%*@Wh!ph0aK?BbuV`u#T{7%^dXbrMbB>CWIU~g_=<U}Ip
z477I!g3RO{?Tnp`OdMaoWn*Mz1b}3}ec8Q=XyEzz;Z1CfU!!38S2QF{%5L^1BupX(
zKm%(#GZH3w12YpR5>}9cBFVoEIR9brs`skxXkrWeQ=Lh~&K3yz=tRQx+nh<!#L3Rt
z(Fimi+kcJss_<7iq)lwifaWADtSlTPOyU;SKu`uu;?|(zh?*GL8JoPu?F0nX%m&^y
zV@5+Jnz9bHy`$qt!!Xh0ZH<kd3#0@c6x^qLg10a%i2ftqu)(l29sDX+%1vcvAxDVR
zFvt%lyx%H{mQ_qzGJg!k{5U+D+8&~Me%LxG8_Nz!Z;|?G7}@!<vfQe49-)R3s@Hbw
z|H0$->em$~Y49l{O>4eS1rtr2@a==xb9^>d*~97XqLiyGqA-mXVQ#)nq5MsK51u9A
z;@5oaXVeu!j7<(F<mY6XXQl_K+h?Qbw+*UMA3|@4z(;r6KNGF!9uQ~^!LMrTmT-Np
zA93i%B!+KlVCWGMI?t8AZr%n`Q}bh3(rFEZ_#3L?c%>_4zm$LB+2BYEoyyR%mHIYL
zgW*1lmFzEEB=7X|eqAUR#bFTQt?gYm_PC6Oj246Kpzm@5<`;$K1n?FbDW}EsbQXP1
z69;QA7bR%=K&ypo&tbS6k1tEK9fnHsdsH^74a!j;Aur=+A5?n*7QV8ilF;{hvO6dr
zV)OVihv?5DPu{=CX6*zim<8>Ie?)mMFzZ_7{5n_fm1e!F$%+5S@GOeL1<<R@O*`C?
zd&c7H_D3Yp-@IMc&Wbm!;Kj~p6v3hUre)&?|3xlJNpeq7_cs0A__P9`%-Z^0I?^+i
z&lec!RIHKUINEf9JFUyN2rTx#pEqH8YwW_Z41N+op2kW=#)b2YWX*HVC=cFqevkmz
zS{w+5D-u^Kf+frMfH_o{&V@dHlRf9OBFh;Rl#b;6I4-LUpO|j2WQn<9MgFmByV1lF
z&0~eOPPlY;5(%y-W5Ir)V1`~xYy<x%4k`3juXUJK<`C!GT2!s=K`9%#)+>!-z!A%?
zLtoyW;D)~xyVEu5M2deBz_+_J@nSarGy2aXux>hK6JPZagHLaYKbEP16-!FT!{gUD
zX(vAGvX!{M@ggkNqJp%%knmuUs{KTI>ar~=0~N+J`7M>{la5jc9Tj0t8r_iW&ssN3
zu$j_*rpYgNOrN64CW)u#c*gX_q2>kG{Xf>ucVqP%Y>3%Qpu)xo2B@Z)PkxMN8J6{P
z#Z)28ahXaY-2SMi3-z7OuoT03dei~xNtz)$JVR3>Rd^NpUiH1-fc4&TA{GmSM*-cM
zAUB0A3+2d?*EGmJh!z@e&Tm0w)Zjw1T%DdWe-QdpR@r0`yIrzrY$8AO-ujTNFMRcR
zmzkqV^;SC=obncXTCJ@O{d@2;cojs%6~|&iJ{rDQcPQV|ZZ?{|lFZH0%ZbPRFzDtO
zB;fN%f|-bEwzsn%RE!`SO$`H06ZH_XV5=}D<n+u#^*a|QP}VBl$~5!EuXb+nl(QQ9
zhhVzc<+KSQ#)3_07EK>D{N8-2qwb@~%vd&sk)Mm0^!w31!gxWSHrJ|~(EFFEg-)?L
zkwP<vJa;EcIb}x3=OJo2iw-b?;&C43bYu+8T&|X$QQy~eSw3*Ikh>M?Lcw}{E>WyR
zZ9!KXPd+9%0Kj?Sr&O1JQOVCmIC_KDzy!7~Y9${nrOCG3^Ub_8IlmNVPm9sCk79n9
z4osP%P%;teh*zyZPO5X*w#Eoj`L;VODTcI0&wSF{4_5>f^tM?&g-zVDB}`>KFGmb_
z^}>Mc&sKO}+CS$=<$kH|y?_;MP6cG<RLvQUR5<Ie39$q)@^MB<%EXspWI@Cno?4dv
z%!Q967qnodp(T^0iPW!4qo@FEfSuh9m~+OWV-6>TLc~^3`#m975Mtt2!1IWrW=hP7
zb&+6ecquZuH;gIv9bt+Z_*r-gEfjv8leSQ{9BoU&wdigYxVtWYpY(u;B@if`2W83x
z&m~{5UDQW061QBbM53D;?niHiUnk}0saydes$wKuYOTN9MC+vqU#KkGd$cH+oJq_;
z`flJl?5o9J37<%5HW8*1{;qTJ#f*ts)7FSpU9dht1)jrlw~^NCpg@@{A+6l69a%an
zUXm}xvQUXBCzDpz84NA-o^z(;Bb8d^0(huev0=`T2Ii-POAodRF!{t3%HOf4dhd{N
zd5w%;;EjeEkn(>HIXn0g^6rx6JuRUjU*Oss;}3Bnt)s|X7W3Y#R(pj64>xt`pWvH?
zdXU|Y)iV8U!K+~E2y>MxWw1WES&eOD0PzF7rhm?GK`m437gb*)$qFSK9ziV)wqgLQ
zx-45}+eJE)(7&2lIXl4rth|!J0mA=8n}}Bn0bPOPS3kgIItTH*OjMeb#(xNHf8%k4
zPfIIC8d5GgD@Jkm(Jbimrhh80;ChzYe2OlhLE)!9&()G`E&+zD!bV+0Ud0f#Y#=`=
zL&bNHbGh_jHl}%oJh#(j#f;IrlIgBDLl(9nZ>nM@le+@DhKq>qQ49{$0nb>n{!-58
zu%}(;>wNY=PdFtN;xUm%_M6r~50?ueR}N0SgdbMVuKZ24WgBnY;_fDHHqSo_48ed?
zo&HKqa8Q&!Z0#Nj7t<4<ENZ<#%{sSjbng6RkN`MN$1N)4GUCQjsT!8stos1iVM-NP
zf%<x<+PV^ba9TYkus)yhqJwvZnjL+)eEu%-QTaJR&Z~L%UF1#Pm|P<pr$9!<#0c@O
z6Y#?7?gjGOTYlZzv5F_8A0os$0(bNw7~eWjXv*;s?!aZ;26xwW_`9_`kkqa;8oj;l
zF7F2#7;auNn-cqP@cK@}7I~G>Qwp%kFrY8Odn2B_mtv&+chVRZjTol%hoqT8!V?*$
z@&bO;k}sTw@1y6<&Nl{CTSTF>LCpEq%oBGIuc&(2&!9Kx&b0+~^u(Y%lj+Q16=_9^
zwA(VsCgX&7>S$)1FeBeJ*Kp|=R6g#zxOX_=0Yk*Ov+o@^{qG&s-cbrHCp4y`#%D81
zF7<>C$lZ`X3%%%DoRxiiRQW8RDL~W3I)g5-R_*B-ocyFa-e=}|tWju^%Ac`LY)96S
z*Llt<z%%E}&!V_jQPJ-HLO)LCbY*qBv1<3?_Ha>S_v32lK{@N~Y|njJUW{?YP(wlO
zP{Q}@)%HAdzx&6T20QOZ6!j%<Y~5C$mx>PFyiKRZW(1#y^OIz|)#WAr{)1>~*(OeI
z|2njNVD`>S#8z_h3j`t&+XV>J{fUqMgm>KkLRWvmyZ;fb0hs@PaP7CoKkyni;D3YH
zq@s1$(N<RV)=jFv-+Fo0gG1AxNyH?X`l$zIk%T%11Vf+S?3N|VW#;F$7_dhyjqbWH
zny4*Mg(Q=PG^}?d+hv`dPPN!|c%s^#GIAFt_`kHZcvWnT=`M?SUoFhBH3tbWu6hwG
z-i}d(B;Vd#$D!l=P{V%i65wOM*`4D|KKsE7Ti44WGU70n`<X}yp?su~^DM|G<R&D7
z=A_`wixlH?;v>c5tvoEvwR>O;Kdex~uBPF^y1T_gWE+mRi(3Jr$rxKZ4hAt|dDr@r
zzPvP|Q{GG1i<+2CAy<#!%8mN{7djv!X8Uz54hB!Ks^bze5`0_O6m57{TlWD3+8)HH
z{KkfYC-XEz9}6;ckM*}1u0hmT!Ckq1hEW>;^(AREsuQ0u)}F{pz4|U{wP0V2>(dic
z!DL@v9H|yRm0Qjo4y0RUUvuYoJad|QEcefC;KpfeuF+>0OK4WnMtakz5A|DbLRC#?
z{BVbC9i1Z>(pABSX*m*0Pr@>8uC9MU3~S|u4V26O(pv^gA`2;>(%hLfT<qd7|AiSt
zI=sDt<=m^0lmRh2_?0Kd2ww_$pF3F<dbpG$sB|SX<3?p^wISsLwO23g$>Mj}&Ufqf
z&<U~BuGY!gh%vbtio-3Rx4H)^;t;9~*JootSI-$&(MwRNF!-j-R4UF;VDWO(1Ga`K
z7uOOy?os*QI~u?&Qdh0F{wVA@pfMbiB<&QdNRW|IMLoo_#of<%Ly6#18C<9pV;^Bt
zP^L%g@k`67Q<Ajj;5!;sSdX~<TRDwZ=-NJR|6?B^Ym|Y&Y28a%jSgeXD%<s0)Gc{F
zV~?eW-1J!cvftK+PqXUc1)g*kZ}$&A)He#5vJFnhQ5JXxYm!%!)x*s6Ux=-mBhYJz
zebs7o&R1Gt(nygAu@(POUZ42QyX5PQmQgAf4nS#Mn_Y^Y!=IpM@XODfC>84w|HRCA
zG)>gSbZQQZ373p8rvW^}<zbt53`w8*$<BL00qm(RN~(1z&pt3zjR$fJV@P_W5)9!j
zd=&Y38(-OnW<g!-@60NCT`XMZJ&s@_7SG<@9HjN(t($JO%~{R?R`-&n*8Wa%f7auy
zx@9=zd7?}+G^hVpKQfGOQB5@O=QM%%2~+;_M3(eu`V;sR-{Vw%+e5<`RJVpD{YM5u
zZVPj?+RLz}%?2FuDU{_>g>}=K-$IX_S5hKI@tp?X(UyOdzvE4(j7YV$*vBGYGF;p{
zC>eC`V8<dS`1Nf|64(3cUGb;efawpTNs+Ch$Oc;S;BULH<<DX+BHE|C)vLqR8|)7y
z*862Rrka|8X74+b|5bc>C+bt?dZ}3LcYCp<0hTH~b7AjHcDr;?7s`sA%NHfY2)Cxq
z7J_~ckc;j)V9L{(Y%5L>v7p*~-&1YYL6H=m7SKx9f40P(#O^N=Jjtmrs_1SiX_Ij<
z^VTE8y`dmL*FxaW2{g`|w5y&6r+)f#`K0TSy)GZ25uufBLRK0NpU;Tb$XD>(Qm{@+
z3S&ELpd3z%V{qG?NLl|}0deDT<wu$OrsS{XKS=&vgz0yK@2w^I;^$!TPDaI&9kIMb
zON17IV%8`#(T*#7@izejw0NH=D`H23dW`YeIiy@_7(HI&FQ{Iz`wYbn<Cc7pZ(d>h
z1|egO@IXZ-Qr?vA05^>nG{6-yw;ga=_y{Img}ga+QJg~)A|C2TW)~0vZTYWRphyxQ
z30cCg&M>kvF8WlPt9v!N{7D=?)SaW0w0drjBN^F!9a_W$NoEn;>dgmL-D>CYv@~0x
zVbBt`^$t540ELi_1tZiMXdzH9t`EDCt1H7Z+8Dz=UBiu3#W&sJMF5-*{nSq@IUBtI
zi<N&Wg@xcxm-FX3*hn6DT?f;7k*R4t2hVsV$}RYD><X;YTO&D2vUG+X36*y_?+C9&
zR28b4JH#8^G^)c6?k6+9o-qBn{L8Cyc5+0OqU(n+T|ZktLW|fKng$wmF>Zkt?*XB$
zOGLZUfcBp&D9^Ii7+)o?;BU8Yp*s9X_dZlV<D6s5TuVK_U*Oq+VS`!||Dnn3#IY&$
zjZzym&wcaA!19=_epc5Bz%j4d<MVFiqTiCzQT-bC*WIL(Rd@ToDjtpflJc4>f>~n8
zBFm>QyOrOh?gb8i@P=HmOsZ(uZE}AViP7!hf$pn1_yz=14dDw>VJ#_BUFHp$fh|#0
zW#tXYY*mCjBL75i$QvI&Gs-z;Koygr_Q}!}dC|iWeqQ&IAZenyfxe5dFONv}#V&O=
z{`d*S(dX$PBUbED8*tl`XvJJl%BMFTK{v7j8NISgwJv>Q^^19EQ(RoYXU!#YKk^bs
z24i!$(FxOL=?K2pD`)+(Uwriqb}mY4tng{0f5do0ZpyKWTL7i>6*P{A_}M&;oyNKN
z^V!|=ZpyXI6FP6Lp|lxFN_OK2vA*(X*Xn4?VQA4sO(!bTS-&m6Z^QBkcDdfCi-nMp
z?BR~Dr!95-L}N#{ehT!mmtTQd0`4-YkwpgE$#;;W!l(Xs2^;RsG%N6Z755HaIFWtI
zUoVteHfS%BrHT1uD?s>g_QRc@)K$t@gNV;2%vu+}{lKWeR`gcU{pu~bzy{>EsUEqz
zrWO$H8yO+qc5-@!hw+|^Zm!84jpx~)PCTMtF7=9ktnYlt=s#3Kdx0GdOKzO-(1+)o
zJNwFq!fD5(x5T@TKJWE1+CC#noDsDeZ}+gk9vn0WL4o>&eU~`QH2F|T(uQjQ_pw?)
z6^%V2o1039OM<27fwGIRJZL+5H9C4Wx;>cUAnkZ6=0(0+&m^U}9Wj$1gYfdoVDB<B
zbWiTZ`i`?}+!X%irInC4r!;HsgH+z-)GqbS-N0T^g~G{U{Dd5Fo@dSr^T6I!+I^|$
z?peHJP#XIQp$v-vkI>=mcXRLBnfh}$r_dx1*30$dTh1tW_p9#x;bk^ypPQpay(}1D
z+}Ja$-i!Cs`bfoC2MQ;X0B{@5G_zbT1EEQ>Uwe^er$YDT@_za7rsMJKLg)5zw`Eng
zqr>YWMZ<LvaT>+cx5_)Dcbo0>r9a+|@ddIOPciWyfDiER?Zy8c;(OfyR3!NyHs}EV
z-Xer&`fG=d9l-YgM|{T`Qg%8Gs364Gsk1GH(>}I7Ym5F}5WW^`Tr4IihdNAK52{WZ
zW3S@W^Mb_9cp)k6%ib3R!sY0THG3`s<bgfimc5ncj?CxF$c~EBm&ao%7cZZ1{=l`_
z)2FM<IK4Kvz{i~#;)&)apN#gGPN<trn&!F7n`;b6K;x<GW3&K!!uIMzdc21PR%6Cf
zl1^0n2Q|%eEf)&x)lVBvJ`OHZ-D@6_wogC0#|E;O<n-`_&)COE17SqoZUe}()8A)H
zu6oG?uV${`bhSUEU_*DtliJ{3BSH7Xlg=T>!vt&$$V*-o$U^t%hmzhZPD7am6tL^H
zFkf!n@Z!E<M6{^LTCzZd#}00_X_tWxLJ(Q#H*#R6*0$ja437-!RZjvSXtzp=NDswe
zGU*na2y3>|x}IBFhlcE{bVjIIP#&_>gm=nE$W4{_&z>jjEN%e(%H*7G40v?b50uCf
zdx{{`vI~!HKv(SG?+UHnRV4hjiCn824kbKl=;dU|{g*h!0SlFPjZKovP!KIBGZVKr
zBj4H8%TAr;@gWAcFe*VnF*tkx(6nKMtBK28c1@#PC%&bsRUyX|#O;!Nz1KjwV+9$+
zc_^FjQRKFx>jcQ*h^_hDd`A%c8zcUF>v>4o;4B?viTVBOWRc<wUaXc4=7pAT1IvK*
zhL=^b_o7aOy>;QUrW3<m3Sv%Pk{5h{q=qxGBfL6Yjityn14`3Y189#!U-ek2*z$9K
zd42euNtpy|75Uf5i=V<q1zgK%W9qeR$q8{ZMlp~WI;=M*Lmd0dyF4bpabD7%QF?V}
zk=Fs}clu|;<P2k?&M}b^`1O45E-${KAm<jAqoy<}R?L)nIOd-6cjf{`uq{$5!P*05
z`!3{dT5L}BJWZ%Q<s7+!Qa7@OpFCI)lO|L<wg`Us)I~P1gII9DBWw)v+H&ggAVE}%
zX{3+CuSL^=UHUHdeXF#>#M~qY!S4a1RP&u<S^0t2Xj4iSqn|MnE<*NWzR(OLZ$Cab
z<h7C*3!7zp@mp@FuZbGzJ&4-`iFx@;Ef$zLPT~BCrL&<+mSyhQzW|@v)jF*3CWCi8
zc8>g_IAr0_>LKRk2$3hD(-z{SN5QPxl47=k{Z(uod99~-G7P=oRmxGDW>ncbq=e#v
z@4S?$(cBXnI3%axgfuuH=XeyH<!kcAW0hip?|#0m0#DceK5D9oyC=~kXUlA*exiSP
z=3O)%0e98oX3^goWW5$yb%s^l#RByc6_e(tp#CYPnQt8}9p{@USWvt;A*#>^UMHYg
zLK-5pR=J>|7*|*7bQEAh!c2{r#VlA(z(d3SsPscoe<}moIRrXh>APmC3O>uAV!Tvw
zp)?)j@_dtAZYld~l^gWo2RZ0qZDJ&SxsBrS{1}vSQeoiaOYn=t6P1{Me<8S?FD_qR
z<3G+H^hbfx<ieurMKtl#+`i^kffdYWf;&nH`YY$)rpY3U%e0Hev&xI_Oau?(if+e=
zqglx$xqbaBO-(Q{M~NdTaE$eDg*_7F+~?>Sl=@hL;@fC_TNrGMCFQx6f;9woCWVQe
z^uCe5-5sho!?umCa5Ev{?t&d8mc~ybL}@bD;L8i-PR2niQ_XB4f|fz^k{~{i$nnT!
z3a{6+_+Qf@V6YQ=Doi}%v}R2`!@9jF+7qlgHsY%?<|dC{-7!F;U#`mcVbo2N&ulbU
zlaCWplg}x8qI8!TbmM(E1JdqOv(_o_IYF9ePRRIuEGgBffS;BFer2ZJZt?1;fmJQB
z6xo}%Q(-*fiY!(I$K%~oDlz;#SY98s)>$jMD*h8j)>70E>2o|knjul4OB7Wv25w1l
z<g=+49{YE{7PHjJ1;$yb4hD9yyf7rPn)xCV*$?)KKZ<fMCMylH7jJuW$KArOxzmJ7
z7cg_M)5O9J@=W<Zmi}?S%PVB3YIl;Sy!M1%jsY57jxbsV&*!&MR&krCOMG1Z?{9W2
z8s5K{-Dex)HJC=PQhGjbR2(JD$y0EC4QC~uyY)?yr*gC&7nc)|T8V`m4HeZ~faM=_
z9zF|-7}%`S$#wEfC(MH1-d9EHFw$;#ck<ERFr(@Ma1wqv+KBH3vl1iVhcJt29~E9R
zUV`7bM4asJzrRLc@zlF)Z>Ak@h{tu-?ghQ6Sm_NRqPGn>bcB)!)IB2Qy)|Gmprzb#
z-J0g^TfCLDybrKPi0FVVd+E+*NPFpi33R!R9vuruc??`=4?QIA)IK|RhO^)J)#AW|
z`Sm*V1(qt;tSgrp`RJ67EY_g4G#e8<`5R6jPXhc%<*kccaMy-ikCAsmcD<jxYU}2S
z*#T9?;iQ1=*rXwGPegaoPrch)fxZxg`u50<GjE=!08bt47t*Z54ZVmLy{+e`A&a%`
zZTF|=d7I6-)p6eDq3k_%)^87Y?{G>;emH^CZT~uA97J%ryS@w=Vrt{^Jiv3(xj{ik
zsn`8^G8t^npIv<NK^65DIVPiCI9w<=#fb--+(1@?<aonZqon2gpfVPO@e2pPeRA!l
zc7IGOQc8f}FCKj<_@E|Ezy81yDOxvzB$BXQBUqI;QdxwkOZ67wVF9?}Xu}T#5S^+!
zU~nfk$Q|Tx=FZ4M=ZTHTYUUdH=(KuchxGJs=Vp7QE2*7Oaa0nGu@WkmV#!sPkiM5P
zx?$X*%TGakLJ=H+=ok0>h=nvu&^l4`T@hCj{v8eurAp!h@T+h-l!@%UAizMIt!Oz8
z_eQF+)&jFry=WwCXan%2c7T4*HNM>A9NMs2PLE;);^1S%?#s%T!qR#X3&BexnIip&
zdv3ev=qEy(A#f(I3y;cjVi_YGmv;NDE<D5s+EfGiN0W~U)YB*ty+z@Lh{Q$R5vX6*
z-Qu`SMG?uf7tBrw@p<JsjaVzmr@A*T9vm=lCKij9H1_eA?+yDZ)T9aJfFal|ABmZO
zAyYr{qhXr6JeKpLZ<@NC2=~T(ww^hcNtk@p4;Niw5_;-*A@T0cXd7E0@Se`78+l(Z
zl#Q+6pbL3p^E&6Z^udrDj2MY*%+5#8c<h1uOV}N&IU3yTB@5YK@SI}VOIEqg$Uevu
zMMd9MGdRu(CaV|iF-3HIe9n(YmPa12qYm;S%3K;=^7UqPmLB(XpY9L!d~SXy1i3k!
zdi!|1tRI}*@AROY{JIi(pKUN@*sm7hmEd8Iw&zI-5z>0UeQ>gi@8YrB-t@HA-qO_k
z^6-#ebi2p>Jb?1pr9OANQq|@)iQ#f}F>}296w+t<SmWN{pLrv&bsPOG`vP&+{gD1o
zu=)CM_>axN{~g%;<KgfB1U6qa{%d=Ym6?_Ge*rc%{|z>6a=6zYs_Le2m}4Ny5sgue
zb7{jF)Z!F|#QdJ`N8>Z{v^A=adUb&7R$duYAzStmS!=q@Yb(t@1kb;=U${Hm#sqrI
zE><6$Zt`0AJv}bwu*a4bMo0zoT&{0+qzA_2qGIH-UU0S2ALiDJvNK&e2>d<ZURH=8
zQXkgs8}<V{I0&UJtOg1yKZE(~mUASZcRb@`KXDiWud#I=Qyx2?-1V>@RVNSpJ7%C%
zi9DASoz+5=AClTOhPuL7F582Qad4P|=h5ytao;cnD?>Mf;3QLU#!%9b;oC^Wyf1*L
z?}<6uhpN!AA~ss)@r>il-O~*?lmUWWF^}_xsfZBeVavSwA-RwXz|dSc+*@K9pEv0-
ze%!<|{W$3{LCtjr!oBl>h`}W(2I!H9Gj`5bDY3^W@&^Z%g7_@*Dt+*{B17ZDwXKTg
zMrjwbv&FD%h^y^fVFgWdE{r4CRNy69nnz4>hgjEiQ%#!nZyZLLg?Ll6k!&uO{kCJ^
zF4yg&$f3;o;n8_yEkkBFx8KC58WaRq2QhE|N_w1@^i}VnRO-P`s%g(6%p!H_!s#Qb
zS6%S3kZ4rxe0#}t9T}e=WrG1_UL3@LXk+dWp}UO{qiUE4ErF=8{c~PYz8o>X$Xuy-
z18Lyc)8xisS;$rr3OcPzKnt%va0j3aM`=Xu8z*E=U)!8_Q;wq{l%L*IK#>zcn{KgQ
zYAZf`%%JmRxRANJl#D@ghp1ljwoKY<A_=bne3-`25&msztPsYLvzR%P!F@zL|KsdK
z0%mt}n3ksDe0X+sWpYA6I;um2O>Tr9S<6NG1DZi8IIMxSg}h!Zg+CdN?w${%2ViqL
z`fW)PW-g!RPv?oUM38-|mgYgy&jh4_15?o@L#@37nxT`&ByMO~jJ*f*`4PIsNd=M6
zpq+x7Mr7f9DEiPHp+M0EEH8cCjB0-gPA4e@2E|yUcXoNDpzR!)dMceqG>yZCNQ7^E
z<H?Ytv2=?NO{5oC@GPmU(pK}-p@wtZX~H2h$PvmWm{6IqKOsEyXp}fGr3XvXW23l#
zs`5NgRKM?VKK2NQKD(i{($nNIZIOo@*%!s7oF2Z1RX8G1*s0vhOez_KCtDt7U<eET
zU;=s1;|y~=RLcN7mohG2WnpAKfTBw)_DW@d3J88+`#K8?N4@qf(;*n)%REItmD?n%
ziY9FtwWVGyl#^X{bxPdi^hl_4&7Xzy_$<}Oj#1^?RL$Tv6sCv<4W&lV+qR~R2BgF#
zC{6oFjA4<8DcYQxu9f2%*BV|;l;oP-j0@mg1jSxZJ!`iw<xhY-;ah-V;xmDT!<Plt
zSZ`qPg%n3$`ysNQ2zw-6m*#|)E)>YxP%F5;nQT38sMib!Xqd-pzNG+}A4`3Qz%Q~L
zZ<j!z7rUf#{Qe#b*}Oto>?d~GurO1ksOmVqx~Wv5jkpQTuIw@^mE58*c8BQfJ%g8@
zDA=PbTP!b!mU&WehIyWU-7@Kk{CRhN=#daaULMkvffDo#C~qsGbcV{V8484mV#eh_
z^?!NZN07DI^C(fr(pg9ncOtq`3)R`54ux`47QGhE+p4&XLU8)C=v$o$+sAs&=!Ujc
zf7ZN=NU9Qm__L(KTB;Y1K#eYcj?26ODQ0F<zMEa?3wJq?r@i<q+2qL5<c0&Ra=jCQ
zsEeKh+M=*phx(^o&KruUSY#K9v92#jxECUY-%0ntc%VhWe)^}O3_06U6Xnr5!b{Lq
zTZ(ie<~oA6DU`n%od>TTP(Bad{#v29K!s@FFXECosEbFI)SM`s22s~J=s!)i7Ks9`
zz6)Qs2`B%0n={Pif-6#UkvPWa7<|AiM+{B0flq8h<vuJMCa2aS!KnOIYE3p7!<=t`
z7^z7k@i~h_3%b@wC>7p}tRTsTJXYx~79G3bQi<rhrlm+#Xf~gDO3qMD`i2V?1P0n?
zR^yPY{P7K}VG3Swk~w2_$y&P_4smTBq4pGpDr}zp@UJAZpzP5A7sd4CoY{r8;+p2U
zHU$wy^eq?o5B?u^TB5rjq~Kfk7xdT~wXl&!GkK5v4_>5S;-z`;cc9ip3VV*<55M6?
z;iMWz#D&<A2vc_J+`%!48W)q%OFC(lUe;sZ2#)G`n}vbD3??+Qf}w|oY!oe&Wskj-
z+JW%{5`ntA4uLVUWwPp97~Ff$k9;CGU`6COa+AJR7#U;+633CO#50C!8L)GgsE9o(
z$Swvnqv2H{LEFsF`OUWzdf!LVpTUS`F0pCz#$|ZHiYzx0e})tHqE~iDZ)TetHD+uY
zISd37Ge@BqPNMWxF<w%q$R_B^Ro!--2M9n<wHeMURtoST^o|f$?%AnRkVJ{uuhlGs
zhQC{%yh+)la*Z1Mv3%-Ke-9p_f7In^zI=-vKP-1nY>!ju!~7x+c!F6f4MQ27)87cW
z;7++cc=YMAgTASFFsq1S7+7WWZWetaazE02g2&wO{JLf#&Sf&@B7iwJTt;_8F5q0T
zLiB7bv*gD25OY*B2G>v{_jOA=>T@<W+3F@-J2~RR26E9T7lQ3DF{D9?b*lZKKl<hu
zU@IKVQrz0;B>>i?e$EYfinJZ|n=W%ZV)rO_0P{Vo{lU`b^O)6>A3^rWy+p`Tfak58
z$0IxB)}1;5Uz#EdAMhLN4RENSXi8mQvq$DES3R6JQs8Jb{4D%feWtf%MgEjZiCvA%
znT{A6)YmjbCd@uDx9fCo9~Lo+gc@BtbJ0zoe7*8-#}Cnl<sTsd+i$uvH*sq=tTz~3
zqoU6sN1dYsW>7Yl)2gD?{1$Gpx4EUCLhh!^_KDZ<WM1;Nq(7sKu9|LO?;}=ilwG`7
zZH218`^3?Oi)U>Z6{;KGz01G0vTge77NXhP#Wt!f`OP{xl80Juiy81dy_oS!*aHP*
z#alEa6wIDX!|@=r$B-}2Q{?CJFj~Zt-d{%E>}Hb(I9St1k?%@a&Yqs5cYR~KOu^V?
zEpofu;#A82rY8Ha*R*#SJ4Gah<<;i!lWI&%PcGxS2Yo-AWu#q#>lT%LlJw5<=DVe?
za$>88K08XO<Z4?*+RgOGhOx38Zk7F(ITidUrur{ywFp|$`E&efW?F7%ap6!kYjF&V
z0cz0vDT0`7#*1CYy=muuHUk^ue9Bw~Iz(S9_s@+k=HE>+=o$*Qp+JREae8xkIk~&J
zKb_cNbDp$Bq@xtgneFy1xhEj>o^bXOwp#yaVBn$2MIT$u(Cy%N!Ja#Lh(I|BxtL1#
zX%{%?=~G<`7>y}G%5c%g8<RCSOC|)zM({cLrHANsfNQ6D#*I}fc64*NySh9xlqddy
zX!KbzT9l56D#s`mzA`*>$;I3Y#bK3G>X$)xTjI?z6S*~q6XPZ|io0Efr(FfZXXI*b
z==+2G2|^p$%K~4LX7G`W&Y&ze%&_18G?pMjv}wlWIi^ek{``sFNbWy*?6s=uLxx*f
z`9Zr$S_0*gRKZK$zM-V>E43q0lff^5V2^@S%$<Cy60kuswOBKWAP^D?8ct?u$$#!d
zS;VY;g6NV)K8-k78P~SH9C$5d;^OTy8WV_mmfs^X=}>lmAhH2~GhCkxM}%uIz6YZS
ziEKiP*h(;+$2pGy{!|Ur-AC0lfvU3INWzj2=UgW{-0#$JX5An)`2B<QY4>Iq5zd=%
z!bJ4;<}gYEFDR?kF?DrELaHD=8*!l3Nayq%kwyW%dnyY#Lb9&{X7CluxQRAGt%eOH
zTm^(HO^k$o2()pE06k&z=r|7boK^kwHEG~>U+qhD5sRJijVg=Y5{I;8D88CON5eVn
zu0#2T_f6K60t@f2oX18>%H>j1H0)ZS;Q??6cI6S+8tzmC+B^Ja2S&`mCgLFXe0ynj
zbFh6(cAuw;PeEFJzwvng$VdI_WMF^0r$b}=j_KsC#l`2PO9{s6Y=C7?)@~_;bN9rU
zMRt)GX?b1130|=CA3ewaz2$HsVUn^ib|TRtdF7`nlDzUB6iI$>hJ!AzHi{%VAa1Lj
zGl-)K+Pao8F}5%ewsR%XdSyJYGqb(^>HLl6`u9Er#HKY;G69lkF@gBJphm~U6{rKx
zBo3mulJN2}iL;OZIDR_-u~V7ET!9~yfUk63US1NW4@v+Mz#kET-h&kXq{D*VvvZOF
z{)qDR8+I-dmOr=`uY-u#S=%`(*&7)BF^ZGo2VoNSzg+yAiu)>h?Ul5cK8lD48#tL5
zll%@u=Pz>a>!-ikGoe>Dt?g@IOrj>AER0NE2W6545qAH*sn7j?r`^7m<zKYh%p@0O
zr3I{zu{O>&`M9^?3FB7giVzY=Z4$AJv;erblHX&%$T&om?L`%03t?HP*L?FtVyno)
zgTNK44&UI4PCM9R;3s`FF2`MMbvqqNUnRcrw#sfBTQ)oOnr#M)c?^L1A*~MQi<qm9
zvA)|A`Z>Z(9|j!p3=Gm5Y|Xx=Mi}c>93mCMb3W2K75x^wtr?l*?tyknAhHweB|D)>
zgA4#BpAxzI>vlmK7khw<Pc4I@L2e+Wo=s+yS#p$2!=FVp-su;;mgjY<2@3ZGOCT@q
zkkjaf$=TKt)#hhzwD!nYeO&Ba*xM$Sc1oiF;KyLWUc0R-Tu)gJFzUD^#vLamV`~PZ
zIBFJuhD|rHHHiVevq>iP;W6|fl%COHO-3~GqEK8)ZKdGP@#qojV`DPA5k14q$ef?w
zZ2;y|aWx3m9l=b;9yKuBY;dc|-p_M!b1~a6*tAnm`{+ET$K8BN(hc-E!4)OG?YdU;
zPDKK%EHtTGxDT<hzy#Ad1eb4l3K`L70t>HhcSiPdAZqs1RqWoNH=`?CTGEyR&$gmr
z+B?0^M)ncW^BhxbUq0P7aJ`7<iUR^fbUAHJ-3vv#`$sqI?+|ygLvD#-o}G!XP#?+L
zLzxo9M28BcNu-nNwEEFQp;Co|1#h#{Zc)0P^B0<u!hYVOjwh)@6p-QxdJ}<B)1s2V
zc^K1;pL!wmet^6iiI_-D$$S{|fA&&qlz#5VBKL_h`XI{Eh^8=ll=i3z!A1fmCJ1>T
zh`$Gk{{2n;chp8Olr=~QA?naK4*GyKe~cJ#Q4(wn-?|)}1pjh<s$yv898|Ud4tpG+
zpV1l{7Yw7{kp7z`a7AY1EMJ;4_+i1ffkFLmF%?6ABqIF;)j{ep_*;-rBB=5N39u9*
z!V(aef`>VZ#RR3n2f3zuT6@$$yk;;Yk%k-vbF_O*cbIp~aD|<km8jNGHTrT~uu*~7
zodahqT*Qq24c`;(AwA!o1j=<Kuc4~L3-}eG`o2Zw`4Yj5iuFTWIM$UM4wbYl1}z`L
zLhxq{Q9f}~xI}DGzZM3xm>57_Issvj7&B<&OH!_CE=OLe0*riSG5@s6G_@uC6{?K5
z)t8Zp(LFFAJr8bX1f#fGPWm|EUd<lrCD%PBbU?0ABp?dMN1A{Vvj(CDw+4SIiaQ9m
zH`1W6O8qO(EWu$oQjcS$p*5f$-s7zYng`VTD6fEb8$kPdE<J5t*-*M}l{HO!Dtk41
z^hVrPw3WA~ktjVEXG7;Ex8j~C-mn5u*`(Cyk#OGvC4^}W7;@vr#kr&!rCf!`(XK*e
zyY0XDHU_JI<Q3(WBcN!fP>?Vs=b-T;C$%9xrphE+CTk;!9}qW2apcvIUJ<sX%qCr>
z>QKm`Ri}MJeTJb*H%$$s)mCsV3{_)OzEsx!no^;j>niFXabAL^J0KBO-mRcnpqvL(
z-VdIvWl>oW=T_zxajS8R_yx870fW9LbVFoaA)NYw+Kk#rIj~gDM5$@od0PBH<A4iE
zgtZWAT+yPS*3LYjmIgCCJX5k(GE4IFpaUiiVJ`jz<lg=z%e}$9^Zk1j1{FD#_5;cp
zjv2brikY2ax8m~|Nf6Ofe0rsXb=qrsYPzAsrDSzV&?3bm&}z%_Xo9Qk9EfZaHAHZd
z5yu^}yEGffKEn2nP0jpdieNlrc5?RD)Zk?4Bv0Oza{9O9+|1H8k%z$Rd}NhAnGZQN
zBD?vcM#1o1@v)=D#+QIXAiMLEGez@Nb3ikzd#}5KdyxCF`^_cXJ^d5;6XAW^y%IDJ
zG$}L`bUAbb(s?j_7e<#(@KbQb=g7~7ksLVAxO^<CtjNrxoVM)kR>B+(oEoN+`i>(Q
zMuRLetjjj;3y2&woKk=?06t5|2;nf)PRmFL!#I{KmT&}g1im!dRD5cewMmrbBtrlL
z0z(UZwg#?xkCsi%LBqYdja{bgvbk4Zupzp+lKGU;MOk?&!(4su^-t|z+N(`vyPw9q
zY|U)@8WxRQrm`%F4Lds9=bi%s@&fvEM+_z>z5^zzX3wD`df|<1Ofy}C-@9>_f8yD>
zoL%VDIP#po!^SR@N|TBh%5&?zxV^yJ{oHiW<fwh9Q(%kOAllGotLUQR66F!_0Eyy-
zlJz0JFR-Y@taMjs%uIJym#hN0V$Ry^TJoC4u6uE;rQ2gXWCML$>DXyiZH#lFVa8=}
zY{f2UmBy#mC+3CY1s_}!e93pmw?we-d!nD0P<{@S(2Nj@|A@a>Cm*T3gns{ZOv?q*
zF2b%$kpH;v_;K#dkI1g1ATQ!5I242mI3FAXbcMGQ9lp)_1O7=+6=ZE>Lu99-siGgh
zxaVQ!p?^uw)9Z8XOX!m@DmPT?)eWBtMela*Guq7COy4NNzNcSmvghsO&*fUbC)*Tr
zi(sUmRnKU@y+!Mx$e{RuHW3C0V+cFt&dq=xy`1zdsFNycv_4@xLC=D*gRzS{9c0>-
zBX6Qiq~eqOB(svmo4S#u82d~4JC-s+R?&UYY(cv$g3My}H4PtS8y{qHY>I=j1ZpAN
z14|w21xrZHa7}nk{!m8pu{H+<JQ}-+$U;PQOAp(^*tBoAXtgd&+Gdhc`~@W^&y7Tu
zV-~ZWe$&D)nIWoQUKbs@ewv*+6IwZ1mG%0-g5AM;mGRPoGIzDZ3Pw%7_FlUYzq-IW
z3o9S$T1I*%&n8Aw9l2Q5iL@_OSKfzO5t!j(J)9hk*2RmVKdI-r8ZAx=OJ#_qTq&@n
z%BAQeSM7Fe+iw#`)Drc)e#X_BR1W7_kEsr;uFoDWE;SHRt1*r-5ib%|w>iNWeHe&J
zexH2QDsP~s&=656RV}57tG%FOB41Z`;N>iNCGlBJ#-rTyxbWv(i(5;@(fRz8YLSYX
z=8bKI{^^K&$NeaDUU0ZqiWm0-5r;VEwa>_XbJNGXlG^<)?{)8ezfef8k*JXohBO9~
za{ID;=U@6iI)_$WeWsRst;_6+8j^`1LdZjie9~`^ueyXnJfbV3YQ4;Eb!=ej^77FS
zv{$%i5+^B+qYjDq@my;g`U;9RJ%&zih^P2eowT=I7vO$`eL>Hb-O0}2Pjo)m-`|=t
zYLjw5%?RUjDrai7JF^~NxY+d^%iO;?5^!6&daObi#g5W7);Vh4vM+J^?wzv<^FYi4
zO$hxygdn7u>AlVgZWgEKM)*VGlv-OwRtuK>@nrvF?E0a?oZX5<BWc?}%iW&pefe&K
zlkM1&h7-0|!&cU0TlhoT4Zb(qQ<csE2_OFfAJ3V7XTW=-r79``lb1{H);-cYgPy=q
z(TjJBANnZ0sF(zFZ(Sch&c>$}^$5_tXg>^ZEx2^JKC&=lzGr{dzI>djNY^Xwn17i(
zF*ya!Up}efR+#E0w<`EJtuY+hytJgaptCleR-Mubq<HV$;ar}d6h+Eqdqds?UOfB*
zvIG7Cs;@`h{y1bNW(zvcWnl|CW%i0&|NTkL!dA%E$>J{oJd>z}si_I*WZCP{F)b1{
zb|xDOTW2Q{oxdGn`;ET;g1LVT{~Mb#DLESge;;Xk6@nPbvIaIF^b5l4zn%V}BVl6w
z$pmO&WB`Khf5r4W-oKxD`wij$dMfc}0)OFr5Ly=^VP<4!X8Bw9@Ehm<PX+)+W@hHU
z8~mH7{jU=J+k=t8A0Gae^zUxr|LFl_@Y}=R3|t*eOyQYHSOD-K&wsv1IM~@)*-1=E
z{v~4p(T84NB)0#O0XRVI(|^bS0A_9wDd@js0B+DU|CfxF9hBF9%l=`@%*^siGWria
zknEp&EF9dRCh0%*SU5TV5i>grD69X{V`qJx#s9SB<OVf<|B!Kk+NS@|1F$i3{kJWE
z4Z!u^F$36Gxc)mn02?dIKWf0n$_Z)&|H}s(%inwe9Stn3O&mdoL`j&GEZj{%^G3q-
s(bm+C<dxs`%Ga{91CqQlvwp8rC!m2N@b?S@*w{EZ;3+7?<iz3s4{gE0i2wiq


From 471224c65e6910c513e13d71c707b93d4afcedd4 Mon Sep 17 00:00:00 2001
From: Shu Fai Cheung <shufai.cheung@gmail.com>
Date: Thu, 5 Oct 2023 23:08:49 +0800
Subject: [PATCH 6/6] Update to 0.1.5.6

---
 DESCRIPTION | 2 +-
 NEWS.md     | 4 +++-
 README.md   | 2 +-
 3 files changed, 5 insertions(+), 3 deletions(-)

diff --git a/DESCRIPTION b/DESCRIPTION
index 77a3c15..bc8124e 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: semfindr
 Title: Influential Cases in Structural Equation Modeling
-Version: 0.1.5.5
+Version: 0.1.5.6
 Authors@R: c(
     person(given = "Shu Fai",
            family = "Cheung",
diff --git a/NEWS.md b/NEWS.md
index be00b21..bec6f8b 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -1,4 +1,4 @@
-# semfindr 0.1.5.5
+# semfindr 0.1.5.6
 
 - Fixed typos in the article on
   multiple-group models. (0.1.5.1)
@@ -19,6 +19,8 @@
 - Add `index_plot()` for generating
   an index plot for an arbitrary
   statistic. (0.1.5.5)
+- Make plot functions to work with NAs.
+  (0.1.5.6)
 
 
 # semfindr 0.1.5
diff --git a/README.md b/README.md
index e99e0d2..3c0ea66 100644
--- a/README.md
+++ b/README.md
@@ -9,7 +9,7 @@
 [![R-CMD-check](https://github.com/sfcheung/semfindr/actions/workflows/R-CMD-check.yaml/badge.svg)](https://github.com/sfcheung/semfindr/actions/workflows/R-CMD-check.yaml)
 <!-- badges: end -->
 
-(Version 0.1.5.5, updated on 2023-10-05, [release history](https://sfcheung.github.io/semfindr/news/index.html))
+(Version 0.1.5.6, updated on 2023-10-05, [release history](https://sfcheung.github.io/semfindr/news/index.html))
 
 # semfindr: Finding influential cases in SEM <img src="man/figures/logo.png" align="right" height="150" />