GLIM 3.77 update 3 (copyright)1985 Royal Statistical Society, London $TRANS I O $ $C ************************************************** $C THIS CODE ESTIMATE THE BASE LINE MODELS OF HGR $C 1966 MIGRATION PAPER $C $C US MIGRATION DATA $c 1975-80 8% SAMPLE, AND 1990 5% SAMPLE EXCLUDING AK, HI, DC. $C LAT LON (IN DECIMAL DEGREES) 1980. $C IMPLEMENTING DISTANCE MEASURED USING SPHERICAL $C COORDINATES. $C The original coding was kindly provided by Grusky, David $C $C Differs from xie.glm in coding CA and NY dispersistence effects $C $C Jan 29, 1996 $C GE LIN $C ************************************************** $CYCLE 25 0 1.0E-4 $ $ACCURACY 6 $ $UNIT 2304 $Echo -- $data list abolished -- $data list abolished -- invalid function/operator argument(s) -- invalid function/operator argument(s) -- $data list abolished $C *********Factorizing crossing effects************************* $CAL LCDIS=LDIS*cross$ $CAL cont=%if(cross>1,0,1) $ ! contiguity !*********6th neighbor was not significant in the model $CAL cnt2=%if(cross>2,0,cross)+1 $ !get to 5th neighbors $Fac cnt2 3 $ $CAL cnt3=%if(cross>3,0,cross)+1 $ !get to 5th neighbors $Fac cnt3 4 $ $CAL cnt5=%if(cross>5,0,cross)+1 $ !get to 5th neighbors $Fac cnt5 6 $ $C ************************************************** $CAL DIAG=%IF(%EQ(O,D),O,0):DIAG=DIAG+1$FAC DIAG %Z$ $C ******************************************** $C CREATING UNIFORM PERSISTENCE VARIABLE $C ******************************************** $CAL UNIF=%IF(D==O,1,0) $ $C TPRINT UNIF 48;48 $ $C ******************************************** $C CREATING REGIONAL EFFECTS $C ******************************************** $CAL REG=0$ $CAL NE=%IF(O<10&D<10,1,0) $ $CAL NCEN=%IF(O>9&O<22&D>9&D<22,1,0) $ $CAL S=%IF(O>21&O<38&D>21&D<38,1,0) $ $CAL W=%IF(O>37&D>37,1,0) $ $CAL REG=(NE+(NCEN*2)+(S*3)+(W*4))+1 $ $FAC REG 5$ $C DELETE NE NCEN S W $ $C TPRINT REG 48;48 $ $C ******************************************** $C CREATING SUBREGIONAL EFFECTS $C ******************************************** $CAL SUBR=0$ $CAL NEW=%IF(O<7&D<7,1,0) $ $CAL MA=%IF(O>6&O<10&D>6&D<10,1,0) $ $CAL ENC=%IF(O>9&O<15&D>9&D<15,1,0) $ $CAL WNC=%IF(O>14&O<22&D>14&D<22,1,0) $ $CAL SA=%IF(O>21&O<30&D>21&D<30,1,0) $ $CAL ESC=%IF(O>29&O<34&D>29&D<34,1,0) $ $CAL WSC=%IF(O>33&O<38&D>33&D<38,1,0) $ $CAL NM=%IF(O>37&O<41&D>37&D<41,1,0) $ $CAL SM=%IF(O>40&O<46&D>40&D<46,1,0) $ $C CAL NSM=%IF(O>37&O<46&D>37&D<46,1,0) $ !put SM to NM $CAL P=%IF(O>45&D>45,1,0) $ $CAL SUBR=(NEW+(MA*2)+(ENC*3)+(WNC*4)+(SA*5)+(ESC*6)+(WSC*7)+ (NM*8)+(SM*9)+(P*10))+1 $c CAL SUBR=(NEW+(MA*2)+(ENC*3)+(WNC*4)+(SA*5)+(ESC*6)+(WSC*7)+ (NSM*8)+(P*9))+1 !combine NM and SM $FAC SUBR 11$ $C DELETE NEW MA ENC WNC SA ESC WSC M NM SM P $ $C TPRINT SUBR 48;48 $ $C ******************************************** $C CREATING WITHIN SUBREGION CONT VARIABLE $C ******************************************** $CAL RSUB=0 $ $CAL RSUB=NEW+MA+ENC+WNC+SA+ESC+WSC+NM+SM+P $ $tab the rsub total $ 258.000 $C CAL RSUB=NEW+MA+ENC+WNC+SA+ESC+WSC+NSM+P $ !new NSM subregion $CAL CONS=CONT*RSUB $ $DELETE NEW MA WNC SA ESC WSC NM SM P $ $C ******************************************** $C CREATING WITHIN REGION CONT VARIABLE $C ******************************************** $CAL RREG=0 $ $CAL RREG=NE+NCEN+S+W $ $tab the rreg total $ 602.000 $CAL CONW=CONT*RREG $ $CAL CONR=CONT*(RREG-RSUB) $ $DELE NE NCEN S W $ $C ************************************ $C TRANSFORMING CONTIGUOUS STATE VARIABLE $C ************************************ $CAL CONB=CONT-CONR-CONS $ $C TPRINT CONT 48;48 $ $C TPRINT CONS 48;48 $ $C TPRINT CONR 48;48 $ $C TPRINT CONB 48;48 $ $C TPRINT CONW 48;48 $ $C ************************************ $C CREATING regional restructuring dummies $C ************************************ $Cal NtoS=%if((O<15&(D>25&D<30)) ,1,0) !north to south $Cal NtoW=%if((O<15&(D>40&D/=44)),1,0) !north to west $C ************************************ $C CREATING AFFINITY/DISAFFINITY MATRICES $C ************************************ $C FLORIDA BY NORTHEAST AFFINITY $CAL A1=%IF((D==29& O<10),1,0) $ $CAL A2=%IF((O==29&(D<10)),1,0) $ $C ARIZONA BY NORTHEAST AND NORTH CENTRAL AFFINITY $CAL B1=%IF((D==43&O<22),1,0) $ $CAL B2=%IF((O==43&D<22),1,0) $ $C UTAH BY NON-WEST DISAFFINITY $CAL F1=%IF((D==44&O<38),1,0) $ $CAL F2=%IF((O==44&D<38),1,0) $ $CAL F=F1+F2 $ $C CALIFORNIA DISPERSISTENCE EFFECT $CAL L1=%IF((O<38&D==48),1,0) $ $CAL L2=%IF((O==48&D<38),1,0) $ $C FLORIDA DISPERSISTENCE EFFECT $CAL M1=%IF((O>21&O<29&D==29)?(O>29&O<38&D==29),1,0) $ $CAL M2=%IF((O==29&D>21&D<29)?(O==29&D>29&D<38),1,0) $ $C NY DISPERSISTENCE EFFECT $CAL NY1=%IF((O==7&D>9),1,0) $ $CAL NY2=%IF((D==7&O>9),1,0) $ $C TPRINT cnt2 48;48 $ $C *************************************************************** $C get population and other association terms into interaction terms $C *************************************************************** $Cal intpop=%log(pop(D))*%log(pop(O))$ !pop term $Cal intarea=%log(area(D))*%log(area(O))$ !area term $CAL grprat=(grow(d))/(grow(o)) $CAL gr90=(gr86(d))/(gr86(o)) $C AL gr90=(gr85(d))/(gr85(o)) $C AL grprat=(grow(d))/(%sqrt(grow(o))) ! better off $C AL grprat=(grow(d)*grow(d))/(%sqrt(grow(o))) !better off $Cal indp=Ldis*intpop $ !interaction between pop and distance $cal wgt=1-(o==d) $tab the cons total $ 164.000 $tab the conr total $ 52.0000 $tab the conb total $ 46.0000 $tab the freq total $ 18600398. $tab the fr90 total $ 11115874. $Cal FREQ=(Freq)*wgt $Cal FR90=(Fr90+0.5)*wgt $tab the freq total $ 1745327. $tab the fr90 total $ 1036468. $CAL %n=%CU(freq*wgt) $ ! $CAL %m=%CU(fr90*wgt) $ ! $PRINT 'SAMPLE SIZE N ' %N $ SAMPLE SIZE N 1745327. $CAL %L=%LOG(%N) $PRINT 'log(n) ' %L $ log(n) 14.3725 $PRINT 'SAMPLE SIZE N ' %N $ SAMPLE SIZE N 1745327. $CAL %L=%LOG(%m) $PRINT 'log(m) ' %L $ log(m) 13.8513 $C ************************************ $C DEFINING DEPENDENT VARIABLE AND ERROR DIST $C ************************************ $YVAR FREQ $ $ERR P $ $C ************************************ $C MODEL 2 vs 3. $C ************************************ $weight wgt $ $tab the freq total $yvar freq 1745327. $error p $link L $macro wtsta! $PRINT 'PEARSONS CHISQUARE ' %X2! $PRINT 'LOGLIKELIHOOD RATIO ' %DV! $PRINT 'Degrees of Freedom ' %df! $CAL %U=%CU(%SQRT((%FV-freq)**2)*wgt ) $C AL %s=%CU(%FV*wgt ) $CAL %U=50*(%U/%n) $PRINT 'INDEX OF DISSIMILARITY ' %U! $calc %b= %DV - %df*%log(%n) $PRINT 'BIC Statistic' %b ! $endmac ! $c *********Fitting three final models ********************** ! grusk's final model $FIT O+D+REG+SUBR+A1+A2+B1+B2+F+L1+L2+M1+LDIS+CONB+CONR+CONS $ scaled deviance = 167997.5 at cycle 4 d.f. = 2135 from 2256 observations $dis E estimate s.e. parameter 1 9.66914 0.0236748 1 2 -0.249474 0.0143992 O(2) 3 -0.820104 0.0171643 O(3) 4 1.31397 0.0113825 O(4) 5 -0.674585 0.0154582 O(5) 6 0.780153 0.0119294 O(6) 7 2.49638 0.0109615 O(7) 8 1.69384 0.0113754 O(8) 9 1.69003 0.0112957 O(9) 10 2.01570 0.0113622 O(10) 11 1.27276 0.0119345 O(11) 12 2.22961 0.0112926 O(12) 13 1.82635 0.0115545 O(13) 14 1.09262 0.0122514 O(14) 15 1.03955 0.0122967 O(15) 16 0.579628 0.0126426 O(16) 17 1.07244 0.0119568 O(17) 18 -0.136653 0.0152476 O(18) 19 -0.309732 0.0151814 O(19) 20 0.194708 0.0134092 O(20) 21 0.660465 0.0125712 O(21) 22 -0.655771 0.0160977 O(22) 23 1.07317 0.0121376 O(23) 24 1.50030 0.0117964 O(24) 25 -0.0850426 0.0139479 O(25) 26 1.11469 0.0121413 O(26) 27 0.586477 0.0129132 O(27) 28 1.08113 0.0121502 O(28) 29 2.35625 0.0114655 O(29) 30 0.515619 0.0126917 O(30) 31 0.796813 0.0124799 O(31) 32 0.574156 0.0129465 O(32) 33 0.419750 0.0133212 O(33) 34 0.287394 0.0133102 O(34) 35 0.706161 0.0128616 O(35) 36 0.542997 0.0128457 O(36) 37 2.20627 0.0114685 O(37) 38 0.239054 0.0152416 O(38) 39 0.159844 0.0149718 O(39) 40 -0.411962 0.0166718 O(40) 41 1.50518 0.0120477 O(41) 42 0.486653 0.0136997 O(42) 43 1.21479 0.0133454 O(43) 44 0.577219 0.0155772 O(44) 45 -0.0193941 0.0151921 O(45) 46 1.30905 0.0129955 O(46) 47 0.609719 0.0140412 O(47) 48 2.85811 0.0134089 O(48) 49 -0.255981 0.0116468 D(2) 50 -0.694773 0.0131620 D(3) 51 0.335891 0.0105177 D(4) 52 -1.56821 0.0157225 D(5) 53 -0.193825 0.0111381 D(6) 54 1.03128 0.0103009 D(7) 55 0.663438 0.0107330 D(8) 56 0.970669 0.0102502 D(9) 57 1.29056 0.0104020 D(10) 58 0.908048 0.0108205 D(11) 59 1.48366 0.0103543 D(12) 60 1.17429 0.0106590 D(13) 61 0.949427 0.0109056 D(14) 62 0.943540 0.0109740 D(15) 63 0.199309 0.0117982 D(16) 64 0.831058 0.0108156 D(17) 65 -0.0500422 0.0134155 D(18) 66 -0.352216 0.0138548 D(19) 67 0.105948 0.0122153 D(20) 68 0.795741 0.0109979 D(21) 69 -1.40177 0.0159725 D(22) 70 0.174532 0.0116111 D(23) 71 1.09500 0.0106132 D(24) 72 -0.487067 0.0128174 D(25) 73 0.943409 0.0107766 D(26) 74 0.471014 0.0114732 D(27) 75 1.04403 0.0106989 D(28) 76 2.53046 0.0104454 D(29) 77 0.198181 0.0114386 D(30) 78 0.739090 0.0109847 D(31) 79 0.506607 0.0114703 D(32) 80 0.245549 0.0121189 D(33) 81 0.435726 0.0115025 D(34) 82 0.636226 0.0114686 D(35) 83 0.847940 0.0109382 D(36) 84 2.46575 0.00986630 D(37) 85 0.0960517 0.0141088 D(38) 86 0.483612 0.0126373 D(39) 87 -0.0459008 0.0136881 D(40) 88 1.56542 0.0105260 D(41) 89 0.442504 0.0123881 D(42) 90 1.03142 0.0123910 D(43) 91 0.926576 0.0135969 D(44) 92 0.193177 0.0127788 D(45) 93 1.73672 0.0109170 D(46) 94 1.03334 0.0118787 D(47) 95 2.30226 0.0127636 D(48) 96 0.00662815 0.00769033 REG(2) 97 0.0271154 0.00557933 REG(3) 98 0.233819 0.00466816 REG(4) 99 0.242012 0.00868075 REG(5) 100 0.888388 0.00923248 SUBR(2) 101 0.0231399 0.00892070 SUBR(3) 102 -0.288024 0.00712268 SUBR(4) 103 0.905310 0.00745232 SUBR(5) 104 0.246064 0.00596702 SUBR(6) 105 0.371046 0.0101744 SUBR(7) 106 0.185172 0.00805296 SUBR(8) 107 0.335135 0.0213051 SUBR(9) 108 -0.370561 0.0115279 SUBR(10) 109 0.540699 0.00860009 SUBR(11) 110 1.05733 0.00693042 A1 111 0.716404 0.0102341 A2 112 1.18551 0.0105238 B1 113 0.509188 0.0137426 B2 114 -0.551608 0.0130616 F 115 1.05884 0.00909767 L1 116 0.575649 0.00870186 L2 117 -0.482051 0.00806132 M1 118 -0.827379 0.00251797 LDIS 119 0.781755 0.00495824 CONB 120 0.619411 0.00442742 CONR 121 0.470091 0.00435462 CONS scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 176745. LOGLIKELIHOOD RATIO 167998. Degrees of Freedom 2135.00 INDEX OF DISSIMILARITY 11.9497 BIC Statistic 137312. !delete reg and subr $FIT O+D+A1+A2+B1+B2+F+L1+L2+M1+LDIS+CONB+CONR+CONS $ scaled deviance = 208540.8 at cycle 4 d.f. = 2149 from 2256 observations $dis E estimate s.e. parameter 1 10.5376 0.0185536 1 2 -0.219219 0.0143387 O(2) 3 -0.810556 0.0171518 O(3) 4 1.31741 0.0113415 O(4) 5 -0.632325 0.0153862 O(5) 6 0.754186 0.0119026 O(6) 7 2.29236 0.0107329 O(7) 8 1.47757 0.0111238 O(8) 9 1.51422 0.0110730 O(9) 10 1.78137 0.0110232 O(10) 11 0.968373 0.0115658 O(11) 12 1.97833 0.0109318 O(12) 13 1.55906 0.0112069 O(13) 14 0.830368 0.0119212 O(14) 15 1.07137 0.0118710 O(15) 16 0.641288 0.0121255 O(16) 17 1.06758 0.0115426 O(17) 18 -0.0936981 0.0149286 O(18) 19 -0.253723 0.0147974 O(19) 20 0.249039 0.0129629 O(20) 21 0.660096 0.0121751 O(21) 22 -0.684563 0.0156958 O(22) 23 1.02954 0.0116075 O(23) 24 1.47323 0.0113065 O(24) 25 -0.100220 0.0135622 O(25) 26 1.08886 0.0116475 O(26) 27 0.561889 0.0124239 O(27) 28 1.03893 0.0116908 O(28) 29 2.37721 0.0110906 O(29) 30 0.498552 0.0122402 O(30) 31 0.764197 0.0119349 O(31) 32 0.541333 0.0123876 O(32) 33 0.385909 0.0127863 O(33) 34 0.220549 0.0128249 O(34) 35 0.656129 0.0123202 O(35) 36 0.469194 0.0123659 O(36) 37 2.13846 0.0110774 O(37) 38 0.185617 0.0146276 O(38) 39 0.129305 0.0143359 O(39) 40 -0.436492 0.0161593 O(40) 41 1.33811 0.0117126 O(41) 42 0.306008 0.0133484 O(42) 43 1.00781 0.0127555 O(43) 44 0.332140 0.0146152 O(44) 45 -0.210531 0.0147850 O(45) 46 1.40060 0.0121649 O(46) 47 0.701874 0.0128483 O(47) 48 3.11434 0.0115735 O(48) 49 -0.230720 0.0115576 D(2) 50 -0.700666 0.0131409 D(3) 51 0.304107 0.0104450 D(4) 52 -1.53694 0.0156091 D(5) 53 -0.267994 0.0110769 D(6) 54 0.736958 0.00975696 D(7) 55 0.334358 0.00999334 D(8) 56 0.709499 0.00963737 D(9) 57 0.994804 0.00964600 D(10) 58 0.515342 0.0100028 D(11) 59 1.17809 0.00958799 D(12) 60 0.835775 0.00989201 D(13) 61 0.609207 0.0101661 D(14) 62 0.933997 0.0101821 D(15) 63 0.213570 0.0109412 D(16) 64 0.789972 0.0100174 D(17) 65 -0.0368497 0.0128337 D(18) 66 -0.327279 0.0132164 D(19) 67 0.129117 0.0114663 D(20) 68 0.779340 0.0102296 D(21) 69 -1.51980 0.0152657 D(22) 70 0.0509575 0.0106599 D(23) 71 1.00478 0.00967124 D(24) 72 -0.560517 0.0121239 D(25) 73 0.874389 0.00987841 D(26) 74 0.414327 0.0106216 D(27) 75 0.979195 0.00984259 D(28) 76 2.43405 0.0100145 D(29) 77 0.136600 0.0107204 D(30) 78 0.666278 0.0101167 D(31) 79 0.446324 0.0106129 D(32) 80 0.179758 0.0113377 D(33) 81 0.335674 0.0108039 D(34) 82 0.555755 0.0107544 D(35) 83 0.746866 0.0102446 D(36) 84 2.37565 0.00919573 D(37) 85 -0.0205239 0.0135369 D(38) 86 0.380743 0.0118985 D(39) 87 -0.137009 0.0130352 D(40) 88 1.36023 0.00994851 D(41) 89 0.223149 0.0118433 D(42) 90 0.767340 0.0115559 D(43) 91 0.621545 0.0123239 D(44) 92 -0.0675588 0.0120644 D(45) 93 1.77883 0.0100019 D(46) 94 1.06062 0.0106673 D(47) 95 2.49120 0.0109607 D(48) 96 1.03526 0.00672473 A1 97 0.619064 0.00990225 A2 98 1.36735 0.0100141 B1 99 0.646779 0.0133313 B2 100 -0.402282 0.0122463 F 101 0.756061 0.00761513 L1 102 0.222561 0.00646620 L2 103 -0.311389 0.00754186 M1 104 -0.914383 0.00193884 LDIS 105 0.666458 0.00451148 CONB 106 0.531751 0.00399015 CONR 107 0.628855 0.00317053 CONS scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 225542. LOGLIKELIHOOD RATIO 208541. Degrees of Freedom 2149.00 INDEX OF DISSIMILARITY 13.4280 BIC Statistic 177654. $FIT O+D+REG+SUBR+A1+A2+B1+B2+F+L1+L2+M1+ny1+ny2+LDIS+CONB+CONR+CONS $ scaled deviance = 159278.3 at cycle 4 d.f. = 2133 from 2256 observations $dis E estimate s.e. parameter 1 9.76387 0.0237295 1 2 -0.252044 0.0143984 O(2) 3 -0.814244 0.0171642 O(3) 4 1.31777 0.0113839 O(4) 5 -0.681901 0.0154602 O(5) 6 0.784856 0.0119325 O(6) 7 1.86563 0.0131686 O(7) 8 1.61065 0.0114566 O(8) 9 1.62920 0.0113298 O(9) 10 2.01595 0.0113711 O(10) 11 1.27192 0.0119438 O(11) 12 2.22356 0.0113029 O(12) 13 1.82770 0.0115634 O(13) 14 1.08747 0.0122612 O(14) 15 1.03649 0.0123069 O(15) 16 0.569364 0.0126534 O(16) 17 1.06341 0.0119676 O(17) 18 -0.136347 0.0152558 O(18) 19 -0.311924 0.0151907 O(19) 20 0.190462 0.0134197 O(20) 21 0.654449 0.0125825 O(21) 22 -0.652381 0.0161027 O(22) 23 1.07694 0.0121442 O(23) 24 1.49826 0.0118040 O(24) 25 -0.0887341 0.0139555 O(25) 26 1.11454 0.0121497 O(26) 27 0.588206 0.0129227 O(27) 28 1.07646 0.0121598 O(28) 29 2.36347 0.0114732 O(29) 30 0.502124 0.0127052 O(30) 31 0.781413 0.0124923 O(31) 32 0.556670 0.0129646 O(32) 33 0.410830 0.0133355 O(33) 34 0.282709 0.0133183 O(34) 35 0.705245 0.0128724 O(35) 36 0.537937 0.0128574 O(36) 37 2.21146 0.0114777 O(37) 38 0.256573 0.0152498 O(38) 39 0.154245 0.0149920 O(39) 40 -0.415477 0.0166856 O(40) 41 1.50842 0.0120559 O(41) 42 0.487122 0.0137099 O(42) 43 1.20658 0.0133682 O(43) 44 0.584162 0.0155913 O(44) 45 -0.0402370 0.0152244 O(45) 46 1.32140 0.0130064 O(46) 47 0.621495 0.0140530 O(47) 48 2.87090 0.0134311 O(48) 49 -0.257042 0.0116464 D(2) 50 -0.682853 0.0131601 D(3) 51 0.344664 0.0105185 D(4) 52 -1.57439 0.0157249 D(5) 53 -0.182764 0.0111436 D(6) 54 0.337281 0.0130204 D(7) 55 0.553195 0.0108545 D(8) 56 0.870514 0.0103355 D(9) 57 1.28553 0.0104218 D(10) 58 0.904160 0.0108419 D(11) 59 1.47480 0.0103768 D(12) 60 1.17183 0.0106806 D(13) 61 0.943019 0.0109261 D(14) 62 0.941683 0.0109953 D(15) 63 0.186395 0.0118191 D(16) 64 0.819966 0.0108367 D(17) 65 -0.0478080 0.0134323 D(18) 66 -0.352075 0.0138719 D(19) 67 0.104280 0.0122343 D(20) 68 0.792628 0.0110190 D(21) 69 -1.40955 0.0159855 D(22) 70 0.167931 0.0116308 D(23) 71 1.08731 0.0106347 D(24) 72 -0.498654 0.0128351 D(25) 73 0.939216 0.0107984 D(26) 74 0.469026 0.0114958 D(27) 75 1.03781 0.0107206 D(28) 76 2.53569 0.0104642 D(29) 77 0.179202 0.0114629 D(30) 78 0.722547 0.0110081 D(31) 79 0.487988 0.0114981 D(32) 80 0.237091 0.0121434 D(33) 81 0.430340 0.0115220 D(34) 82 0.634961 0.0114902 D(35) 83 0.843972 0.0109603 D(36) 84 2.47067 0.00988852 D(37) 85 0.113055 0.0141220 D(38) 86 0.486004 0.0126608 D(39) 87 -0.0459071 0.0137078 D(40) 88 1.57066 0.0105441 D(41) 89 0.445759 0.0124053 D(42) 90 1.01284 0.0124240 D(43) 91 0.936142 0.0136253 D(44) 92 0.164828 0.0128241 D(45) 93 1.75011 0.0109345 D(46) 94 1.04725 0.0118929 D(47) 95 2.31085 0.0127954 D(48) 96 0.373570 0.00835593 REG(2) 97 0.00200745 0.00560695 REG(3) 98 0.218096 0.00468637 REG(4) 99 0.195939 0.00872913 REG(5) 100 0.490288 0.00994926 SUBR(2) 101 0.254731 0.00920012 SUBR(3) 102 -0.269203 0.00714025 SUBR(4) 103 0.933088 0.00747204 SUBR(5) 104 0.254706 0.00597939 SUBR(6) 105 0.404859 0.0101933 SUBR(7) 106 0.196560 0.00806294 SUBR(8) 107 0.364614 0.0213182 SUBR(9) 108 -0.322311 0.0115656 SUBR(10) 109 0.559701 0.00862423 SUBR(11) 110 1.05186 0.00693148 A1 111 0.708328 0.0102334 A2 112 1.21344 0.0105556 B1 113 0.527310 0.0137617 B2 114 -0.552837 0.0130727 F 115 1.06151 0.00912644 L1 116 0.574406 0.00872780 L2 117 -0.482428 0.00806444 M1 118 0.700145 0.00863125 NY1 119 0.860712 0.0113005 NY2 120 -0.841382 0.00252230 LDIS 121 0.800788 0.00496462 CONB 122 0.686886 0.00446014 CONR 123 0.460784 0.00435572 CONS scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 165925. LOGLIKELIHOOD RATIO 159278. Degrees of Freedom 2133.00 INDEX OF DISSIMILARITY 11.5992 BIC Statistic 128622. $FIT O+D+REG+SUBR+A1+A2+B1+B2+F+L1+L2+M1+ny1+ny2+LDIS+CONB+CONR+CONS +cnt5+grprat+NtoS+NtoW $ scaled deviance = 137067.0 at cycle 3 d.f. = 2126 from 2256 observations $dis E estimate s.e. parameter 1 8.75456 0.0297305 1 2 -0.131608 0.0145539 O(2) 3 -0.748321 0.0172410 O(3) 4 1.35901 0.0114243 O(4) 5 -0.626377 0.0154530 O(5) 6 0.853524 0.0119885 O(6) 7 1.88712 0.0131754 O(7) 8 1.66314 0.0114960 O(8) 9 1.58376 0.0113644 O(9) 10 1.86172 0.0115357 O(10) 11 1.10927 0.0121465 O(11) 12 2.06322 0.0114899 O(12) 13 1.35350 0.0142543 O(13) 14 1.03323 0.0123197 O(14) 15 1.22609 0.0124744 O(15) 16 0.670088 0.0128234 O(16) 17 1.12903 0.0121452 O(17) 18 0.0929937 0.0154947 O(18) 19 -0.130658 0.0153716 O(19) 20 0.315984 0.0135884 O(20) 21 0.791047 0.0127409 O(21) 22 -0.675915 0.0165350 O(22) 23 1.18615 0.0122825 O(23) 24 1.66268 0.0119519 O(24) 25 -0.230780 0.0149903 O(25) 26 1.27940 0.0123030 O(26) 27 0.834495 0.0131363 O(27) 28 1.26174 0.0123547 O(28) 29 2.59073 0.0117844 O(29) 30 0.594845 0.0128416 O(30) 31 0.941835 0.0126924 O(31) 32 0.747963 0.0131209 O(32) 33 0.555661 0.0134901 O(33) 34 0.440298 0.0134956 O(34) 35 0.795170 0.0131334 O(35) 36 0.704277 0.0130609 O(36) 37 2.37640 0.0117724 O(37) 38 0.314834 0.0153761 O(38) 39 0.202962 0.0151853 O(39) 40 -0.257687 0.0168690 O(40) 41 1.63089 0.0123009 O(41) 42 0.531503 0.0139044 O(42) 43 1.35131 0.0137488 O(43) 44 0.638319 0.0158308 O(44) 45 0.151504 0.0154114 O(45) 46 1.47875 0.0131912 O(46) 47 0.750322 0.0141779 O(47) 48 2.95194 0.0136883 O(48) 49 -0.390953 0.0119454 D(2) 50 -0.821060 0.0132914 D(3) 51 0.273546 0.0105503 D(4) 52 -1.53786 0.0157205 D(5) 53 -0.231994 0.0111731 D(6) 54 0.329228 0.0130290 D(7) 55 0.494229 0.0108920 D(8) 56 0.822084 0.0103939 D(9) 57 1.19142 0.0105184 D(10) 58 0.816826 0.0109159 D(11) 59 1.36235 0.0104781 D(12) 60 1.12088 0.0107639 D(13) 61 0.790819 0.0110854 D(14) 62 0.816067 0.0111509 D(15) 63 0.102671 0.0119247 D(16) 64 0.697449 0.0109804 D(17) 65 -0.286873 0.0138920 D(18) 66 -0.546696 0.0141064 D(19) 67 -0.0345793 0.0123613 D(20) 68 0.644681 0.0111681 D(21) 69 -1.37671 0.0160216 D(22) 70 0.123428 0.0116902 D(23) 71 0.960302 0.0108080 D(24) 72 -0.498885 0.0129375 D(25) 73 0.601430 0.0112847 D(26) 74 0.105646 0.0120324 D(27) 75 0.672730 0.0112197 D(28) 76 2.01750 0.0115869 D(29) 77 0.132278 0.0115636 D(30) 78 0.547966 0.0113091 D(31) 79 0.383398 0.0116624 D(32) 80 0.121275 0.0122949 D(33) 81 0.302385 0.0117088 D(34) 82 0.579296 0.0116096 D(35) 83 0.685426 0.0112018 D(36) 84 2.32349 0.0101755 D(37) 85 -0.0787563 0.0143277 D(38) 86 0.242818 0.0130033 D(39) 87 -0.360684 0.0142587 D(40) 88 1.15067 0.0113740 D(41) 89 0.103384 0.0128285 D(42) 90 0.529275 0.0135905 D(43) 91 0.545867 0.0141194 D(44) 92 -0.302994 0.0139127 D(45) 93 1.33381 0.0116604 D(46) 94 0.682184 0.0123827 D(47) 95 1.94383 0.0132904 D(48) 96 0.361768 0.00849506 REG(2) 97 0.0486315 0.00570881 REG(3) 98 0.129843 0.00504644 REG(4) 99 0.327750 0.00894998 REG(5) 100 0.514511 0.00997708 SUBR(2) 101 0.257786 0.00919837 SUBR(3) 102 -0.113659 0.00734225 SUBR(4) 103 0.766787 0.00756757 SUBR(5) 104 0.325526 0.00616883 SUBR(6) 105 0.361874 0.0102505 SUBR(7) 106 0.0961740 0.00819962 SUBR(8) 107 0.228635 0.0212989 SUBR(9) 108 -0.340047 0.0115890 SUBR(10) 109 0.502593 0.00873689 SUBR(11) 110 0.865811 0.00767311 A1 111 0.660925 0.0104060 A2 112 0.995620 0.0110196 B1 113 0.510291 0.0138261 B2 114 -0.406216 0.0131569 F 115 1.08598 0.00920758 L1 116 0.641344 0.00878717 L2 117 -0.268847 0.00835688 M1 118 0.665538 0.00862672 NY1 119 0.831757 0.0113034 NY2 120 -0.721922 0.00357172 LDIS 121 1.12575 0.00773591 CONB 122 1.00378 0.00759005 CONR 123 0.822111 0.00785853 CONS 124 0.00000 aliased CNT5(2) 125 0.313433 0.00541947 CNT5(3) 126 0.140951 0.00435387 CNT5(4) 127 0.131536 0.00398089 CNT5(5) 128 0.0534446 0.00408779 CNT5(6) 129 0.111254 0.00219583 GRPR 130 0.610759 0.00623665 NTOS 131 0.409726 0.00558071 NTOW scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 141289. LOGLIKELIHOOD RATIO 137067. Degrees of Freedom 2126.00 INDEX OF DISSIMILARITY 10.8709 BIC Statistic 106511. $fit -subr $ scaled deviance = 156762.7 (change = +19695.7) at cycle 4 d.f. = 2136 (change = +10 ) from 2256 observations $dis E $ estimate s.e. parameter 1 9.12595 0.0278628 1 2 -0.0959029 0.0145256 O(2) 3 -0.718075 0.0172467 O(3) 4 1.36278 0.0114093 O(4) 5 -0.628102 0.0154504 O(5) 6 0.847876 0.0119796 O(6) 7 1.75451 0.0120513 O(7) 8 1.59159 0.0111876 O(8) 9 1.49489 0.0111352 O(9) 10 1.70960 0.0113991 O(10) 11 0.909735 0.0119830 O(11) 12 1.90015 0.0113459 O(12) 13 1.05021 0.0139367 O(13) 14 0.894327 0.0122150 O(14) 15 1.34903 0.0122040 O(15) 16 0.790128 0.0125324 O(16) 17 1.19501 0.0119380 O(17) 18 0.246267 0.0152890 O(18) 19 0.0218669 0.0151310 O(19) 20 0.440385 0.0133275 O(20) 21 0.880049 0.0125243 O(21) 22 -0.746479 0.0164392 O(22) 23 1.16389 0.0121407 O(23) 24 1.67942 0.0118029 O(24) 25 -0.329680 0.0148614 O(25) 26 1.29397 0.0121403 O(26) 27 0.872417 0.0129568 O(27) 28 1.27534 0.0121999 O(28) 29 2.65094 0.0116673 O(29) 30 0.566263 0.0126491 O(30) 31 0.929146 0.0124310 O(31) 32 0.724984 0.0128523 O(32) 33 0.534353 0.0132525 O(33) 34 0.364184 0.0132956 O(34) 35 0.687428 0.0128900 O(35) 36 0.634918 0.0128392 O(36) 37 2.33009 0.0116176 O(37) 38 0.303506 0.0150779 O(38) 39 0.228727 0.0149139 O(39) 40 -0.188430 0.0166230 O(40) 41 1.57417 0.0122368 O(41) 42 0.436304 0.0138249 O(42) 43 1.27504 0.0136083 O(43) 44 0.466324 0.0155929 O(44) 45 0.100978 0.0152999 O(45) 46 1.64901 0.0127989 O(46) 47 0.911313 0.0135601 O(47) 48 3.22011 0.0129467 O(48) 49 -0.466000 0.0119088 D(2) 50 -0.872120 0.0132849 D(3) 51 0.219349 0.0105140 D(4) 52 -1.56799 0.0157033 D(5) 53 -0.306894 0.0111368 D(6) 54 0.177470 0.0114753 D(7) 55 0.363342 0.0100905 D(8) 56 0.727665 0.00974353 D(9) 57 1.05246 0.0102693 D(10) 58 0.613773 0.0106427 D(11) 59 1.20664 0.0102197 D(12) 60 0.960463 0.0105276 D(13) 61 0.587350 0.0108215 D(14) 62 0.835591 0.0108407 D(15) 63 0.166067 0.0115708 D(16) 64 0.715562 0.0106744 D(17) 65 -0.292505 0.0137207 D(18) 66 -0.514792 0.0138560 D(19) 67 0.0160464 0.0120600 D(20) 68 0.663203 0.0108775 D(21) 69 -1.41373 0.0158140 D(22) 70 0.0631242 0.0113897 D(23) 71 0.902800 0.0104811 D(24) 72 -0.506014 0.0126790 D(25) 73 0.546269 0.0108999 D(26) 74 0.0408934 0.0116500 D(27) 75 0.601519 0.0108450 D(28) 76 1.91611 0.0113300 D(29) 77 0.0764995 0.0112602 D(30) 78 0.450218 0.0109198 D(31) 79 0.304913 0.0112954 D(32) 80 0.0438404 0.0119600 D(33) 81 0.167250 0.0114664 D(34) 82 0.471060 0.0113686 D(35) 83 0.545969 0.0109329 D(36) 84 2.20850 0.00993302 D(37) 85 -0.175396 0.0140690 D(38) 86 0.158395 0.0126552 D(39) 87 -0.464855 0.0139135 D(40) 88 0.934435 0.0110746 D(41) 89 -0.0944048 0.0126082 D(42) 90 0.250053 0.0131124 D(43) 91 0.219217 0.0135513 D(44) 92 -0.578184 0.0134512 D(45) 93 1.32669 0.0112003 D(46) 94 0.689191 0.0117518 D(47) 95 2.06888 0.0125089 D(48) 96 0.464230 0.00732226 REG(2) 97 0.0975239 0.00467803 REG(3) 98 0.157564 0.00466573 REG(4) 99 0.244469 0.00799004 REG(5) 100 0.851091 0.00760954 A1 101 0.611313 0.0103802 A2 102 1.04106 0.0107409 B1 103 0.601647 0.0136937 B2 104 -0.195317 0.0127991 F 105 0.868731 0.00850540 L1 106 0.395081 0.00781163 L2 107 -0.193276 0.00824151 M1 108 0.738466 0.00757341 NY1 109 0.910376 0.0104320 NY2 110 -0.768694 0.00338299 LDIS 111 1.11370 0.00758324 CONB 112 0.911683 0.00729291 CONR 113 0.987520 0.00718835 CONS 114 0.00000 aliased CNT5(2) 115 0.336190 0.00529052 CNT5(3) 116 0.159186 0.00430447 CNT5(4) 117 0.142044 0.00393317 CNT5(5) 118 0.0551515 0.00407791 CNT5(6) 119 0.155558 0.00210492 GRPR 120 0.624358 0.00598295 NTOS 121 0.502342 0.00546874 NTOW scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 166076. LOGLIKELIHOOD RATIO 156763. Degrees of Freedom 2136.00 INDEX OF DISSIMILARITY 11.5524 BIC Statistic 126063. $fit -Reg $ scaled deviance = 162672.5 (change = +5909.8) at cycle 4 d.f. = 2140 (change = +4 ) from 2256 observations $dis E $ estimate s.e. parameter 1 9.73284 0.0255411 1 2 -0.0858978 0.0145017 O(2) 3 -0.710598 0.0172369 O(3) 4 1.35895 0.0113967 O(4) 5 -0.637729 0.0154406 O(5) 6 0.843515 0.0119692 O(6) 7 1.84831 0.0118958 O(7) 8 1.56622 0.0111785 O(8) 9 1.44265 0.0111120 O(9) 10 1.58699 0.0111266 O(10) 11 0.778966 0.0116696 O(11) 12 1.77766 0.0110469 O(12) 13 0.890969 0.0134210 O(13) 14 0.790375 0.0119970 O(14) 15 1.25584 0.0120144 O(15) 16 0.668789 0.0122696 O(16) 17 1.07934 0.0117166 O(17) 18 0.160654 0.0151676 O(18) 19 -0.0824674 0.0149795 O(19) 20 0.324876 0.0131213 O(20) 21 0.763549 0.0123335 O(21) 22 -0.887456 0.0160345 O(22) 23 1.04074 0.0116912 O(23) 24 1.58447 0.0114191 O(24) 25 -0.449901 0.0144744 O(25) 26 1.21080 0.0117607 O(26) 27 0.797524 0.0125821 O(27) 28 1.19242 0.0118535 O(28) 29 2.63046 0.0113664 O(29) 30 0.492409 0.0123510 O(30) 31 0.847759 0.0120976 O(31) 32 0.640700 0.0125025 O(32) 33 0.453043 0.0128957 O(33) 34 0.282475 0.0129698 O(34) 35 0.598481 0.0125171 O(35) 36 0.557118 0.0125332 O(36) 37 2.25141 0.0113347 O(37) 38 0.270647 0.0147816 O(38) 39 0.194776 0.0145168 O(39) 40 -0.224949 0.0163748 O(40) 41 1.52810 0.0120116 O(41) 42 0.396766 0.0135727 O(42) 43 1.24260 0.0132077 O(43) 44 0.463783 0.0149287 O(44) 45 0.0501823 0.0150515 O(45) 46 1.62922 0.0124379 O(46) 47 0.869458 0.0131017 O(47) 48 3.25684 0.0119601 O(48) 49 -0.498525 0.0119013 D(2) 50 -0.897510 0.0132757 D(3) 51 0.189681 0.0104972 D(4) 52 -1.59270 0.0156854 D(5) 53 -0.343249 0.0111255 D(6) 54 0.242513 0.0113249 D(7) 55 0.295443 0.0100607 D(8) 56 0.656540 0.00970217 D(9) 57 0.878902 0.00980737 D(10) 58 0.435275 0.0101175 D(11) 59 1.04127 0.00971954 D(12) 60 0.795878 0.0100654 D(13) 61 0.420973 0.0102712 D(14) 62 0.685957 0.0103306 D(15) 63 0.00263084 0.0110880 D(16) 64 0.557535 0.0102126 D(17) 65 -0.445801 0.0132555 D(18) 66 -0.675858 0.0134349 D(19) 67 -0.142491 0.0116303 D(20) 68 0.501370 0.0104103 D(21) 69 -1.59923 0.0153465 D(22) 70 -0.126736 0.0107393 D(23) 71 0.733765 0.00983969 D(24) 72 -0.644378 0.0122575 D(25) 73 0.410563 0.0103163 D(26) 74 -0.0909374 0.0110555 D(27) 75 0.469224 0.0102832 D(28) 76 1.80541 0.0110031 D(29) 77 -0.0450440 0.0108418 D(30) 78 0.304075 0.0103989 D(31) 79 0.170311 0.0107865 D(32) 80 -0.0862615 0.0114880 D(33) 81 0.0352994 0.0110075 D(34) 82 0.353075 0.0109145 D(35) 83 0.415979 0.0104915 D(36) 84 2.08110 0.00948721 D(37) 85 -0.256492 0.0136809 D(38) 86 0.0724047 0.0121068 D(39) 87 -0.576499 0.0134273 D(40) 88 0.813140 0.0105384 D(41) 89 -0.192221 0.0121209 D(42) 90 0.123949 0.0124005 D(43) 91 0.144104 0.0125859 D(44) 92 -0.732885 0.0128782 D(45) 93 1.23287 0.0105575 D(46) 94 0.583290 0.0110560 D(47) 95 2.03896 0.0113380 D(48) 96 0.840659 0.00750299 A1 97 0.506335 0.0101633 A2 98 1.01622 0.0106446 B1 99 0.575145 0.0133875 B2 100 -0.266328 0.0122805 F 101 0.794147 0.00782767 L1 102 0.303429 0.00662842 L2 103 -0.148967 0.00808170 M1 104 0.560353 0.00686664 NY1 105 0.694440 0.00965008 NY2 106 -0.828184 0.00315307 LDIS 107 1.04221 0.00740520 CONB 108 0.985132 0.00702733 CONR 109 1.07981 0.00687190 CONS 110 0.00000 aliased CNT5(2) 111 0.394476 0.00515476 CNT5(3) 112 0.193642 0.00424451 CNT5(4) 113 0.176555 0.00384181 CNT5(5) 114 0.0678323 0.00403811 CNT5(6) 115 0.170244 0.00205109 GRPR 116 0.564161 0.00584305 NTOS 117 0.499845 0.00515167 NTOW scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 172240. LOGLIKELIHOOD RATIO 162673. Degrees of Freedom 2140.00 INDEX OF DISSIMILARITY 11.8145 BIC Statistic 131915. $c **********************90 data***************************** $c **********************90 data***************************** $YVAR FR90 $ -- model changed $macro wts90! $PRINT 'PEARSONS CHISQUARE ' %X2! $PRINT 'LOGLIKELIHOOD RATIO ' %DV! $PRINT 'Degrees of Freedom ' %df! $CAL %U=%CU(%SQRT((%FV-fr90)**2)*wgt ) $C AL %s=%CU(%FV*wgt ) $CAL %U=50*(%U/%m) $PRINT 'INDEX OF DISSIMILARITY ' %U! $calc %b= %DV - %df*%log(%m) $PRINT 'BIC Statistic' %b ! $endmac ! $c *********Fitting three final models ********************** ! grusk's final model $FIT O+D+REG+SUBR+A1+A2+B1+B2+F+L1+L2+M1+LDIS+CONB+CONR+CONS $ scaled deviance = 99703.6 at cycle 4 d.f. = 2135 from 2256 observations $dis E estimate s.e. parameter 1 8.68059 0.0308942 1 2 -0.0846254 0.0190670 O(2) 3 -0.827876 0.0236632 O(3) 4 1.37811 0.0155879 O(4) 5 -0.611203 0.0207285 O(5) 6 0.819150 0.0162773 O(6) 7 2.46107 0.0150719 O(7) 8 1.71827 0.0155786 O(8) 9 1.66756 0.0155220 O(9) 10 1.95633 0.0157307 O(10) 11 1.34587 0.0164454 O(11) 12 2.32070 0.0155219 O(12) 13 1.86411 0.0159405 O(13) 14 1.27301 0.0167535 O(14) 15 1.19107 0.0169709 O(15) 16 0.831986 0.0172746 O(16) 17 1.23349 0.0164801 O(17) 18 0.226110 0.0202097 O(18) 19 -0.114916 0.0211599 O(19) 20 0.455132 0.0183782 O(20) 21 0.947884 0.0171152 O(21) 22 -0.771650 0.0229700 O(22) 23 1.07389 0.0167995 O(23) 24 1.63491 0.0161874 O(24) 25 0.338573 0.0181033 O(25) 26 1.29747 0.0165655 O(26) 27 0.805560 0.0173930 O(27) 28 1.36802 0.0164653 O(28) 29 2.58153 0.0157514 O(29) 30 0.730392 0.0172730 O(30) 31 1.01151 0.0169569 O(31) 32 0.803537 0.0174333 O(32) 33 0.629428 0.0181156 O(33) 34 0.637572 0.0179280 O(34) 35 1.56374 0.0164528 O(35) 36 1.30230 0.0166541 O(36) 37 3.02780 0.0153605 O(37) 38 0.614761 0.0196854 O(38) 39 0.566621 0.0192298 O(39) 40 0.265490 0.0201101 O(40) 41 1.90770 0.0162250 O(41) 42 0.838918 0.0181170 O(42) 43 1.47638 0.0176525 O(43) 44 1.10084 0.0195178 O(44) 45 0.285593 0.0195338 O(45) 46 1.62474 0.0173950 O(46) 47 0.843680 0.0187649 O(47) 48 2.80448 0.0181282 O(48) 49 -0.0614781 0.0162114 D(2) 50 -0.896244 0.0203949 D(3) 51 0.840344 0.0141658 D(4) 52 -0.859086 0.0186930 D(5) 53 0.233305 0.0150041 D(6) 54 1.36700 0.0140181 D(7) 55 0.969859 0.0145138 D(8) 56 1.16981 0.0140979 D(9) 57 1.32603 0.0143687 D(10) 58 0.891111 0.0149866 D(11) 59 1.46742 0.0143253 D(12) 60 1.16925 0.0147508 D(13) 61 0.666238 0.0156040 D(14) 62 0.711527 0.0155807 D(15) 63 -0.0827117 0.0170892 D(16) 64 0.791653 0.0149838 D(17) 65 -0.877252 0.0231767 D(18) 66 -0.866557 0.0217174 D(19) 67 -0.247872 0.0179224 D(20) 68 0.397841 0.0158947 D(21) 69 -0.902431 0.0199834 D(22) 70 0.812555 0.0150478 D(23) 71 1.55863 0.0142538 D(24) 72 -0.545554 0.0184794 D(25) 73 1.46604 0.0143400 D(26) 74 0.855974 0.0151769 D(27) 75 1.53551 0.0142885 D(28) 76 2.89025 0.0140818 D(29) 77 0.244082 0.0158973 D(30) 78 0.943689 0.0149319 D(31) 79 0.563148 0.0157256 D(32) 80 0.117666 0.0169242 D(33) 81 0.274934 0.0162368 D(34) 82 0.319987 0.0166012 D(35) 83 0.445562 0.0159178 D(36) 84 2.33043 0.0138374 D(37) 85 -0.236572 0.0206694 D(38) 86 0.0328609 0.0186758 D(39) 87 -0.886411 0.0230170 D(40) 88 1.34145 0.0147460 D(41) 89 0.350269 0.0168171 D(42) 90 1.41741 0.0156975 D(43) 91 0.556536 0.0188061 D(44) 92 0.702714 0.0160436 D(45) 93 1.73404 0.0149522 D(46) 94 0.885483 0.0163729 D(47) 95 2.52432 0.0165817 D(48) 96 -0.00725737 0.00992681 REG(2) 97 0.0806835 0.00806448 REG(3) 98 0.200095 0.00596498 REG(4) 99 0.447482 0.0110732 REG(5) 100 0.872553 0.0121537 SUBR(2) 101 0.238163 0.0114531 SUBR(3) 102 -0.191887 0.0101341 SUBR(4) 103 0.967812 0.0109085 SUBR(5) 104 0.207513 0.00723239 SUBR(6) 105 0.533159 0.0131308 SUBR(7) 106 0.135930 0.0101692 SUBR(8) 107 0.436123 0.0294413 SUBR(9) 108 -0.252850 0.0137534 SUBR(10) 109 0.584675 0.0110381 SUBR(11) 110 1.14717 0.00846004 A1 111 0.705013 0.0123089 A2 112 0.860131 0.0129823 B1 113 0.767270 0.0164520 B2 114 -0.419314 0.0166362 F 115 0.913448 0.0109373 L1 116 0.844724 0.0112876 L2 117 -0.380347 0.00897809 M1 118 -0.804537 0.00316665 LDIS 119 0.726360 0.00674399 CONB 120 0.581311 0.00575117 CONR 121 0.396462 0.00565037 CONS scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 105115. LOGLIKELIHOOD RATIO 99703.6 Degrees of Freedom 2135.00 INDEX OF DISSIMILARITY 22.5500 BIC Statistic 69018.4 $FIT O+D+REG+SUBR+A1+A2+B1+B2+F+L1+L2+M1+ny1+ny2+LDIS+CONB+CONR+CONS $ scaled deviance = 95190.02 at cycle 4 d.f. = 2133 from 2256 observations $dis E estimate s.e. parameter 1 8.75926 0.0309533 1 2 -0.0855577 0.0190662 O(2) 3 -0.819661 0.0236629 O(3) 4 1.38457 0.0155889 O(4) 5 -0.615409 0.0207301 O(5) 6 0.828859 0.0162796 O(6) 7 1.91670 0.0177879 O(7) 8 1.64681 0.0156710 O(8) 9 1.60715 0.0155676 O(9) 10 1.95132 0.0157375 O(10) 11 1.34021 0.0164528 O(11) 12 2.31140 0.0155301 O(12) 13 1.86004 0.0159479 O(13) 14 1.26551 0.0167609 O(14) 15 1.18612 0.0169790 O(15) 16 0.820506 0.0172834 O(16) 17 1.22222 0.0164887 O(17) 18 0.223649 0.0202165 O(18) 19 -0.119012 0.0211669 O(19) 20 0.449325 0.0183861 O(20) 21 0.940630 0.0171237 O(21) 22 -0.774859 0.0229734 O(22) 23 1.07119 0.0168044 O(23) 24 1.62683 0.0161935 O(24) 25 0.330148 0.0181097 O(25) 26 1.29104 0.0165725 O(26) 27 0.800472 0.0174020 O(27) 28 1.35841 0.0164727 O(28) 29 2.58221 0.0157582 O(29) 30 0.712154 0.0172858 O(30) 31 0.988618 0.0169696 O(31) 32 0.777556 0.0174516 O(32) 33 0.618512 0.0181285 O(33) 34 0.629677 0.0179350 O(34) 35 1.55999 0.0164622 O(35) 36 1.29570 0.0166643 O(36) 37 3.02853 0.0153687 O(37) 38 0.628649 0.0196919 O(38) 39 0.558856 0.0192493 O(39) 40 0.262859 0.0201233 O(40) 41 1.91029 0.0162310 O(41) 42 0.837834 0.0181256 O(42) 43 1.46236 0.0176731 O(43) 44 1.10683 0.0195358 O(44) 45 0.259962 0.0195641 O(45) 46 1.63036 0.0174073 O(46) 47 0.846512 0.0187792 O(47) 48 2.79995 0.0181471 O(48) 49 -0.0642297 0.0162115 D(2) 50 -0.891304 0.0203948 D(3) 51 0.842257 0.0141684 D(4) 52 -0.866726 0.0186963 D(5) 53 0.234465 0.0150138 D(6) 54 0.678875 0.0175609 D(7) 55 0.868050 0.0146784 D(8) 56 1.07821 0.0142180 D(9) 57 1.32931 0.0143949 D(10) 58 0.895528 0.0150132 D(11) 59 1.46722 0.0143540 D(12) 60 1.17523 0.0147772 D(13) 61 0.668052 0.0156295 D(14) 62 0.716158 0.0156074 D(15) 63 -0.0868344 0.0171157 D(16) 64 0.788588 0.0150114 D(17) 65 -0.869617 0.0231942 D(18) 66 -0.861069 0.0217372 D(19) 67 -0.244277 0.0179458 D(20) 68 0.399434 0.0159211 D(21) 69 -0.901341 0.0200030 D(22) 70 0.815631 0.0150747 D(23) 71 1.56068 0.0142812 D(24) 72 -0.546843 0.0185014 D(25) 73 1.47075 0.0143676 D(26) 74 0.862417 0.0152048 D(27) 75 1.53805 0.0143167 D(28) 76 2.90110 0.0141060 D(29) 77 0.235512 0.0159280 D(30) 78 0.936291 0.0149627 D(31) 79 0.553708 0.0157609 D(32) 80 0.113389 0.0169558 D(33) 81 0.277095 0.0162605 D(34) 82 0.325144 0.0166269 D(35) 83 0.447448 0.0159448 D(36) 84 2.34047 0.0138649 D(37) 85 -0.213344 0.0206858 D(38) 86 0.0404705 0.0187046 D(39) 87 -0.882454 0.0230393 D(40) 88 1.35172 0.0147691 D(41) 89 0.356971 0.0168411 D(42) 90 1.41619 0.0157306 D(43) 91 0.571667 0.0188364 D(44) 92 0.692874 0.0160886 D(45) 93 1.75325 0.0149730 D(46) 94 0.904046 0.0163927 D(47) 95 2.53834 0.0166155 D(48) 96 0.331730 0.0108708 REG(2) 97 0.0574565 0.00809681 REG(3) 98 0.188592 0.00598244 REG(4) 99 0.407510 0.0111269 REG(5) 100 0.510388 0.0130898 SUBR(2) 101 0.453992 0.0118562 SUBR(3) 102 -0.174477 0.0101549 SUBR(4) 103 0.993223 0.0109311 SUBR(5) 104 0.213497 0.00724334 SUBR(6) 105 0.564683 0.0131563 SUBR(7) 106 0.143526 0.0101780 SUBR(8) 107 0.459546 0.0294536 SUBR(9) 108 -0.212243 0.0137984 SUBR(10) 109 0.607036 0.0110685 SUBR(11) 110 1.14388 0.00846136 A1 111 0.699167 0.0123086 A2 112 0.876009 0.0130045 B1 113 0.785429 0.0164744 B2 114 -0.422232 0.0166477 F 115 0.915193 0.0109651 L1 116 0.855257 0.0113169 L2 117 -0.379290 0.00898177 M1 118 0.594287 0.0114372 NY1 119 0.859208 0.0141280 NY2 120 -0.816510 0.00317122 LDIS 121 0.745341 0.00675250 CONB 122 0.641709 0.00579612 CONR 123 0.388129 0.00565207 CONS scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 99538.0 LOGLIKELIHOOD RATIO 95190.0 Degrees of Freedom 2133.00 INDEX OF DISSIMILARITY 22.5456 BIC Statistic 64533.6 $FIT O+D+REG+SUBR+A1+A2+B1+B2+F+L1+L2+M1+ny1+ny2+LDIS+CONB+CONR+CONS +cnt5+gr90+NtoS+NtoW $ scaled deviance = 85053.09 at cycle 4 d.f. = 2126 from 2256 observations $dis E estimate s.e. parameter 1 8.16532 0.0387030 1 2 -0.0962995 0.0191228 O(2) 3 -0.838220 0.0236737 O(3) 4 1.36214 0.0156365 O(4) 5 -0.609069 0.0208220 O(5) 6 0.820410 0.0163138 O(6) 7 1.89667 0.0180819 O(7) 8 1.63900 0.0157019 O(8) 9 1.57845 0.0157259 O(9) 10 1.85066 0.0160816 O(10) 11 1.25863 0.0166179 O(11) 12 2.21834 0.0158320 O(12) 13 1.74579 0.0168097 O(13) 14 1.18989 0.0169487 O(14) 15 1.33950 0.0172283 O(15) 16 0.941255 0.0176809 O(16) 17 1.30717 0.0170967 O(17) 18 0.258217 0.0241591 O(18) 19 0.00666612 0.0215439 O(19) 20 0.568605 0.0188033 O(20) 21 1.04142 0.0178398 O(21) 22 -0.646332 0.0231497 O(22) 23 1.17731 0.0170080 O(23) 24 1.74831 0.0164141 O(24) 25 0.368172 0.0194680 O(25) 26 1.40763 0.0168036 O(26) 27 0.946148 0.0176424 O(27) 28 1.45807 0.0167991 O(28) 29 2.72725 0.0160784 O(29) 30 0.813057 0.0176308 O(30) 31 1.10235 0.0173175 O(31) 32 0.916216 0.0178395 O(32) 33 0.722093 0.0185461 O(33) 34 0.740966 0.0183764 O(34) 35 1.66547 0.0175664 O(35) 36 1.36803 0.0181968 O(36) 37 3.14243 0.0161125 O(37) 38 0.718382 0.0203912 O(38) 39 0.678613 0.0195292 O(39) 40 -0.174638 0.0591646 O(40) 41 1.99610 0.0169124 O(41) 42 0.896246 0.0189653 O(42) 43 1.56538 0.0181404 O(43) 44 1.15350 0.0201253 O(44) 45 0.424697 0.0197503 O(45) 46 1.76845 0.0176029 O(46) 47 0.981180 0.0189763 O(47) 48 2.91632 0.0184475 O(48) 49 -0.0480988 0.0163110 D(2) 50 -0.918720 0.0204223 D(3) 51 0.836350 0.0142587 D(4) 52 -0.829810 0.0188164 D(5) 53 0.239509 0.0150621 D(6) 54 0.702810 0.0177558 D(7) 55 0.863859 0.0147262 D(8) 56 1.06804 0.0144868 D(9) 57 1.29401 0.0148930 D(10) 58 0.851406 0.0152811 D(11) 59 1.41470 0.0147954 D(12) 60 1.16705 0.0155824 D(13) 61 0.621439 0.0159397 D(14) 62 0.700538 0.0157393 D(15) 63 -0.0990709 0.0174497 D(16) 64 0.756506 0.0156703 D(17) 65 -0.839609 0.0244031 D(18) 66 -0.892010 0.0221054 D(19) 67 -0.283412 0.0183447 D(20) 68 0.370368 0.0166507 D(21) 69 -0.935337 0.0201975 D(22) 70 0.789251 0.0151556 D(23) 71 1.51541 0.0144141 D(24) 72 -0.537558 0.0192436 D(25) 73 1.21274 0.0148287 D(26) 74 0.628037 0.0156012 D(27) 75 1.31146 0.0148068 D(28) 76 2.62020 0.0145759 D(29) 77 0.229015 0.0162783 D(30) 78 0.872565 0.0152977 D(31) 79 0.518602 0.0161259 D(32) 80 0.0622258 0.0173766 D(33) 81 0.244339 0.0167020 D(34) 82 0.326040 0.0174869 D(35) 83 0.434313 0.0170938 D(36) 84 2.31517 0.0146767 D(37) 85 -0.254231 0.0215347 D(38) 86 -0.0808026 0.0190315 D(39) 87 -0.884277 0.0250967 D(40) 88 1.22695 0.0159795 D(41) 89 0.221612 0.0179853 D(42) 90 1.27728 0.0164334 D(43) 91 0.450599 0.0196900 D(44) 92 0.503571 0.0167921 D(45) 93 1.59614 0.0153794 D(46) 94 0.754577 0.0168146 D(47) 95 2.37466 0.0170653 D(48) 96 0.383454 0.0110317 REG(2) 97 0.0951498 0.00823816 REG(3) 98 0.176459 0.00654422 REG(4) 99 0.431098 0.0113750 REG(5) 100 0.522717 0.0131506 SUBR(2) 101 0.448872 0.0118772 SUBR(3) 102 -0.0757464 0.0102922 SUBR(4) 103 0.851532 0.0110768 SUBR(5) 104 0.266992 0.00744558 SUBR(6) 105 0.504271 0.0132520 SUBR(7) 106 0.0352513 0.0104043 SUBR(8) 107 0.371854 0.0294515 SUBR(9) 108 -0.236558 0.0138373 SUBR(10) 109 0.562721 0.0112574 SUBR(11) 110 0.879438 0.00939880 A1 111 0.671848 0.0125164 A2 112 0.808615 0.0133220 B1 113 0.751945 0.0165682 B2 114 -0.332630 0.0167452 F 115 0.921401 0.0110615 L1 116 0.888364 0.0113822 L2 117 -0.224415 0.00920247 M1 118 0.576700 0.0114458 NY1 119 0.845686 0.0141546 NY2 120 -0.757456 0.00451356 LDIS 121 0.934552 0.0101807 CONB 122 0.820651 0.00978704 CONR 123 0.593299 0.0101650 CONS 124 0.00000 aliased CNT5(2) 125 0.190917 0.00701522 CNT5(3) 126 0.0679790 0.00566782 CNT5(4) 127 0.0546472 0.00520573 CNT5(5) 128 0.0276821 0.00537963 CNT5(6) 129 0.0771518 0.00735025 GR90 130 0.659239 0.00725840 NTOS 131 0.369273 0.00773339 NTOW scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 88745.3 LOGLIKELIHOOD RATIO 85053.1 Degrees of Freedom 2126.00 INDEX OF DISSIMILARITY 22.2928 BIC Statistic 54497.3 $fit -subr $ scaled deviance = 97213.5 (change = +12160.4) at cycle 4 d.f. = 2136 (change = +10 ) from 2256 observations $dis E $ estimate s.e. parameter 1 8.54219 0.0366348 1 2 -0.103421 0.0191089 O(2) 3 -0.851204 0.0236656 O(3) 4 1.33664 0.0156175 O(4) 5 -0.616577 0.0208055 O(5) 6 0.779462 0.0163036 O(6) 7 1.84770 0.0166145 O(7) 8 1.57980 0.0152639 O(8) 9 1.53633 0.0154079 O(9) 10 1.75963 0.0159541 O(10) 11 1.12058 0.0164560 O(11) 12 2.12331 0.0156881 O(12) 13 1.65135 0.0167265 O(13) 14 1.08051 0.0168155 O(14) 15 1.45167 0.0169332 O(15) 16 1.08784 0.0173133 O(16) 17 1.40763 0.0168722 O(17) 18 0.453071 0.0242612 O(18) 19 0.152816 0.0212843 O(19) 20 0.711046 0.0185017 O(20) 21 1.16470 0.0176243 O(21) 22 -0.640349 0.0230026 O(22) 23 1.18295 0.0168026 O(23) 24 1.78125 0.0162011 O(24) 25 0.435653 0.0193668 O(25) 26 1.43541 0.0165716 O(26) 27 0.987432 0.0173958 O(27) 28 1.48587 0.0166014 O(28) 29 2.79348 0.0159261 O(29) 30 0.861111 0.0173746 O(30) 31 1.14984 0.0169904 O(31) 32 0.975116 0.0175107 O(32) 33 0.773097 0.0182459 O(33) 34 0.699075 0.0182149 O(34) 35 1.63767 0.0174231 O(35) 36 1.34636 0.0180982 O(36) 37 3.13191 0.0159936 O(37) 38 0.759989 0.0201684 O(38) 39 0.735329 0.0192180 O(39) 40 0.200718 0.0609258 O(40) 41 1.95648 0.0168869 O(41) 42 0.851903 0.0189356 O(42) 43 1.51111 0.0179994 O(43) 44 1.00678 0.0198648 O(44) 45 0.359488 0.0195759 O(45) 46 1.96629 0.0170275 O(46) 47 1.17963 0.0180396 O(47) 48 3.21798 0.0174392 O(48) 49 -0.0813812 0.0162781 D(2) 50 -0.938209 0.0204098 D(3) 51 0.785824 0.0142125 D(4) 52 -0.871238 0.0187711 D(5) 53 0.170298 0.0150162 D(6) 54 0.640360 0.0159643 D(7) 55 0.813589 0.0139089 D(8) 56 1.01921 0.0138518 D(9) 57 1.15998 0.0146844 D(10) 58 0.666499 0.0150294 D(11) 59 1.27934 0.0145803 D(12) 60 1.00397 0.0154147 D(13) 61 0.465897 0.0157376 D(14) 62 0.790395 0.0153427 D(15) 63 0.00531552 0.0170323 D(16) 64 0.803946 0.0153459 D(17) 65 -0.777043 0.0242638 D(18) 66 -0.793295 0.0217975 D(19) 67 -0.190037 0.0179849 D(20) 68 0.422568 0.0163623 D(21) 69 -0.937617 0.0200007 D(22) 70 0.761034 0.0148839 D(23) 71 1.50118 0.0141389 D(24) 72 -0.570501 0.0190968 D(25) 73 1.17341 0.0145219 D(26) 74 0.601542 0.0152866 D(27) 75 1.26496 0.0145103 D(28) 76 2.59609 0.0143816 D(29) 77 0.216995 0.0159454 D(30) 78 0.862884 0.0148554 D(31) 79 0.508516 0.0157051 D(32) 80 0.0460731 0.0170027 D(33) 81 0.122041 0.0164235 D(34) 82 0.194970 0.0172453 D(35) 83 0.297866 0.0168682 D(36) 84 2.21592 0.0144642 D(37) 85 -0.312833 0.0211181 D(38) 86 -0.0837231 0.0185324 D(39) 87 -0.915948 0.0248543 D(40) 88 1.08265 0.0157878 D(41) 89 0.0664506 0.0177958 D(42) 90 1.11806 0.0161161 D(43) 91 0.202046 0.0191795 D(44) 92 0.390410 0.0164973 D(45) 93 1.68772 0.0148900 D(46) 94 0.834319 0.0160855 D(47) 95 2.56821 0.0160771 D(48) 96 0.577389 0.00947528 REG(2) 97 0.165838 0.00645866 REG(3) 98 0.189701 0.00612337 REG(4) 99 0.347589 0.0100781 REG(5) 100 0.851971 0.00931573 A1 101 0.616866 0.0124667 A2 102 0.885284 0.0130125 B1 103 0.836832 0.0163974 B2 104 -0.137882 0.0163370 F 105 0.687915 0.00993540 L1 106 0.624429 0.0100404 L2 107 -0.177101 0.00910230 M1 108 0.594257 0.0100747 NY1 109 0.852797 0.0129126 NY2 110 -0.801215 0.00429025 LDIS 111 0.911375 0.00998392 CONB 112 0.685493 0.00943069 CONR 113 0.765161 0.00919869 CONS 114 0.00000 aliased CNT5(2) 115 0.202635 0.00682196 CNT5(3) 116 0.0822210 0.00558901 CNT5(4) 117 0.0680699 0.00513461 CNT5(5) 118 0.0319356 0.00535697 CNT5(6) 119 0.0348361 0.00762820 GR90 120 0.725457 0.00704034 NTOS 121 0.485606 0.00755748 NTOW scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 103901. LOGLIKELIHOOD RATIO 97213.5 Degrees of Freedom 2136.00 INDEX OF DISSIMILARITY 22.3692 BIC Statistic 66513.9 $fit -Reg $ scaled deviance = 103157.1 (change = +5943.7) at cycle 4 d.f. = 2140 (change = +4 ) from 2256 observations $dis E $ estimate s.e. parameter 1 9.29757 0.0339611 1 2 -0.118007 0.0190962 O(2) 3 -0.860891 0.0236608 O(3) 4 1.31035 0.0156054 O(4) 5 -0.641318 0.0207914 O(5) 6 0.754436 0.0162944 O(6) 7 1.95016 0.0164302 O(7) 8 1.52449 0.0152503 O(8) 9 1.45773 0.0153759 O(9) 10 1.59822 0.0156423 O(10) 11 0.961506 0.0161099 O(11) 12 1.97387 0.0153573 O(12) 13 1.49452 0.0164084 O(13) 14 0.945946 0.0165107 O(14) 15 1.33009 0.0166617 O(15) 16 0.937959 0.0169827 O(16) 17 1.25464 0.0165933 O(17) 18 0.288260 0.0241021 O(18) 19 0.00939153 0.0210583 O(19) 20 0.560948 0.0182315 O(20) 21 1.00186 0.0173512 O(21) 22 -0.809109 0.0225110 O(22) 23 1.00581 0.0161425 O(23) 24 1.62718 0.0155917 O(24) 25 0.280109 0.0188704 O(25) 26 1.29578 0.0159742 O(26) 27 0.850590 0.0167926 O(27) 28 1.33873 0.0160642 O(28) 29 2.73462 0.0154448 O(29) 30 0.746222 0.0168973 O(30) 31 1.01659 0.0164747 O(31) 32 0.837079 0.0169889 O(32) 33 0.645578 0.0177286 O(33) 34 0.575109 0.0177374 O(34) 35 1.50936 0.0168901 O(35) 36 1.21195 0.0176347 O(36) 37 3.01032 0.0155498 O(37) 38 0.712623 0.0197772 O(38) 39 0.700870 0.0186894 O(39) 40 -0.0454366 0.0609869 O(40) 41 1.88042 0.0165040 O(41) 42 0.787912 0.0185542 O(42) 43 1.45018 0.0174397 O(43) 44 0.991190 0.0189887 O(44) 45 0.279471 0.0191759 O(45) 46 1.93410 0.0165149 O(46) 47 1.11918 0.0174111 O(47) 48 3.25921 0.0162806 O(48) 49 -0.0876785 0.0162707 D(2) 50 -0.952757 0.0203994 D(3) 51 0.763196 0.0142009 D(4) 52 -0.888465 0.0187605 D(5) 53 0.141541 0.0150110 D(6) 54 0.745430 0.0157695 D(7) 55 0.744486 0.0138956 D(8) 56 0.946821 0.0138275 D(9) 57 0.986230 0.0142139 D(10) 58 0.481574 0.0144566 D(11) 59 1.11200 0.0140645 D(12) 60 0.843625 0.0149270 D(13) 61 0.305294 0.0152269 D(14) 62 0.645707 0.0149083 D(15) 63 -0.158948 0.0165697 D(16) 64 0.647434 0.0149167 D(17) 65 -0.898012 0.0240326 D(18) 66 -0.948150 0.0214948 D(19) 67 -0.350335 0.0176116 D(20) 68 0.261342 0.0159885 D(21) 69 -1.14499 0.0194405 D(22) 70 0.564963 0.0141829 D(23) 71 1.33505 0.0134937 D(24) 72 -0.706705 0.0186845 D(25) 73 1.03731 0.0138939 D(26) 74 0.479873 0.0146511 D(27) 75 1.14100 0.0139286 D(28) 76 2.49459 0.0140800 D(29) 77 0.0969255 0.0155002 D(30) 78 0.721515 0.0143197 D(31) 79 0.377941 0.0151673 D(32) 80 -0.0857796 0.0164995 D(33) 81 -0.0152176 0.0159504 D(34) 82 0.0736453 0.0167908 D(35) 83 0.174743 0.0164695 D(36) 84 2.08742 0.0140483 D(37) 85 -0.361756 0.0207328 D(38) 86 -0.153299 0.0179767 D(39) 87 -0.967223 0.0245317 D(40) 88 0.993914 0.0153421 D(41) 89 -0.000652459 0.0173395 D(42) 90 1.04027 0.0153970 D(43) 91 0.179935 0.0182260 D(44) 92 0.261465 0.0159440 D(45) 93 1.62760 0.0142712 D(46) 94 0.753528 0.0153738 D(47) 95 2.58988 0.0148084 D(48) 96 0.850163 0.00909970 A1 97 0.508067 0.0122114 A2 98 0.847504 0.0127723 B1 99 0.792012 0.0159734 B2 100 -0.247940 0.0157025 F 101 0.558600 0.00891927 L1 102 0.502152 0.00884210 L2 103 -0.110182 0.00879375 M1 104 0.363685 0.00908427 NY1 105 0.595722 0.0119085 NY2 106 -0.878406 0.00402261 LDIS 107 0.809854 0.00978718 CONB 108 0.779036 0.00910881 CONR 109 0.883601 0.00880025 CONS 110 0.00000 aliased CNT5(2) 111 0.276817 0.00664933 CNT5(3) 112 0.125861 0.00549674 CNT5(4) 113 0.114358 0.00499493 CNT5(5) 114 0.0504366 0.00528661 CNT5(6) 115 0.0614588 0.00762975 GR90 116 0.653712 0.00682714 NTOS 117 0.495358 0.00704605 NTOW scale parameter taken as 1.00000 $use wtsta $ PEARSONS CHISQUARE 110705. LOGLIKELIHOOD RATIO 103157. Degrees of Freedom 2140.00 INDEX OF DISSIMILARITY 22.4099 BIC Statistic 72400.1 $stop