Kiểm định ảnh hưởng của biến động lên chênh lệch sở hữu nước ngoài

Một phần của tài liệu (LUẬN văn THẠC sĩ) sở hữu nước ngoài và biến động giá chứng khoán, bằng chứng thực nghiệm tại việt nam (Trang 43)

CHƯƠNG 4 NỘI DUNG VÀ KẾT QUẢ NGHIÊN CỨU

4.2 Kết quả nghiên cứu

4.2.1. Kiểm định ảnh hưởng của biến động lên chênh lệch sở hữu nước ngoài

Tác giả ước lượng mơ hình bình phương tối thiểu thơng thường (OLS)

_cons -1.50852 .1881937 -8.02 0.000 -1.877384 -1.139656 r1 .1738834 .0681248 2.55 0.011 .0403572 .3074096 r .0952601 .0668735 1.42 0.154 -.0358136 .2263339 toveri .0075517 .015846 0.48 0.634 -.0235069 .0386103 lnmacp .1466138 .0180061 8.14 0.000 .1113214 .1819062 lnst .0217067 .0272906 0.80 0.426 -.0317836 .075197 fh1 -.0183043 .0010291 -17.79 0.000 -.0203213 -.0162874 deltafh Coef. Std. Err. t P>|t| [95% Conf. Interval] Total 188166.406 39223 4.79734864 Root MSE = 2.1813 Adj R-squared = 0.0082 Residual 186593.435 39217 4.7579732 R-squared = 0.0084 Model 1572.97043 6 262.161739 Prob > F = 0.0000 F( 6, 39217) = 55.10 Source SS df MS Number of obs = 39224 . reg deltafh fh1 lnst lnmacp toveri r r1

Kiểm định VIF cho kết quả là 1.07 nên vấn đề đa cộng tuyến khơng ảnh hưởng đến kết quả mơ hình.

Kiểm tra tự tương quan cho thấy hiện tượng tự tương quan

Kiểm định phương sai thay đổi

Tác giả kiểm tra hiện tượng phương sai không đồng nhất với giả thiết :

Ho: Khơng có hiện tượng phương sai không đồng nhất Giả thiết đối : có hiện tượng phương sai khơng đồng nhất

Nguyên tắc bác bỏ là nếu P-value của chi-sq thấp hơn mức tin cậy của nghiên cứu là 5% thì bác bỏ Ho.

Kết quả kiểm định Prob = 0.000 nên bác bỏ Ho, tức mơ hình có hiện tượng phương sai thay đổi.

Prob > F = 0.0354 F( 1, 708) = 4.441 H0: no first-order autocorrelation

Wooldridge test for autocorrelation in panel data . xtserial deltafh fh1 lnst lnmacp toveri r r1

Prob > chi2 = 0.0000 chi2(1) = 6389.38

Variables: fitted values of fh1 Ho: Constant variance

Breusch-Pagan / Cook-Weisberg test for heteroskedasticity . hettest Mean VIF 1.07 toveri 1.00 0.996686 r 1.00 0.996484 r1 1.02 0.977948 lnmacp 1.14 0.875572 lnst 1.16 0.858455 Variable VIF 1/VIF . vif

Tác giả ước lượng mơ hình bình phương tối thiểu thơng thường (OLS) khơng có hiện tương đa cộng tuyến và tự tương quan nhưng vi phạm giả thiết phương sai do đó tác giả sử dụng mơ hình bình phương tối thiểu tổng qt (GLS).

a) Kiểm định các mơ hình hồi quy FEM – REM và GLS

Ước lượng mơ hình (1) theo phương pháp GLS kết quả như sau :

Ước lượng mơ hình (1) theo phương pháp FEM có robust kết quả như sau : _cons -1.50852 .1881769 -8.02 0.000 -1.87734 -1.1397 r1 .1738834 .0681187 2.55 0.011 .0403732 .3073936 r .0952601 .0668675 1.42 0.154 -.0357978 .2263181 toveri .0075517 .0158446 0.48 0.634 -.0235031 .0386066 lnmacp .1466138 .0180045 8.14 0.000 .1113257 .181902 lnst .0217067 .0272882 0.80 0.426 -.0317772 .0751906 fh1 -.0183043 .001029 -17.79 0.000 -.0203211 -.0162876 deltafh Coef. Std. Err. z P>|z| [95% Conf. Interval] Log likelihood = -86244.17 Prob > chi2 = 0.0000 Wald chi2(6) = 330.66 max = 65 avg = 51.95232 Estimated coefficients = 7 Obs per group: min = 2 Estimated autocorrelations = 0 Number of groups = 755 Estimated covariances = 1 Number of obs = 39224 Correlation: no autocorrelation

