5. Thảo luận kết quả về mối quan hệ giữa giới tính CEO và tỷ suất sinh lợi của
5.3 Giới hạn của bài nghiên cứu
Trên thực tế, bài nghiên cứu vẫn còn tồn tại một số hạn chế như sau:
(v) Kết quả nghiên cứu của bài được xem xét đối với tỷ suất sinh lợi vượt trội của cổ phiếu công ty trong khoản thời gian từ 2007 đến 2014, khoản thời gian này chưa đủ dài để nhiệm kỳ của CEO có biến động. Bên cạnh đó, đặc thù các cơng ty tại Việt Nam thì tình trạng kiêm nhiệm chức vụ của Chủ tịch hội đồng quản trị và CEO của công ty là khá phổ biến, làm cho nhiệm kỳ của CEO càng ít biến động hơn.
(vi) Mẫu quan sát thu thập gần như đầy đủ các công ty đủ điều kiện niêm yết trên thị trường từ năm 2007 đến 2014. Do đặc điểm thị trường nên mẫu thu thập theo dạng panel data khơng cân bằng (do có nhiều giá trị khơng đủ
trong các giai đoạn) và bỏ qua sự mất cân bẳng trong phân bố của cơng ty có CEO nữ trong các giai đoạn xem xét.
(vii) Về mặt mơ hình áp dụng trong bài, mơ hình bốn nhân tố được xây dựng dựa trên mơ hình ba nhân tố của Fama - French (1993) và nhân tố thứ tư MOM theo Carhart (1997), mơ hình có những hạn chế nhất định, theo các kết quả nghiên cứu của nhiều tác giả trên thế giới đã chỉ ra rằng mơ hình này chỉ tập trung vào nguồn gốc của lợi nhuận hơn là tổng rủi ro của nó. Bên cạnh đó, mơ hình dự báo này chỉ được được vận hành tốt khi các nhà đầu tư có thơng tin như nhau, thơng tin khơng được rị rỉ và phải minh bạch. Tuy nhiên, mức độ hiệu quả tại thị trường chứng khốn của Việt Nam cịn rất kém, và đã tạo cơ hội kiếm được lợi nhuận cao cho những nhà đầu cơ là những người có khả năng tiếp cận với các nguồn thông tin sớm nhất.
(viii) Bài nghiên cứu đánh giá chung cho tổng thể thị trường nhưng trên thực tế sẽ có một số ngành là phù hợp với nữ hơn so với nam giới và ngược lại. Bên cạnh đó, các yếu tố mang tính vĩ mơ (khủng hoảng kinh tế, suy thối, lạm phát, luật pháp, chính trị), yếu tố con người (tâm lý, mối quan hệ, thông tin hành lang), các thông tin (tỷ giá hối đoái đồng Việt Nam so với đồng đô la Mỹ, tỷ lệ lãi suất ngắn hạn trong nước) cũng chưa được đưa vào trong mơ hình để đánh giá tác động tổng thể.
(ix) Dữ liệu của mẫu được lấy từ các công ty niêm yết trên thị trường chứng khoán Việt Nam tại hai sàn giao dịch tập trung là HOSE và HXN nhưng lại bỏ qua các công ty hoạt động trên thị trường phi tập trung OTC và UpCom mà hoạt động của thị trường này được coi là khá sôi nổi và chứa nhiều rủi ro.
Phụ lục 1: Thống kê mơ tả các biến trong mơ hình
Phụ lục 2: Ma trận tương quan và các nhân tố VIF
r_MOM 40601 .0376806 .0639304 -.1775545 .308794 r_HML 40601 -.044125 .0612258 -.2657585 .129827 r_SMB 40601 .0056694 .0459912 -.1741283 .370661 r_Mkt 40601 .0012838 .0795143 -.240092 .2798557 female 40601 .0681018 .2519236 0 1 Ri 40601 .010608 .1729851 -.7621622 3.125 Variable Obs Mean Std. Dev. Min Max . sum Ri female r_Mkt r_SMB r_HML r_MOM
femaler_MOM -0.0208 0.4944 -0.0314 -0.1141 0.0690 0.3212 -0.2060 -0.9290 0.2269 1.0000 r_MOM -0.0883 -0.0000 -0.1258 -0.0364 0.2947 0.0790 -0.9049 -0.1942 1.0000 femaler_HML 0.0201 -0.5711 0.0316 0.1137 -0.0837 -0.3863 0.2146 1.0000 r_HML 0.0856 -0.0002 0.1326 0.0388 -0.3822 -0.1016 1.0000 femaler_SMB -0.0190 0.1175 -0.0660 -0.2542 0.2590 1.0000 r_SMB 0.0081 -0.0002 -0.2555 -0.0672 1.0000 femaler_Mkt 0.0948 0.0144 0.2576 1.0000 r_Mkt 0.4450 -0.0004 1.0000 female 0.0021 1.0000 Ri 1.0000 Ri female r_Mkt female~t r_SMB female~B r_HML female~L r_MOM female~M (obs=40601)
Phụ lục 3: Kiểm định phương sai sai số thay đổi Mean VIF 4.10 r_Mkt 1.15 0.869767 femaler_Mkt 1.15 0.867626 r_SMB 1.35 0.742230 femaler_SMB 1.37 0.727645 female 1.59 0.630421 r_MOM 6.03 0.165710 r_HML 6.44 0.155193 femaler_MOM 8.09 0.123567 femaler_HML 9.71 0.103014 Total 1873.52 39 0.0000 Kurtosis 14.74 1 0.0001 Skewness 163.61 9 0.0000 Heteroskedasticity 1695.16 29 0.0000 Source chi2 df p Cameron & Trivedi's decomposition of IM-test
Prob > chi2 = 0.0000 chi2(29) = 1695.16
against Ha: unrestricted heteroskedasticity White's test for Ho: homoskedasticity
Hồi quy Return với Female, MktRf, Female* MktRf
_cons .0105114 .0008893 11.82 0.000 .0087683 .0122545 female .0014184 .0034079 0.42 0.677 -.005261 .0080979 Ri Coef. Std. Err. t P>|t| [95% Conf. Interval] Total 1214.90821 40600 .029923847 Root MSE = .17299 Adj R-squared = -0.0000 Residual 1214.90302 40599 .029924457 R-squared = 0.0000 Model .005184181 1 .005184181 Prob > F = 0.6772 F( 1, 40599) = 0.17 Source SS df MS Number of obs = 40601 . reg Ri female _cons .0092451 .0007963 11.61 0.000 .0076843 .0108059 femaler_Mkt -.1796808 .0388364 -4.63 0.000 -.255801 -.1035607 r_Mkt .9800837 .0100043 97.97 0.000 .9604751 .9996923 female .0017481 .0030514 0.57 0.567 -.0042327 .007729 Ri Coef. Std. Err. t P>|t| [95% Conf. Interval] Total 1214.90821 40600 .029923847 Root MSE = .15488 Adj R-squared = 0.1984 Residual 973.781081 40597 .023986528 R-squared = 0.1985 Model 241.127125 3 80.3757083 Prob > F = 0.0000 F( 3, 40597) = 3350.87 Source SS df MS Number of obs = 40601 . reg Ri female r_Mkt femaler_Mkt
Hồi quy Return với Female, MktRf, Female* MktRf, SMB, Female* SMB, HML, Female* HML _cons .0061975 .0007948 7.80 0.000 .0046396 .0077554 femaler_SMB -.4412849 .0679776 -6.49 0.000 -.5745225 -.3080473 r_SMB .5198329 .0177187 29.34 0.000 .4851039 .5545619 femaler_Mkt -.244533 .039773 -6.15 0.000 -.322489 -.166577 r_Mkt 1.056789 .0102394 103.21 0.000 1.03672 1.076859 female .0043393 .0030453 1.42 0.154 -.0016297 .0103082 Ri Coef. Std. Err. t P>|t| [95% Conf. Interval] Total 1214.90821 40600 .029923847 Root MSE = .15326 Adj R-squared = 0.2150 Residual 953.529957 40595 .023488852 R-squared = 0.2151 Model 261.378249 5 52.2756498 Prob > F = 0.0000 F( 5, 40595) = 2225.55 Source SS df MS Number of obs = 40601
_cons .0162521 .0009779 16.62 0.000 .0143355 .0181687 femaler_HML -.0671468 .0533846 -1.26 0.208 -.1717818 .0374883 r_HML .2431764 .0138795 17.52 0.000 .2159722 .2703805 femaler_SMB -.4729352 .072865 -6.49 0.000 -.6157522 -.3301182 r_SMB .6403225 .0189424 33.80 0.000 .6031949 .67745 femaler_Mkt -.2447464 .0396739 -6.17 0.000 -.3225082 -.1669846 r_Mkt 1.049989 .0102071 102.87 0.000 1.029982 1.069995 female .0015689 .0037526 0.42 0.676 -.0057862 .0089241 Ri Coef. Std. Err. t P>|t| [95% Conf. Interval] Total 1214.90821 40600 .029923847 Root MSE = .15267 Adj R-squared = 0.2211 Residual 946.103619 40593 .023307063 R-squared = 0.2213 Model 268.804587 7 38.4006553 Prob > F = 0.0000 F( 7, 40593) = 1647.60 Source SS df MS Number of obs = 40601 . reg Ri female r_Mkt femaler_Mkt r_SMB femaler_SMB r_HML femaler_HML
Phụ lục 5: Kết quả chạy hồi quy GLS cho tổng mẫu
Hồi quy Return với Female
_cons .0162777 .0009866 16.50 0.000 .014344 .0182115 femaler_MOM -.0386319 .1122592 -0.34 0.731 -.2586625 .1813988 r_MOM .0056971 .0291145 0.20 0.845 -.0513679 .0627621 femaler_HML -.1045659 .1211401 -0.86 0.388 -.3420033 .1328715 r_HML .2486908 .0314138 7.92 0.000 .1871191 .3102625 femaler_SMB -.4763514 .0735377 -6.48 0.000 -.620487 -.3322158 r_SMB .6408343 .0191226 33.51 0.000 .6033536 .678315 femaler_Mkt -.2453157 .0397089 -6.18 0.000 -.3231461 -.1674854 r_Mkt 1.050078 .0102175 102.77 0.000 1.030051 1.070104 female .0013917 .003788 0.37 0.713 -.0060328 .0088163 Ri Coef. Std. Err. t P>|t| [95% Conf. Interval] Total 1214.90821 40600 .029923847 Root MSE = .15267 Adj R-squared = 0.2211 Residual 946.100575 40591 .023308137 R-squared = 0.2213 Model 268.80763 9 29.8675145 Prob > F = 0.0000 F( 9, 40591) = 1281.42 Source SS df MS Number of obs = 40601
. reg Ri female r_Mkt femaler_Mkt r_SMB femaler_SMB r_HML femaler_HML r_MOM femaler_MOM
_cons .0105114 .0008893 11.82 0.000 .0087684 .0122544 female .0014184 .0034078 0.42 0.677 -.0052607 .0080975 Ri Coef. Std. Err. z P>|z| [95% Conf. Interval] Log likelihood = 13626.74 Prob > chi2 = 0.6772 Wald chi2(1) = 0.17 max = 96 avg = 66.77796 Estimated coefficients = 2 Obs per group: min = 1 Estimated autocorrelations = 0 Number of groups = 608 Estimated covariances = 1 Number of obs = 40601 Correlation: no autocorrelation
Panels: homoskedastic
Coefficients: generalized least squares Cross-sectional time-series FGLS regression . xtgls Ri female
Hồi quy Return với Female, MktRf, Female* MktRf, SMB, Female* SMB _cons .0092451 .0007963 11.61 0.000 .0076844 .0108058 femaler_Mkt -.1796808 .0388345 -4.63 0.000 -.255795 -.1035667 r_Mkt .9800837 .0100038 97.97 0.000 .9604767 .9996908 female .0017481 .0030513 0.57 0.567 -.0042322 .0077285 Ri Coef. Std. Err. z P>|z| [95% Conf. Interval] Log likelihood = 18117.88 Prob > chi2 = 0.0000 Wald chi2(3) = 10053.60 max = 96 avg = 66.77796 Estimated coefficients = 4 Obs per group: min = 1 Estimated autocorrelations = 0 Number of groups = 608 Estimated covariances = 1 Number of obs = 40601 Correlation: no autocorrelation
Panels: homoskedastic
Coefficients: generalized least squares Cross-sectional time-series FGLS regression
_cons .0061975 .0007948 7.80 0.000 .0046398 .0077552 femaler_SMB -.4412849 .0679725 -6.49 0.000 -.5745086 -.3080611 r_SMB .5198329 .0177174 29.34 0.000 .4851075 .5545583 femaler_Mkt -.244533 .0397701 -6.15 0.000 -.3224809 -.1665851 r_Mkt 1.056789 .0102386 103.22 0.000 1.036722 1.076857 female .0043393 .0030451 1.42 0.154 -.001629 .0103076 Ri Coef. Std. Err. z P>|z| [95% Conf. Interval] Log likelihood = 18544.51 Prob > chi2 = 0.0000 Wald chi2(5) = 11129.40 max = 96 avg = 66.77796 Estimated coefficients = 6 Obs per group: min = 1 Estimated autocorrelations = 0 Number of groups = 608 Estimated covariances = 1 Number of obs = 40601 Correlation: no autocorrelation
Panels: homoskedastic
Coefficients: generalized least squares Cross-sectional time-series FGLS regression
_cons .0162521 .0009778 16.62 0.000 .0143357 .0181685 femaler_HML -.0671468 .0533794 -1.26 0.208 -.1717684 .0374749 r_HML .2431764 .0138781 17.52 0.000 .2159757 .270377 femaler_SMB -.4729352 .0728578 -6.49 0.000 -.6157338 -.3301366 r_SMB .6403225 .0189406 33.81 0.000 .6031996 .6774453 femaler_Mkt -.2447464 .03967 -6.17 0.000 -.3224982 -.1669946 r_Mkt 1.049989 .010206 102.88 0.000 1.029985 1.069992 female .0015689 .0037522 0.42 0.676 -.0057853 .0089232 Ri Coef. Std. Err. z P>|z| [95% Conf. Interval] Log likelihood = 18703.23 Prob > chi2 = 0.0000 Wald chi2(7) = 11535.45 max = 96 avg = 66.77796 Estimated coefficients = 8 Obs per group: min = 1 Estimated autocorrelations = 0 Number of groups = 608 Estimated covariances = 1 Number of obs = 40601 Correlation: no autocorrelation
Panels: homoskedastic
Coefficients: generalized least squares Cross-sectional time-series FGLS regression
Phụ lục 6: Thống kê mơ tả các biến trong mơ hình (nhóm các cơng ty thương mại) _cons .0162777 .0009865 16.50 0.000 .0143443 .0182112 femaler_MOM -.0386319 .1122454 -0.34 0.731 -.2586288 .1813651 r_MOM .0056971 .0291109 0.20 0.845 -.0513592 .0627533 femaler_HML -.1045659 .1211252 -0.86 0.388 -.341967 .1328351 r_HML .2486908 .0314099 7.92 0.000 .1871285 .3102531 femaler_SMB -.4763514 .0735287 -6.48 0.000 -.6204649 -.3322378 r_SMB .6408343 .0191202 33.52 0.000 .6033593 .6783093 femaler_Mkt -.2453157 .039704 -6.18 0.000 -.3231342 -.1674973 r_Mkt 1.050078 .0102162 102.79 0.000 1.030054 1.070101 female .0013917 .0037875 0.37 0.713 -.0060316 .0088151 Ri Coef. Std. Err. z P>|z| [95% Conf. Interval] Log likelihood = 18703.3 Prob > chi2 = 0.0000 Wald chi2(9) = 11535.62 max = 96 avg = 66.77796 Estimated coefficients = 10 Obs per group: min = 1 Estimated autocorrelations = 0 Number of groups = 608 Estimated covariances = 1 Number of obs = 40601 Correlation: no autocorrelation
Panels: homoskedastic
Coefficients: generalized least squares Cross-sectional time-series FGLS regression
MOM 26342 .0372333 .0637792 -.1775545 .308794 HML 26342 -.0438095 .0609975 -.2657585 .129827 SMB 26342 .005329 .0440159 -.1741283 .370661 RmRf 26342 -.0882543 .0878728 -.4745579 .2931591 Female 26342 .0679144 .2516037 0 1 RiRf 26342 -.0813185 .1682819 -.8280522 .9351542 Variable Obs Mean Std. Dev. Min Max . sum RiRf Female RmRf SMB HML MOM
Hồi quy Return với Female, MktRf, Female* MktRf
_cons -.0817205 .0102554 -7.97 0.000 -.1020801 -.061361 Female .0059201 .0038067 1.56 0.123 -.0016371 .0134773 RiRf Coef. Std. Err. t P>|t| [95% Conf. Interval] Robust
(Std. Err. adjusted for 96 clusters in time) Root MSE = .16828 R-squared = 0.0001 Prob > F = 0.1232 F( 1, 95) = 2.42 Linear regression Number of obs = 26342 . reg RiRf Female , vce(cluster time)
_cons .0006065 .0087369 0.07 0.945 -.0167385 .0179514 FemaleRmRf -.1060317 .0490849 -2.16 0.033 -.2034775 -.0085858 RmRf .9322742 .0567159 16.44 0.000 .8196789 1.04487 Female -.0040931 .0060254 -0.68 0.499 -.0160551 .0078689 RiRf Coef. Std. Err. t P>|t| [95% Conf. Interval] Robust
(Std. Err. adjusted for 96 clusters in time) Root MSE = .14731 R-squared = 0.2338 Prob > F = 0.0000 F( 3, 95) = 98.14 Linear reg ression Number of obs = 26342 . reg RiRf Female RmRf FemaleRmRf , vce(cluster time)
Hồi quy Return với Female, MktRf, Female* MktRf, SMB, Female* SMB, HML, Female* HML _cons .0027976 .0088247 0.32 0.752 -.0147216 .0203169 FemaleSMB -.2132292 .0956458 -2.23 0.028 -.4031102 -.0233483 SMB .3848257 .1115562 3.45 0.001 .1633586 .6062928 FemaleRmRf -.1302262 .0537529 -2.42 0.017 -.2369393 -.0235131 RmRf .9803975 .0575338 17.04 0.000 .8661784 1.094617 Female -.0050571 .0061987 -0.82 0.417 -.0173631 .0072488 RiRf Coef. Std. Err. t P>|t| [95% Conf. Interval] Robust
(Std. Err. adjusted for 96 clusters in time) Root MSE = .14645 R-squared = 0.2428 Prob > F = 0.0000 F( 5, 95) = 63.33 Linear regression Number of obs = 26342
_cons .0172342 .0079063 2.18 0.032 .0015381 .0329303 FemaleHML -.0566657 .079089 -0.72 0.475 -.2136772 .1003458 HML .3985172 .1114071 3.58 0.001 .1773461 .6196883 FemaleSMB -.2377079 .1098745 -2.16 0.033 -.4558363 -.0195795 SMB .5881292 .1212185 4.85 0.000 .3474801 .8287784 FemaleRmRf -.1327929 .0534869 -2.48 0.015 -.2389778 -.026608 RmRf .9583894 .050642 18.92 0.000 .8578524 1.058926 Female -.0076745 .0074495 -1.03 0.306 -.0224636 .0071146 RiRf Coef. Std. Err. t P>|t| [95% Conf. Interval] Robust
(Std. Err. adjusted for 96 clusters in time) Root MSE = .14477 R-squared = 0.2601 Prob > F = 0.0000 F( 7, 95) = 57.78 Linear regression Number of obs = 26342 . reg RiRf Female RmRf FemaleRmRf SMB FemaleSMB HML FemaleHML , vce(cluster time)
_cons .0159561 .0077206 2.07 0.041 .0006289 .0312834 FemaleMOM .1215404 .2118156 0.57 0.567 -.2989666 .5420475 MOM -.2221319 .2109242 -1.05 0.295 -.6408693 .1966055 FemaleHML .0612852 .218542 0.28 0.780 -.3725754 .4951459 HML .1817862 .2515557 0.72 0.472 -.3176151 .6811874 FemaleSMB -.2259303 .1197981 -1.89 0.062 -.4637597 .0118991 SMB .5649591 .1305977 4.33 0.000 .3056898 .8242284 FemaleRmRf -.1321766 .0543652 -2.43 0.017 -.2401053 -.024248 RmRf .9563521 .0508009 18.83 0.000 .8554995 1.057205 Female -.007058 .0076724 -0.92 0.360 -.0222895 .0081736 RiRf Coef. Std. Err. t P>|t| [95% Conf. Interval] Robust
(Std. Err. adjusted for 96 clusters in time) Root MSE = .14467 R-squared = 0.2612 Prob > F = 0.0000 F( 9, 95) = 45.34 Linear reg ression Number of obs = 26342
Hồi quy Return với Female
Hồi quy Return với Female, MktRf, Female* MktRf
rho .01072468 (fraction of variance due to u_i)
sigma_e .16875994 sigma_u .01757126 _cons -.0811294 .0134872 -6.02 0.000 -.1075638 -.054695 Female .0070041 .0049666 1.41 0.158 -.0027304 .0167385 RiRf Coef. Std. Err. z P>|z| [95% Conf. Interval] Robust
(Std. Err. adjusted for 96 clusters in time) corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.1585 Wald chi2(1) = 1.99 overall = 0.0001 max = 96 between = 0.0032 avg = 64.1 R-sq: within = 0.0002 Obs per group: min = 1 Group variable: congty Number of groups = 411 Random-effects GLS regression Number of obs = 26342 . xtreg RiRf Female , vce(cluster time) nonest
rho .02232761 (fraction of variance due to u_i) sigma_e .14771917 sigma_u .02232345 _cons .0002032 .0119014 0.02 0.986 -.0231231 .0235295 FemaleRmRf -.1059731 .0494868 -2.14 0.032 -.2029654 -.0089808 RmRf .9332731 .0566231 16.48 0.000 .8222939 1.044252 Female -.0028307 .0074131 -0.38 0.703 -.0173601 .0116986 RiRf Coef. Std. Err. z P>|z| [95% Conf. Interval] Robust
(Std. Err. adjusted for 96 clusters in time) corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 Wald chi2(3) = 297.01 overall = 0.2338 max = 96 between = 0.0138 avg = 64.1 R-sq: wit hin = 0.2340 Obs per group: min = 1 Group variable: congty Number of groups = 411 Random-effects GLS regression Number of obs = 26342 . xtreg RiRf Female RmRf FemaleRmRf , vce(cluster time) nonest
Hồi quy Return với Female, MktRf, Female* MktRf, SMB, Female* SMB, HML, Female* HML
rho .01458373 (fraction of variance due to u_i)
sigma_e .14686302 sigma_u .01786639 _cons .0024909 .0110155 0.23 0.821 -.019099 .0240809 FemaleSMB -.2161178 .0960517 -2.25 0.024 -.4043756 -.02786 SMB .3830954 .1113212 3.44 0.001 .1649099 .6012809 FemaleRmRf -.1305942 .0541367 -2.41 0.016 -.2367002 -.0244883 RmRf .9809646 .0574135 17.09 0.000 .8684363 1.093493 Female -.0043573 .0071373 -0.61 0.542 -.0183462 .0096316 RiRf Coef. Std. Err. z P>|z| [95% Conf. Interval] Robust
(Std. Err. adjusted for 96 clusters in time) corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 Wald chi2(5) = 319.40 overall = 0.2428 max = 96 between = 0.0097 avg = 64.1 R-sq: wit hin = 0.2429 Obs per group: min = 1 Group variable: congty Number of groups = 411
rho .00698458 (fraction of variance due to u_i)
sigma_e .14516912 sigma_u .01217493 _cons .0170084 .0085957 1.98 0.048 .0001612 .0338556 FemaleHML -.0576351 .0793341 -0.73 0.468 -.2131271 .0978569 HML .3988759 .1114896 3.58 0.000 .1803603 .6173916 FemaleSMB -.2400158 .1101836 -2.18 0.029 -.4559717 -.0240599 SMB .5871311 .1208905 4.86 0.000 .3501901 .8240721 FemaleRmRf -.132855 .0537446 -2.47 0.013 -.2381924 -.0275175 RmRf .9588393 .0506848 18.92 0.000 .8594989 1.05818 Female -.0074051 .0078807 -0.94 0.347 -.022851 .0080407 RiRf Coef. Std. Err. z P>|z| [95% Conf. Interval] Robust
(Std. Err. adjusted for 96 clusters in time) corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 Wald chi2(7) = 404.67 overall = 0.2601 max = 96 between = 0.0225 avg = 64.1 R-sq: within = 0.2603 Obs per group: min = 1 Group variable: congty Number of groups = 411 Random-effects GLS regression Number of obs = 26342