BÀI 3 (tiếp theo)
MÔ HÌNH NHIỀU PHƯƠNG TRÌNH
Sau đây ta sẽ nghiên cứu một số phương pháp theo cách tiếp cận thứ nhất.
1. MÔ HÌNH ĐỆ QUY VÀ PHƯƠNG PHÁP OLS.
Mô hình đệ quy là mô hình trong đó ta có thể sắp xếp các phương trình theo
một trình tự đặc biệt. Biến nội sinh trong phương trình thứ nhất chỉ xác định theo
các biến ngoại sinh. Biến nội sinh trong phương trình thứ hai thì xác định theo các
biến nội sinh của phương trình thứ nhất và các biến ngoại sinh - nhưng không xác
định theo biến nội sinh nào khác. Biến nội sinh trong phương trình thứ ba được
xác định theo biến nội sinh của phương trình thứ nhất và phương trình thứ hai
và các biến ngoại sinh, chứ không theo bất cứ một biến nội sinh nào khác đưa
vào phương trình như một biến hồi quy. Hệ ba phương trình sau đây là hệ đệ
quy:
Y
1t
= α
1
+ α
2
X
1t
+ α
3
X
2t
+ u
1t
Y
2t
= β
1
+ β
2
Y
1t
+ β
3
X
2t
+ u
2t
Y
3t
= γ
1
+ γ
2
Y
1t
+ γ
3
Y
2t
+ γ
4
X
3t
+ u
3t
Mô hình đệ quy như trên có hai đặc trưng quan trọng là:
Thứ nhất, không có tác động ngược trở lại của biến nội sinh ở phương trình sau
đối với biến nội sinh ở phương trình trước theo quan hệ dây chuyền.
Thứ hai, sai số ngẫu nhiên (nhiễu) được giả thiết là độc lập, theo nghĩa sau đây:
cov (u
1t
, u
2t
) = cov (u
2t
, u
3t
) = cov (u
1t
, u
3t
) = 0
Đối với hệ phương trình đệ quy có mối quan hệ nhân quả đa chiều, phương
pháp OLS áp dụng đối với phương trình cấu trúc sẽ cho các kết quả ước lượng
vững, tiệm cận hiệu quả.
Tuy nhiên, trong một vài trường hợp tác động qua lại của các biến nội sinh có
tính trễ và chính điều đó cho phép áp dụng phương pháp OLS ước lượng một hệ
đệ quy hoàn toàn có tính kỹ thuật thuần tuý. Hãy xét mô hình tiền công sau đây (
thuộc nhóm mô hình Phillip):
Phương trình giá: p
t
= b
01
+ b
11
w
t-1
+ b
12
r
t
+ b
13
m
t
+ b
14
l
t
+ u
1t
Phương trình tiền công: w
t
= b
20
+ b
21
un
t
+ b
22
p
t
+ u
2t
Trong đó:
p : tỷ lệ thay đổi giá sản phẩm,
w : tỷ lệ thay đổi tiền công,
r : tỷ lệ thay đổi giá vốn,
m : tỷ lệ thay đổi giá nhập khẩu,
l : tỷ lệ thay đổi của năng suất lao động
un : tỷ lệ thất nghiệp.
Mô hình trên về mặt hình thức không phải là một mô hình đệ quy, tuy vậy có
thể sử dụng OLS để ước lượng từng phương trình với một vài giả thiết khá rộng
rãi đối với biến tiền công (w). Quan hệ tác động qua lại của hai biến nội sinh p và
w hoàn toàn có thể biến đổi một cách hình thức bằng cách coi w
t-1
là một biến
khác w
t
, cách quan niệm như vậy không làm mất các tính chất cần thiết của các
ước lượng khi bản thân quá trình w nội sinh theo thời gian - có thể thấy rằng, giả
thiết này hầu như đạt được với biến w là tỷ lệ thay đổi tiền công theo thời gian.
2. PHƯƠNG TRÌNH ĐỊNH DẠNG ĐÚNG VÀ PHƯƠNG PHÁP BÌNH
PHƯƠNG NHỎ NHẤT GIÁN TIẾP (ILS)
Xét mô hình hai phương trình biểu diễn tiêu dùng trong một nền kinhtế đóng:
C
t
= β
1
+ β
2
Y
t
+ u
t
0 < β
2
< 1
Y
t
= C
t
+ I
t
Thay Yt hoặc Ct có phương trình rút gọn dạng
Yt = BETA1/(1-BETA2) +1/(1-BETA2)*It +Ut
Pi1 Pi2
Trong mô hình này hàm tiêu dùng là hàm định dạng đúng . Điều đó có nghĩa là
với giá trị các tham số ước lượng trong hàm rút gọn chỉ có một cách duy nhất để
xác định các tham số cấu trúc ban đầu.
Phương pháp ILS gồm có ba bước :
1. Biến đổi thành những phương trình rút gọn;
2. Ước lượng từng phương trình rút gọn bằng OLS;
3. Xác định các ước lượng vững cho các tham số cấu trúc từ các tham số của các
phương trình rút gọn.
Trên thực tế, phương pháp ước lượng ILS có hai hạn chế cơ bản.
* Nhiều phương trình thay vì định dạng đúng thì
l lại vô định .
• Ước lượng sai số chuẩn của các hệ số của phương trình cấu trúc không
hoàn toàn đơn giản phản ánh sai lệch của thống kê ước lượng các hệ số đó.
Do đó, ILS không hay được sử dụng.
4. PHƯƠNG TRÌNH VÔ ĐỊNH VÀ PHƯƠNG PHÁP 2SLS
Phương pháp 2SLS có thể áp dụng để thu được các ước lượng vững và
tiệm cận hiệu quả. Đối với các phương trình định dạng đúng, phương pháp
cho kết quả giống như ILS. Tuy nhiên nó còn có thể áp dụng được với cả các
mô hình định dạng đúng và vô định.
Với các mô hình nhiều phương trình có các phương trình vô định, phương
pháp 2SLS cho phép ước lượng các tham số hiệu quả hơn vì chính các bước
ước lượng của phương pháp này sẽ làm cho các phương trình chỉ còn là các
phương trình địng dạng đúng.
Để làm rõ nội dung của phương pháp, ta xét mô hình sau:
Hàm thu nhập: Y
1t
= β
10
+ β
11
Y
2t
+ α
11
X
1t
+ α
12
X
2t
+ u
1t
Hàm cung tiền: Y
2t
= β
20
+ β
21
Y
1t
+ u
2t
Trong đó: Y
t
- thu nhập; Y
2t
- dự trữ tiền; X
1
- đầu tư; X
2
- chi tiêu của chính
phủ về hàng hóa và dịch vụ. X
1
và X
2
là các biến ngoại sinh.
Theo điều kiện cần và điều kiện đủ của định dạng, phương trình thu nhập
không định dạng được, phương trình cung là vô định . Do đó, không thể áp dụng
ILS. Nếu như áp dụng OLS cho hàm cung tiền thì các ước lượng sẽ không có tính
vững. Để giải quyết vấn đề này, ngưới ta dùng phương pháp biến công cụ.
Giả thiết rằng ta tìm được một biến xấp xỉ biến giải thích Y
1
ở hàm cung tiền
song lại không tương quan với u
2
. Một biến như vậy được gọi là biến công cụ.
Dùng biến này để ước lượng hàm cung tiền. Các ước lượng nhận được có tính chất
vững. Để tìm biến công cụ gười ta dùng phương pháp bình phương nhỏ nhất hai
giai đoạn (2SLS).
Giai đoạn 1: Để loại bỏ tính tự tương quan giữa Y
1
và u
2
, trước hết hồi quy Y
1
phụ thuộc vào tất cả các biến ngoại sinh trong hệ phương trình. Trong ví dụ này:
121 32
ˆˆ ˆ
tt
X
X
π
ππ
+
+
Từ đó, ta có:
Y
1t
= + e
t
Y
1
ˆ
t
Giai đoạn 2: Phương trình cung tiền bây giờ có thể viết lại như sau:
Y
2t
= β
20
+ β
21
Y
1t
+ u
2t
= β
20
+ β
21
( + e
t
Y
1
ˆ
t
) + u
2t
= β
20
+ β
21
t
Y
1
ˆ
+ β
21
e
t
+ u
2t
= β
20
+ β
21
t
Y
1
ˆ
+ u
*
t
Phương trình trên giống hàm cung tiền xuất phát trong đó Y
1
đã được thay
bằng . Như vậy, dù Y
1
ˆ
Y
1
tương quan với u
2
nhưng không tương quan với u
1
ˆ
Y
2
nên ước lượng của β
20
và β
21
là các ước lượng vững.
Thực chất của 2SLS là lọc khỏi Y
1
ảnh hưởng của yếu tố ngẫu nhiên u
2
.
Việc này được thực hiện khi hồi quy phương trình rút gọn. Các ước lượng tìm
được là các ước lượng vững, chúng hội tụ về giá trị thực khi kích thước của mẫu
tăng vô hạn.
Bây giờ ta sửa đổi hàm thu nhập và hàm dự trữ tiền như sau:
Y
1t
= β
10
+ β
12
Y
2t
+ α
11
X
1t
+ α
12
X
2t
+ u
1t
Y
2t
= β
20
+ β
21
Y
1t
+ α
23
X
3t
+ α
24
X
24
+ u
2t
Trong đó X
3
- thu nhập, X
4
- mức cung tiền. Hai biến này là các biến ngoại
sinh.
Ta có thể thấy rằng, cả hai phương trình đều là vô định.
Thủ tục 2SLS như sau:
Giai đoạn 1: Ước lượng các hàm rút gọn:
Y
1t
= π
10
+ π
11
X
1t
+ π
12
X
2t
+ π
13
X
3t
+ π
14
X
4t
+ v
1t
Y
2t
= π
20
+ π
21
X
1t
+ π
22
X
2t
+ π
23
X
3t
+ π
24
X
4t
+ v
2t
Giai đoạn 2: Ước lượng mô hình xuất phát bằng cách thay Y
1
và Y
2
ở vế phải
của các phương trình bằng và nhận được ở giai đoạn 1.
1
Y
ˆ
2
Y
ˆ
Các ưu điểm của 2SLS:
• Có thể áp dụng cho từng phương trình riêng rẽ, không cần chú ý đến các
phương trình khác. Điều này là thuận lợi khi ước lượng một hệ gồm nhiều
phương trình.
• ILS đưa ra các ước lượng của các hệ số của phương trình thu gọn, để tìm được
ước lượng của các hệ số ban đầu, ta phải thực hiện một vài tính toán, 2SLS cho
ngay ước lượng của từng hệ số.
• Dễ áp dụng vì chỉ cần biết tổng số các biến ngoại sinh.
• Khi áp dụng cho các phương trình định dạng đúng thì kết quả không khác so
với kết quả của ILS.
• Nếu R
2
ở giai đoạn 1 khá cao thì có nghĩa là chúng ta đã tìm được một xấp xỉ
khá tốt của Y
1
và Y
2
.
• 2SLS cho Se của các ước lượng, trong khi đó ILS không cho.
• Hạn chế : 2SLS chỉ dùng trong trường hợp mẫu lớn.
5. PHƯƠNG PHÁP BÌNH PHƯƠNG NHỎ NHẤT BA GIAI ĐOẠN- 3SLS.
Là sự mở rộng của 2SLS và dùng để ứơc lượng các phương trình của một hệ
phương trình. Phương pháp này gồm ba giai đoạn sau:
Giai đoạn 1. Ứơc lượng các phương trình rút gọn.
Giai đoạn 2. Ước lượng các phương trình cấu trúc bằng cách thay các biến nội sinh ở
vế phải của các phương trình bằng các ước lượng thu được ở giai đoạn 1.
Giai đoạn 3. Gồm các công việc sau:
+ Tính các phần dư e
1t
, e
2t
, . . . e
mt
thu được ở giai đoạn 2 tương ứng với các phương
trình cấu trúc.
+ Tính ma trận hiệp phương sai của các sai số ngẫu nhiên đối với từng phương trình
của hệ.
+ Biến đổi các biến số theo phương pháp bình phương nhỏ nhất tổng quát.
Các ước lượng nhận được theo phương pháp 3SLS là các ước lượng chệch nhưng
vững. Phương pháp 3SLS hiệu quả hơn so với 2SLS.
Ví dụ 1: Sử dụng tệp sốliệu ch10bt14 hãy ước lượng mô hình sau:
LM: R
t
= α
1
+ α
2
M
t
+ α
3
Y
t
+ α
4
M
t-1
+ u
1t
IS: Y
t
= β
1
+ β
2
R
t
+ β
3
I
t
+ u
2t
Kết quả hồi quy bằng 2SLS như sau:
Đặt trên cửa sổ lệnh SYStem Quy
inst c i m m(-1)
r=c(1)+c(2)*m+c(3)*y+c(4)*m(-1)
y=c(5)+c(6)*r+c(7)*i
Estimate Two least square
System: MYSYS
Estimation Method: Two-Stage Least Squares
Date: 11/23/08 Time: 21:36
Sample: 1960 1990
Included observations: 31
Total system (balanced) observations 62
Instruments: C M M(-1) I
Coefficien
t
Std. Error t-Statistic Prob.
C(1) 10.11926 0.970995 10.42154 0.0000
C(2) -0.045620 0.018716 -2.437471 0.0181
C(3) 0.012456 0.001912 6.513333 0.0000
C(4) -0.037661 0.016342 -2.304583 0.0250
C(5) 671.9805 212.1567 3.167379 0.0025
C(6) -143.2998 40.59820 -3.529709 0.0008
C(7) 7.039408 0.366547 19.20465 0.0000
Determinant residual
covariance
122035.4
Equation: R=C(1)+C(2)*M+C(3)*Y+C(4)*M(-1)
Observations: 31
R-squared 0.771641 Mean dependent
var
7.38548
4
Adjusted R-
squared
0.746268 S.D. dependent var 2.82187
2
S.E. of regression 1.421429 Sum squared resid 54.5524
0
Durbin-Watson
stat
0.730862
Equation: Y=C(5)+C(6)*R+C(7)*I
Observations: 31
R-squared 0.966006 Mean dependent
var
2185.54
6
Adjusted R-
squared
0.963578 S.D. dependent var 1578.48
8
S.E. of regression 301.2472 Sum squared resid 2540996
.
Durbin-Watson
stat
0.706409
Kết quả hồi quy bằng 3SLS như sau:
System: MYSYS
Estimation Method: Iterative Three-Stage Least Squares
Date: 11/23/08 Time: 21:40
Sample: 1960 1990
Included observations: 31
Total system (balanced) observations 62
Instruments: C M M(-1) I
Simultaneous weighting matrix & coefficient iteration
Convergence achieved after: 3 weight matrices, 4 total coef
iterations
Coefficien
t
Std. Error t-Statistic Prob.
C(1) 10.07889 0.903808 11.15158 0.0000
C(2) -0.041006 0.016238 -2.525321 0.0145
C(3) 0.012356 0.001778 6.949073 0.0000
C(4) -0.041822 0.014104 -2.965193 0.0045
C(5) 671.9805 201.6299 3.332743 0.0015
C(6) -143.2998 38.58379 -3.713990 0.0005
C(7) 7.039408 0.348360 20.20730 0.0000
Determinant residual
covariance
120405.1
Equation: R=C(1)+C(2)*M+C(3)*Y+C(4)*M(-1)
Observations: 31
R-squared 0.771966 Mean dependent
var
7.38548
4
Adjusted R-
squared
0.746628 S.D. dependent var 2.82187
2
S.E. of regression 1.420418 Sum squared resid 54.4748
8
Durbin-Watson
stat
0.764417
Equation: Y=C(5)+C(6)*R+C(7)*I
Observations: 31
R-squared 0.966006 Mean dependent
var
2185.54
6
Adjusted R-
squared
0.963578 S.D. dependent var 1578.48
8
S.E. of regression 301.2472 Sum squared resid 2540996
.
Durbin-Watson
stat
0.706409
Ví dụ 2: Xét mô hình kinhtế vĩ mô gồm hai phương trình:
Y
1t
= α
1
+ α
2
Y
2t
+ α
3
X
1t
+ α
4
X
2t
+ u
1t
và Y
2t
= β
1
+ β
2
Y
1t
+ u
2t
Trong mô hình này, Y
1t
là tổng thu nhập thực của nền kinh tế, Y
2t
là tổng khối lượng
tiền mặt và cả hai đều là biến nội sinh. Các biến ngoại sinh là đầu tư X
1t
và chi tiêu của
chính phủ X
2t
.
Với các sốliệu của Mỹ( tr.691 Gujarati) hãy ứơc lượng mô hình trên.
. t-Statistic Prob.
C(1) 10.0 788 9 0.90 380 8 11.151 58 0.0000
C(2) -0 .041006 0.0162 38 -2 .525321 0.0145
C(3) 0.012356 0.0017 78 6.949073 0.0000
C(4) -0 .04 182 2. 0.0000
C(4) -0 .037661 0.016342 -2 .304 583 0.0250
C(5) 671. 980 5 212.1567 3.167379 0.0025
C(6) -1 43.29 98 40.5 982 0 -3 .529709 0.00 08
C(7) 7.0394 08 0.366547 19.20465