Bài giảng Kinh tế lượng - Phần 4

CHƯƠNG 4: HỒI QUI ĐA BIẾN

Mô hình hồi quy đơn đã trình bày ở các chương 2 và 3 là khá hữu dụng cho rất nhiều

trường hợp khác nhau. Mặc dù vậy, nó trở nên không còn phù hợp nữa khi có nhiều hơn

một yếu tố tác động đến biến cần được giải thích. Hồi quy đa biến cho phép chúng ta

nghiên cứu những trường hợp như vậy. Hãy xét các ví dụ sau:

4.1 Giới thiệu về hồi quy đa biến

Ví dụ 4.1: Rất nhiều các nghiên cứu trên thế giới quan tâm tới mối quan hệ giữa thu nhập

và trình độ học vấn. Chúng ta kỳ vọng rằng, ít ra về trung bình mà nói, học vấn càng cao,

thì thu nhập càng cao. Vì vậy, chúng ta có thể lập phương trình hồi quy sau:

pdf12 trang | Chia sẻ: thanhthanh29 | Lượt xem: 716 | Lượt tải: 0download
Bạn đang xem nội dung tài liệu Bài giảng Kinh tế lượng - Phần 4, để tải tài liệu về máy bạn click vào nút DOWNLOAD ở trên
Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-1 CHƯƠNG 4: HỒI QUI ĐA BIẾN Mô hình hồi quy đơn đã trình bày ở các chương 2 và 3 là khá hữu dụng cho rất nhiều trường hợp khác nhau. Mặc dù vậy, nó trở nên không còn phù hợp nữa khi có nhiều hơn một yếu tố tác động đến biến cần được giải thích. Hồi quy đa biến cho phép chúng ta nghiên cứu những trường hợp như vậy. Hãy xét các ví dụ sau: 4.1 Giới thiệu về hồi quy đa biến Ví dụ 4.1: Rất nhiều các nghiên cứu trên thế giới quan tâm tới mối quan hệ giữa thu nhập và trình độ học vấn. Chúng ta kỳ vọng rằng, ít ra về trung bình mà nói, học vấn càng cao, thì thu nhập càng cao. Vì vậy, chúng ta có thể lập phương trình hồi quy sau: Thu nhập = 21 ββ + Học vấn ε+ Tuy nhiên, mô hình này đã bỏ qua một yếu tố khá quan trọng là mọi người thường có mức thu nhập cao hơn khi họ làm việc lâu năm hơn, bất kể trình độ học vấn của họ thế nào. Vậy nên, mô hình tốt hơn cho mục đích nghiên cứu của chúng ta sẽ là: Thu nhập = 21 ββ + Học vấn 3β+ Tuổi ε+ Nhưng người ta cũng thường quan sát thấy, thu nhập có xu hướng tăng chậm dần khi người ta càng nhiều tuổi hơn so với thời trẻ. Để thể hiện điều đó, chúng ta mở rộng mô hình như sau: Thu nhập = 21 ββ + Học vấn 3β+ Tuổi 4β+ Tuổ 2i ε+ Và chúng ta sẽ kỳ vọng rằng, 3β mang dấu dương, và 4β mang dấu âm. Như vậy, chúng ta đã rời bỏ thế giới của hồi quy đơn và bước sang hồi quy đa biến. Ví dụ 4.2: Nghiên cứu về nhu cầu đầu tư ở Mỹ trong khoảng thời gian từ năm 1968 – 1982. Ở Mỹ, thời kỳ này mang dấu ấn lịch sử là cuộc chiến tranh Việt Nam kéo dài, dẫn đến bội chi ngân sách và lạm phát. Một năm sau khi chiến tranh kết thúc, lạm phát ở Mỹ đã đạt tới mức kỷ lục là 9.31% vào năm1976. Điều đó dẫn đến việc ngân hàng trung ương phải áp dụng mạnh mẽ chính sách tiền tệ chặt, vốn đã được áp dụng trong vài năm trước, và đưa Trần Thiện Trúc Phượng Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-2 mức lãi suất lên tới mức cao kỷ lục là 7.83%. Khi sự dính líu của Mỹ về quân sự tại Việt Nam đã hoàn toàn chấm dứt, nguồn nhân lực trước đây phục vụ cho chiến tranh nay chuyển ào ạt sang khu vực thương mại. Và điều này lại lại làm dấy lên một đợt lạm phát mới, đạt tới 9.44% vào năm 1981, sau đó được đưa về mức 5.99% vào năm 1982 nhờ vào việc nâng lãi suất lên tới 13.42%. Như vậy, lịch sử kinh tế Mỹ trong thời kỳ này được đặc trưng bởi chính sách tiền tệ chặt, kéo theo xu hướng cắt giảm liên tục về đầu tư qua các năm. Chính vì vậy, các nhà nghiên cứu Mỹ đã đề xuất mô hình nghiên cứu sau về cầu đầu tư vào giai đoạn này: INV = 21 ββ + T 3β+ G 4β+ INT ε+ Trong đó, INV và G lần lượt là cầu về đầu tư và GNP thực tế, đơn vị trillions dollars; INT là lãi suất; và T là biến xu thế, tính theo thời gian đã trôi qua, kể từ năm 1968. Từ lý luận kinh tế vĩ mô, chúng ta kỳ vọng rằng, 3β mang dấu dương, và 4β mang dấu âm. Và vì đây là thời kỳ đầu tư ở Mỹ có xu thế bị co hẹp, chúng ta cũng kỳ vọng rằng 2β mang dấu âm. Sử dụng dữ liệu thống kê vĩ mô của nền kinh tế Mỹ, từ năm 1968 - 1982 [xem bảng dữ liệu 4.2 phía dưới], kết quả ước lượng của mô hình hồi quy này như sau: Bảng Error! No text of specified style in document..1: Bảng kết xuất mô hình hồi qui các yếu tố ảnh hưởng đến cầu về đầu tư của Mỹ trong giai đoạn từ 1968 - 1982 Dependent Variable: INV Method: Least Squares Date: 04/09/07 Time: 16:14 Sample: 1 15 Included observations: 15 Variable Coefficient Std. Error t-Statistic Prob. C -0.509237 0.052526 -9.694973 0.0000 T -0.016583 0.001880 -8.819528 0.0000 G 0.670266 0.052426 12.78506 0.0000 INT -0.002365 0.001034 -2.287282 0.0430 R-squared 0.972420 Mean dependent var 0.203333 Adjusted R-squared 0.964898 S.D. dependent var 0.034177 S.E. of regression 0.006403 Akaike info criterion -7.040816 Sum squared resid 0.000451 Schwarz criterion -6.852003 Log likelihood 56.80612 F-statistic 129.2784 Durbin-Watson stat 1.958353 Prob(F-statistic) 0.000000 Trần Thiện Trúc Phượng Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-3 Dưới dạng báo cáo, kết quả đó có thể được viết tóm tắt như dưới đây: INV = -0.5092 - 0.0165T + 0.67G - 0.0023 INT (0.0525) (0.0018) (0.052) (0.001) , N= 15, ESS = 0.00045 972.02 =R Nếu viết dưới dạng sai phân, ta có: INV = - 0.0165Δ Δ T + 0.67ΔG - 0.0023Δ INT Nói khác đi, nếu các yếu tố khác được giữ không đổi, cứ sau mỗi một năm, kể từ năm 1968 (tức là T = 1), nhu cầu đầu tư sẽ bị giảm là -0.0165 trillions dollars. Cũng như vậy, nếu bỏ qua yếu tố xu thế và lãi suất, tác động riêng phần của việc tăng GNP lên 0.1 trillions dollars ( G = 0.1), sẽ làm cầu về đầu tư tăng lên thêm 0.067 trillions; và nếu đẩy lãi suất lên thêm 1% ( INT = 1), trong khi giữ nguyên các yếu tố còn lại, thì sẽ làm đầu tư giảm đi là -0.0023 trillions dollars. Δ Δ Δ Những tính toán trên đây cho thấy có sự tương đồng rõ rệt về cách diễn giải ý nghĩa của các hệ số ước lượng trong mô hình hồi quy đa biến so với trường hợp đơn biến. Điều đó gợi ý rằng, về mặt bản chất, mô hình hồi quy đa biến sẽ chỉ là sự mở rộng của hồi quy đơn biến. Ta sẽ thấy rõ hơn điều đó ở các phần sau. 4.2 Biểu diễn đại số của mô hình hồi quy đa biến Chúng ta hãy đưa ra bảng so sánh về dạng hàm của mô hình hồi quy đa biến so với trường hợp đơn biến: Hồi quy đơn biến Hồi quy đa biến Ví dụ εββ ++= INCCONS 21 INV = 21 ββ + T 3β+ G 4β+ INT ε+ Dạng mô hình εββ ++= XY 21 Y = 221 Xββ + 33 Xβ+ 4β+ 4X ε+ Với mỗi quan sát nnn xy εββ ++= 21 nnnnn xxxy εββββ ++++= 4433221 Như vậy, hồi quy đa biến là một sự mở rộng tự nhiên của trường hợp đơn biến, khi số biến giải thích lớn hơn 2, kể cả hằng số. Để cho tiện lợi, chúng ta sẽ đưa vào các ký hiệu vector: Trần Thiện Trúc Phượng Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-4 Gọi vector hàng là vector các quan sát thứ ),1( 4,3,2 ' nnnn xxxx = Nn ...,2,1= của các biến giải thích. [Lưu ý, dấu phẩy ở bên phải, phía trên vector là dấu chuyển vị. Như vậy, theo mặc định, mọi vector (mà không có dấu chuyển vị) đều được coi là vector cột]. nx Từng “cặp” quan sát dữ liệu do vậy, sẽ là . Nnnn xy 1 ' },{ = Để minh họa, trong ví dụ 4.2 về cầu về đầu tư ở Mỹ (1968 – 82), những cặp và được tô màu: ),( '55 xy ),( '1111 xy Bảng Error! No text of specified style in document..2: Dữ liệu vĩ mô về đầu tư và các biến giải thích của nền kinh tế Mỹ (1968 – 82). Obs INV C T G INT (n) (Y) (X1) (X2) (X3) (X4) 1 0.161 1 1 1.058 5.16 2 0.172 1 2 1.088 5.87 3 0.158 1 3 1.086 5.95 4 0.173 1 4 1.122 4.88 5 0.195 1 5 Trần Thiện Trúc Phượng 1.186 4.5 6 0.217 1 6 1.254 6.44 7 0.199 1 7 ),( '55 xy 1.246 7.83 8 0.163 1 8 1.232 6.25 9 0.195 1 9 1.298 5.5 10 0.231 1 10 1.37 5.46 11 0.257 1 11 1.439 7.46 12 0.259 1 12 1.479 10.28 13 0.225 1 13 ),( '1111 xy 1.474 11.77 14 0.241 1 14 1.503 13.42 15 0.204 1 15 1.475 11.02 Nguồn: Economic Report of the President. Government, Printing Office, Washington D.C., 1983. Tiếp theo, ta gọi vector cột là vector các tham số tổng thể, cần được ước lượng Lưu ý rằng, tích vô hướng giữa hai vector và ⎟⎟ ⎟⎟ ⎟ ⎠ ⎞ ⎜⎜ ⎜⎜ ⎜ ⎝ ⎛ = 4 3 2 1 β β β β β ' nx β sẽ tạo lại phần xu thế trong vế phải của phương trình hồi quy (4.2): Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-5 = ×× 1441 ' βnx ),1( 4,3,2 nnn xxx ⎟⎟ ⎟⎟ ⎟ ⎠ ⎞ ⎜⎜ ⎜⎜ ⎜ ⎝ ⎛ 4 3 2 1 β β β β = Nnxxx nnn ..2,1,4433221 =+++ ββββ Vì vậy, ứng với từng “cặp” quan sát , ta có thể viết lại phương trình hồi quy đó như sau: N nnn yx 1 ' },{ = (4.3) Nnxy nnn ,..,2,1 ' =+= εβ Như vậy, mọi ký hiệu ta đã sử dụng trong ước lượng mô hình hồi quy đơn biến, nay có thể được sử dụng lại cho mô hình hồi quy đa biến. Cụ thể là: (4.3) Nnxy nn ,..,2,1 ^ ' ^ == β Và sai số ước lượng hay số dư (residual) sẽ có dạng: nnn yye ˆ−= (4.4) Việc tiến hành ước lượng các tham số của mô hình bằng phương pháp bình phương cực tiểu tương đương với việc giải bài toán sau: ^min)()( 2 ^ '2 ^ β ββ →−== ∑ ∑n nnn xyeS (4.5) Tương tự như trong hồi quy đơn, ở đây, ta sử dụng điều kiện cực trị, (first order condition, FOC), để tìm các tham số ước lượng . Nói khác đi, ta đi giải hệ phương trình sau: .4,3,2,1, ^ =kkβ 0) ˆ( 1 ^ =∂ ∂ β βS 0) ˆ( 2 ^ =∂ ∂ β βS 0ˆ )ˆ( 3 =∂ ∂ β βS Trần Thiện Trúc Phượng Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-6 0) ˆ( 4 ^ =∂ ∂ β βS (4.6) Đây là hệ gồm 4 phương trình với 4 ẩn số, mà việc giải nó cho chúng ta tham số ước lượng . Sử dụng phần mềm Eviews, kết quả tính toán các tham số này đã được nêu trong bảng báo cáo 4.1 ở trên. ,( 1 ^^ ββ = '4 ^ 3 ^ 2 ^ ),, βββ Mặc dù dạng biểu diễn giải tích của vector là khá phức tạp. Tuy nhiên, về bản chất chúng vẫn không khác gì trường hợp đơn biến. Cụ thể là, tương tự như , phương trình đầu tiên của hệ (4.6) để ước lượng dẫn đến cái điều là, đường hồi quy đi qua điểm trung bình . Và vì vậy, ta cũng có thể nói đến tiêu chuẩn đo lường độ phù hợp của đường hồi quy ^ 14× β ^α 1 ^β ),( ' −− nn xy 2R . Cụ thể là từ mối quan hệ (4.4): nnn eyy += ^ Hay cũng hệt như thế: nnn eyyyy +−=− −− ^ )( Người ta có thể viết lại nó như sau: nnn exxyy +−=− −− ^ '' )()( β Tức là, sự giao động so với trung bình của biến Y được giải thích một phần bởi mô hình, và phần còn lại là sai số , chưa được giải thích bởi mô hình. Sử dụng các điều kiện tìm cực trị (4.6), ta cũng có thể viết lại quan hệ đó như sau: ne 22 ^ 2 )()( nnnnnn eyyyy ∑∑∑ +−=− −− Hay cũng vậy, TSS = RSS + ESS Vì thế, chúng ta có thể đưa ra định nghĩa: Trần Thiện Trúc Phượng Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-7 TSS ESSR −= 12 ( ). 10 2 ≤≤ R và sử dụng nó làm thước đo mức độ phù hợp của đường hồi quy với dữ liệu quan sát. Phần tiếp sau sẽ đề cập đến bản chất thống kê của mô hình hồi quy đa biến. 4.3 Bản chất thống kê của mô hình hồi quy đa biến Từ bây giờ, chúng ta sẽ sử dụng dạng tổng quát của mô hình hồi quy đa biến: Y = 221 Xββ + kk Xβ++ .... Kβ++ ... KX ε+ (4.7) nnKKnkknn xxxy εββββ ++++++= .........221 Nnx nn ,...,3,2,1, ' =+= εβ Việc hồi quy mô hình (4.7) sẽ cho ta biểu diễn sau: Y = 22 ^ 1 ^ Xββ + kk X ^ .... β++ K ^ ... β++ KX e+ (4.8) nnKKnkknn exxxy ++++++= ^^ 22 ^ 1 ^ ......... ββββ Nnex nn ,...,3,2,1, ^ ' =+= β Trong đó, N là số quan sát, và K là số biến giải thích. Ta phát biểu định lý sau1: Định lý 4.1: Phương pháp bình phương cực tiểu, áp dụng cho mô hình hồi quy đa biến, sẽ cho ta các tham số ước lượng dưới dạng sau: (4.9) ∑ =+= n nknkk Kkc ,..,2,1,ˆ εββ Cũng như trường hợp đơn biến, phương trình (4.9) chỉ ra rằng: bị tác động bởi các yếu tố ngẫu nhiên kβˆ nε , làm giá trị của nó không trùng khít với kβ tổng thể. Và vì bị tác động Trần Thiện Trúc Phượng 1 Xem chứng minh chi tiết ở chương 8, phần Maximum likelihood. Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-8 bởi các yếu tố ngẫu nhiên, cũng là một biến ngẫu nhiên. Do đó, độ tốt của ước lượng sẽ phụ thuộc trực tiếp vào bản chất của các quá trình ngẫu nhiên . kβˆ N nn 1}{ =ε Điều này dẫn đến việc cần phải khắc họa bản chất thống kê của mô hình hồi quy, như chúng ta đã làm cho trường hợp đơn biến. Ta sẽ tiếp tục sử dụng các giả thuyết đã đưa ra về nε . Cụ thể là: A1 ,0=nEε với mọi n. A2 với mọi n. ,2σε =nVar A3 , với mọi n. Và: ),0(~ 2σε Niidn A4 , với mọi n. β'' )|( nnn xxyE = Đối với trường hợp đa biến, chúng ta đưa thêm đòi hỏi sau: Gọi là ma trận tạo bởi các vector cột của K biến giải thích [xem lại ví dụ minh họa về ma trận X ở bảng 4.2 về dữ liệu của mô hình đầu tư]. Khi đó, ta đòi hỏi rằng: ],...,,...,,[ 21 KkKN XXXXX =× A5 Các cột là độc lập tuyến tính. Hay cũng vậy, rank X = K. },...,,...,,{ 21 Kk XXXX Về mặt hình học, giả thuyết này có ý nghĩa như sau. Hãy xét trường hợp K = 2, phương pháp bình phương cực tiểu có thể được biểu diễn bởi lược đồ dưới đây: Trần Thiện Trúc Phượng Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-9 Đồ thị Error! No text of specified style in document..1: Biểu diễn hình học của hồi quy 1X 2X ^ y 2 ^β 1 ^β Việc ước lượng tham số cũng giống như là tìm các hệ số sao cho . Để làm được điều đó, điều kiện cần là các vector không được trùng khít với nhau. Hay cũng vậy, phải độc lập tuyến tính. Đây được gọi là điều kiện xác định (identification condition). Trong trường hợp tổng quát, khi ^β '^2 ^ 1 ^ 12 ),( βββ = × 2 ^ 21 ^ 1 ^ XXy ββ += 21 , XX 21 , XX 2≥K , điều kiện đó được phát biểu dưới dạng giả thuyết A5. Chúng ta sẽ sử dụng giả thuyết này khi bàn tới vấn đề đa cộng tuyến (multicollinearity) trong chương 7. 4.4 Kiểm định các giả thuyết thống kê Bây giờ hãy chỉ chú ý đến giả thuyết đầu tiên A1 – A3, và sử dụng chúng để đánh giá tính tốt của ước lượng theo các tiêu chuẩn thống kê. Từ phương trình (4.9), ta đã có: ∑+= nknkk c εββˆ . Bây giờ, hãy áp dụng toán tử kỳ vọng vào hai vế của (4.9). Ta có: ∑+= )(ˆ nknkk cEE εββ ∑+= nknk Ec εβ kβ= (4.10) Trần Thiện Trúc Phượng Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-10 [ở đây, ta sử dụng giả thiết A1: 0=nEε ]. Do vậy, là ước lượng không chệch của kβˆ kβ . Tiếp theo, sử dụng lại công thức: )()( ExxVarxVar −= [xem chương 1, phần ôn tập], và lưu ý (4.9), (4.10), ta có: )(ˆ ^ kkk VarVar βββ −= = ∑ )( nkncVar ε Sử dụng giả thiết A3 về tính độc lập của các yếu tố ngẫu nhiên, cuối cùng ta nhận được: ∑= nknk VarcVar εβ 2ˆ = , hay ∑ 22 kncσ Kk S Var kk k ,..,2,1, 2^ == σβ (4.11) (ở đây, mặc dù ta không đưa ra được tính toán trực tiếp; nhưng về cơ bản cũng là phương sai mẫu của biến , tương tự như trong trường hợp đơn biến). kkS kX XXS Định Lý 4.2 [Gauss – Markov]: Phương pháp bình phương cực tiểu có sai số ước lượng, đo lường bởi , là nhỏ nhất trong lớp tất cả các ước lượng tuyến tính và không chệch. KkVar k ,..,2,1, ^ =β Ta cũng nên nhấn mạnh lại rằng, chúng ta có được những tính chất rất tốt: không chệch và hiệu quả của ước lượng bình phương cực tiểu, mà chỉ đòi hỏi có trung bình bằng zero, tính độc lập, và phương sai giống nhau của các yếu tố ngẫu nhiên – tức là giả thiết A3. Sử dụng (4.9) – (4.11), chúng ta đi đến kết luận rằng: ),(~ˆ 2 kk kk S N σββ . Điều đó có nghĩa là, sau khi chuẩn hóa, )1,0(~ ˆ 2 N S Z kk kk k σ ββ −= . Thay thế bởi ước lượng không chệch 2σ Trần Thiện Trúc Phượng Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-11 của nó là ∑−= n neKNs 22 1 , ta có thống kê )(~)( ˆ ^ ^ 2 KNt seSs t k kk kk kk k −−=−= β ββββ . Chúng ta bây giờ có thể xây dựng khoảng tin cậy cho Kkk ,..,2,1, =β , và tiến hành kiểm định các giả thuyết thống kê về tham số của tổng thể. Chẳng hạn như nếu muốn kiểm định tính có ý nghĩa của biến giải thích , chúng ta lập giả thuyết sau: KkX k ,..,2,1, = 0:0 =kH β .vs. 0:1 ≠kH β Việc kiểm định bao gồm các bước sau: Bước 1: Xác định thống kê – t (t-stat): )(~ )( ^ ^ KNt se t k k k −= β β Bước 2: Tra bảng thống kê t-student với (N-K) bậc tự do )( KNt − để tìm giá trị t-tra bảng (t-critical) )( 2 KNt −λ , ứng với mỗi mức ý nghĩa (significance) λ [Chẳng hạn, 0.05 (5%); hay 0.1 (10%)] Bước 3: Bác bỏ giả thuyết (viết tắt là ), nếu 0H 0RH )( )( || 2^ ^ KNt se t k k k −≥= λβ β , và không bác bỏ giả thuyết đó ( ), nếu 0DNRH )( )( 2 ^ ^ KNt se k k −≤ λβ β . Cũng như trong trường hợp đơn biến, người ta thường hay sử dụng p-value, hơn là phải tính toán và tra bảng qua các bước 1 đến 3 như trên. Cụ thể, ứng với từng biến giải thích KkX k ,..,2,1, = , ta cũng đặt: |)||)({|Pr ktKNtobvaluep ≥−=− Trần Thiện Trúc Phượng Khoa Kinh tế Kinh tế lượng ©2007 ĐHQG TP.HCM Lê Hồng Nhật 4-12 Cũng hệt như ở đồ thị 3-9, chúng ta sẽ bác bỏ giả thuyết 0:0 =kH β , nếu: λ≤− valuep , [trong trường hợp đó, ta nói là có ý nghĩa ở mức kX %λ ]. Và chúng ta sẽ không bác bỏ giả thuyết đó, nếu λ≥− valuep . Trong ví dụ 4.1 về đầu tư ở Mỹ (1968-82), p-value của cả 3 biến giải thích: T, G, INT , đều nhỏ hơn 5%. Vì vậy, ta nói rằng tất cả các biến này là có ý nghĩa ở mức . %5=λ Trần Thiện Trúc Phượng

File đính kèm:

  • pdfBai Giang Kinh Te Luong. Phan 4.pdf
Giáo án liên quan