Panels: homoskedastic

Coefficients: generalized least squares Cross-sectional time-series FGLS regression . xtgls deltafh fh1 lnst lnmacp toveri r r1

Ước lượng mơ hình (1) theo phương pháp REM có robust kết quả như sau :

.

rho .15287882 (fraction of variance due to u_i)

sigma_e 2.1533123 sigma_u .91476157 _cons -1.389836 .6953743 -2.00 0.046 -2.754936 -.0247358 r1 .2002042 .0792528 2.53 0.012 .0446219 .3557865 r .0826925 .0707864 1.17 0.243 -.0562693 .2216543 toveri .0051506 .0089392 0.58 0.565 -.0123981 .0226993 lnmacp .2201661 .0673874 3.27 0.001 .0878769 .3524553 lnst -.1287629 .0363234 -3.54 0.000 -.2000699 -.0574559 fh1 -.0873258 .007353 -11.88 0.000 -.1017605 -.072891 deltafh Coef. Std. Err. t P>|t| [95% Conf. Interval] Robust

(Std. Err. adjusted for 755 clusters in code) corr(u_i, Xb) = -0.8664 Prob > F = 0.0000 F(6,754) = 26.37 overall = 0.0069 max = 65 between = 0.0001 avg = 52.0 R-sq: within = 0.0406 Obs per group: min = 2 Group variable: code Number of groups = 755 Fixed-effects (within) regression Number of obs = 39224 . xtreg deltafh fh1 lnst lnmacp toveri r r1, fe robust

rho 0 (fraction of variance due to u_i)

sigma_e 2.1533123 sigma_u 0 _cons -1.50852 .2416906 -6.24 0.000 -1.982225 -1.034815 r1 .1738834 .0783282 2.22 0.026 .0203629 .3274038 r .0952601 .0699677 1.36 0.173 -.0418739 .2323942 toveri .0075517 .0169636 0.45 0.656 -.0256963 .0407998 lnmacp .1466138 .0231916 6.32 0.000 .1011592 .1920685 lnst .0217067 .0260181 0.83 0.404 -.0292879 .0727013 fh1 -.0183043 .0022917 -7.99 0.000 -.022796 -.0138126 deltafh Coef. Std. Err. z P>|z| [95% Conf. Interval] Robust

(Std. Err. adjusted for 755 clusters in code) corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 Random effects u_i ~ Gaussian Wald chi2(6) = 85.77 overall = 0.0084 max = 65 between = 0.0016 avg = 52.0 R-sq: within = 0.0323 Obs per group: min = 2 Group variable: code Number of groups = 755 Random-effects GLS regression Number of obs = 39224 . xtreg deltafh fh1 lnst lnmacp toveri r r1, re robust

Bảng 4 : Kiểm định ảnh hưởng của biến động giá chứng khoán lên chênh lệch sở hữu của nhà đầu tư nước ngồi theo mơ hình (1)

Variable FEM REM GLS

fh1 -0.08732576*** -0.01830434*** -0.01830434*** lnst - 0.12876287*** 0.02170668 0.02170668 lnmacp 0.22016612** 0.14661383*** 0.14661383*** toveri 0.00515061 0.00755173 0.00755173 r 0.08269249 0.09526015 0.09526015 r1 0 .2002042* 0.1738834* 0.1738834* _cons -1.3898357* -1.5085201*** -1.5085201*** N 39223 39223 39223 r2 .04063244 r2_a .04048566

( Nguồn : từ số liệu tác tính tốn với mức ý nghĩa * p<.05; ** p<.01; *** p<.001)

Kiểm định phương trình (1) với 3 mơ hình, mơ hình REM và mơ hình GLS có kết quả ước lượng giống nhau. Từ kết quả trên tác giả tiếp tục thực hiện các bước kiểm tra tiếp theo.

b) Kiểm định Hausman

