------------------------------------------------------------------------------- name: log: c:\users\bchavez2\desktop\week3_lecture.txt log type: text opened on: 3 Feb 2016, 13:59:27 r; t=0.00 13:59:27 . end of do-file r; t=0.02 13:59:27 . doed r; t=0.14 14:13:28 . use AApredation.dta, clear r; t=0.05 14:15:25 . ed r; t=0.28 14:15:28 . summ year Variable | Obs Mean Std. Dev. Min Max -------------+--------------------------------------------------------- year | 1,760 1995 1.414615 1993 1997 r; t=0.02 14:15:46 . ed r; t=0.36 14:17:19 . encode destination, gen(edest) r; t=0.00 14:19:35 . do "C:\Users\bchavez2\AppData\Local\Temp\STD0o000000.tmp" . twoway (line NJfare period if destin=="MCI", sort lcolor(blue) lwidth(medthic > k)cmissing(n)) /// > (line AAfare period if destin=="MCI", sort lcolor(red) lwidth(medthick)) /// > (line DLfare period if destin=="MCI", sort lcolor(green) lwidth(medthick) lpa > ttern(longdash)cmissing(n)), /// > ytitle(Fare price) xtitle(Time) title(DFW/Kansas City Market(MCI)) legend(siz > e(small)) r; t=8.03 14:19:53 . end of do-file r; t=8.03 14:19:53 . do "C:\Users\bchavez2\AppData\Local\Temp\STD0o000000.tmp" . twoway (line NJfare period if destin=="MCI", sort lcolor(blue) lwidth(medthic > k)cmissing(n)) /// > (line AAfare period if destin=="MCI", sort lcolor(red) lwidth(medthick)) /// > (line DLfare period if destin=="MCI", sort lcolor(green) lwidth(medthick) lpa > ttern(longdash)cmissing(n)), /// > ytitle(Fare price) xtitle(Time) title(DFW/Kansas City Market(MCI)) legend(siz > e(small)) r; t=1.12 14:19:54 . end of do-file r; t=1.12 14:19:54 . do "C:\Users\bchavez2\AppData\Local\Temp\STD0o000000.tmp" . /*Fixed effects: Heterogeneity across entities */ . bysort edest: egen aamean=mean(AAfare) r; t=0.09 14:30:54 . twoway scatter AAfare edest, msymbol(circle_hollow) || connected aamean edest > , r; t=1.11 14:30:55 . msymbol(diamond) || command msymbol is unrecognized r(199); t=0.03 14:30:55 end of do-file r(199); t=1.25 14:30:55 . do "C:\Users\bchavez2\AppData\Local\Temp\STD0o000000.tmp" . /*Fixed effects: Heterogeneity across time */ . bysort period: egen aamean_p=mean(AAfare) r; t=0.02 14:31:53 . twoway scatter AAfare period, msymbol(circle_hollow) || connected aamean_p pe > riod, // r; t=0.62 14:31:53 . msymbol(diamond) || command msymbol is unrecognized r(199); t=0.02 14:31:53 end of do-file r(199); t=0.67 14:31:53 . reg AAfare dist Source | SS df MS Number of obs = 1,729 -------------+---------------------------------- F(1, 1727) = 3490.67 Model | 70621771.9 1 70621771.9 Prob > F = 0.0000 Residual | 34939940.2 1,727 20231.5809 R-squared = 0.6690 -------------+---------------------------------- Adj R-squared = 0.6688 Total | 105561712 1,728 61088.9538 Root MSE = 142.24 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | .2318662 .0039245 59.08 0.000 .2241689 .2395634 _cons | 173.8493 7.18437 24.20 0.000 159.7583 187.9403 ------------------------------------------------------------------------------ r; t=1.17 14:33:16 . lfit AAfare dist too many variables specified r(103); t=0.09 14:35:33 . twoway scatter AAfare dist, mlabel(dest) || lfit AAfare dist, clstyle(p2) r; t=1.69 14:36:06 . reg AAfare dist i.edest Source | SS df MS Number of obs = 1,729 -------------+---------------------------------- F(88, 1640) = 115.68 Model | 90914902.6 88 1033123.89 Prob > F = 0.0000 Residual | 14646809.6 1,640 8930.98145 R-squared = 0.8612 -------------+---------------------------------- Adj R-squared = 0.8538 Total | 105561712 1,728 61088.9538 Root MSE = 94.504 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | 9.134217 5.649705 1.62 0.106 -1.947178 20.21561 | edest | ACT | 8723.954 5423.799 1.61 0.108 -1914.347 19362.26 AMA | 4584.856 2892.803 1.58 0.113 -1089.122 10258.83 ATL | -2819.999 1842.046 -1.53 0.126 -6433.01 793.0112 AUS | 6984.426 4361.674 1.60 0.109 -1570.613 15539.46 BDL | -15892.31 10192.11 -1.56 0.119 -35883.24 4098.611 BHM | -145.3567 317.7917 -0.46 0.647 -768.6771 477.9637 BNA | -833.9304 701.2005 -1.19 0.234 -2209.273 541.4123 BOS | -17467.41 11220.35 -1.56 0.120 -39475.14 4540.317 BPT | 5461.001 3378.656 1.62 0.106 -1165.932 12087.94 BTR | 3415.882 2101.948 1.63 0.104 -706.903 7538.667 BUR | -11689.45 7480.269 -1.56 0.118 -26361.34 2982.431 BWI | -11325.26 7322.078 -1.55 0.122 -25686.87 3036.348 CLE | -7698.048 5107.42 -1.51 0.132 -17715.8 2319.705 CLL | 7381.852 4576.358 1.61 0.107 -1594.27 16357.97 CLT | -6292.523 4158.29 -1.51 0.130 -14448.64 1863.594 CMH | -6099.748 4045.299 -1.51 0.132 -14034.24 1834.748 COS | -280.3255 261.599 -1.07 0.284 -793.4288 232.7778 CRP | 3871.218 2429.557 1.59 0.111 -894.1423 8636.579 CVG | -4059.369 2745.973 -1.48 0.140 -9445.353 1326.615 DAY | -4915.255 3299.563 -1.49 0.137 -11387.06 1556.546 DCA | -10850.15 7039.595 -1.54 0.123 -24657.69 2957.394 DEN | -1188.765 831.0441 -1.43 0.153 -2818.784 441.2545 DSM | -719.5364 622.1856 -1.16 0.248 -1939.898 500.8257 DTW | -7218.988 4723.248 -1.53 0.127 -16483.22 2045.245 ELP | 328.2468 205.5732 1.60 0.111 -74.96678 731.4604 EWR | -14151.52 9073.475 -1.56 0.119 -31948.34 3645.299 FAT | -13168.57 8406.814 -1.57 0.117 -29657.79 3320.652 FLL | -9664.642 6214.747 -1.56 0.120 -21854.32 2525.034 FYV | 5567.697 3378.656 1.65 0.100 -1059.236 12194.63 GSO | -7487.097 4858.838 -1.54 0.124 -17017.28 2043.083 HOU | 5792.695 3638.532 1.59 0.112 -1343.965 12929.35 HRL | 1982.608 1220.702 1.62 0.105 -411.6912 4376.907 HSV | -189.9168 385.3405 -0.49 0.622 -945.7281 565.8945 IAD | -10444.56 6813.609 -1.53 0.125 -23808.85 2919.73 IAH | 6226.797 3898.411 1.60 0.110 -1419.59 13873.19 ICT | 4383.489 2723.322 1.61 0.108 -958.0651 9725.044 IND | -3106.097 2180.991 -1.42 0.155 -7383.917 1171.723 JAN | 3219.581 1819.481 1.77 0.077 -349.1701 6788.332 JAX | -6150.915 3954.906 -1.56 0.120 -13908.11 1606.283 JFK | -14599.54 9288.162 -1.57 0.116 -32817.45 3618.369 LAS | -8628.749 5491.594 -1.57 0.116 -19400.03 2142.527 LAW | 7793.091 4847.539 1.61 0.108 -1714.927 17301.11 LAX | -11656.36 7525.466 -1.55 0.122 -26416.9 3104.174 LBB | 5155.208 3243.068 1.59 0.112 -1205.783 11516.2 LGA | -14258.21 9265.564 -1.54 0.124 -32431.79 3915.375 LIT | 4760.071 2994.493 1.59 0.112 -1113.361 10633.5 MAF | 4661.987 2937.998 1.59 0.113 -1100.636 10424.61 MCI | 2019.196 1231.998 1.64 0.101 -397.2589 4435.652 MCO | -7333.776 4689.35 -1.56 0.118 -16531.52 1863.969 MEM | 2780.146 1548.307 1.80 0.073 -256.7221 5817.014 MFE | 1891.149 1141.632 1.66 0.098 -348.0604 4130.358 MIA | -9714.339 6237.345 -1.56 0.120 -21948.34 2519.662 MKE | -4923.658 3209.301 -1.53 0.125 -11218.42 1371.101 MSP | -4749.191 3197.872 -1.49 0.138 -11021.54 1523.152 MSY | 2198.744 1367.555 1.61 0.108 -483.5942 4881.082 OAK | -15763.55 10022.62 -1.57 0.116 -35422.03 3894.932 OKC | 7127.469 4452.068 1.60 0.110 -1604.867 15859.81 OMA | 149.2239 160.9898 0.93 0.354 -166.5434 464.9912 ONT | -10868.83 7005.697 -1.55 0.121 -24609.88 2872.229 ORD | -3854.233 2632.932 -1.46 0.143 -9018.496 1310.03 ORF | -11376.93 7276.881 -1.56 0.118 -25649.89 2896.028 PBI | -9423.664 6033.958 -1.56 0.119 -21258.74 2411.412 PDX | -18636.78 11830.52 -1.58 0.115 -41841.29 4567.742 PHL | -12800.14 8282.521 -1.55 0.122 -29045.57 3445.297 PHX | -5178.207 3378.656 -1.53 0.126 -11805.14 1448.726 PIT | -8596.69 5638.484 -1.52 0.128 -19656.08 2462.698 PSP | -9805.087 6293.842 -1.56 0.119 -22149.9 2539.727 RDU | -8456.338 5570.689 -1.52 0.129 -19382.75 2470.076 RIC | -10263.65 6655.419 -1.54 0.123 -23317.66 2790.37 RNO | -13892.56 8768.392 -1.58 0.113 -31090.99 3305.861 SAN | -10571.89 6802.31 -1.55 0.120 -23914.02 2770.236 SAT | 5803.489 3638.532 1.60 0.111 -1333.171 12940.15 SDF | -2769.318 1853.344 -1.49 0.135 -6404.489 865.8522 SEA | -19454.48 12327.69 -1.58 0.115 -43634.16 4725.195 SFO | -15876.59 10113.02 -1.57 0.117 -35712.38 3959.191 SGF | 3828.7 2316.572 1.65 0.099 -715.0503 8372.45 SHV | 6939.565 4282.58 1.62 0.105 -1460.338 15339.47 SJC | -15227.89 9819.232 -1.55 0.121 -34487.45 4031.663 SJT | 6255.739 3853.214 1.62 0.105 -1302.001 13813.48 SLC | -7248.765 4734.547 -1.53 0.126 -16535.16 2037.63 SMF | -15282.49 9740.137 -1.57 0.117 -34386.91 3821.927 SNA | -11120.09 7186.486 -1.55 0.122 -25215.75 2975.564 SPS | 8313.63 5152.617 1.61 0.107 -1792.773 18420.03 STL | 675.4977 206.1322 3.28 0.001 271.1877 1079.808 TPA | -6259.084 4067.897 -1.54 0.124 -14237.9 1719.736 TUL | 5991.441 3751.523 1.60 0.110 -1366.839 13349.72 TUS | -4143.191 2757.218 -1.50 0.133 -9551.229 1264.848 | _cons | -10118.4 6429.399 -1.57 0.116 -22729.1 2492.297 ------------------------------------------------------------------------------ r; t=0.16 14:37:21 . predict yhat (option xb assumed; fitted values) r; t=0.05 14:38:35 . separate yhat, by(edest) storage display value variable name type format label variable label ------------------------------------------------------------------------------- yhat1 double %10.0g yhat, edest == ABQ yhat2 double %10.0g yhat, edest == ACT yhat3 double %10.0g yhat, edest == AMA yhat4 double %10.0g yhat, edest == ATL yhat5 double %10.0g yhat, edest == AUS yhat6 double %10.0g yhat, edest == BDL yhat7 double %10.0g yhat, edest == BHM yhat8 double %10.0g yhat, edest == BNA yhat9 double %10.0g yhat, edest == BOS yhat10 double %10.0g yhat, edest == BPT yhat11 double %10.0g yhat, edest == BTR yhat12 double %10.0g yhat, edest == BUR yhat13 double %10.0g yhat, edest == BWI yhat14 double %10.0g yhat, edest == CLE yhat15 double %10.0g yhat, edest == CLL yhat16 double %10.0g yhat, edest == CLT yhat17 double %10.0g yhat, edest == CMH yhat18 double %10.0g yhat, edest == COS yhat19 double %10.0g yhat, edest == CRP yhat20 double %10.0g yhat, edest == CVG yhat21 double %10.0g yhat, edest == DAY yhat22 double %10.0g yhat, edest == DCA yhat23 double %10.0g yhat, edest == DEN yhat24 double %10.0g yhat, edest == DSM yhat25 double %10.0g yhat, edest == DTW yhat26 double %10.0g yhat, edest == ELP yhat27 double %10.0g yhat, edest == EWR yhat28 double %10.0g yhat, edest == FAT yhat29 double %10.0g yhat, edest == FLL yhat30 double %10.0g yhat, edest == FYV yhat31 double %10.0g yhat, edest == GSO yhat32 double %10.0g yhat, edest == HOU yhat33 double %10.0g yhat, edest == HRL yhat34 double %10.0g yhat, edest == HSV yhat35 double %10.0g yhat, edest == IAD yhat36 double %10.0g yhat, edest == IAH yhat37 double %10.0g yhat, edest == ICT yhat38 double %10.0g yhat, edest == IND yhat39 double %10.0g yhat, edest == JAN yhat40 double %10.0g yhat, edest == JAX yhat41 double %10.0g yhat, edest == JFK yhat42 double %10.0g yhat, edest == LAS yhat43 double %10.0g yhat, edest == LAW yhat44 double %10.0g yhat, edest == LAX yhat45 double %10.0g yhat, edest == LBB yhat46 double %10.0g yhat, edest == LGA yhat47 double %10.0g yhat, edest == LIT yhat48 double %10.0g yhat, edest == MAF yhat49 double %10.0g yhat, edest == MCI yhat50 double %10.0g yhat, edest == MCO yhat51 double %10.0g yhat, edest == MEM yhat52 double %10.0g yhat, edest == MFE yhat53 double %10.0g yhat, edest == MIA yhat54 double %10.0g yhat, edest == MKE yhat55 double %10.0g yhat, edest == MSP yhat56 double %10.0g yhat, edest == MSY yhat57 double %10.0g yhat, edest == OAK yhat58 double %10.0g yhat, edest == OKC yhat59 double %10.0g yhat, edest == OMA yhat60 double %10.0g yhat, edest == ONT yhat61 double %10.0g yhat, edest == ORD yhat62 double %10.0g yhat, edest == ORF yhat63 double %10.0g yhat, edest == PBI yhat64 double %10.0g yhat, edest == PDX yhat65 double %10.0g yhat, edest == PHL yhat66 double %10.0g yhat, edest == PHX yhat67 double %10.0g yhat, edest == PIT yhat68 double %10.0g yhat, edest == PSP yhat69 double %10.0g yhat, edest == RDU yhat70 double %10.0g yhat, edest == RIC yhat71 double %10.0g yhat, edest == RNO yhat72 double %10.0g yhat, edest == SAN yhat73 double %10.0g yhat, edest == SAT yhat74 double %10.0g yhat, edest == SDF yhat75 double %10.0g yhat, edest == SEA yhat76 double %10.0g yhat, edest == SFO yhat77 double %10.0g yhat, edest == SGF yhat78 double %10.0g yhat, edest == SHV yhat79 double %10.0g yhat, edest == SJC yhat80 double %10.0g yhat, edest == SJT yhat81 double %10.0g yhat, edest == SLC yhat82 double %10.0g yhat, edest == SMF yhat83 double %10.0g yhat, edest == SNA yhat84 double %10.0g yhat, edest == SPS yhat85 double %10.0g yhat, edest == STL yhat86 double %10.0g yhat, edest == TPA yhat87 double %10.0g yhat, edest == TUL yhat88 double %10.0g yhat, edest == TUS r; t=0.28 14:38:59 . ed r; t=0.39 14:39:03 . twoway connected yhat1-yhat88 dist || lfit AAfare dist, clwidth(thick) clcolo > r(black) r; t=2.15 14:41:03 . reg AAfare dist Source | SS df MS Number of obs = 1,729 -------------+---------------------------------- F(1, 1727) = 3490.67 Model | 70621771.9 1 70621771.9 Prob > F = 0.0000 Residual | 34939940.2 1,727 20231.5809 R-squared = 0.6690 -------------+---------------------------------- Adj R-squared = 0.6688 Total | 105561712 1,728 61088.9538 Root MSE = 142.24 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | .2318662 .0039245 59.08 0.000 .2241689 .2395634 _cons | 173.8493 7.18437 24.20 0.000 159.7583 187.9403 ------------------------------------------------------------------------------ r; t=0.05 14:44:03 . estimates store ols r; t=0.06 14:44:14 . qui reg AAfare dist i.edest r; t=0.09 14:44:33 . estimates store ols_dummies r; t=0.03 14:44:46 . estimates table ols ols_dummies, star stats(N) ---------------------------------------------- Variable | ols ols_dummies -------------+-------------------------------- distance | .23186619*** 9.1342173 | edest | ACT | 8723.9541 AMA | 4584.8558 ATL | -2819.9993 AUS | 6984.4257 BDL | -15892.313 BHM | -145.35667 BNA | -833.93044 BOS | -17467.413 BPT | 5461.0014 BTR | 3415.8822 BUR | -11689.454 BWI | -11325.261 CLE | -7698.0484 CLL | 7381.852 CLT | -6292.5234 CMH | -6099.7476 COS | -280.3255 CRP | 3871.2184 CVG | -4059.3691 DAY | -4915.2549 DCA | -10850.149 DEN | -1188.7649 DSM | -719.5364 DTW | -7218.9877 ELP | 328.24682 EWR | -14151.519 FAT | -13168.569 FLL | -9664.642 FYV | 5567.6974 GSO | -7487.0974 HOU | 5792.6949 HRL | 1982.6079 HSV | -189.91678 IAD | -10444.562 IAH | 6226.7974 ICT | 4383.4892 IND | -3106.0969 JAN | 3219.581 JAX | -6150.9151 JFK | -14599.54 LAS | -8628.7492 LAW | 7793.0914 LAX | -11656.362 LBB | 5155.2082 LGA | -14258.209 LIT | 4760.0712 MAF | 4661.9875 MCI | 2019.1964 MCO | -7333.7758 MEM | 2780.146 MFE | 1891.1489 MIA | -9714.3394 MKE | -4923.6578 MSP | -4749.1915 MSY | 2198.7441 OAK | -15763.551 OKC | 7127.4692 OMA | 149.22392 ONT | -10868.827 ORD | -3854.2328 ORF | -11376.93 PBI | -9423.6641 PDX | -18636.775 PHL | -12800.135 PHX | -5178.2074 PIT | -8596.6899 PSP | -9805.0866 RDU | -8456.3377 RIC | -10263.645 RNO | -13892.565 SAN | -10571.894 SAT | 5803.4889 SDF | -2769.3183 SEA | -19454.482 SFO | -15876.594 SGF | 3828.7001 SHV | 6939.5647 SJC | -15227.892 SJT | 6255.7387 SLC | -7248.7646 SMF | -15282.489 SNA | -11120.093 SPS | 8313.6297 STL | 675.49774** TPA | -6259.0844 TUL | 5991.4408 TUS | -4143.1905 | _cons | 173.84929*** -10118.4 -------------+-------------------------------- N | 1729 1729 ---------------------------------------------- legend: * p<0.05; ** p<0.01; *** p<0.001 r; t=0.41 14:45:05 . xtset edest period panel variable: edest (strongly balanced) time variable: period, 1 to 20 delta: 1 unit r; t=0.33 14:46:00 . xtdes edest: 1, 2, ..., 88 n = 88 period: 1, 2, ..., 20 T = 20 Delta(period) = 1 unit Span(period) = 20 periods (edest*period uniquely identifies each observation) Distribution of T_i: min 5% 25% 50% 75% 95% max 20 20 20 20 20 20 20 Freq. Percent Cum. | Pattern ---------------------------+---------------------- 88 100.00 100.00 | 11111111111111111111 ---------------------------+---------------------- 88 100.00 | XXXXXXXXXXXXXXXXXXXX r; t=0.36 14:46:20 . xtline AAfare r; t=26.97 14:48:18 . xtline AAfare, overlay r; t=4.59 14:48:58 . xtreg AAfare distance, fe Fixed-effects (within) regression Number of obs = 1,729 Group variable: edest Number of groups = 88 R-sq: Obs per group: within = 0.0016 min = 7 between = 0.7752 avg = 19.6 overall = 0.6690 max = 20 F(1,1640) = 2.61 corr(u_i, Xb) = -0.9999 Prob > F = 0.1061 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | 9.134217 5.649705 1.62 0.106 -1.947178 20.21561 _cons | -14157.38 9095.038 -1.56 0.120 -31996.49 3681.731 -------------+---------------------------------------------------------------- sigma_u | 7756.2678 sigma_e | 94.50387 rho | .99985157 (fraction of variance due to u_i) ------------------------------------------------------------------------------ F test that all u_i=0: F(87, 1640) = 26.12 Prob > F = 0.0000 r; t=0.42 14:49:27 . xtreg AAfare distance, re Random-effects GLS regression Number of obs = 1,729 Group variable: edest Number of groups = 88 R-sq: Obs per group: within = 0.0016 min = 7 between = 0.7752 avg = 19.6 overall = 0.6690 max = 20 Wald chi2(1) = 296.84 corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | .2316695 .0134464 17.23 0.000 .205315 .258024 _cons | 174.6588 24.53939 7.12 0.000 126.5625 222.7551 -------------+---------------------------------------------------------------- sigma_u | 107.13235 sigma_e | 94.50387 rho | .56238531 (fraction of variance due to u_i) ------------------------------------------------------------------------------ r; t=0.34 14:50:17 . help areg r; t=0.62 14:53:30 . areg AAfare distance, absorb(edest) Linear regression, absorbing indicators Number of obs = 1,729 F( 1, 1640) = 2.61 Prob > F = 0.1061 R-squared = 0.8612 Adj R-squared = 0.8538 Root MSE = 94.5039 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | 9.134217 5.649705 1.62 0.106 -1.947178 20.21561 _cons | -14157.38 9095.038 -1.56 0.120 -31996.49 3681.731 -------------+---------------------------------------------------------------- edest | F(87, 1640) = 26.117 0.000 (88 categories) r; t=0.41 14:54:01 . xtreg AAfare distance, fe Fixed-effects (within) regression Number of obs = 1,729 Group variable: edest Number of groups = 88 R-sq: Obs per group: within = 0.0016 min = 7 between = 0.7752 avg = 19.6 overall = 0.6690 max = 20 F(1,1640) = 2.61 corr(u_i, Xb) = -0.9999 Prob > F = 0.1061 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | 9.134217 5.649705 1.62 0.106 -1.947178 20.21561 _cons | -14157.38 9095.038 -1.56 0.120 -31996.49 3681.731 -------------+---------------------------------------------------------------- sigma_u | 7756.2678 sigma_e | 94.50387 rho | .99985157 (fraction of variance due to u_i) ------------------------------------------------------------------------------ F test that all u_i=0: F(87, 1640) = 26.12 Prob > F = 0.0000 r; t=0.17 14:59:09 . est store fixed r; t=0.11 14:59:14 . xtreg AAfare distance, re Random-effects GLS regression Number of obs = 1,729 Group variable: edest Number of groups = 88 R-sq: Obs per group: within = 0.0016 min = 7 between = 0.7752 avg = 19.6 overall = 0.6690 max = 20 Wald chi2(1) = 296.84 corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | .2316695 .0134464 17.23 0.000 .205315 .258024 _cons | 174.6588 24.53939 7.12 0.000 126.5625 222.7551 -------------+---------------------------------------------------------------- sigma_u | 107.13235 sigma_e | 94.50387 rho | .56238531 (fraction of variance due to u_i) ------------------------------------------------------------------------------ r; t=0.14 14:59:22 . est store random r; t=0.02 14:59:26 . hausman fixed random ---- Coefficients ---- | (b) (B) (b-B) sqrt(diag(V_b-V_B)) | fixed random Difference S.E. -------------+---------------------------------------------------------------- distance | 9.134217 .2316695 8.902548 5.649689 ------------------------------------------------------------------------------ b = consistent under Ho and Ha; obtained from xtreg B = inconsistent under Ha, efficient under Ho; obtained from xtreg Test: Ho: difference in coefficients not systematic chi2(1) = (b-B)'[(V_b-V_B)^(-1)](b-B) = 2.48 Prob>chi2 = 0.1151 r; t=0.20 14:59:39 . hausman fixed random, sigmamore ---- Coefficients ---- | (b) (B) (b-B) sqrt(diag(V_b-V_B)) | fixed random Difference S.E. -------------+---------------------------------------------------------------- distance | 9.134217 .2316695 8.902548 5.652207 ------------------------------------------------------------------------------ b = consistent under Ho and Ha; obtained from xtreg B = inconsistent under Ha, efficient under Ho; obtained from xtreg Test: Ho: difference in coefficients not systematic chi2(1) = (b-B)'[(V_b-V_B)^(-1)](b-B) = 2.48 Prob>chi2 = 0.1152 r; t=0.06 15:02:22 . xtreg AAfare distance i.period Random-effects GLS regression Number of obs = 1,729 Group variable: edest Number of groups = 88 R-sq: Obs per group: within = 0.1825 min = 7 between = 0.7755 avg = 19.6 overall = 0.6948 max = 20 Wald chi2(20) = 660.41 corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | .2317864 .0134083 17.29 0.000 .2055066 .2580662 | period | 2 | -64.69106 12.96495 -4.99 0.000 -90.10191 -39.28022 3 | -53.45679 12.96495 -4.12 0.000 -78.86763 -28.04594 4 | -82.22591 13.0047 -6.32 0.000 -107.7146 -56.73717 5 | 18.73503 12.96495 1.45 0.148 -6.675812 44.14588 6 | -3.804172 12.96495 -0.29 0.769 -29.21502 21.60667 7 | -38.64099 12.96495 -2.98 0.003 -64.05183 -13.23015 8 | -57.5612 13.08928 -4.40 0.000 -83.21572 -31.90667 9 | -24.44524 13.04796 -1.87 0.061 -50.01877 1.128284 10 | -36.88024 13.04796 -2.83 0.005 -62.45377 -11.30672 11 | -65.4575 13.08928 -5.00 0.000 -91.11202 -39.80298 12 | -95.47245 13.04796 -7.32 0.000 -121.046 -69.89893 13 | -103.5668 13.04796 -7.94 0.000 -129.1403 -77.99323 14 | -122.2619 13.08898 -9.34 0.000 -147.9159 -96.60799 15 | -90.4904 13.08898 -6.91 0.000 -116.1443 -64.83646 16 | -101.2924 13.04724 -7.76 0.000 -126.8645 -75.72029 17 | -85.44977 13.04754 -6.55 0.000 -111.0225 -59.87706 18 | -79.66233 13.04754 -6.11 0.000 -105.235 -54.08962 19 | -110.6708 13.00702 -8.51 0.000 -136.1641 -85.17756 20 | -108.8706 13.08857 -8.32 0.000 -134.5238 -83.2175 | _cons | 239.1718 26.04539 9.18 0.000 188.1237 290.2198 -------------+---------------------------------------------------------------- sigma_u | 107.23725 sigma_e | 86.007861 rho | .60854731 (fraction of variance due to u_i) ------------------------------------------------------------------------------ r; t=0.30 15:06:49 . testparm i.period ( 1) 2.period = 0 ( 2) 3.period = 0 ( 3) 4.period = 0 ( 4) 5.period = 0 ( 5) 6.period = 0 ( 6) 7.period = 0 ( 7) 8.period = 0 ( 8) 9.period = 0 ( 9) 10.period = 0 (10) 11.period = 0 (11) 12.period = 0 (12) 13.period = 0 (13) 14.period = 0 (14) 15.period = 0 (15) 16.period = 0 (16) 17.period = 0 (17) 18.period = 0 (18) 19.period = 0 (19) 20.period = 0 chi2( 19) = 361.85 Prob > chi2 = 0.0000 r; t=0.09 15:07:17 . xtreg AAfare distance i.edest Random-effects GLS regression Number of obs = 1,729 Group variable: edest Number of groups = 88 R-sq: Obs per group: within = 0.0016 min = 7 between = 1.0000 avg = 19.6 overall = 0.8612 max = 20 Wald chi2(88) = 10179.72 corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | 9.134217 5.649705 1.62 0.106 -1.939 20.20743 | edest | ACT | 8723.954 5423.799 1.61 0.108 -1906.496 19354.4 AMA | 4584.856 2892.803 1.58 0.113 -1084.934 10254.65 ATL | -2819.999 1842.046 -1.53 0.126 -6430.343 790.3447 AUS | 6984.426 4361.674 1.60 0.109 -1564.299 15533.15 BDL | -15892.31 10192.11 -1.56 0.119 -35868.48 4083.857 BHM | -145.3567 317.7917 -0.46 0.647 -768.217 477.5037 BNA | -833.9304 701.2005 -1.19 0.234 -2208.258 540.3973 BOS | -17467.41 11220.35 -1.56 0.120 -39458.9 4524.075 BPT | 5461.001 3378.656 1.62 0.106 -1161.042 12083.04 BTR | 3415.882 2101.948 1.63 0.104 -703.8603 7535.625 BUR | -11689.45 7480.269 -1.56 0.118 -26350.51 2971.603 BWI | -11325.26 7322.078 -1.55 0.122 -25676.27 3025.749 CLE | -7698.048 5107.42 -1.51 0.132 -17708.41 2312.312 CLL | 7381.852 4576.358 1.61 0.107 -1587.645 16351.35 CLT | -6292.523 4158.29 -1.51 0.130 -14442.62 1857.575 CMH | -6099.748 4045.299 -1.51 0.132 -14028.39 1828.892 COS | -280.3255 261.599 -1.07 0.284 -793.0501 232.3992 CRP | 3871.218 2429.557 1.59 0.111 -890.6253 8633.062 CVG | -4059.369 2745.973 -1.48 0.139 -9441.378 1322.64 DAY | -4915.255 3299.563 -1.49 0.136 -11382.28 1551.769 DCA | -10850.15 7039.595 -1.54 0.123 -24647.5 2947.204 DEN | -1188.765 831.0441 -1.43 0.153 -2817.581 440.0515 DSM | -719.5364 622.1856 -1.16 0.247 -1938.998 499.925 DTW | -7218.988 4723.248 -1.53 0.126 -16476.38 2038.407 ELP | 328.2468 205.5732 1.60 0.110 -74.6692 731.1628 EWR | -14151.52 9073.475 -1.56 0.119 -31935.2 3632.165 FAT | -13168.57 8406.814 -1.57 0.117 -29645.62 3308.482 FLL | -9664.642 6214.747 -1.56 0.120 -21845.32 2516.038 FYV | 5567.697 3378.656 1.65 0.099 -1054.346 12189.74 GSO | -7487.097 4858.838 -1.54 0.123 -17010.24 2036.05 HOU | 5792.695 3638.532 1.59 0.111 -1338.698 12924.09 HRL | 1982.608 1220.702 1.62 0.104 -409.9241 4375.14 HSV | -189.9168 385.3405 -0.49 0.622 -945.1703 565.3367 IAD | -10444.56 6813.609 -1.53 0.125 -23798.99 2909.867 IAH | 6226.797 3898.411 1.60 0.110 -1413.947 13867.54 ICT | 4383.489 2723.322 1.61 0.107 -954.123 9721.101 IND | -3106.097 2180.991 -1.42 0.154 -7380.76 1168.566 JAN | 3219.581 1819.481 1.77 0.077 -346.5362 6785.698 JAX | -6150.915 3954.906 -1.56 0.120 -13902.39 1600.558 JFK | -14599.54 9288.162 -1.57 0.116 -32804 3604.924 LAS | -8628.749 5491.594 -1.57 0.116 -19392.08 2134.578 LAW | 7793.091 4847.539 1.61 0.108 -1707.91 17294.09 LAX | -11656.36 7525.466 -1.55 0.121 -26406 3093.28 LBB | 5155.208 3243.068 1.59 0.112 -1201.088 11511.5 LGA | -14258.21 9265.564 -1.54 0.124 -32418.38 3901.962 LIT | 4760.071 2994.493 1.59 0.112 -1109.026 10629.17 MAF | 4661.988 2937.998 1.59 0.113 -1096.383 10420.36 MCI | 2019.196 1231.998 1.64 0.101 -395.4755 4433.868 MCO | -7333.776 4689.35 -1.56 0.118 -16524.73 1857.181 MEM | 2780.146 1548.307 1.80 0.073 -254.4809 5814.773 MFE | 1891.149 1141.632 1.66 0.098 -346.4078 4128.706 MIA | -9714.339 6237.345 -1.56 0.119 -21939.31 2510.633 MKE | -4923.658 3209.301 -1.53 0.125 -11213.77 1366.456 MSP | -4749.192 3197.872 -1.49 0.138 -11016.91 1518.523 MSY | 2198.744 1367.555 1.61 0.108 -481.6146 4879.103 OAK | -15763.55 10022.62 -1.57 0.116 -35407.53 3880.424 OKC | 7127.469 4452.068 1.60 0.109 -1598.423 15853.36 OMA | 149.2239 160.9898 0.93 0.354 -166.3103 464.7582 ONT | -10868.83 7005.697 -1.55 0.121 -24599.74 2862.088 ORD | -3854.233 2632.932 -1.46 0.143 -9014.685 1306.219 ORF | -11376.93 7276.881 -1.56 0.118 -25639.35 2885.494 PBI | -9423.664 6033.958 -1.56 0.118 -21250.01 2402.677 PDX | -18636.78 11830.52 -1.58 0.115 -41824.17 4550.616 PHL | -12800.14 8282.521 -1.55 0.122 -29033.58 3433.307 PHX | -5178.207 3378.656 -1.53 0.125 -11800.25 1443.836 PIT | -8596.69 5638.484 -1.52 0.127 -19647.92 2454.536 PSP | -9805.087 6293.842 -1.56 0.119 -22140.79 2530.617 RDU | -8456.338 5570.689 -1.52 0.129 -19374.69 2462.012 RIC | -10263.65 6655.419 -1.54 0.123 -23308.03 2780.736 RNO | -13892.56 8768.392 -1.58 0.113 -31078.3 3293.169 SAN | -10571.89 6802.31 -1.55 0.120 -23904.18 2760.389 SAT | 5803.489 3638.532 1.60 0.111 -1327.904 12934.88 SDF | -2769.318 1853.344 -1.49 0.135 -6401.806 863.1693 SEA | -19454.48 12327.69 -1.58 0.115 -43616.31 4707.35 SFO | -15876.59 10113.02 -1.57 0.116 -35697.74 3944.552 SGF | 3828.7 2316.572 1.65 0.098 -711.6969 8369.097 SHV | 6939.565 4282.58 1.62 0.105 -1454.138 15333.27 SJC | -15227.89 9819.232 -1.55 0.121 -34473.23 4017.449 SJT | 6255.739 3853.214 1.62 0.104 -1296.423 13807.9 SLC | -7248.765 4734.547 -1.53 0.126 -16528.31 2030.776 SMF | -15282.49 9740.137 -1.57 0.117 -34372.81 3807.828 SNA | -11120.09 7186.486 -1.55 0.122 -25205.35 2965.161 SPS | 8313.63 5152.617 1.61 0.107 -1785.315 18412.57 STL | 675.4977 206.1322 3.28 0.001 271.4861 1079.509 TPA | -6259.084 4067.897 -1.54 0.124 -14232.02 1713.847 TUL | 5991.441 3751.523 1.60 0.110 -1361.409 13344.29 TUS | -4143.191 2757.218 -1.50 0.133 -9547.238 1260.857 | _cons | -10118.4 6429.399 -1.57 0.116 -22719.79 2482.99 -------------+---------------------------------------------------------------- sigma_u | 0 sigma_e | 94.50387 rho | 0 (fraction of variance due to u_i) ------------------------------------------------------------------------------ r; t=0.56 15:08:57 . testparm i.edest ( 1) 2.edest = 0 ( 2) 3.edest = 0 ( 3) 4.edest = 0 ( 4) 5.edest = 0 ( 5) 6.edest = 0 ( 6) 7.edest = 0 ( 7) 8.edest = 0 ( 8) 9.edest = 0 ( 9) 10.edest = 0 (10) 11.edest = 0 (11) 12.edest = 0 (12) 13.edest = 0 (13) 14.edest = 0 (14) 15.edest = 0 (15) 16.edest = 0 (16) 17.edest = 0 (17) 18.edest = 0 (18) 19.edest = 0 (19) 20.edest = 0 (20) 21.edest = 0 (21) 22.edest = 0 (22) 23.edest = 0 (23) 24.edest = 0 (24) 25.edest = 0 (25) 26.edest = 0 (26) 27.edest = 0 (27) 28.edest = 0 (28) 29.edest = 0 (29) 30.edest = 0 (30) 31.edest = 0 (31) 32.edest = 0 (32) 33.edest = 0 (33) 34.edest = 0 (34) 35.edest = 0 (35) 36.edest = 0 (36) 37.edest = 0 (37) 38.edest = 0 (38) 39.edest = 0 (39) 40.edest = 0 (40) 41.edest = 0 (41) 42.edest = 0 (42) 43.edest = 0 (43) 44.edest = 0 (44) 45.edest = 0 (45) 46.edest = 0 (46) 47.edest = 0 (47) 48.edest = 0 (48) 49.edest = 0 (49) 50.edest = 0 (50) 51.edest = 0 (51) 52.edest = 0 (52) 53.edest = 0 (53) 54.edest = 0 (54) 55.edest = 0 (55) 56.edest = 0 (56) 57.edest = 0 (57) 58.edest = 0 (58) 59.edest = 0 (59) 60.edest = 0 (60) 61.edest = 0 (61) 62.edest = 0 (62) 63.edest = 0 (63) 64.edest = 0 (64) 65.edest = 0 (65) 66.edest = 0 (66) 67.edest = 0 (67) 68.edest = 0 (68) 69.edest = 0 (69) 70.edest = 0 (70) 71.edest = 0 (71) 72.edest = 0 (72) 73.edest = 0 (73) 74.edest = 0 (74) 75.edest = 0 (75) 76.edest = 0 (76) 77.edest = 0 (77) 78.edest = 0 (78) 79.edest = 0 (79) 80.edest = 0 (80) 81.edest = 0 (81) 82.edest = 0 (82) 83.edest = 0 (83) 84.edest = 0 (84) 85.edest = 0 (85) 86.edest = 0 (86) 87.edest = 0 (87) 88.edest = 0 chi2( 87) = 2272.22 Prob > chi2 = 0.0000 r; t=0.11 15:09:01 . xtreg AAfare distance, re Random-effects GLS regression Number of obs = 1,729 Group variable: edest Number of groups = 88 R-sq: Obs per group: within = 0.0016 min = 7 between = 0.7752 avg = 19.6 overall = 0.6690 max = 20 Wald chi2(1) = 296.84 corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | .2316695 .0134464 17.23 0.000 .205315 .258024 _cons | 174.6588 24.53939 7.12 0.000 126.5625 222.7551 -------------+---------------------------------------------------------------- sigma_u | 107.13235 sigma_e | 94.50387 rho | .56238531 (fraction of variance due to u_i) ------------------------------------------------------------------------------ r; t=0.20 15:09:46 . xttest0 Breusch and Pagan Lagrangian multiplier test for random effects AAfare[edest,t] = Xb + u[edest] + e[edest,t] Estimated results: | Var sd = sqrt(Var) ---------+----------------------------- AAfare | 61088.95 247.1618 e | 8930.981 94.50387 u | 11477.34 107.1323 Test: Var(u) = 0 chibar2(01) = 5063.53 Prob > chibar2 = 0.0000 r; t=0.19 15:09:58 . ssc install xtscc checking xtscc consistency and verifying not already installed... installing into c:\ado\plus\... installation complete. r; t=5.52 15:16:32 . xtreg AAfare dist, fe Fixed-effects (within) regression Number of obs = 1,729 Group variable: edest Number of groups = 88 R-sq: Obs per group: within = 0.0016 min = 7 between = 0.7752 avg = 19.6 overall = 0.6690 max = 20 F(1,1640) = 2.61 corr(u_i, Xb) = -0.9999 Prob > F = 0.1061 ------------------------------------------------------------------------------ AAfare | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | 9.134217 5.649705 1.62 0.106 -1.947178 20.21561 _cons | -14157.38 9095.038 -1.56 0.120 -31996.49 3681.731 -------------+---------------------------------------------------------------- sigma_u | 7756.2678 sigma_e | 94.50387 rho | .99985157 (fraction of variance due to u_i) ------------------------------------------------------------------------------ F test that all u_i=0: F(87, 1640) = 26.12 Prob > F = 0.0000 r; t=0.28 15:16:51 . ssc install xttest3 checking xttest3 consistency and verifying not already installed... installing into c:\ado\plus\... installation complete. r; t=2.82 15:17:50 . xttest3 Modified Wald test for groupwise heteroskedasticity in fixed effect regression model H0: sigma(i)^2 = sigma^2 for all i chi2 (88) = 2.0e+05 Prob>chi2 = 0.0000 r; t=0.17 15:18:04 . xtreg AAfare dist, fe robust Fixed-effects (within) regression Number of obs = 1,729 Group variable: edest Number of groups = 88 R-sq: Obs per group: within = 0.0016 min = 7 between = 0.7752 avg = 19.6 overall = 0.6690 max = 20 F(1,87) = 111.51 corr(u_i, Xb) = -0.9999 Prob > F = 0.0000 (Std. Err. adjusted for 88 clusters in edest) ------------------------------------------------------------------------------ | Robust AAfare | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- distance | 9.134217 .8649972 10.56 0.000 7.414942 10.85349 _cons | -14157.38 1392.494 -10.17 0.000 -16925.11 -11389.65 -------------+---------------------------------------------------------------- sigma_u | 7756.2678 sigma_e | 94.50387 rho | .99985157 (fraction of variance due to u_i) ------------------------------------------------------------------------------ r; t=0.31 15:18:37 . log close name: log: c:\users\bchavez2\desktop\week3_lecture.txt log type: text closed on: 4 Feb 2016, 11:52:56 -------------------------------------------------------------------------------