From e554cc04c079291103b911446bb2893808ca2aa9 Mon Sep 17 00:00:00 2001
From: bcjaeger
oobag_fun_brier(y_mat = pbc_orsf[,c('time', 'status')],
s_vec = fit$pred_oobag)
-#> [1] 0.113472
Second, you can pass your function into orsf()
, and it
will be used in place of Harrell’s C-statistic:
@@ -327,12 +327,12 @@User-supplied function importance = 'negate') fit_tdep_cstat$importance -#> bili copper stage sex albumin age -#> 0.11431938 0.04283192 0.03011800 0.02978767 0.02434698 0.02379080 -#> protime chol ascites ast spiders edema -#> 0.02289864 0.01823123 0.01264641 0.00986446 0.00913007 0.00769056 -#> hepato trig alk.phos trt platelet -#> 0.00636780 0.00619816 0.00258887 0.00046192 -0.00201958
vPjs<<9wa|=fD4u1@2GDpqjaUtu* zZCo7a0w|_{@~cO;u8Nv_h-qLo!fOxccEy>_Fi-cXrd&jR(*^&xvpEQ_cUQ)9RCVAB zhZ3EWEKhI>_5=DI2knn+T^m*#Sk1h5wd2}G|Jnld9uT5!Ya)=NX(7^QXk*~<76r z<3P@S
Gaq5DCn^C+G{C$H`vOB>(T>Hx1T0Kjoer>NfUCkHO z@rAn6ajHMqubDFa$aPzKWuN!wznds<2Hb1v-5Zn)hY-*dDkn~KE`XT5G+(Wj-o`$V z^FdWy-)*;27GG+oP~QF|r_yK9FG}4n^J;3C48jR5n6GPED^DCohl7c1tIT(??#}3Q zR#-aukLsSr&@9)~SvyEX!mJ>q$0*Q@$~A89qH`==d^vQm#$quyN%nbXl2vH*E+-iw z)tU;vA+uc%W}1|_qp8w}5Tj{E{!5S5)Y(I__yyq64$;q=MagWtv> UHc6ms zM%#LS)pWvGhJA?y3!c}iV1{IYHH>P%5v1p~wB+WTXrl>1TOKFo!%WdSdAdb9>cp(` zwaP>hOga!AT}=D??_(YIRU3O#rd?#)KR<_LrUo`BXSusuR5-MmvvktoqfOGelhrrQ z*TOoqzWkZ6vbG6!?fsfCI&7}xIw?zpT1DV=ayL457zyOLc1v0ldP0@n7~w !d=c@67U}U*8hR%PjgiZ99laUXR$llKzLGHe
i#@M<`sUD4ML}2^tb%!TiG;Y$euF|$2gU91uH{DUQP|w+SXQ`hBMUQK|9#Sr! zvK!F8OcZYai;TDQFN<=(P%;%DvQ>EL3cTP)*M*pDC7iA&N>)VA0u9VZg*kn2tFbF=m5@@gv@{a@3bp%OgII@hYA zmNpkH%^$t@PZBhLW1CU6Ny-_@AH18k5fU0mQBocrln>i^-GD5(TRZoh){_4hxBo%R zSwzZUQ57}9swl{K>~N$yQrs*^rQ$qsTXx|OiwRwl>ER={sfylAQoqC7OHlE04kf~1 z+{+hQc%{YiOCd!=y+iJyujP|_PZPPQwk~QQbRJg>_-JaTcX_N7BV+_$T*hsl+Ie5Z z6p_8s2kEuKtg_W0nFe =n&7SOVRj+~zi)?aZ_&Up^XS=4&%>=(zI-XvhApbM5x?PYVqW95tf#dtbJ3 zH|=Qc8IX30h)!spS8X(eB?Z3M@i{g(0H90yo6lp -52Lx#FR_n z>}R7L5--Kvtxu>HIEqd2@;2BtZr!-S9I{#{1=rh|=mh*#)*|miI3UUxvEK*HbZpy} zQ1r$vleLH%mszG1eOY!*F?|qSox`%XFFf@}IE=NM8#C8oq>V2=;A@qtr=#uBG*Nt) zu#i8P;6E;5`Q!x|;@xw{9i5dsLD$#FJrreuGKS9$T)Udc;risQ=q*&(s>vL$`)UT+ z_385PM+eeMXd5qkxq5RXl8|hBA1gsD-BAV+7{&ZODM2hB8?0|;TH8#Po-172mEGU< zb5{4PdGKm-|EZ~hU~_M=Z=UXX^tww6u7UYr;WW2e;k=0&Rlc3t^D7|g#Vslgo+Ce- zR#le&KHrF3fF`z2!>oL;;wiIHqoHQbNymTE)FxFP6*fW8vwP&M6#Lr)Zlc@cYP)hg z9oh!H>B)`mQBnrRO&F`oH7SoPWgEMP3V(1*e>3~!Z0Avx`0#EP&u)q9Gdo9 <)l5~QVI5{X5ku_^9 zDqb<35J-5moxG?zr!Fg8;*XgPm0chsF}gWP`6^ap()#Kxo#?B+c|yMNOfR*81N+a9 zn;gDxI4QU0czNnYt)j9k;|F5@+IWIlc-2brLY}NzBz7Kk)%5m)a56VZbmh-5e8O ?JCc~vxH>FBjS(1WENfA^K2L|t;g9uPW5?C_pPDP@p;SR@WzJKjKCkAs#i9} zUYkuWRVbZ7x4+%_VW+5?YV>G-#Hsdaj+IY$EbUk@6g$4GgZHd$jGA4+$d02H=a&*i zS8OJvyRos5F)dHF5S!x%L7c9FRY@w<=OvtrUv=!oT}&+YZHFO^ZXAA^SkFz~F1~np ztiw*M{d(DHd~@&DiT*)qX>l2fdKSUC)d64bM^GXYB6n053ImuNv%*~oB_XAoqj}>w z%Pml=jv24~i5jJ^a%)g$1EEfGf8QXdQ@-;pr2BhtJLDSp4CrBwWQo#s2oawcE}63s zQSRd&)x1eFncTX~&8n-362T^%_*4&WXq>lzkmxCpgjcdEC=}6gqIE9Vi1B<`>*S*K zq;94orrJX9!@SO5 t=P$%?VA zF5bvVoijV5Qj9%3Tg4>RP^v!GTROY#`gRB6&nTmumFM1%22(=|s+Anz&DMw3s|F_#~ai_MWbMviB;)m!B|lj|70lw^tEAOj!M zCXVlr1*x#h`Qm)jMUsa4OyYj%oUao`T#Wd2pW*n$nRbpumCdHXd+7|O#~Q*18YOsT zDYl-vGzQzbYV@VN=V&4R(2zC4V3p}p9}b(@^%$Bna+KD?ceBSn5I5}Zf*u+hnSPq- zpK@$7^~*Q+T3aUsbzgou^@GFeCACficb?%Ko`Jq$d=l9Yd8R?Wv0{8;7s`s9$FKrg zL!(dLT)LPzj%H 2X&;=$`)3U~D;Js4(o#JH8HNU{`vZ9{ zX@G-xRj1BX8Oz=hR0#1nQ2pHYZBMCLkMyCuWedy{GAd^{CyO-up14w } +1yNWR*mgZg2js3SRt<6gDOq(8nycOhAHzn7awg|zPA zGk j$_PGEhKNOsZ^pz
tX@PgJe%)_$d<#MeeQ3;jh_72NZ5x_+7Wop5*mT!HWx z1yZl=sB1{3whAL7Iz5TyTaxE@vyh|v*R>SW<_4h{TY?6Hdl?^?3z6eA3_Z)R+%r*r zdhFl j zt))z~@=ULe&)v8hP@X{ENRbu(pyB1aLWJcsK-UIS8Oy21D^@~0x_69MB!&xd0%1h) zf!mCIhTbiQRi-e75dWU2N?1_qzNK(fThrsQ*RzKp1OV %@D5U9F*)6rBhtxt96eBz2#QEcxmLgK0CxP8joJyjC)X~L z(+DVsPk|n48N+q{#zsPuF_zjt)Wa!cpGOc+F&K*qAtZg4Y-_m05KcO<-*9O)dRQ2O zFNwV!kumCd8 X0KNadiYOIrg2!dp&jF=c9t)1Xvzke}=aY~rl3&@fuhYp^rPx2q zMlUzDeLJV5+1I_^{4N@)cGwhJUw4?4(>UyUs@rG>YUYzpty7 &7maZ-Ek*816hD_LV=>xY(`wSOf z-Vu3%*=LBcgzgk1Wa@GCW(;05dUY%X81E|I(f7RbSt;`dqOrI{W5^q~qts70lT(uM zr-hIv}{s?Gos zbh`)lSLGAZK0+&3GiD>Fw)V;;>H+rM%DdDVUK!vgk5=#83W$h?zU2&L=8ru(S>x|$ zUFlkPT;6h9KZc!zNRYF!9p@1zuZ^n7b_xB@4_xqOwA$3gn4dt!VW+F;k-Y?4Aa`1w zfdoG4uDEiAsLo&G9aXXs?Nfrx7EscbK5{}2du}bsF{xrr1)hZHd$a?^Y2WVO`@Eb% zyKRzj&L@u(-mxKqbS2GE{m;p}4rAsu$Nh&>y~n{?k7D6!P<$-Gb$6l6L&Sp}cQ=)) zWrAs+9s}cxQkb=}AKP=QjFM=&%cQJ1ioapJWZ80RY!xbP$XZn8Sk+_rTO@hWmd(9ELOVJ%n?+S;KBV}iL zeN9Cxotcm*2S2j6&bqhpWJ~T9;c&{-SPinf0_u~8l6fi>S(!3~l?MGVGm7qaW7=8y z1nHh5Dp-w`7W(v3)R8lGd91Gm*nwxHqSZYpD7t;vBy01<@H;9_&XNPBD;A sm0!Y zbx2%c%I>&mPOT#G@a3IXiE|Wf57vBSFGc$+dnwhsx+G;Y?q%;da0qn*^xShtnv@|X zuk;r#GfMyLjD{}39i0R)dn4L>LQMA+Nesnqk3NPa$+%wmjg4I;vV48=8kI>RLa>f6 z>TbOy$=oT?iK@|HS-gwsyNeGnv1B?`b!>tVQ&c$ffSIaFamlI+>aMs5D^9Vsc?w&3 zJvZ;2e%TX}*1jz7ls-GNSD$=|gkO`$8r7g)jtQ?s{C&SQxUkj*U(3+ =Ze8t@Ve84)MqjI?@5@SlcYrX3%wLtKiM*8JJlO;}eZ zNM3e<Pw^~r_ h37Fmf;M$0e5WM3&0jLcs6~?CQ zX2u$Sff=`iS&kiY3@yIj3zTzp9!LLjgJB)!aigIpDim`SQiGy%RKES{PGqIdMW6%E zj<)S>P_Kqrx>$);m`n=Yw;Izt%EYHIAt$%%TsrpqlLFiA hZ7c1+#tAD^hi LbZRe6zj}i*g@6IZndW9R( zGy5mcUE=ZF4$U!r9?uP!TGW` ?#;@nOI-b^*6EyjJl#oqYV! 28~ZSpj5WRVc6gJmX^mczgZeQjN4@ zuPE$H839sRiNv%XjZ_a1_rD3hHH}T6_iQ#Is!ePA>L1+HhHG4fLAP|Eg%IE4`<6 z;P_F8{ewff6`g8EyZX9qq50jgfo&i5h #;ko*51#55i( zhr3|G-7r`u{P;!o=Kp%Ey?b)FrhLYL<#Xl;z z=i3!*?sU47#-q;$s@{nnCjG0IbUGu)=r_998Ey7nC$^#5u24y_m}4M4z4s(fCw(Bt z+3^lc&XX2)vC|VK5t*qQ+V@ipS4sgR%yGbETsN(81@ukzs&nn-67_$rG-|(LyJ)8F zuc5$LoZ%m7wa^a6wCZB)Sk%Kb?_#4E7a=j*GMuG<_yI<2_4&EVCmL0$@88aC#7b># zcY?`KO{o|HUOj(b5l(ylbgMUHSC;|BLXF@p+^6B{#HT6<>Qk5zGRz3-?aNlg=Eb)x z OF zH2>kl!IGP_-xD&iH`HP-6O0@^xaT3Oucxp4dA`X^%iLMHT-