Lựa chọn mơ hình hồi quy thích hợp REM hay FEM, tác giả thực hiện kiểm định Hausman test với giả thiết :

Ho: REM là mơ hình thích hợp hơn FEM

Giả thiết đối : FEM là mơ hình thích hợp hơn REM

Nguyên tắc bác bỏ là nếu P-value của chi-sq thấp hơn mức tin cậy của nghiên cứu là 5% thì bác bỏ Ho.

Kết quả kiểm định Prob = 0.000 nên bác bỏ Ho, tức mơ hình FEM là mơ hình phù hợp hơn. Do mơ hình GLS và REM có kết quả ước lượng giống nhau do đó mơ hình phù hợp nhất của 3 mơ hình REM, FEM và GLS là mơ hình FEM.

c) Kiểm định phương sai (heteroskedastic)

Để có mơ hình FEM phù hợp nhất tác giả kiểm tra hiện tượng phương sai không đồng nhất với giả thiết :

Ho: Khơng có hiện tượng phương sai khơng đồng nhất Giả thiết đối : có hiện tượng phương sai khơng đồng nhất

Nguyên tắc bác bỏ là nếu P-value của chi-sq thấp hơn mức tin cậy của nghiên cứu là 5% thì bác bỏ Ho.

Kết quả kiểm định Prob = 0.000 nên bác bỏ Ho, tức mơ hình FEM là mơ hình có hiện tượng phương sai thay đổi. Để kiểm sốt tình trạng phương sai không đồng nhất tác giả đưa thêm kiểm sốt robust vào trong mơ hình FEM.

(V_b-V_B is not positive definite) Prob>chi2 = 0.0000

= 1308.17

chi2(6) = (b-B)'[(V_b-V_B)^(-1)](b-B) Test: Ho: difference in coefficients not systematic

B = inconsistent under Ha, efficient under Ho; obtained from xtreg b = consistent under Ho and Ha; obtained from xtreg r1 .2002042 .1738834 .0263208 . r .0826925 .0952601 -.0125677 . toveri .0051506 .0075517 -.0024011 .0046879 lnmacp .2201661 .1466138 .0735523 .0464191 lnst -.1287629 .0217067 -.1504695 .0198602 fh1 -.0873258 -.0183043 -.0690214 .0019118 fixed random Difference S.E.

(b) (B) (b-B) sqrt(diag(V_b-V_B)) Coefficients

. hausman fixed random .

d) Kiểm định đa cộng tuyến

Với kết quả kiểm tra cho thấy mơ hình khơng xảy ra hiện tượng đa cộng tuyến và các biến giải thích khơng tương quan nhau. Do đó mơ hình là phù hợp

Kết luận :

Tác giả ước tính tỷ lệ sở hữu nước ngoài và biến động bị ảnh hưởng bởi các yếu tố khác như vốn chủ sở hữu nước ngồi, lợi nhuận chứng khốn, doanh thu chứng

khốn có thể ảnh hưởng đến sở hữu nước ngồi. Biến kiểm sốt ln (MACPi,t),

TOVERi,t,ri,t,ri,t-1 đưa vào phương trình kiểm sốt các yếu tố này. Độ trễ biến phụ thuộc

FHi,t-1 được đưa vào để kiểm sốt tính bền vững trong quyền sở hữu nước ngoài.

Sau khi ước lượng mơ hình theo các phương pháp hồi quy REM – FEM – GLS và kiểm tra các hiện tượng phương sai, tự tương quan và đa cộng tuyến, kết quả ước lượng mơ hình theo phương pháp FEM ( có tuỳ chọn robust để xử lý phương sai thay đổi), kết q của mơ hình (1) có mối quan hệ giữa biến động chứng khoán và sở hữu nước ngồi do đó tác giả bác bỏ giả thiết Ho.

Prob>chi2 = 0.0000 chi2 (755) = 3.7e+09

H0: sigma(i)^2 = sigma^2 for all i in fixed effect regression model

Modified Wald test for groupwise heteroskedasticity . xttest3 Mean VIF 1.07 toveri 1.00 0.996686 r 1.00 0.996484 r1 1.02 0.977948 lnmacp 1.14 0.875572 lnst 1.16 0.858455 Variable VIF 1/VIF . vif

Mô hình ước lượng (1)

