Var Obs Mean Stand.Dev Min Max
SIZE 3,728 27.1527 1.5087 22.7801 32.9960 15.13% -1.98% 37.15% 60.11% 38.29% 6.90% -16.42% 25.34% 15.20% -1.42% 36.92% 53.12% 37.32% 7.48% -16.44% 24.74% 2010 2011 2012 2013 2014 2015 2016 2017 CTA2 CTA4
CAPINT2 3,728 .9307 .0438 .6702 .9989
ROA 3,728 .0744 .0798 -.6246 .7499
D/E 3,728 1.7858 2.357 0.01 33.03
STATE 3,728 .2611 .2477 0 .9672
FOREI 3,728 .0918 .1304 0 .6462
Nguồn: Kết quả tổng hợp từ phần mềm Stata - Biến FOREI có mức trung bình là 9.18% thể hiện số cổ phần mà nhà đầu tư nước nắm giữ trong 466 công ty niêm yết trên sàn chứng khốn Việt Nam. Ngồi ra có những cơng ty có mức nắm giữa bởi các cá nhân và tổ chức nước ngoài đạt trên 50% và cao nhất là 64.62%, do ở đây chỉ xét sự tăng lên của mức sở hữu nước ngoài trong khoảng thời gian ngắn từ 2015-2017 và chính phủ chưa thực sự thả lỏng hoàn tồn mà cịn nắm giữ các ngành kinh doanh trọng điểm nên tác nhân sở hữu chưa thực sự nổi trội do đó chưa cần phải tính riêng các trường hợp có mức sở hữu trên 50% hoặc phải thêm biến mới.
- Biến STATE có mức trung bình là 26.11% cho thấy tỷ lệ trung bình mà nhà nước sở hữu trong các công ty niêm yết trên sàn chứng khoán Việt Nam, tỉ lệ sở hữu lớn nhất mà chính phủ giữ tại một doanh nghiệp có thể lên tới 96.72%, tuy nhiên trong thực tế có thể lên tới 100% nhưng do trong giai đoạn xem xét thì chưa thấy xuất hiện tỉ lệ này. Vấn đề này sẽ trở nên quan trọng hơn trong thời gian tới khi chủ trương của nước ta là thoái vốn nhà nước trong giai đoạn 2019-2022 các doanh nghiệp sẽ IPO, niêm yết và tiến hành thoái vốn một cách đáng kể trên thị trường chứng khốn. - Biến SIZE có giá trị trung bình là 27.15 với thang đo logarit tự nhiên điều này cho thấy quy mô doanh nghiệp là lớn.
- Biến ROA có giá trị trung bình là 7.44% cho thấy mức độ tạo ra lợi nhuận của doanh nghiệp đây chính là biến số tác động quyết định đến mức thuế TNDN phải đóng của doanh nghiệp. Giá trị nhỏ nhất là lớn nhất của biến này lần lượt là -62.46%
và 74.99% với giá trị trung bình trên cho thấy mẫu nghiên cứu phân bổ đều cho các loại doanh nghiệp có lợi nhuận thuộc 3 loại cao, trung bình, thấp.
- Biến DE là biến đại diện cho tỷ lệ nợ trên vốn chủ sở hữu có giá trị trung bình là 1.78, cho thấy tỉ lệ nợ trung bình của các doanh nghiệp cao hơn 1.78 lần lượng vốn chủ sở hữu, giá trị cao nhất là 33.03 – đây là do đặc thù của các doanh nghiệp Việt Nam còn đang trong sự bảo trợ của nhà nước và sự yếu kém của doanh nghiệp trong điều hành, giá trị thấp nhất là 0.01 tượng trưng cho doanh nghiệp hầu như khơng vay nợ mà chỉ có các khoản phải thu.
- Biến Cường độ vốn CAPINT2 có giá trị trung bình là 93.07% với cấu trúc hơi đặc biệt là tỷ số của logarit tự nhiên của tài sản cố định chia cho logarit tự nhiên của tổng tài sản, cho thấy mức độ lien quan của biến số này trong do các ưu đãi về thuế đối với chi phí khấu hao tài sản cố định và mức độ đầu tư cho tăng trưởng của doanh nghiệp đến gánh nặng thuế.
4.2. Phân tích tương quan:
Hình 4.4. Phân tích tương quan các biến
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA
4.2.1. Hệ số tương quan Pearson:
Trước khi thực hiện hồi quy, tác giả tiến hành kiểm tra mối tương quan tuyến tính giữa biến phụ thuộc với các biến độc lập. Điều này giúp chọn lựa được những nhân tố thực sự có lien quan để đưa vào mơ hình. Mức dao động của hệ số tương quan là từ -1 đến 1 tương ứng với hai chiều tương quan thuận và nghịch, khi mức tương quan
của hai biến là 0 nghĩa là khơng có sự tương quan giữa hai biến đó. Nếu có mối tương quan chặt và mạnh giữa các biến độc lập thì sẽ có hiện tượng đa cộng tuyến xảy ra.
Hình 4.5. Kết quả tính tốn hệ số tương quan Pearson
Hình 4.6. Kết quả tính tốn hệ số tương quan Pearson (tiếp theo)
Nguồn: Theo tính tốn của tác giả bằng phần mềm STATA Bảng kết quả trên thể hiện các giá trị của hệ số tương quan Pearson, các giá trị này có giá trị từ -1 đến 1, nếu hệ số này bằng 0 (hoặc nằm gần 0) nghĩa là hai biến số khơng có tương quan gì với nhau, cịn nếu gần 1 thì có mối lien hệ với nhau. Theo Guajarati và Porter (2009) thì hệ số tương quan phải nhỏ hơn 0.8 để tránh trường hợp đa cộng tuyến trong mơ hình mà trong thống kê trên các hệ số tương quan đều có giá trị nhỏ hơn 0.8 nên mơ hình ít có khả năng xảy ra đa cộng tuyến.
Từ kiểm định tương quan Pearson cho thấy khơng có sự tương quan trong các biến CTA2 và CTA4 với các biến độc lập STATE và FOREI nên ta khơng đưa ra mơ hình nghiên cứu đối với 2 biến này.
4.3. Phân tích hồi quy và kiểm định:
Để phân tích hồi quy, nghiên cứu này lần lượt thông qua các phương pháp (Pool Ordinary Least Square, Pool OLS), (Fixed Effects Models, FEM) và (Random Effects Model, REM). Sau đó để lựa chọn ra mơ hình phù hợp nhất, tác giả tiến hành kiểm định Hausman.
4.3.1. Phân tích hồi quy mơ hình 1 với biến CTA1:
4.3.1.1. Mơ hình hồi quy Pool OLS:
Hình 4.7. Kết quả hồi quy OLS
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA
Kiểm định tự tương quan:
Sử dụng kiểm định Wooldridge cho mơ hình với giả thiết H0: Khơng có hiện tượng tự tương quan bậc nhất.
Hình 4.8. Kết quả kiểm định tự tương quan
Nguồn: Tổng hợp từ phần mềm STATA của tác giả Kết quả trên cho thấy p-value = 0.8852 lớn hơn 5% nên khơng có hiện tượng tự tương quan bậc nhất trong mơ hình.
4.3.1.2. Mơ hình hồi quy FEM:
Hình 4.9. Bảng hồi quy mơ hình FEM
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Để lựa chọn giữa mơ hình FEM và OLS, từ F = 1.18 > 1.13 nên bác bỏ giả thiết H0 (H0: nên chọn theo mơ hình OLS) điều này cho thấy mơ hình FEM tốt hơn mơ hình Pool OLS.
4.3.1.3. Mơ hình hồi quy REM:
Hình 4.10. Kết quả hồi quy của mơ hình REM cho biến CTA1
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Tiếp theo tác giả tiến hành kiểm định Hausman để lựa chọn giữa mơ hình FEM và REM.
Hình 4.11. Kết quả kiểm định Hausman của biến CTA1
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Kết quả từ bảng trên Prob>chi2 = 0.2051 > 5% nên chấp nhận giải thiết H0 (H0: nên chọn mơ hình REM). Như vậy phương pháp ước lượng phù hợp là REM. Như vậy tác giả sẽ chọn mơ hình tác động ngẫu nhiên để nghiên cứu và tiếp tục tiến hành kiểm định các khuyết tật của mơ hình như là tự tương quan, đa cộng tuyến, phương sai thay đổi và phần dư không phân phối chuẩn.
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA
Kiểm định phương sai thay đổi: tác giả sử dụng kiểm định Breusch – Pagan cho mơ hình REM với giả thiết H0: khơng có hiện tượng phương sai thay đổi, nếu Prob Chi2 < 0.05 thì có hiện tượng phương sai thay đổi.
Hình 4.12. Kết quả kiểm định phương sai thay đổi của CTA1
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Kết quả cho thấy Prob Chi2 <0.05 nên bác bỏ giả thiết H, vậy có hiện tượng phương sai sai số thay đổi trong mơ hình.
Khi bị hiện tượng phương sai thay đổi, ước lượng tính bằng phương pháp REM khơng cịn là ước lượng hiệu quả nữa, các kiểm định về hệ số hồi quy được tính tốn cũng khơng cịn đáng tin cậy. Để xử lý hiện tượng phương sai thay đổi này bằng cách dung ma trận ước lượng vững của hiệp phương sai hoặc tìm ra ước lượng hiệu quả hơn bằng phương pháp GLS (còn gọi là phương pháp WLS)
Kiểm định hiện tượng đa cộng tuyến:
Hình 4.13. Kết quả kiểm định đa cộng tuyến của CTA1
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Kết quả trong bảng trên cho thấy trị trung bình Mean VIF = 1.16 < 10 nên khơng có hiện tượng đa cộng tuyến trong mơ hình.
Kiểm định phân bố chuẩn của phần dư:
Có nhiều nguyên nhân dẫn đến phân phối của phần dư không phải là phân phối chuẩn như sai dạng hàm, số lượng phần dư không đủ để tiến hành phân tích,…
Hình 4.14. Kiểm định phân bố chuẩn của CTA1
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Kết quả kiểm định theo CTA1 cho thấy các p-value của kiểm định Skewness và Kurtosis đều có p-value <0.05, suy ra phần dư của mơ hình khơng có phân phối chuẩn. (Chỉ khi p-value của Skewness>0.05 và Kurtosis p<0.05 thì biến kiểm định mới có phân phối chuẩn
Hồi quy Mơ hình theo phương pháp GLS:
Từ các kết quả kiểm tra các khuyết tật của mơ hình như kiểm định tự tương quan, kiểm định phương sai thay đổi, kiểm định đa cộng tuyến và kiểm định phân phối chuẩn của phần dư cho thấy mơ hình có hiện tượng phương sai thay đổi và phần dư khơng có phân phối chuẩn. Theo Wooldridge (2002) thì cách khắc phục khi xảy ra hai hiện tượng này là dùng ước lượng hồi quy theo phương pháp bình phương nhỏ nhất tổng quát (GLS Generalize Least Square). Theo đó các biến mới được biến đổi từ mơ hình vi phạm các giả thuyết Gauss-Markov sang mơ hình mới thoả các giả thiết này dẫn đến các tham số ước lượng thu được sẽ tin cậy hơn.
Hình 4.15. Kết quả ước lượng GLS của CTA1
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA
Hình 4.16. Bảng tổng hợp hệ số của các phương pháp GLS
Hình 4.17. Bảng mơ tả hệ số ước lượng của 4 phương pháp
(Ghi chú: *, ** và *** tương ứng ở mức 10%, 5% và 1%)
Kết quả mơ hình hồi quy:
CTA1 = 0.384 + 0.0241*STATE - 0.0201*FOREI
- 0.0025*SIZE + 0.0033*DE – 0.1101*CAPINT2 – 0.1177*ROA
Trong mơ hình 1, tỷ lệ sở hữu nhà nước và tỷ lệ nợ trên vốn chủ sở hữu tác động cùng chiều với tỷ lệ đóng thuế TNDN, ngược lại tỷ lệ sở hữu nước ngồi, quy mơ doanh nghiệp và cường độ vốn và tỷ lệ lợi nhuận trên tổng tài sản tác động ngược chiều đến tỷ lệ đóng thuế TNDN.
4.3.2. Phân tích hồi quy mơ hình 2 với biến CTA3:
4.3.2.1. Mơ hình Pool OLS:
Hình 4.18. Kết quả hồi quy mơ hình OLS biến CTA3
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Tác giả cũng tiến hành ước lượng bằng phương pháp Pool OLS.
Hình 4.19. Kết quả kiểm định tự tương quan
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Từ kiểm định Wooldridge trên cho thấy Prob > 0.05 khơng bác bỏ giả thiết H0 (H0: Khơng có hiện tượng tự tương quan trong mơ hình).
4.3.2.2. Mơ hình FEM:
Hình 4.20. Kết quả ước lượng mơ hình FEM biến CTA3
Để lựa chọn giữa mơ hình FEM và OLS, từ F = 1.24 > 1.13 nên bác bỏ giả thiết H0 (H0: nên chọn theo mơ hình OLS) điều này cho thấy mơ hình FEM tốt hơn mơ hình Pool OLS.
4.3.2.3. Mơ hình REM:
Hình 4.21. Kết quả hồi quy mơ hình REM
Kiểm định Hausman:
Hình 4.22. Kết quả kiểm định Hausman
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Kết quả từ bảng trên Prob > chi2 = 0.2051 > 5% nên chấp nhận giải thiết H0 (H0: nên chọn mơ hình REM). Như vậy phương pháp ước lượng phù hợp là REM. Như vậy tác giả sẽ chọn mơ hình tác động ngẫu nhiên để nghiên cứu và tiếp tục tiến hành kiểm định các khuyết tật của mơ hình như là tự tương quan, đa cộng tuyến, phương sai thay đổi và phần dư không phân phối chuẩn.
Kiểm định phương thay đổi:
Hình 4.23. Kết quả kiểm định phương sai
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA
Kiểm định phương sai thay đổi: tác giả sử dụng kiểm định Breusch – Pagan cho mơ hình REM với giả thiết H0: khơng có hiện tượng phương sai thay đổi, nếu Prob Chi2 < 0.05 thì có hiện tượng phương sai thay đổi. Vậy theo kết quả trên cho kết luận cũng có hiện tượng phương sai sai số thay đổi trong mơ hình như với biến CTA1.
Kiểm định phân bố chuẩn phần dư:
Hình 4.24. Kết quả kiểm định phân phối chuẩn của phần dư
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA
Từ kết quả kiểm định cho thấy p-value của kiểm định Skewness và Kurtosis đều cho kết quả p < 0.05 do đó phần dư của biến CTA3 có phân phối khơng chuẩn. (Chỉ khi p-value của Skewness > 0.05 và Kurtosis p < 0.05 thì biến kiểm định mới có phân phối chuẩn)
Kiểm định đa cộng tuyến:
Hình 4.25. Kết quả kiểm định đa cộng tuyến biến CTA3
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Mean VIF = 1.16 < 10 khơng có hiện tượng đa cộng tuyến
4.3.2.4. Mơ hình GLS:
Hình 4.26. Kết quả ước lượng theo GLS của biến CTA3
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA
Hình 4.27. Bảng tổng kết các phương pháp Pool OLS, FEM, REM và GLS
Hình 4.28. Kết quả đánh giá theo CTA3
Nguồn: Theo tính tốn của tác giả từ phần mềm STATA Kết quả mơ hình hồi quy:
CTA3 = 0.420 + 0.0317*STATE - 0.0241*FOREI
- 0.00417*SIZE + 0.00224*DE – 0.102*CAPINT2 – 0.118*ROA
Trong mơ hình 2, mối lien hệ giữa các biến cũng xảy ra như với mơ hình 1 là tỷ lệ sở hữu nhà nước và tỷ lệ nợ trên vốn chủ sở hữu tác động cùng chiều với tỷ lệ đóng thuế TNDN, ngược lại tỷ lệ sở hữu nước ngồi, quy mơ doanh nghiệp và cường độ vốn và tỷ lệ lợi nhuận trên tổng tài sản tác động ngược chiều đến tỷ lệ đóng thuế TNDN.
4.4. Thảo luận kết quả nghiên cứu:
Vì dạng dấu của 2 phương trình hồi quy của 2 biến CTA1 và CTA3 là như nhau chỉ khác nhau về giá trị nên ở đây ta chỉ làm cho biến CTA1 cịn biến CTA3 thì hồn tồn tương tự.