From 611f14f8a317db6d07794c8fa13ff4f450479da5 Mon Sep 17 00:00:00 2001
From: merliseclyde Posterior Distributions of Coef
Density intervals from the summaries from
coef
.
confint(coef.ZS)
## 2.5% 97.5% beta
-## Intercept 6.668044e+00 6.783787269 6.72493620
-## M 0.000000e+00 2.192320648 1.14359433
-## So -5.751335e-02 0.307426845 0.03547522
-## Ed 6.476458e-01 3.227313256 1.85848834
-## Po1 0.000000e+00 1.462613496 0.60067372
-## Po2 -2.798292e-01 1.425802768 0.31841766
-## LF -5.305761e-01 0.989284449 0.05933737
-## M.F -2.390702e+00 1.728841837 -0.02702786
-## Pop -1.287358e-01 0.005151215 -0.02248283
-## NW -2.381958e-05 0.167769999 0.06668437
-## U1 -4.897097e-01 0.366114580 -0.02456854
-## U2 0.000000e+00 0.669002969 0.20702927
-## GDP -3.073416e-02 1.216292967 0.20625063
-## Ineq 6.537824e-01 2.105661202 1.39012647
-## Prob -4.100642e-01 0.000000000 -0.21536203
-## Time -5.249789e-01 0.052986444 -0.08433479
+## 2.5% 97.5% beta
+## Intercept 6.670469697 6.785152320 6.72493620
+## M 0.000000000 2.202308895 1.14359433
+## So -0.056058236 0.301798257 0.03547522
+## Ed 0.665575996 3.213773443 1.85848834
+## Po1 0.000000000 1.450391161 0.60067372
+## Po2 -0.221171926 1.458348086 0.31841766
+## LF -0.403621472 1.124187338 0.05933737
+## M.F -2.192806188 2.172851229 -0.02702786
+## Pop -0.124776006 0.008026921 -0.02248283
+## NW 0.000000000 0.167554995 0.06668437
+## U1 -0.493080018 0.390273539 -0.02456854
+## U2 -0.003165238 0.659330844 0.20702927
+## GDP -0.062907223 1.218823142 0.20625063
+## Ineq 0.695530392 2.143256295 1.39012647
+## Prob -0.414228765 0.000000000 -0.21536203
+## Time -0.522037075 0.058013400 -0.08433479
## attr(,"Probability")
## [1] 0.95
## attr(,"class")
@@ -505,7 +505,7 @@ Alternative algorithms }
)
## user system elapsed
-## 1.398 0.000 1.399
+## 1.375 0.003 1.378
system.time(
for (i in 1:10) {
@@ -517,7 +517,7 @@ Alternative algorithms }
)
## user system elapsed
-## 1.325 0.003 1.328
+## 1.278 0.008 1.285
which is faster for enumeration than the default method=“BAS”.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
zUk!n&I03={3eTTPDXCbd7Fjp=;CG7*)pZss3gnuKK->n;^?L7Naw`+gQfC5F!w>{RQ0+H$yWc&bwe5@;-fe_P2-golKeLfhRce z9G12kenJ~GD%l~ZErIq|hv-e_gYTxc&_vx&ryIAhte6tPKG#bY?{?X(d61&NmaV$@ zwAoBjJ;9o^V$p24F*5ZSC&`{HK&wCIvvhdu-D&UeJPROMOY;vc{jSBNtPi{F`c2%T z-pICbmM@ExVX0^^CVaZCF;l!Zo8X%Il+jJsLb1Tq^Q>rNspEZuF)w|{#2D9dDf(9W z-lK%auwc9k@vXU_SMyT&6FM=hh^hqEEZf$rV>toO_9t=i&xABUq+xj*vpzm*vzOXD zba&^sSNF&lmoo%|m0NaUe5xW&^*GNrFJd#2mVIut^2)|?fv+`lXh?gz4uHU?x4~Kr zH0~vXgA&JVYE>0}*T$rzh%g_-;Cd~)RI4bqBQAYpZGQKvXyXg+CCh;$1ON$wyOk=m z(1H=bc>1)0-;{r=2>kE=i8%8A@+NYr6Qi6pJOQ#g_sI QhZEJ!<>kvyeiP4PUfC};WOhCgDtkD$ z(ksSoa95ziCdqV$RDqVxTH9@i=yC&BE!SR0Al_O!L3=7Thvpd9f+BvQ1x=SdIVr*$ zZwX37oI0G7{9f?KY&tH6yzP8~zv8Y7H8?CbC@Op==iS(C5LO%C_Vnp+v95N}n7b9* zKlb{ZJm0jm(e3pxMh~WfCix(SB{Jf9Q=YTqMsEQ~?J^Q6-B#$v(c8h7mK(Jtk09a% zG*^%aqxcP8y6ilRLT@ENMKaJFnx5WJ2c|xHiXXH2^IdyK#Ja$Fo8 M5u41XKTOHrL{5yGUSrxBWy%$wduavSex6kC^z4a)b#X(_|8oSDcVN# zK9a$eB~_uHa!Yefuarp>A~2m8&%h5$$+=V84y3rxQGr4P{fL%8NdR_^+QmQQcs EMj_iT-7|iK8)Hx&ObjSsQ}cHJ>g!%DLHm _$pGVq;qJHP#0tiu#w{qXRaL4g*^sf0v0SWeBvC z{+xvZKRn$`ehtwsrw7%^DYnVatjP2=NUJdqHjJg$6)bv=+L7Gb=`r=m(Cy`n*t4}u zcQ==d`13|GBsLg6bm=dpPE|8bkK&-TB|GYsgXp2bjj!aNE|K}jJ<_R0y0P%XNBybD zE@qW6r!&|&8|SYVJ2X(9()%u`#(=#-f?z1yRPI{A{jE-I-%@oQ{ANsmYV73jC0p+M z1xvfXI29E0S{@iMZA{=qOn=+Tmy_O{e^oa~xuP^FJ)Ea(i~&DgV`-A{${PH#)h>5c z&Qnj|h;^oYDkxV1r?`Q2Ej1r1C}DzxP!S+Y?Vy_FcnmUyoKhaNS9&g%{w)D@l}rL!+IC=g0|1M>$m zu4>##jtJ0yYLnpf96 hKb^NP=HPa`$;0=n^gbFnf@F%qy&y`{+`lmmfxV6NdN zC>q>V9OcfRpj;7}Uu>zkGVN`uTw0bln7$j{xqoHN@cy>=J(do;xi$G(wSB|1ASGw> zm=w*4W$7R@uKT8{o)4tGFypc5E#uWb3#uup11)ATkdRsDasDED3?o>9Q;oUt5t?z( z&SBj%ld 0s^sS=7e9pJ8_na$7EwYtqZ;kDbDF^%JLO4ExqGU zmaEqpIDY;INVLW{ZTX3haX01@Bs@MZdSNhV!GnWs^HxXm`|o2ZSG1$ISV+fUeD4}7 z;NbxMWB$yC*K?op|Jn^^uKoP6#IQ!Mb=#F`t-st=E1qkl &ynKl^`0w!R9mFyD@uH5uf{_vz-^d< zZAF6%c1r6m+CG2m^W@tslp8lVgC8KVbK_PmGE8L0-F1xqIaoNOC<*IB)IyEZtux9>tu zo+}N6jp!{`C-KNe2i^tb1TpVb_}QgK)@3}68`!Xu%^Ryu(ku}Rz=Zj##VnhAJG5Y$ z<|2fEczWsBVB4?W>{y%X%;JVTHRp;+t=8eD*|c9h$MSaeWf|wg&<%+qd^^H(a-zON zLHRNt4tC{ZJiJrA4Gv!=@5mtv!4J?=AQ+ZYmN_vo+}+j57K#qu2Fz+>0)+#P0ZhV= z6gVcu7H8{YXER+%nDt~d3BcLs<@lz5_+Zwj9%vQ6D(>1CzV$G|RBwN8!1cZeqOjHI zJrZH%P6Oz9+Zx=^#@1(Ux6$YraP*F*Z2#FTOEpaWM#)C;_*~)e;DSkWfZnE|h?+07 zcRhhHmZLkc7<8vKw_lt$mGAW7X;Q(vph)%gs^$YeC@}qL3qwtpew%Dzi_c+Tq{qwb zl2-)_zlsp}aFItWU`Zxb;}9>e*vd?U&BLJ9`PvV~heTwngPj+Yntt{1>K5+htBo3g zrOPgw#=71?5KTtuFFSK~FMvjmTyh+&+`B(Hn~Y&k=~TNB9uTvlfT`~uI_yPrEUbBe z(@SKHw2Ynl3e4=oRBycov#RW_77B!DD}x0B``zz)?1URYbK`Km#RYN2GAkHQj7`JC zPes-y>g^po{@t{#9DNfYcdt_e36@*#=t1jOQ&9z{o}`GqjR%&8q@UdqcF)bD%wG5K zbg8em1Std$1ug4Fucr)f`E*^~r^sOGkKe7!{|rN)c#(n3+@k8ru@tD#xHAyy(z`p- zAv8L?wQp;aTf(uDi>nBU9db17S}xjqs>fc!UuWL#q@h||Rwlnc5aP~tP$?m-^|8qr z6HgfFl)X|UIR(rMu)M!_ (`|D{1x?XArAiyEpYd=Js7&eBip_;m-wxcs7~a$8rrwAP=AG1obP zdg4@`@x##%>Z_gU+f>c^4x5M3XADwJod<2Bi(OUM!^Kh$K>|W$xT_7EK$Tc(q$|s7 z7(#yg(y@4{4b$0H$1Kt2QZRq#KKMB{>;7+%=2vWe5|7vRP_7K4EmhrG^vA%@PH8}O zQKJL5YC3OY+RGBkDOWD_> wC@8kse2?>@N5U#9R9?`Cw?yHw>9rE%*@># z>4=2nKz%pjEq_~;fGzi|ihNT&zyNhmc=j%j@M8L1Q=-0aRA_G6SJmaWxWKkXQN8Ne zn3gQANZMBVS;6l+6Ynqa?WQ`iQXz6>-~jV;k(BJ}zvx3gf)?bgwO)WZ5VL(A7TJC+ zdBt(hEOqKDXKdZi6m%xnH-6i^UQr$B{)p}gn~B@&%jigrb)r}-)0x Ln!d_ z^j*UdU1%9ZAv-wAshS?sMt75>l;Y*p)U@`G=5jM$QMCeuWsE}XoVq KT;U+3!I$g;%(CbCxH;cMhjdQa`7%t(2D<3{TiC0xViGq;nlv~ z!02i47#m+`!d#1~!Vlk--OMdIUud#t+v-=t-Ps8LYRj1y{YgagsF|zd75zitE0hu$ z(#xs3UZ(KG_nRFxSGlLf)*j=$vA+3L{6cNxNg@LRk>AKQ2<4Jt=2|aamK2efEH~|J zS|5*}hUBVe$^nMt=rRU$&*}<4+TUy}0(;h1trfbAK)6VMAr~s)$^ZZH&tT>dIrRHu U4o?#<2c{UKqIn}<$?Bi~2Wpd8!vFvP literal 37824 zcmeFYXH-*L)HVvDB4Xeu3QDygSSd=8eiRi^DM|?~6ghx&5JTtzOHgSF3MfrPDG4O> z5=c>m9H}CNp3np&KmuY2B!t`@IPZJEG4AjC@5>mPkiGU=bIvvAGoLxv&T~t16VP7i zy#fLPAk!OH?+OU)3KI|za@(^T_y)40*BN-*eeagZRp1?Xts6JOfj5z$8}=ar0&?;E ze}XB~%N_y(M+Ho;UbYU)UZB7-&v`?ZxC|W`(_87~JI cRc5T(9hh2%ULx z!+}XAEv*$FR%kwORZZa06DjMfS9e{#^Vii6ikB?~FT0BD5Eaep=9n ?O$-D&wrfkA|NGzn9~LN=9ze-k8lWRQn47dmv?f{(IS-tIq&w%- zemE?+O4~(#I`UA}!NW4>?;(Y^KaawV#C}&+>2sDNYZpA;dBn88^K@RFVm*oe`PD{P zz<2TqEeq*_4$$|^PBUo^aH7C&@uX}CcW~9%hjHpVI3rw?(Yt#aC& OXi(5BOdYFzoKNPh}T2@-L;5D*qPPSH_5nks6zs!i1&rS*TGK%R=$K9^# z2(<6*Xph;*El-hkh6 0Fs_c*x+{at@*@_pnD zCW5vbCAQPmoPLMv$b4oad_j0B;BWD=3&kggEk+(_i?LUC+ku4_qrPrVNk5LZc$D#Y z!30IZ= ty2Qne=q=85ZjTtTKw5B`GWM^MG4HvWO8(cY3d$OLI=x%&DDIJ1rn zvjL-WE6Z;P?l$UO+B?)we@*KLkqS(k%p08jZq=^TM6a?YEt|!E1ic*6mL#n|3wP+p zW82&xP3`wx Td^%hY$HvoSWTbDeJ!D#t^E5TIVALHmI;urr77`!U2O;R@A6wo`A-BsMSmRa2 zhL52qIKeE0|5_Z)4K(7T7mw7l5iiE0vGn`zGHfl-&O?*63}k=L2fQOVAnd{BEtD!U zzy?y`rVwSroTPQj^fS(>V|=0)`#+8RX)KSpWh|5er>r6575bM7$RDzMd{Vr!Jnct* zCrZ{q%n!d=CzCTY7Ul `VrsQnWQ{xuY~yY5Cp5-|`XhPXZh$r4#i* zRS1uUQ0SXGI=tTU^VGhe(zPN4qqC{g{cR@oIObOWy4^lCZDnnzdf6kr2Jp`-1b1e# z$chnOuLwJPaGLOb(NQ?M=l J-0W=y+%f_tnREW{Gz8JuRHSsrXeBBWhiJrn0d?nUVTtTyCn)@63Gon3c(>n znG}`Y5E+$CnER`FUrF&yHAYUrhX8fA;S>G}EKyP-2iBMf+fS+rqz}U8Cq+CUKi0U@ zts}2FKh>vltJ*zs`YY?EZ4(T;M+p<2?se0Fi}VO5`-EYtqkmh2dr@U*<-d#gXCXf| zRx6DbGeQQcN%O`@V_q=<@q>8=o>=gL8bv<^ ;%-Vrk%fpj~h zKBD7baNuGW3hyTsY~~P{F+XS}FajNcBvOkf?nP^><1xp+aZnPiL*I}CJxzE=M>%AS zPpil9>C^Z`re*?e^hu5`U^DWp|HaNy5$ZMhZoG&1wESW~rCEc8%tiW#+&1Uc+S2K4 zub2q&X@$9}Kv~HeoG+8q;@vGbUvH#xPsghwgX>MmZuB%=O3vyo2(?2P$z7xvS!8-r zNr QBx9} z>|z4;3%ll5XA_!D{cEbGFgc){2G!_h^xInaZNU6hZE<)limq(;nqT{%iZ&`HSbSQw z7s*<&?~<#sV{X=$PCGp@(}R^&LG`stBgTqbqDflhbEE|5P;$3lv8P8}t|R3^>-u8r z90 a|EU`S;$>0Y`YA@^_BY|}&UjXW-|Y-}nir^Dt&N7hRcg7MY)cj~>J z2-tQluFH>Ai?^}?vwLSqkoVjW=Ge#|qbhd4dhbSb^V3#Skj0JB`iXm9da41J)XnG3 z7n>Hv $bo?mqe1Z#D-n00a>-sC%DIsT zM@BrH9_tdumJSYl1I@-s$i)21c|?nsaW;0UwPDOwW$SrX(73_&?G1afLh0{{Zp}DR zb>(Gjae;mJ@b2we-2kK{OU+lxF|?Z>!hdv;2q_x+!QmuYh;A(=+cKwr338bAb7NLp zKgS^9nu^4TV=MHVu6%-^BrUa-kgA5>ABU?4!X?VNd45A#LpI1F@X{WO2