∆FHi,t = – 1.3898– 0.87FHi,t-1– 0.128ln(i,t) + 0.22ln (MACPi,t) + 0.005TOVERi,t +

0.08ri,t + 0.2ri,t-1 + εi,t.

Với R2 = 4.06% (FE)

 Biến độc lập FHi,t-1 , ln(i,t) có dấu ngược chiều với biến phụ thuộc ∆FHi,t với

mức ý nghĩa là 10%. Tức biến FHi,t-1 và ln(i,t) làm giảm chênh lệch mức sở hữu của

nhà đầu tư nước ngồi.

 Biến kiểm sốt ln (MACPi,t),TOVERi,t,ri,t,ri,t-1đều có dấu cùng chiều với biến

phụ thuộc. Biến ln (MACPi,t), ri,t, đều có ý nghĩa thống kê với mức 5% và 1% còn 2

biến cịn lại khơng có ý nghĩa thống kê. Các biến kiểm sốt đồng thời có ảnh hưởng đến biến phụ thuộc cho thấy quyết định đầu tư chịu ảnh hưởng của đặc điểm riêng của từng công ty. Mức ảnh hưởng lớn nhất là lợi nhuận chứng khoán trong quá khứ và mức vốn hoá thị trường.

 Vốn hố thị trường có ảnh hưởng lớn nhất trong phương trình hồi quy điều này

cho thấy đúng với thống kê mô tả.

 R2 là 4.06% cho thấy mức độ giải thích của mơ hình là 4.06%, nghĩa là những

thay đổi trong các biến trên chỉ giải thích được 4.06% sự thay đổi của sở hữu nước ngoài.

4.2.2. Kiểm định sở thích của nhà đầu tư nước ngoài đối với tính biến động của cổ phiếu.

Kiểm định độ ảnh hưởng của các biến trễ ln(i,t-1), ln(i,t-2),ln (MACPi,t-1), ri,t-

1, FHi,t-1 đến biến động ln(i,t). Tác giả ước lượng mơ hình (2) hồi quy theo mơ hình FEM, REM và GLS để lựa chọn mơ hình thích hợp nhất.

Ước lượng hệ số hồi quy theo phương pháp FEM ( có tùy chọn robust), kết quả như sau:

Ước lượng hệ số hồi quy theo phương pháp REM ( có tùy chọn robust), kết quả như sau :

lnmacp1 -.1096337 .0041085 -26.68 0.000 -.1176862 -.1015813 lnst2 .0584645 .0037567 15.56 0.000 .0511016 .0658275 lnst1 .3846798 .0142371 27.02 0.000 .3567756 .412584 lnst Coef. Std. Err. z P>|z| [95% Conf. Interval] Robust

(Std. Err. adjusted for 755 clusters in code) corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 Random effects u_i ~ Gaussian Wald chi2(5) = 1481.62 overall = 0.3071 max = 65 between = 0.5297 avg = 52.0 R-sq: within = 0.2285 Obs per group: min = 2 Group variable: code Number of groups = 755 Random-effects GLS regression Number of obs = 39224 . xtreg lnst lnst1 lnst2 lnmacp1 r1 fh1, re robust

.

. estimates store fixed .

rho .25764315 (fraction of variance due to u_i)

sigma_e .32456342 sigma_u .19120644 _cons 2.820181 .015805 178.44 0.000 2.789154 2.851208 fh1 -.0011587 .0006403 -1.81 0.071 -.0024156 .0000983 r1 .3195607 .0114028 28.02 0.000 .2971758 .3419457 lnmacp1 -.1078551 .0041936 -25.72 0.000 -.1160877 -.0996226 lnst2 .0545401 .0037971 14.36 0.000 .0470859 .0619944 lnst1 .3732622 .0143431 26.02 0.000 .3451051 .4014194 lnst Coef. Std. Err. t P>|t| [95% Conf. Interval] Robust

(Std. Err. adjusted for 755 clusters in code) corr(u_i, Xb) = 0.1618 Prob > F = 0.0000 F(5,754) = 285.28 overall = 0.2788 max = 65 between = 0.4206 avg = 52.0 R-sq: within = 0.2301 Obs per group: min = 2 Group variable: code Number of groups = 755 Fixed-effects (within) regression Number of obs = 39224 . xtreg lnst lnst1 lnst2 lnmacp1 r1 fh1, fe robust

