2.2.1.1. Xác định độ trễ tối ƣu từ mơ hình VAR
Kết quả của việc xác định độ trễ tối ưu theo các tiêu chuẩn LR, FPE, AIC và HQ đều là 2, riêng tiêu chuẩn SC là 1, cịn tiêu chuẩn LogL khơng xác định. Do đó, độ trễ tối ưu chọn theo các tiêu chuẩn là p =2.
Bên cạnh đó, tác giả cũng tiến hành kiểm tra hiện tượng tương quan phần dư theo kiểm tra LM và kiểm tra hiện tượng tự tương quan của phần dư theo kiểm tra Portmanteau đều cho thấy độ trễ 2 là phù hợp vì các phần dư khơng tương quan với nhau và không tự tương quan ở mức ý nghĩa 1%. Mơ hình VAR ổn định với độ trễ p =2 khi tất cả nghiệm đơn vị đều nằm trong vịng trịn đơn vị. Ngồi ra, tác giả cũng kiểm tra đồng liên kết - Cointegration cho thấy có 4 mối quan hệ đồng liên kết, vì vậy các biến trong mơ hình đều là chuỗi dừng (hỗ trợ cho kết luận kiểm tra nghiệm đơn vị trên).
Các bước xác định được độ trễ tối ưu p =2 của các biến đều phù hợp với nghiên cứu của Abdul Aleem, Amine Lahiani (2014).
2.2.1.2. Xác định số lƣợng và giá trị ngƣỡng, tham số độ trễ của biến ngƣỡng tỷ lệ lạm phát từ các trƣờng hợp mơ hình TVAR
- Mơ hình TVAR có 1 giá trị ngưỡng r (nthreshold =1) và tham số độ trễ (thdelay) d của biến ngưỡng tỷ lệ lạm phát, độ trễ p=2:
yt = (c1 + A1,1yt-1+ A1,2yt-2) (I(inft-d ≤ r)) + (c2 + A2,1yt-1+ A2,2yt-2) (1-I(inft-d ≤ r)) +t
Với yt là các biến nội sinh dt_neer_sa, Output gap, inf, dt_lir. I()là một hàm
số chỉ dẫn bằng 1 khi inft-d ≤ r và bằng 0 khi inft-d > r.
- Mơ hình TVAR có 2 giá trị ngưỡng r1, r2 (nthreshold =2) và tham số độ trễ (thdelay) d của biến ngưỡng tỷ lệ lạm phát, độ trễ p=2:
yt = (c1 + A1,1yt-1+ A1,2yt-2) (I1(inft-d ≤ r1)) +
(c2 + A2,1yt-1+ A2,2yt-2) (1-I1(inft-d ≤ r1)) I2(inft-d ≤ r2) + (c3 + A3,1yt-1+ A3,2yt-2) (1-I2(inft-d ≤ r2)) +t
Với: yt là các biến nội sinh dt_neer_sa,Output gap, inf, dt_lir. I1(),I2()lần lượt là các hàm số chỉ dẫn bằng 1 khi inft-d ≤ r1, inft-d ≤ r2 và bằng 0 khi inft-d > r1, inft-d > r2.
Như đã trình bày ở phần mở đầu của mục 2.2, tác giả trình bày các chỉ tiêu SSR, AIC và BIC lần lượt của các kết quả hồi quy như sau:
Bảng 2.2: Chỉ tiêu AIC, BIC và SSR từ mơ hình TVAR.
Trường hợp AIC BIC SSR
TVAR(nthreshold=1,thdelay = d =1) -2764.307 -2552.222 26.26031 TVAR(nthreshold=1,thdelay = d =2) -2778.912 -2566.827* 24.01217 TVAR(nthreshold=2,thdelay = d =1) -2790.528* -2470.948 23.69683 TVAR(nthreshold=2,thdelay = d =2) -2778.497 -2458.917 22.54563*
Nguồn: Kết quả tính tốn từ phần mềm R-packages tác giả thực hiện.
Theo Josip và Petra (2009), chọn lựa mơ hình tốt nhất căn cứ vào giá trị nhỏ nhất của các chỉ tiêu AIC, BIC, SSR. Vì vậy, bảng kết quả trên có giá trị gắn dấu” * “ thể hiện mơ hình được lựa chọn. Trong bài nghiên cứu này, tác giả chọn mơ hình TVAR theo chỉ tiêu AIC, giống như sự lựa chọn mơ hình TVAR của Abdul Aleem, Amine Lahiani (2014), do đó kết quả chọn mơ hình TVAR (thdelay =d =1, nthreshold = 2).
2.2.1.3. Kiểm tra tính tuyến tính
Kiểm tra tính tuyến tính của mơ hình TVAR (thdelay =d =1) theo phương pháp mở rộng đa biến đề nghị bởi Lo và Zivot (2001) trên cơ sở kiểm tra ngưỡng của Hansen (1999) với phân phối tự khởi động (bootstrap distribution) với nboot =1000.
Bảng kết quả nhƣ sau:
Bảng A: Giả thuyết H0 : Mơ hình VAR tuyến tính.
H1 : Mơ hình TVAR có 2 chế độ hay có 1 giá trị ngưỡng. LR statistic = 65.13662
p_value = 0.07600 < 10%
Giá trị ngưỡng được ước lượng : r = 0.849
Bác bỏ giả thuyết H0: Mơ hình VAR tuyến tính hay chọn mơ hình TVAR có 2 chế
độ với mức ý nghĩa 10% .
Bảng B: Giả thuyết H0 : Mơ hình VAR tuyến tính.
H1 : Mơ hình TVAR có 3 chế độ hay có 2 giá trị ngưỡng. LR statistic = 165.3577
p_value = 0.0000 < 1%
Giá trị ngưỡng được ước lượng: r = 0.226 và r = 0.849
Bác bỏ giả thuyết H0: Mơ hình VAR tuyến tính hay chọn mơ hình TVAR có 3 chế
độ với mức ý nghĩa 1%.
Bảng C: Giả thuyết H0 : Mơ hình TVAR có 2 chế độ hay có 1 giá trị ngưỡng. H1 : Mơ hình TVAR có 3 chế độ hay có 2 giá trị ngưỡng. LR statistic = 100.2211
p_value = 0.0030 < 5%
Giá trị ngưỡng được ước lượng: 0.226 và 0.849
Bác bỏ giả thuyết H0: Mơ hình TVAR có 2 chế độ hay chọn mơ hình TVAR có 3
chế độ với mức ý nghĩa 5%.
2.2.1.4. Kết quả hồi quy từ mơ hình TVAR đã chọn
Mơ hình TVAR có 3 chế độ với biến ngưỡng là tỷ lệ lạm phát có độ trễ 1 được chọn từ mục 2.2.1.2 và 2.2.1.3. Kết quả được trình bày cho biểu thức của biến phụ
thuộc là tỷ lệ lạm phát (inf) thể hiện trong bảng 2.3, và cho thấy rằng: ở chế độ
1(Lower) inft-1 ≤ 0.226 thì các hệ số của dt_lneer_sa lần lượt tại độ trễ bằng 1 và 2 có ý nghĩa thống kê lần lượt với mức ý nghĩa 10% và 5%, còn các hệ số khác tại các độ trễ của dt_lneer_sa khơng có ý nghĩa thống kê ở các chế độ 2 và 3. Các kết quả này khác với kết quả của Abdul Aleem, Amine Lahiani (2014) ở hệ số của NEER chỉ có ý nghĩa thống kê tại độ trễ 3 ở chế độ 3. Mặc khác, bài nghiên cứu cũng hướng vào mục tiêu phân tích mức độ truyền dẫn của tỷ giá hối đoái theo các chế độ của tỷ lệ lạm phát nên chú trọng vào phân tích phản ứng xung.
Bảng 2.3: Kết quả hồi quy của biến phụ thuộc tỷ lệ lạm phát từ mơ hình
TVAR có 3 chế độ. Kết quả hệ số hồi quy Chế độ 1(Lower): inft-1 ≤ 0.226 Chế độ 2(Middle): 0.226 < inft-1 ≤ 0.849 Chế độ 3(Upper): 0.849 < inft-1 Phần trăm quan sát: 23.7% Phần trăm quan sát: 50.4% Phần trăm quan sát: 25.9% Intercept 0.0023(0.0917) 0.4257(0.1768)* 0.3077(0.2321) dt_lneer_sa(-1) -11.0086(5.8483). -3.6542(5.6367) -5.6794(4.8255) Output gap(-1) 0.3409(1.2978) -0.0208(0.5759) 5.2094(1.7151)** inf(-1) 0.1096(0.2524) 0.0311(0.3596) 0.3407(0.1490)* dt_lir(-1) -0.8082(1.0521) -3.3828(2.7629) -1.3109(1.3178) dt_lneer_sa(-2) 11.7316(5.8678)* 0.7036(5.6300) -5.1977(5.1609) Output gap(-2) -1.2600(1.0740) 0.3132(0.6899) 2.0076(1.0946). inf(-2) -0.1635(0.2045) 0.2383(0.2038) 0.2638(0.1315)* dt_lir(-2) -0.0236(0.8250) 2.0093(2.5517) 0.2816(1.3795) AIC -2790.528 BIC -2470.948 SSR 23.69683
Nguồn: Kết quả tính tốn từ phần mềm R-packages tác giả thực hiện.
Ghi chú: - Mơ hình TVAR có 3 chế độ:
yt = (c1 + A1,1yt-1+ A1,2yt-2) (I1(inft-1 ≤ r1)) +
(c2 + A2,1yt-1+A2,2yt-2) (1-I1(inft-1 ≤ r1)) I2(inft-1 ≤ r2) + (c3 + A3,1yt-1+ A3,2yt-2) (1-I2(inft-1 ≤ r2)) +t
Với: yt là các biến nội sinh dt_neer_sa,Output gap, inf, dt_lir. I1(),I2() lần lượt là các hàm số chỉ dẫn bằng 1 khi inft-1 ≤ r1, inft-1 ≤ r2 và bằng 0 khi inft-1 >r1, inft-1 > r2 . r1 và r2 là các giá trị ngưỡng của biến inft-1.
- *** thể hiện mức ý nghĩa 0% - ** thể hiện mức ý nghĩa 1% - * thể hiện mức ý nghĩa 5% - . thể hiện mức ý nghĩa 10%
Hình 2.1: Giá trị ngưỡng của biến ngưỡng tỷ lệ lạm phát độ trễ bậc 1 với mơ
hình TVAR 2 giá trị ngưỡng.
Nguồn: Kết quả tính tốn từ phần mềm R-packages tác giả thực hiện.
2.2.2. Phân tích phản ứng xung truyền thống (IRFs) cho mơ hình TVAR có 2 giá trị ngƣỡng. giá trị ngƣỡng.
Trong phần này, chúng ta kì vọng vào mức độ truyền dẫn của tỷ giá hối đoái danh nghĩa hiệu lực sẽ cao hơn ở chế độ lạm phát cao so với ở chế độ lạm phát thấp . Đây cũng là quan điểm của Taylor (2000) và cũng là kết quả của các bài nghiên cứu được review ở trên.
Xây dựng hàm phản ứng xung tổng thể phi tuyến cho mơ hình TVAR có 2 giá trị ngưỡng là rất phức tạp và tốn nhiều thời gian, vì thế tác giả xem xét cú sốc tích cực có độ lớn bằng 1 đơn vị độ lệch chuẩn của tỷ giá hối đoái danh nghĩa hiệu lực xảy ra lần lượt ở chế độ 1 và chế độ 2 với việc đơn giản hóa xây dựng hàm phản ứng xung như trong mơ hình VAR tuyến tính, và xây dựng các khoảng tin cậy cho mức độ truyền dẫn của tỷ giá hối đoái danh nghĩa hiệu lực ở mức tin cậy 95% với số lần lặp lại mô phỏng dữ liệu là replications = 1000 như Abdul Aleem, Amine Lahiani (2014).
Ở đây, các biến được thay đổi kí hiệu như sau:
- inf thành infln.
- Output gap thành outgap
- dt_lneer_sa thành dt_neer
- dt_lir thành dt_ir
Kết quả của phản ứng xung độc lập kết hợp qua các phương ngang được trình bày như sau:(cumulative orthogonal impulse responses)
- Chế độ 1(Lower): inft-1 ≤ 0.226:
Ở chế độ lạm phát thấp, cú sốc tích cực 1 đơn vị độ lệch chuẩn của dt_lneer_sa làm cho tỷ giá hối đoái danh nghĩa hiệu lực tăng lên khơng đáng kể ở tháng 1, sau đó suy giảm dần và gần như hoàn toàn dập tắt ở tháng thứ 10. Trong khi đó, khe hở sản lượng có dao động khơng đáng kể quanh giá trị 0 và hồn tồn dập tắt ở tháng thứ 4. Cịn đối với tỷ lệ lạm phát thì có sự thay đổi luân phiên tăng lên và giảm xuống xung quanh giá trị 0.2% trong 7 tháng đầu và duy trì khơng đổi trong những tháng tiếp theo, điều này ngụ ý rằng trong môi trường lạm phát thấp, tỷ lạm phát luôn ổn định, kỳ vọng về tỷ lệ lạm phát giống nhau giữa công chúng và nhà làm chính sách. Kết quả này cũng phù hợp với kết quả của các nghiên cứu được trình bày ở phần review. Cuối cùng, lãi suất tiền gửi kì hạn 1 năm cũng thay đổi gia tăng lên đến 10% trong 7 tháng và sau đó cũng duy trì khơng đổi, cho thấy rằng chính sách tiền tệ thắt chặt phản ứng tức thì với sự tăng lên của tỷ lệ lạm phát.
Hình 2.2: Chế độ tỷ lệ lạm phát thấp: Phản ứng xung độc lập kết hợp
(cumulative orthogonal impulse responses) của inf, Ouput gap, dt_lneer_sa và dt_lir đối với cú sốc tích cực 1 đơn vị độ lệch chuẩn dt_lneer_sa.
Nguồn: Kết quả tính tốn từ phần mềm R-packages tác giả thực hiện.
- Chế độ 2(Middle): 0.226 < inft-1 ≤ 0.849
Ở chế độ vừa phải của lạm phát, trước cú sốc tích cực 1 đơn vị độ lệch chuẩn của tỷ giá hối đối danh nghĩa hiệu lực thì tỷ giá danh nghĩa hiệu lực tăng lên và duy trì ở mức tăng thấp ổn định 2%, trong khi đó thì thay đổi Output gap là không đáng kể. Tỷ lệ lạm phát giảm xuống và duy trì ở mức giảm ổn định -0.02% bắt đầu từ tháng 9 trở đi, còn lãi suất tiền gửi kì hạn 1 năm thống kê theo tháng thì cũng có một sự
xy$x d t_ n e e r 0 .0 0 .1 0 .2 0 .3 xy$x o u tg a p 0 .0 0 .1 0 .2 0 .3 xy$x in fl n 0 .0 0 .1 0 .2 0 .3 xy$x d t_ ir 0 .0 0 .1 0 .2 0 .3 1 3 5 7 9 11 14 17 20 23 26 29 32 35
Regime Lower-Impulse Response from dt_neer (cumulative)
sụt giảm và duy trì ở mức giảm ổn định ở các tháng là -1% cũng bắt đầu từ tháng 9 trở đi, ngụ ý một chính sách tiền tệ mở rộng và ổn định.
Hình 2.3: Chế độ giữa của tỷ lệ lạm phát: Phản ứng xung độc lập kết hợp
(cumulative orthogonal impulse responses) của inf, Ouput gap , dt_lneer_sa và dt_lir đối với cú sốc tích cực 1 đơn vị độ lệch chuẩn dt_lneer_sa.
Nguồn: Kết quả tính tốn từ phần mềm R-packages tác giả thực hiện.
- Chế độ 3(Upper): 0.849 < inft-1
Ở chế độ lạm phát cao, thì cú sốc tích cực 1 đơn vị độ lệch chuẩn của tỷ giá hối đoái danh nghĩa hiệu lực làm cho tỷ giá hối đoái danh nghĩa hiệu lực tăng lên và duy trì ổn định ở mức không đáng kể, Output gap hầu như không thay đổi.
xy$x d t_ n e e r -0 .2 0 .0 0 .1 xy$x o u tg a p -0 .2 0 .0 0 .1 xy$x in fln -0 .2 0 .0 0 .1 xy$x d t_ ir -0 .2 0 .0 0 .1 1 3 5 7 9 11 14 17 20 23 26 29 32 35
Regime Middle-Impulse Response from dt_neer (cumulative)
Hình 2.4: Chế độ tỷ lệ lạm phát cao: Phản ứng xung độc lập kết hợp (cumulative orthogonal impulse responses) của inf, Ouput gap, dt_lneer_sa và dt_lir đối với cú sốc tích cực 1 đơn vị độ lệch chuẩn dt_lneer_sa.
Nguồn: Kết quả tính tốn từ phần mềm R-packages tác giả thực hiện.
Tỷ lệ lạm phát giảm nhiều hơn qua các tháng, giảm nhiều nhất ở tháng 12, 13 và sau đó có sự tăng lên đến tháng 29 nhưng vẫn còn giữ ở mức giảm và trở lại giảm xuống ở các tháng tiếp theo, điều này phù hợp với các nghiên cứu trước đây là khi tỷ lệ lạm phát cao thì mức độ truyền dẫn của tỷ giá hối đoái đến lạm phát nhiều hơn; trong khi đó thì lãi suất tiền gửi kì hạn 1 năm thống kê theo tháng cũng có một sự sụt giảm ở các tháng, nhưng có một sự phản ứng trễ nhất định trong lãi suất khi tỷ lệ lạm phát có xu hướng tăng lên, điều này thể hiện một chính sách tiền tệ mở rộng phù hợp với sự giảm xuống của tỷ lệ lạm phát và mất thời gian ngắn chính sách tiền tệ thắt chặt mới có hiệu lực khi tỷ lệ lạm phát đang có xu hướng tăng lên.
xy$x d t_ n e e r -3 .0 -2 .0 -1 .0 0 .0 xy$x o u tg a p -3 .0 -2 .0 -1 .0 0 .0 xy$x in fln -3 .0 -2 .0 -1 .0 0 .0 xy$x d t_ ir -3 .0 -2 .0 -1 .0 0 .0 1 3 5 7 9 11 14 17 20 23 26 29 32 35
Regime Uper-Impulse Response from dt_neer (cumulative)
Kết luận chung cho 3 chế độ lạm phát: với cú sốc tích cực 1 đơn vị độ lệch chuẩn của tỷ giá hối danh nghĩa có hiệu lực thì: Tỷ giá hối đối danh nghĩa hiệu lực tăng lên khơng đáng kể, Output gap gần như không thay đổi. Mức độ truyền dẫn của tỷ giá hối đoái danh nghĩa hiệu lực đến tỷ lệ lạm phát sẽ cao hơn ở tỷ lệ lạm phát cao so với tỷ lệ lạm phát ở mức thấp, và có một sự truyền dẫn ổn định đến lạm phát khi tỷ lệ lạm phát ở mức thấp. Lãi suất tiền gửi kì hạn 1 năm thống kê theo tháng cũng có một sự phản ứng tương ứng phù hợp trong từng trường hợp của tỷ lệ lạm phát hay nói cách khác một chính sách tiền tệ linh hoạt phù hợp.
- VAR chuẩn tuyến tính:
Trong trường hợp này, chúng ta khơng biết được phản ứng khác nhau ở các mức độ tỷ lệ lạm phát. Cú sốc tích cực 1 đơn vị độ lệch chuẩn của tỷ giá danh nghĩa hiệu lực làm cho tỷ giá hối đối danh nghĩa hiệu lực tăng lên khơng đáng kể đến tháng thứ 5 và duy trì ở độ đó ở các tháng cịn lại, cịn Output gap thì hầu như khơng thay đổi qua các tháng; như vậy phản ứng của tỷ giá hối đoái danh nghĩa hiệu lực và Output gap giống với trong trường hợp xem xét cụ thể qua 3 mức độ của tỷ lệ lạm phát. Phản ứng của tỷ lệ lạm phát là một sự sụt giảm đến tháng thứ 5 là - 0.5% và duy trì ổn định mức đó ở các tháng cịn lại, hay nói cách khác mức độ truyền dẫn của tỷ giá đến tỷ lệ lạm phát không nhiều và ổn định hơn so với trong trường hợp tỷ lệ lạm phát cao. Phản ứng của lãi suất tiền gửi kì hạn 1 năm thống kê theo tháng là một sự giảm nhưng không nhiều, và thay đổi tương ứng phù hợp với tình hình thay đổi của tỷ lệ lạm phát, ngụ ý rằng một chính sách tiền tệ mở rộng và các ngân hàng đưa ra mức lãi suất tiền gửi thích hợp với diễn biến thay đổi của tỷ lệ lạm phát.
Như vậy, kết hợp với 3 trường hợp xem xét của tỷ lệ lạm phát ở trên thì cho thấy rằng cơng cụ lãi suất có vai trị trong việc điều chỉnh mức độ thay đổi của tỷ lệ lạm phát, điều này cũng phù hợp với thực tế đã xảy ra trong quá khứ ở Việt Nam như những tháng cuối năm 2007 và năm 2008 lãi suất tăng lên rất cao để kìm hãm và điều chỉnh giảm xuống tỷ lệ lạm phát khi tỷ lệ lạm phát đang ở mức cao. Tỷ giá hối