44

Ước lượng hệ số hồi quy theo phương pháp GLS, kết quả như sau

. estimates store random .

rho .11643937 (fraction of variance due to u_i)

sigma_e .32456342 sigma_u .11782334 _cons 2.773011 .0158413 175.05 0.000 2.741963 2.804059 fh1 .0015885 .0004814 3.30 0.001 .0006449 .0025321 r1 .3236646 .011344 28.53 0.000 .3014308 .3458984 lnmacp1 -.1096337 .0041085 -26.68 0.000 -.1176862 -.1015813 lnst2 .0584645 .0037567 15.56 0.000 .0511016 .0658275 lnst1 .3846798 .0142371 27.02 0.000 .3567756 .412584 lnst Coef. Std. Err. z P>|z| [95% Conf. Interval] Robust

(Std. Err. adjusted for 755 clusters in code) corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 Random effects u_i ~ Gaussian Wald chi2(5) = 1481.62 overall = 0.3071 max = 65 between = 0.5297 avg = 52.0 R-sq: within = 0.2285 Obs per group: min = 2 Group variable: code Number of groups = 755 Random-effects GLS regression Number of obs = 39224 . xtreg lnst lnst1 lnst2 lnmacp1 r1 fh1, re robust

.

. estimates store fixed .

rho .25764315 (fraction of variance due to u_i)

sigma_e .32456342 sigma_u .19120644 _cons 2.820181 .015805 178.44 0.000 2.789154 2.851208 fh1 -.0011587 .0006403 -1.81 0.071 -.0024156 .0000983 r1 .3195607 .0114028 28.02 0.000 .2971758 .3419457 lnmacp1 -.1078551 .0041936 -25.72 0.000 -.1160877 -.0996226 lnst2 .0545401 .0037971 14.36 0.000 .0470859 .0619944 lnst1 .3732622 .0143431 26.02 0.000 .3451051 .4014194 lnst Coef. Std. Err. t P>|t| [95% Conf. Interval] Robust

(Std. Err. adjusted for 755 clusters in code) corr(u_i, Xb) = 0.1618 Prob > F = 0.0000 F(5,754) = 285.28 overall = 0.2788 max = 65 between = 0.4206 avg = 52.0 . estimates store gls . _cons 2.553778 .0121842 209.60 0.000 2.529897 2.577659 fh1 .0049467 .0001573 31.45 0.000 .0046384 .005255 r1 .3422167 .0110985 30.83 0.000 .320464 .3639695 lnmacp1 -.1104031 .0015099 -73.12 0.000 -.1133624 -.1074438 lnst2 .0794305 .00328 24.22 0.000 .0730019 .0858591 lnst1 .4403767 .0042636 103.29 0.000 .4320201 .4487332 lnst Coef. Std. Err. z P>|z| [95% Conf. Interval] Log likelihood = -15420.01 Prob > chi2 = 0.0000 Wald chi2(5) = 18692.39 max = 65 avg = 51.95232 Estimated coefficients = 6 Obs per group: min = 2 Estimated autocorrelations = 0 Number of groups = 755 Estimated covariances = 1 Number of obs = 39224 Correlation: no autocorrelation

Panels: homoskedastic

Coefficients: generalized least squares Cross-sectional time-series FGLS regression . xtgls lnst lnst1 lnst2 lnmacp1 r1 fh1

45

Sau khi đó tác giả sử dụng test Breusch và Pagan để lựa chọn giữa mơ hình REM và GLS. Kết quả cho thấy mơ hình REM thích hợp hơn mơ hình GLS.

Test Hausman dùng để lựa chọn giữa hai mơ hình FEM và REM. Kết quả cho thấy mơ hình FEM robust thích hợp nhất

Các bước kiểm định lựa chọn mơ hình hồi quy phù hợp tương tự mơ hình (1). Sau khi thực hiện các bước kiểm định, tác giả lựa chọn mơ hình FEM ( có tùy chọn robust) là mơ hình hồi quy để thực hiện kiểm định. Tác giả đưa thêm hai biến kiểm

soát TOVERi,t và ri,t vào mơ hình đã lựa chọn.

Prob > chi2 = 0.0000 chi2(1) = 22754.06 Test: Var(u) = 0 u .0138823 .1178233 e .1053414 .3245634 lnst .1897807 .4356383 Var sd = sqrt(Var) Estimated results:

lnst[code,t] = Xb + u[code] + e[code,t]

Breusch and Pagan Lagrangian multiplier test for random effects . xttest0

.

rho .11643937 (fraction of variance due to u_i)

sigma_e .32456342 sigma_u .11782334 _cons 2.773011 .0158413 175.05 0.000 2.741963 2.804059 fh1 .0015885 .0004814 3.30 0.001 .0006449 .0025321 r1 .3236646 .011344 28.53 0.000 .3014308 .3458984 lnmacp1 -.1096337 .0041085 -26.68 0.000 -.1176862 -.1015813 lnst2 .0584645 .0037567 15.56 0.000 .0511016 .0658275 lnst1 .3846798 .0142371 27.02 0.000 .3567756 .412584 lnst Coef. Std. Err. z P>|z| [95% Conf. Interval] Robust

(Std. Err. adjusted for 755 clusters in code) corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 Random effects u_i ~ Gaussian Wald chi2(5) = 1481.62 overall = 0.3071 max = 65 between = 0.5297 avg = 52.0

(V_b-V_B is not positive definite) Prob>chi2 = 0.0000

= 50016.45

chi2(5) = (b-B)'[(V_b-V_B)^(-1)](b-B) Test: Ho: difference in coefficients not systematic

B = inconsistent under Ha, efficient under Ho; obtained from xtreg b = consistent under Ho and Ha; obtained from xtreg fh1 -.0011587 .0015885 -.0027472 .0001955 r1 .3195607 .3236646 -.0041039 . lnmacp1 -.1078551 -.1096337 .0017786 .0002483 lnst2 .0545401 .0584645 -.0039244 . lnst1 .3732622 .3846798 -.0114176 .0003516 fixed random Difference S.E.

(b) (B) (b-B) sqrt(diag(V_b-V_B)) Coefficients

Tiếp tục kiểm định độ ảnh hưởng của biến sở hữu nhà đầu tư nước ngồi đến

tính biến động của cổ phiếu tác giả đưa thêm hai biến ∆FHi,t và FTi,t.

legend: * p<.05; ** p<.01; *** p<.001 _cons 2.8173778*** fh1 -.00109301 r .10544278*** toveri .01924434*** r1 .31314205*** lnmacp1 -.1082065*** lnst2 .05387931*** lnst1 .3753602*** Variable fixed

. estimates table fixed, star(.05 .01 .001) . estimates store fixed

r(111);

estimation result fixed not found

. estimates table fixed, star(.05 .01 .001)

rho .25666839 (fraction of variance due to u_i)

sigma_e .32384014 sigma_u .19029421 _cons 2.817378 .0158183 178.11 0.000 2.786325 2.848431 fh1 -.001093 .0006324 -1.73 0.084 -.0023344 .0001484 r .1054428 .0135007 7.81 0.000 .0789393 .1319463 toveri .0192443 .0057038 3.37 0.001 .0080472 .0304415 r1 .3131421 .0112882 27.74 0.000 .2909819 .3353022 lnmacp1 -.1082065 .0042334 -25.56 0.000 -.1165172 -.0998958 lnst2 .0538793 .0038057 14.16 0.000 .0464082 .0613504 lnst1 .3753602 .0145174 25.86 0.000 .346861 .4038594 lnst Coef. Std. Err. t P>|t| [95% Conf. Interval] Robust

(Std. Err. adjusted for 755 clusters in code) corr(u_i, Xb) = 0.1633 Prob > F = 0.0000 F(7,754) = 222.78 overall = 0.2831 max = 65 between = 0.4267 avg = 52.0 R-sq: within = 0.2336 Obs per group: min = 2 Group variable: code Number of groups = 755 Fixed-effects (within) regression Number of obs = 39224 . xtreg lnst lnst1 lnst2 lnmacp1 r1 toveri r fh1, fe robust

Một phần của tài liệu (LUẬN văn THẠC sĩ) sở hữu nước ngoài và biến động giá chứng khoán, bằng chứng thực nghiệm tại việt nam (Trang 43)

Tải bản đầy đủ (PDF)

(91 trang)