Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
CHÖÔNG 3
Moâ Hình Hoài Quy Tuyeán Tính Ñôn
ÔÛ chöông 1 phaùt bieåu raèng böôùc ñaàu tieân trong phaân tích kinh teá löôïng laø vieäc thieát laäp moâ hình moâ taû ñöôïc haønh vi cuûa caùc ñaïi löôïng kinh teá. Tieáp theo ñoù nhaø phaân tích kinh teá/ kinh doanh seõ thu thaäp nhöõng döõ lieäu thích hôïp vaø öôùc löôïc moâ hình nhaèm hoã trôï cho vieäc ra quyeát ñònh. Trong chöông naøy seõ giôùi thieäu moâ hình ñôn giaûn nhaát vaø phaùt trieån caùc phöông phaùp öôùc löôïng, phöông phaùp kieåm ñònh giaû thuyeát vaø phöông phaùp döï baùo. Moâ hình naøy ñeà caäp ñeán bieán ñoäc laäp (Y) vaø moät bieán phuï thuoäc (X). Ñoù chính laø moâ hình hoài quy tuyeán tính ñôn. Maëc duø ñaây laø moät moâ hình ñôn giaûn, vaø vì theá phi thöïc teá, nhöng vieäc hieåu bieát nhöõng vaán ñeà cô baûn trong moâ hình naøy laø neàn taûng cho vieäc tìm hieåu nhöõng moâ hình phöùc taïp hôn. Thöïc teá, moâ hình hoài quy ñôn tuyeán tính coù theå giaûi thích cho nhieàu phöông phaùp kinh teá löôïng. Trong chöông naøy chæ ñöa ra nhöõng keát luaän caên baûn veà moâ hình hoài quy tuyeán tính ñôn bieán. Coøn nhöõng phaàn khaùc vaø phaàn tính toaùn seõ ñöôïc giôùi thieäu ôû phaàn phuï luïc. Vì vaäy, ñoái vôùi ngöôøi ñoïc coù nhöõng kieán thöùc caên baûn veà toaùn hoïc, neáu thích, coù theå ñoïc phaàn phuï luïc ñeå hieåu roõ hôn veà nhöõng keát quaû lyù thuyeát.
3.1 Moâ Hình Cô Baûn
Chöông 1 ñaõ trình baøy ví duï veà moâ hình hoài quy ñôn ñeà caäp ñeán moái lieân heä giöõa giaù cuûa moät ngoâi nhaø vaø dieän tích söû duïng (xem Hình 1.2). Choïn tröôùc moät soá loaïi dieän tích, vaø sau ñoù lieät keâ soá löôïng nhaø coù trong toång theå töông öùng vôùi töøng dieän tích ñaõ choïn. Sau ñoù tính giaù baùn trung bình cuûa moãi loaïi nhaø vaø veõ ñoà thò (quy öôùc caùc ñieåm ñöôïc bieåu thò laø X). Giaû thuyeát cô baûn trong moâ hình hoài quy tuyeán tính ñôn laø caùc trò trung bình naøy seõ naèm treân moät ñöôøng thaúng (bieåu thò baèng α + βSQFT), ñaây laø haøm hoài quy cuûa toång theå vaø laø trung bình coù ñieàu kieän (kyø voïng) cuûa GIAÙ theo SQFT cho tröôùc. Coâng thöùc toång quaùt cuûa moâ hình hoài quy tuyeán tính ñôn döïa treân Giaû thieát 3.1 seõ laø
GIAÛ THIEÁT 3.1 (Tính Tuyeán Tính cuûa Moâ Hình)
(3.1)
Yt = α + βXt + ut
trong ñoù, Xt vaø Yt laø trò quan saùt thöù t (t = 1 ñeán n) cuûa bieán ñoäc laäp vaø bieán phuï thuoäc, tieáp theo α vaø β laø caùc tham soá chöa bieát vaø seõ ñöôïc öôùc löôïng;
Ramu Ramanathan
1
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
vaø ut laø soá haïng sai soá khoâng quan saùt ñöôïc vaø ñöôïc giaû ñònh laø bieán ngaãu nhieân vôùi moät soá ñaëc tính nhaát ñònh maø seõ ñöôïc ñeà caäp kyõ ôû phaàn sau. α vaø β ñöôïc goïi laø heä soá hoài quy. (t theå hieän thôøi ñieåm trong chuoãi thôøi gian hoaëc laø trò quan saùt trong moät chuoãi döõ lieäu cheùo.)
Thuaät ngöõ ñôn trong moâ hình hoài quy tuyeán tính ñôn ñöôïc söû duïng ñeå chæ raèng chæ coù duy nhaát moät bieán giaûi thích (X) ñöôïc söû duïng trong moâ hình. Trong chöông tieáp theo khi noùi veà moâ hoài quy ña bieán seõ boå sung theâm nhieàu bieán giaûi thích khaùc. Thuaät ngöõ hoài quy xuaát phaùt töø Fraccis Galton (1886), ngöôøi ñaët ra moái lieân heä giöõa chieàu cao cuûa nam vôùi chieàu cao cuûa ngöôøi cha vaø quan saùt thöïc nghieäm cho thaáy coù moät xu höôùng giöõa chieàu cao trung bình cuûa nam vôùi chieàu cao cuûa nhöõng ngöôøi cha cuûa hoï ñeå “hoài quy” (hoaëc di chuyeån) cho chieàu cao trung bình cuûa toaøn boä toång theå. α + βXb goïi laø phaàn xaùc ñònh cuûa moâ hình vaø laø trung bình coù ñieàu kieän cuûa Y theo X, ñoù laø E(YtXt) = α + βXt. Thuaät ngöõ tuyeán tính duøng ñeå chæ raèng baûn chaát cuûa caùc thoâng soá cuûa toång theå α vaø β laø tuyeán tính (baäc nhaát) chöù khoâng phaûi laø Xt tuyeán tính. Do ñoù, moâ hình vaãn ñöôïc goïi laø hoài quy quyeán X
Y
u
+
=
+
2βα t
t
t
tính ñôn maëc daàu coù X bình phöông. Sau ñaây laø ví duï veà phöông trình hoài quy phi tuyeán tính Yt = α + Xβ + ut. Trong cuoán saùch naøy seõ khoâng ñeà caäp ñeán moâ hình hoài quy phi tuyeán tính maø chæ taäp trung vaøo nhöõng moâ hình coù tham soá coù tính tuyeán tính maø thoâi. Nhöõng moâ hình tuyeán tính naøy coù theå bao goàm caùc soá haïng phi tuyeán tính ñoái vôùi bieán giaûi thích (Chöông 6). Ñeå nghieân cöùu saâu hôn veà moâ hình hoài quy phi tuyeán tính, coù theå tham khaûo caùc taøi lieäu: Greene (1997), Davidson vaø MacKinnon (1993), vaø Griffths, Hill, vaø Judg (1993).
Soá haïng sai soá ut (hay coøn goïi laø soá haïng ngaãu nhieân) laø thaønh phaàn ngaãu nhieân khoâng quan saùt ñöôïc vaø laø sai bieät giöõa Yt vaø phaàn xaùc ñònh α + βXt. Sau ñaây moät toå hôïp cuûa boán nguyeân nhaân aûnh höôûng khaùc nhau:
1. Bieán boû soùt. Giaû söû moâ hình thöïc söï laø Yt = α + βXt + γZt +vt trong ñoù, Zt laø moät bieán giaûi thích khaùc vaø vt laø soá haïng sai soá thöïc söï, nhöng neáu ta söû duïng moâ hình laø Y = α + βXt +ut thì ut = γZt +vt. Vì theá, ut bao haøm caû aûnh höôûng cuûa bieán Z bò boû soùt. Trong ví duï veà ñòa oác ôû phaàn tröôùc, neáu moâ hình thöïc söï bao goàm caû aûnh höôûng cuûa phoøng nguû vaø phoøng taém vaø chuùng ta ñaõ boû qua hai aûnh höôûng naøy maø chæ xeùt ñeán dieän tích söû duïng thì soá haïng u seõ bao haøm caû aûnh höôûng cuûa phoøng nguû vaø phoøng taém leân giaù baùn nhaø.
u
+=
+
+
2. Phi tuyeán tính. ut coù theå bao goàm aûnh höôûng phi tuyeán tính trong moái quan ,
heä giöõa Y vaø X. Vì theá, neáu moâ hình thöïc söï laø
2γβα X X t t
Y t
t
Ramu Ramanathan
2
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
nhöng laïi ñöôïc giaû ñònh baèng phöông trình Y = α + βXt +ut , thì aûnh höôûng cuûa
seõ ñöôïc bao haøm trong ut.
2 tX
3. Sai soá ño löôøng. Sai soá trong vieäc ño löôøng X vaø Y coù theå ñöôïc theå hieän qua u. Ví duï, giaû söû Yt giaù trò cuûa vieäc xaây döïng môùi vaø ta muoán öôùc löôïng haøm Yt = α + βrt +vt trong ñoù rt laø laõi suaát nôï vay vaø vt laø sai soá thaät söï (ñeå ñôn giaûn, aûnh höôûng cuûa thu nhaäp vaø caùc bieán khaùc leân ñaàu tö ñeàu ñöôïc loaïi boû). Tuy nhieân khi thöïc hieän öôùc löôïng, chuùng ta laïi söû duïng moâ hình Yt = α + βXt +ut trong ñoù Xt = rt +Zt laø laõi suaát caên baûn. Nhö vaäy thì laõi suaát ñöôïc ño löôøng trong sai soá Zt thay rt = Xt – Zt vaøo phöông trình ban ñaàu, ta seõ ñöôïc
Yt = α +β(Xt – Zt) +vt = α + βXt – βZt + vt = α + βXt + ut
Caàn luoân löu yù raèng tính ngaãu nhieân cuûa soá haïng ut bao goàm sai soá khi ño löôøng laõi suaát nôï vay moät caùch chính xaùc.
4. Nhöõng aûnh höôûng khoâng theå döï baùo. Duø laø moät moâ hình kinh teá löôïng toát cuõng coù theå chòu nhöõng aûnh höôûng ngaãu nhieân khoâng theå döï baùo ñöôïc. Nhöõng aûnh höôûng naøy seõ luoân ñöôïc theå hieän qua soá haïng sai soá ut.
Nhö ñaõ ñeà caäp ban ñaàu, vieäc thöïc hieän ñieàu tra toaøn boä toång theå ñeå xaùc ñònh haøm hoài quy cuûa toång theå laø khoâng thöïc teá. Vì vaäy, trong thöïc teá, ngöôøi phaân tích thöôøng choïn moät maãu bao goàm caùc caên nhaø moät caùch ngaãu nhieân vaø ño löôøng caùc ñaëc tính cuûa maãu naøy ñeå thieát laäp haøm hoài quy cho maãu. Baûng 3.1 trình baøy döõ lieäu cuûa moät maãu goàm 14 nhaø baùn trong khu vöïc San Diego. Soá lieäu naøy coù saün trong ñóa meàm vôùi teân taäp tin laø DATA3-1. Trong Hình 3.1, caùc caëp giaù trò (Xt, Yt) ñöôïc veõ treân ñoà thò. Ñoà thò naøy ñöôïc goïi laø ñoà thò phaân taùn cuûa maãu cho caùc döõ lieäu. Hình 3.1 töông töï nhö Hình 1.2, nhöng trong Hình 1.2 lieät keâ toaøn boä caùc giaù trò (Xt, Yt) cuûa toång theå, coøn trong Hình 3.1 chæ lieät keâ döõ lieäu cuûa maãu maø thoâi. Giaû söû, taïi moät thôøi ñieåm, ta bieát ñöôïc giaù trò cuûa α vaø β. Ta coù theå veõ ñöôïc ñöôøng thaúng α + βX treân bieåu ñoà. Ñaây chính laø ñöôøng hoài quy cuûa toång theå. Khoaûng caùch chieáu thaúng xuoáng töø giaù thöïc (Yt) ñeán ñöôøng hoài quy α + βX laø sai soá ngaãu nhieân ut. Ñoä doác cuûa ñöôøng thaúng (β) cuõng laø ∆Y/∆X, laø löôïng thay ñoåi cuûa Y treân moät ñôn vò thay ñoåi cuûa X. Vì vaäy β ñöôïc dieãn dòch laø aûnh höôûng caän bieân cuûa X leân Y. Do ñoù, neáu laø β laø 0.14, ñieàu ñoù coù nghóa laø moät meùt vuoâng dieän tích taêng theâm seõ laøm taêng giaù baùn nhaø leân, ôû möùc trung bình, 0.14 ngaøn ñoâ la (löu yù ñôn vò tính) hay 140 ñoâ la. Moät caùch thöïc teá hôn, khi dieän tích söû duïng nhaø taêng theâm 100 meùt vuoâng thì hy voïng raèng giaù baùn trung bình cuûa ngoâi nhaø seõ taêng theâm $14.000 ñoâ la. Maëc daàu α laø tung ñoä goác vaø laø giaù trò cuûa trò trung bình Y khi X baèng 0, soá haïng naøy vaãn khoâng theå ñöôïc hieåu nhö laø giaù trung bình cuûa moät loâ ñaát troáng. Nguyeân nhaân laø vì α cuõng aån chöùa bieán boû soùt vaø do ñoù khoâng coù caùch giaûi thích cho α (ñieàu naøy ñöôïc ñeà caäp kyõ hôn trong Phaàn 4.5).
Ramu Ramanathan
3
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
BAÛNG 3.1 Giaù trò trung bình öôùc löôïng vaø trung bình thöïc teá cuûa giaù
t
SQFT
nhaø vaø dieän tích söû duïng (meùt vuoâng) Giaù baùn1
1 2 3 4 5 6 7 8 9 10 11 12 13 14
1.065 1.254 1.300 1.577 1.600 1.750 1.800 1.870 1.935 1.948 2.254 2.600 2.800 3.000
Giaù trung bình öôùc löôïng2 200,386 226,657 233,051 271,554 274,751 295,601 302,551 312,281 321,316 323,123 365,657 413,751 441,551 469,351
199,9 288 235 285 239 293 285 365 295 290 385 505 425 415
HÌNH 3.1 Bieåu Ñoà Phaân Taùn Cuûa Maãu Trình Baøy Moái Lieân Heä Giöõa Giaù vaø SQFT
Y
600
(
t YX ,
)t
500
α + βX
tu
400
300
200
tXβα +
α
100
X
1000
1400
1800
2200
2600
0
3000
tX
HÌNH 3.2 Phöông Trình Hoài Quy cuûa Toång Theå vaø cuûa Maãu
1 Ñôn vò tính: 1.000 ñoâ la 2 Phöông phaùp tính giaù trung bình öôùc löôïng seõ ñöôïc trình baøy ôû Phaàn 3.2
4
Ramu Ramanathan
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Y
(Hoài qui maãu)
ˆ +
α
ˆ Xβ
D
(
)t
t YX , tuˆ
tu
C
α + βX
(Hoài qui toång theå)
B
=
βα ˆ ˆ X +
ˆ Y t
t
X + βα
=
|
(
t
XYE t
)t
X
0
A
Muïc tieâu ñaàu tieân cuûa moät nhaø kinh teá löôïng laø laøm sao söû duïng döõ lieäu thu thaäp ñöôïc ñeå öôùc löôïng haøm hoài quy cuûa toång theå, ñoù laø, öôùc löôïng tham soá cuûa toång theå α vaø β. Kyù hieäuαˆ laø öôùc löôïng maãu cuûa α vaø βˆ laø öôùc löôïng ^ = α^ + β^X. Ñaây maãu cuûa β. Khi ñoù moái quan heä trung bình öôùc löôïng laø Y ñöôïc goïi laø haøm hoài quy cuûa maãu. ÖÙng vôùi moät giaù trò quan saùt cho tröôùc t, ta ^ t = α^ + β^Xt. Ñaây laø giaù trò döï baùo cuûa Y vôùi moät giaù trò cho tröôùc laø Xt. seõ coù Y Laáy giaù trò quan saùt ñöôïc Yt tröø cho giaù trò naøy, ta seõ ñöôïc öôùc löôïng cuûa ut ñöôïc goïi laø phaàn dö öôùc löôïng, hoaëc ñôn giaûn laø phaàn dö, vaø kyù hieäu laø tuˆ 1vaø ñöôïc theå hieän trong phöông trình sau:
^ u^t = Yt – Y
t = Yt – α^ – β^Xt
Saép xeáp laïi caùc soá haïng treân, ta coù
Y
ˆ u
ˆ ˆ X βα +
=
+
(3.3)
t
t
t
Vieäc phaân bieät giöõa haøm hoài quy cuûa toång theå Y = α + βX vaø haøm hoài quy laø raát quan troïng. Hình 3.2 trình baøy caû hai ñöôøng vaø
cuûa maãu
βα ˆ ˆ X +
=
ˆ Yt
sai soá vaø phaàn dö (caàn nghieân cöùu kyõ vaán ñeà naøy). Löu yù raèng ut laø kyù hieäu chæ “sai soá”, vaøø
tuˆ laø kyù hieäu chæ “phaàn dö”.
BAØI TAÄP 3.1
Xem xeùt caùc phöông trình sau ñaây:
t.
1 Moät soá taùc giaû vaø giaûng vieân thích söû duïng a thay cho α^, b thay cho β^ vaø et thay cho u^ Chuùng ta söû duïng daáu hieäu ^ theo qui ñònh trong lyù thuyeát thoáng keâ vì noù giuùp phaân bieät roõ raøng giöõa giaù trò thaät vaø giaù trò öôùc löôïng vaø cuõng xaùc ñònh ñöôïc thoâng soá ñang ñöôïc öôùc löôïng.
Ramu Ramanathan
5
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
=
uX +
a.
t
Y t
b.
t
Y t
=
uX +
c.
t
d.
+=
e.
=
ˆ uX βα + t ˆ ˆ ˆ uX + + βα
f.
βα + ˆ ˆ ˆ uX += βα + βα ˆ ˆ Y + t βα+=ˆ Yt X ˆ Y t ˆ Y t
t
Giaûi thích kyõ taïi sao phöông trình (a) vaø (b) ñuùng, nhöng (c), (d), (e) vaø
(f) sai. Hình 3.2 raát coù ích trong vieäc traû lôøi caâu hoûi naøy.
3.2 Öôùc löôïng moâ hình cô baûn baèng phöông phaùp bình phöông toái thieåu thoâng
thöôøng
=
βα ˆ ˆ X −−
Trong phaàn tröôùc, ñaõ neâu roõ moâ hình hoài quy tuyeán tính cô baûn vaø phaân bieät giöõa hoài quy cuûa toång theå vaø hoài quy cuûa maãu. Muïc tieâu tieáp theo seõ laø söû duïng caùc döõ lieäu X vaø Y vaø tìm kieám öôùc löôïng “toát nhaát” cuûa hai tham soá cuûa toång theå laø α vaø β. Trong kinh teá löôïng, thuû tuïc öôùc löôïng ñöôïc duøng phoå bieán nhaát laø phöông phaùp bình phöông toái thieåu. Phöông phaùp naøy thöôøng ñöôïc goïi laø bình phöông toái thieåu thoâng thöôøng, ñeå phaân bieät vôùi nhöõng phöông phaùp bình phöông toái thieåu khaùc seõ ñöôïc thaûo luaän trong caùc chöông sau. Kyù hieäu öôùc löôïng cuûa α vaø β laø αˆ vaø βˆ , phaàn dö öôùc löôïng thì baèng ˆ u . Tieâu chuaån toái öu ñöôïc söû duïng bôûi phöông phaùp bình
Y t
t
t
phöông toái thieåu laø cöïc tieåu hoùa haøm muïc tieâu
nt =
nt =
2
ESS
ˆ u
)
)ˆ,ˆ( βα
=
=
ˆ ˆ X −− βα
t
Y ( t
2 t
∑
∑
t
t
1 =
1 =
vôùi caùc tham soá chöa bieát laø αˆ vaø βˆ . ESS laø toång caùc phaàn dö bình phöông vaø phöông phaùp OLS cöïc tieåu toång caùc phaàn dö bình phöông2. Caàn neân löu yù raèng ESS laø khoaûng caùch bình phöông ñöôïc ño löôøng töø ñöôøng hoài quy. Söû duïng khoaûng caùch ño löôøng naøy, coù theå noùi raèng phöông phaùp OLS laø tìm ñöôøng thaúng “gaàn nhaát” vôùi döõ lieäu treân ñoà thò.
Tröïc quan hôn, giaû söû ta choïn moät taäp hôïp nhöõng giaù trò αˆ vaø βˆ , ñoù laø
Xβα ˆ ˆ −
moät ñöôøng thaúng
. Coù theå tính ñöôïc ñoä leäch cuûa Yt töø ñöôøng thaúng
2 Raát deã nhaàm khi goïi ESS laø toång cuûa caùc phaàn dö bình phöông, nhöng kyù hieäu naøy ñöôïc söû duïng phoå bieán trong nhieàu chöông trình maùy tính noåi tieáng vaø coù töø taøi lieäu veà Phaân tích phöông sai
Ramu Ramanathan
6
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
ˆ u
=
βα ˆ ˆ X −−
ñöôïc choïn theo phaàn dö öôùc löôïng
. Sau ñoù bình phöông giaù
t
Y t
trò naøy vaø coäng taát caû caùc giaù trò bình phöông cuûa toaøn boä maãu quan saùt. Toång caùc phaàn dö bình phöông cuûa caùc trò quan saùt [ñöôïc xem nhö toång bình phöông sai soá (ESS)] do ñoù seõ baèng ∑ 2ˆtu . Töông öùng vôùi moät ñieåm treân ñöôøng thaúng seõ coù moät moät trò toång bình phöông sai soá. Phöông phaùp bình phöông toái thieåu choïn nhöõng giaù trò αˆ vaø βˆ sao cho ESS laø nhoû nhaát.
Vieäc bình phöông sai soá ñaït ñöôïc hai ñieàu sau. Thöù nhaát, bình phöông giuùp loaïi boû daáu cuûa sai soá vaø do ñoù xem sai soá döông vaø sai soá aâm laø nhö nhau. Thöù hai, bình phöông taïo ra söï baát lôïi cho sai soá lôùn moät caùch ñaùng keå. Ví duï, giaû söû phaàn dö cuûa maãu laø 1, 2, –1 vaø –2 cuûa heä soá hoài quy choïn tröôùc trò αˆ vaø βˆ choïn tröôùc. So saùnh caùc giaù trò naøy vôùi moät maãu khaùc coù phaàn dö laø –1, –1, –1 vaø 3. Toång giaù trò sai soá tuyeät ñoái ôû caû hai tröôøng hôïp laø nhö nhau. Maëc duø maãu choïn thöù hai coù sai soá tuyeät ñoái thaáp hôn töø 2 ñeán 1, ñieàu naøy daãn ñeán sai soá lôùn khoâng mong muoán laø 3. Neáu ta tính ESS cho caû hai tröôøng hôïp thì ESS cuûa tröôøng hôïp ñaàu laø 10 (12 + 22+ 12+ 22), ESS cho tröôøng hôïp sau laø 12 (12 + 12+ 12+ 32). Phöông phaùp bình phöông toái thieåu aùp ñaët söï baát lôïi lôùn cho sai soá lôùn vaø do ñoù ñöôøng thaúng trong tröôøng hôïp ñaàu seõ ñöôïc choïn. Phaàn 3.3 seõ tieáp tuïc trình baøy nhöõng ñaëc tính caàn thieát khaùc cuûa phöông phaùp cöïc tieåu ESS.
Phöông Phaùp Thích Hôïp Cöïc Ñaïi
Phaàn naøy chæ ñeà caäp sô veà phöông phaùp thích hôïp cöïc ñaïi. Phöông phaùp naøy seõ ñöôïc trình baøy chi tieát ôû phaàn 2.A.4. Phaàn 3.A.5 seõ trình baøy nguyeân taéc aùp duïng moâ hình hoài quy tuyeán tính ñôn. Maëc duø phöông phaùp thích hôïp cöïc ñaïi döïa treân moät tieâu chuaån toái öu khaùc, nhöng caùc thoâng soá öôùc löôïng vaãn gioáng nhö caùc thoâng soá öôùc löôïng ôû phöông phaùp OLS. Noùi ñôn giaûn, phöông phaùp thích hôïp cöïc ñaïi choïn öôùc löôïng sao cho xaùc suaát xaûy ra cuûa maãu quan saùt laø lôùn nhaát.
Phaàn thaûo luaän tröôùc cho thaáy neáu thöïc hieän hai phöông phaùp öôùc löôïng α vaø β khaùc nhau moät caùch chính xaùc thì ñeàu daãn ñeán cuøng moät keát quaû. Nhö vaäy thì taïi sao caàn phaûi xem xeùt caû hai phöông phaùp? Caâu traû lôøi laø trong caùc chöông sau, ta seõ thaáy raèng khi moät soá giaû thieát cuûa moâ hình ñöôïc giaûm nheï, thì thöïc teá, hai phöông phaùp öôùc löôïng khaùc nhau seõ cho keát quaû khaùc nhau. Moät phöông phaùp khaùc coù theå cho keát quaû khaùc nöõa, ñoù laø phöông phaùp cöïc tieåu toång sai soá tuyeät ñoái ∑ tuˆ . Nhöng phöông phaùp naøy khoâng ñöôïc duøng phoå bieán trong kinh teá löôïng vì khoù tính toaùn.
Ramu Ramanathan
7
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Phöông Trình Chuaån
Trong phaàn 3.A.3 cuûa phuï luïc, phöông phaùp OLS ñöôïc chính thöùc aùp duïng. Phaàn naøy cho thaáy raèng ñieàu kieän ñeå cöïc tieåu ESS vôùi αˆ vaø βˆ seõ theo hai phöông trình sau ñaây, ñöôïc goïi laø phöông trình chuaån (khoâng coù lieân heä gì ñeán phaân phoái chuaån).
ˆ u
X
)
(
0 ==
ˆ ˆ X −− βα
=
ˆ )ˆ n βα −
−
(3.4)
t
Y ( t
t
Y t
t
∑
∑
∑
∑
= 0
( [ )] ˆ ˆ X βα−−
(3.5)
t
∑
ˆ
Trong Phöông trình (3.4), caàn löu yù raèng ∑ = αα ˆ
n bôûi vì moãi soá haïng seõ coù moät αˆ vaø coù n soá haïng. Chuyeån veá caùc soá haïng aâm trong Phöông trình (3.4) sang phaûi vaø chia moïi soá haïng cho n, ta ñöôïc
uX t )ˆ t YX ( t t = ∑
X βα ˆ ˆ +=
(3.6)
t
∑
∑
(1/n)ΣYt laø trung bình maãu cuûa Y, kyù hieäu laø Y , vaø (1/n)ΣYt laø trung bình maãu cuûa X, kyù hieäu laø X . Söû duïng keát quaû naøy thay vaøo Phöông trình (3.6), ta ñöôïc phöông trình sau
Y t 1 n 1 n
Y βαˆ += X
(3.7)
)YX ,
Ñöôøng thaúng α^ +β^ X laø ñöôøng öôùc löôïng vaø laø ñöôøng hoài quy cuûa maãu, hoaëc ñöôøng thaúng thích hôïp. Coù theå thaáy raèng töø Phöông trình (3.7) ñöôøng hoài quy cuûa maãu ñi qua ñieåm trung bình ( . Trong Baøi taäp 3.12c, ta seõ thaáy raèng tính chaát naøy khoâng ñaûm baûo tröø khi soá haïng haèng soá α coù trong moâ hình.
Töø Phöông trình (3.5), coäng taát caû theo töøng soá haïng, vaø ñöa αˆ vaø βˆ ra
laøm thöøa soá chung, ta ñöôïc
t
2 = t
∑
∑
∑
( ) 0 X X − ˆ α − ˆ β YX t t
hay
( ) X X = ˆ α + ˆ β
(3.8)
t
2 t
∑
∑
∑
YX t t
Lôøi Giaûi veà Phöông Trình Chuaån
Ramu Ramanathan
8
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Ñeå thuaän lôïi cho vieäc ñaùp aùn veà hai phöông trình chuaån, caùc tính chaát sau ñaây laø raát caàn thieát. Nhöõng tính chaát naøy ñöôïc chöùng minh trong Phuï luïc Phaàn 3.A.2
TÍNH CHAÁT 3.1
– 2 – n(X )2 = ∑Xt )2 = ∑Xt – Sxx = ∑(Xt – X 1 2 – (∑Xt)2 n
TÍNH CHAÁT 3.2
– Y – Sxy = ∑(Xt – X – ) = (∑XtYt) – n X
– )(Yt – Y = ∑XtYt – [(∑Xt ) – (∑Yt) / n]
Töø Phöông trình (3.7),
Y X ˆ α = ˆ X β − =
(3.9)
t
∑ − Y t
∑
1 n 1ˆ β n
Thay αˆ vaøo (3.8)
t
t
2 t
∑
∑
∑
∑
∑
( ) X X X = − + ˆ β YX t t Y t 1ˆ β n 1 n
2
X
)
Y t
t
( ∑
X
=
−
YX t t
2 t
∑
∑
( )( ∑∑ X t n
n
) +
ˆβ
Nhoùm caùc soá haïng coù thöøa soá βˆ :
)
t
Y t
−
YX t t
∑
ˆβ
=
2
( )( ∑∑ X n X
)
t
( ∑
X
−
2 t
∑
n
Tìm βˆ ta ñöôïc
S
xy
=βˆ
Söû duïng kyù hieäu ñôn giaûn ñaõ ñöôïc giôùi thieäu ôû Tính chaát 3.1 vaø 3.2, coù theå ñöôïc dieãn taû nhö sau
(3.10)
S
xx
Ramu Ramanathan
9
Thuïc Ñoan/Haøo Thi
trong ñoù
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
2
X
)
t
( 2 ∑
S
=
(3.11)
xx
t
∑ − X
n
)
t
vaø
Y t S = −
(3.12)
xy
∑
( )( ∑∑ X n
YX t t
t
t
t
X X Y y x − − = Kyù hieäu Sxx vaø Sxy coù theå ñöôïc nhôù moät caùch tröïc quan nhö sau, ñònh nghóa = , trong ñoù kyù hieäu thanh ngang chæ trung bình cuûa vaø Y t
maãu. Do ñoù xt vaø yt kyù hieäu ñoä leäch giöõa X vaø Y so vôùi giaù trò X vaø Y trung bình. Keát quaû sau ñaây seõ ñöôïc chöùng minh ôû phaàn Phuï luïc Phaàn 2.A.1 vaø 3.A.2.
2
S
X
X
(
)
∑xt = 0
(3.13)
=
=
−
)2
xx
XX − t
2 t
t
∑∑ 2 x = t
∑
( ∑
1 n
(3.14)
S
(
X
)(
)
X
=
=
−
−
=
−
xy
yx t
t
t
YYX t
YX t t
t
∑
∑
∑
]∑∑ [ ( ) )( Y t
1 n
Sxy laø “toång caùc giaù trò cuûa xt nhaân yt “. Töông töï, Sxx “toång caùc giaù trò cuûa xt
0
=
=
nhaân xt , hay toång cuûa xt bình phöông
XX − t
S xx
∑∑ 2 x ( = t
trình (3.10) neáu Phöông trình (3.9) vaø (3.10) laø lôøi giaûi cho phöông trình chuaån [(3.4) vaø (3.5)] vaø cho ta öôùc löôïng αˆ vaø βˆ cuûa maãu cho tham soá α vaø β cuûa toång theå. Caàn löu yù raèng khoâng theå xaùc ñònh ñöôïc öôùc löôïng cuûa β trong Phöông 2 ) . Sxx baèng khoâng khi vaø chæ khi
moïi xt baèng khoâng, coù nghóa laø khi vaø chæ khi moïi Xt baèng nhau. Ñieàu naøy daãn ñeán giaû thuyeát sau ñaây
GIAÛ THIEÁT 3.2 (Caùc Giaù Trò Quan Saùt X Laø Khaùc Nhau) Khoâng phaûi laø taát caû giaù trò Xt laø baèng nhau. Coù ít nhaát moät giaù trò Xt khaùc so vôùi nhöõng giaù trò coøn laïi. Noùi caùch khaùc, phöông sai cuûa maãu
Var
(
X
)
X
2)
=
−
t∑ X (
n
1
1 −
khoâng ñöôïc baèng khoâng.
Ramu Ramanathan
10
Thuïc Ñoan/Haøo Thi
Ñaây laø moät giaû thieát raát quan troïng vaø luoân luoân phaûi tuaân theo bôûi vì neáu khoâng moâ hình khoâng theå öôùc löôïng ñöôïc. Moät caùch tröïc quan, neáu Xt khoâng ñoåi, ta khoâng theå giaûi thích ñöôïc taïi sao Yt thay ñoåi. Hình 3.3 minh hoïa giaû thuyeát treân baèng hình aûnh. Trong ví duï veà ñòa oác, giaû söû thoâng tin thu thaäp chæ taäp trung moät vaøo loaïi nhaø coù dieän tích söû duïng laø 1.500 meùt vuoâng. Ñoà thò phaân taùn cuûa maãu seõ ñöôïc theå hieän nhö ôû Hình 3.3. Töø ñoà thò coù theå thaáy roõ
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
raèng döõ lieäu naøy khoâng ñaày ñuû cho vieäc öôùc löôïng ñöôøng hoài quy toång theå α +βX.
HÌNH 3.3 Ví Duï veà Giaù Trò X Khoâng Ñoåi
Y
X
0
1,500
Ví duï 3.1
,0
351,52
13875351
X
=
+
cuûa maãu) laø
ˆ Y t
t
.
Theo thuaät ngöõ ñöôïïc duøng phoå bieán trong kinh teá löôïng, neáu ta söû duïng döõ lieäu trong Baûng 3.1 vaø thöïc hieän “hoài quy Y (GIAÙ) theo soá haïng haèng soá vaø X (SQFT)”, ta coù theå xaùc ñònh ñöôïc moái quan heä öôùc löôïng (hay haøm hoài quy tYˆ laø giaù öôùc löôïng trung bình (ngaøn ñoâ la) töông öùng vôùi Xt. (xem Baûng 3.1). Heä soá hoài quy cuûa Xt laø aûnh höôûng caän bieân öôùc löôïng cuûa dieän tích söû duïng ñeán giaù nhaø, ôû möùc trung bình. Do vaäy, neáu dieän tích söû duïng taêng leân moät ñôn vò, giaù trung bình öôùc löôïng kyø voïng seõ taêng theâm 0,13875 ngaøn ñoâ la ($138.75). Moät caùch thöïc teá, cöù moãi 100 meùt vuoâng taêng theâm dieän tích söû duïng, giaù baùn öôùc löôïng ñöôïc kyø voïng taêng theâm, möùc trung bình, $ 13.875.
Ramu Ramanathan
11
Thuïc Ñoan/Haøo Thi
Haøm hoài quy cuûa maãu coù theå ñöôïc duøng ñeå öôùc löôïng giaù nhaø trung bình döïa treân dieän tích söû duïng cho tröôùc (Baûng 3.1 coù trình baøy giaù trung bình ôû coät cuoái.) Do ñoù, moät caên nhaø coù dieän tích 1.800 meùt vuoâng thì giaù baùn kyø voïng trung bình laø $302.551[ = 52,351 + (0,139 × 1.800)]. Nhöng giaù baùn thöïc söï cuûa caên nhaø laø $285.000. Moâ hình ñaõ öôùc löôïng giaù baùn vöôït quaù $17.551. Ngöôïc laïi, ñoái vôùi moät caên nhaø coù dieän tích söû duïng laø 2.600 meùt vuoâng, giaù baùn trung bình öôùc löôïng laø $413.751, thaáp hôn giaù baùn thöïc söï $505.000 moät caùch ñaùng keå. Söï khaùc bieät naøy coù theå xaûy ra bôûi vì chuùng ta ñaõ boû qua caùc yeáu toá aûnh höôûng khaùc leân giaù baùn nhaø. Ví duï, moät ngoâi nhaø coù saân vöôøn roäng vaø/ hay hoà bôi, seõ coù giaù cao hôn giaù trung bình. Ñieàu naøy nhaán maïnh taàm quan troïng trong vieäc nhaän dieän ñöôïc caùc bieán giaûi thích coù theå aûnh höôûng ñeán giaù trò cuûa bieán phuï thuoäc vaø ñöa caùc aûnh höôûng naøy vaøo moâ hình ñöôïc thieát laäp. Ngoaøi ra, raát caàn thieát trong vieäc phaân tích ñoä tin caäy cuûa caùc öôùc
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
löôïng cuûa tung ñoä vaø heä soá ñoä doác trong Phöông trình (3.1), vaø möùc ñoä “thích hôïp” cuûa moâ hình ñoái vôùi döõ lieäu thöïc teá.
BAØI TAÄP 3.2
.55
462
.
515
2 =
753
∑ tX
.9
095
.
985
5,
444
9,
2 =∑ tY
.26=∑ tX .4=∑ tY
Sao cheùp hai coät soá lieäu trong Baûng 3.1 vaøo moät baûng môùi. Trong coät ñaàu tieân cuûa baûng tính sao cheùp caùc giaù trò veà Yt (GIAÙ) vaø Xt (SQFT) trong coät thöù hai. Söû duïng maùy tính vaø tính theâm giaù trò cho hai coät khaùc. Bình phöông töøng giaù trò trong coät thöù hai vaø ñieàn giaù trò ñoù vaøo coät thöù ba (x). Nhaân laàn löôït töøng giaù trò ôû coät thöù nhaát vôùi giaù trò töông öùng ôû coät hai vaø ñieàn keát qua vaøo coät thöù tö (XtYt). Tieáp theo, tính toång cuûa töøng coät vaø ñaùnh giaù caùc toång sau ñaây:
Ñeå traùnh tình traïng quaù nhieàu vaø sai soá laøm troøn, caàn söû duïng caøng nhieàu soá thaäp phaân caøng toát. Sau ñoù, tính Sxy töø Phöông trình (3.12) vaø Sxx töø Phöông trình (3.11). Cuoái cuøng, tính βˆ theo (3.10) vaø αˆ theo (3.9) vaø kieåm tra laïi nhöõng giaù trò ñaõ trình baøy ban ñaàu.
3.3 Tính chaát cuûa caùc öôùc löôïng
Maëc duø phöông phaùp bình phöông cho ra keát quaû öôùc löôïng veà moái quan heä tuyeán tính coù theå phuø hôïp vôùi döõ lieäu saün coù, chuùng ta caàn traû lôøi moät soá caâu hoûi sau. Ví duï, Ñaëc tính thoáng keâ cuûa αˆ vaø βˆ ? Thoâng soá naøo ñöôïc duøng ñeå ño ñoä tin caäy cuûa αˆ vaø βˆ ? Baèng caùch naøo ñeå coù theå söû duïng αˆ vaø βˆ ñeå kieåm ñònh giaû thuyeát thoáng keâ vaø thöïc hieän döï baùo? Sau ñaây chuùng ta seõ ñi vaøo thaûo luaän töøng vaán ñeà treân. Seõ raát höõu ích neáu baïn oân laïi Phaàn 2.6, phaàn naøy ñöa ra toùm taét veà nhöõng tính chaát caàn thieát cuûa thoâng soá öôùc löôïng.
Ramu Ramanathan
12
Thuïc Ñoan/Haøo Thi
Tính chaát ñaàu tieân caàn xem xeùt laø ñoä khoâng thieân leäch. Caàn löu yù raèng trong Phaàn 2.4 caùc thoâng soá öôùc löôïng αˆ vaø βˆ ? töï thaân chuùng laø bieán ngaãu nhieân vaø do ñoù tuaân theo phaân phoái thoáng keâ. Nguyeân nhaân laø vì nhöõng laàn thöû khaùc nhau cuûa moät cuoäc nghieân cöùu seõ cho caùc keát quaû öôùc löôïng thoâng soá khaùc nhau . Neáu chuùng ta laëp laïi nghieân cöùu vôùi soá laàn thöû lôùn, ta coù theå ñaït ñöôïc nhieàu giaù trò öôùc löôïng. Sau ñoù chuùng ta coù theå tính tyû soá soá laàn maø nhöõng öôùc löôïng naøy rôi vaøo moät khoaûng giaù trò xaùc ñònh. Keát quaû seõ seõ cho ra phaân phoái cuûa caùc öôùc löôïng cuûa maãu. Phaân phoái naøy coù giaù trò trung bình
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
vaø phöông sai. Neáu trung bình cuûa phaân phoái maãu laø thoâng soá thöïc söï (trong tröôøng hôïp naøy laø α hoaëc β), thì ñaây laø öôùc löôïng khoâng thieân leäch. Ñoä khoâng thieân leäch roõ raøng laø ñieàu luoân ñöôïc mong muoán bôûi vì, ñieàu ñoù coù nghóa laø, ôû möùc trung bình, giaù trò öôùc löôïng seõ baèng vôùi giaù trò thöïc teá, maëc duø trong moät soá tröôøng hôïp caù bieät thì ñieàu naøy coù theå khoâng ñuùng.
Coù theå noùi raèng thoâng soá öôùc löôïng OLS cuûa α vaø β ñöa ra trong Phaàn 3.2 coù tính chaát khoâng thieân leäch. Tuy nhieân, ñeå chöùng minh ñieàu naøy, chuùng ta caàn ñaët ra moät soá giaû thuyeát boå sung veà Xt vaø ut. Caàn nhôù raèng, maëc duø Giaû thieát 3.1 coù theå vaø ñöôïc giaûm nheï ôû phaàn sau, nhöng Giaû thuyeát 3.2 vaø 3.3 laø luoân luoân caàn thieát vaø phaûi tuaân theo. Sau ñaây laø caùc giaû thieát boå sung caàn thieát.
GIAÛ THIEÁT 3.3 (Sai Soá Trung Bình baèng Zero) Moãi laø u moät bieán ngaãu nhieân vôùi E(u) = 0
Trong Hình 3.1 caàn löu yù raèng moät soá ñieåm quan saùt naèm treân ñöôøng α + βX vaø moät soá ñieåm naèm döôùi. Ñieàu naøy coù nghóa laø coù moät giaù trò sai soá mang daáu döông vaø moät soá sai soá mang daáu aâm. Do α + βX laø ñöôøng trung bình, neân coù theå giaû ñònh raèng caùc sai soá ngaãu nhieân treân seõ bò loaïi tröø nhau, ôû möùc trung bình, trong toång theå. Vì theá, giaû ñònh raèng ut laø bieán ngaãu nhieân vôùi giaù trò kyø voïng baèng 0 laø hoaøn toaøn thöïc teá.
GIAÛ THIEÁT 3.4 (Caùc Giaù Trò X Ñöôïc Cho Tröôùc vaø Khoâng Ngaãu Nhieân) Moãi giaù trò Xt ñöôïc cho tröôùc vaø khoâng laø bieán ngaãu nhieân. Ñieàu naøy ngaàm chæ raèng ñoàng phöông sai cuûa toång theå giöõa Xt vaø ut, Cov(Xt, ut) = E(Xt, ut) – E(Xt)E(ut) = XtE(ut) – XtE(ut) = 0. Do ñoù giöõa Xt vaø ut khoâng coù moái töông quan (xem Ñònh nghóa 2.4 vaø 2.5).
Theo tröïc giaùc, neáu X vaø u coù moái töông quan, thì khi X thay ñoåi, u cuõng seõ thay ñoåi. Trong tröôøng hôïp naøy, giaù trò kyø voïng cuûa Y seõ khoâng baèng α + βX. Neáu giaù trò X laø khoâng ngaãu nhieân thì giaù trò kyø voïng coù ñieàu kieän cuûa Y theo giaù trò X seõ baèng α + βX. Keát quaû cuûa vieäc vi phaïm Giaû thieát 3.4 seõ ñöôïc trình baøy trong phaàn sau, ñaëc bieät laø khi nghieân cöùu moâ hình heä phöông trình (Chöông 13). Tính chaát 3.3 phaùt bieåu raèng khi hai giaû thieát ñöôïc boå sung, thoâng soá öôùc löôïng OLS laø khoâng thieân leäch.
TÍNH CHAÁT 3.3 (Ñoä Khoâng Thieân Leäch)
Ramu Ramanathan
13
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
E
) αα =ˆE
(
, vaø ø. Trong hai giaû thieát boå sung 3.3 vaø 3.4, [E(ut) = 0, Cov(Xt, ut) = 0], thoâng soá öôùc löôïng, thoâng soá öôùc löôïng bình phöông toái thieåu αˆ vaø βˆ laø khoâng thieân leäch; nghóa laø
) ββ ˆ ( ˆ =
CHÖÙNG MINH
E
(Neáu ñoäc giaû khoâng quan taâm ñeán chöùng minh, coù theå boû qua phaàn).
( SE
)xx
( ) =βˆ
E
Töø Phöông trình (3.10), . Nhöng theo Giaû thuyeát 3.4, Xt laø
( SE
)xy
( ) ˆ =β
xy S khoâng ngaãu nhieân vaø do ñoù Sxx cuõng khoâng ngaãu nhieân. Ñieàu naøy coù nghóa laø khi tính giaù trò kyø voïng, caùc soá haïng lieân quan ñeán Xt coù theå ñöôïc ñöa ra ngoaøi 1 S
xx
giaù trò kyø voïng. Vì vaäy, ta coù . Trong Phöông trình (3.12),
thay Yt töø Phöông trình (3.1) vaø thay ∑α baèng nα .
)( nX t
t
t
( ∑
∑
X u βα + + S X u = + −
(3.15)
( X βα +
)
xy
t
t
t
∑
∑ n
)
2
X
u
)
t
t
( ∑
X
X
X
=
+
+
−
−
−
α
β
α
t
2 t
uX t
t
t
∑
∑
∑
∑
n
( )( ∑∑ X t n
)
β
2
u
X
)
t
t
( ∑
X
=
−
+
−
t
uX t
t
∑
∑
( )( ∑∑ X t n
n
)
β
xu
S = β S + xx
trong ñoù Sxx ñöôïc cho bôûi Phöông trình (3.13) vaø
)
t
t
u S =
(3.16)
xu
t
∑ − uX t
( )( ∑∑ X n
X
X
=
−
−
(
uX t
t
t
t
) t uX
∑
∑∑ u =
X
X
=
−
=
−
( SE
)
( uXE
)
( uE
)
)
( uE
) 0=
xu
t
t
t
( uEX t
t
t
∑
∑
∑
∑
Ramu Ramanathan
14
Thuïc Ñoan/Haøo Thi
X laø trung bình maãu cuûa X, Xt laø khoâng ngaãu nhieân, X xuaát hieän ôû moïi soá haïng, vaø kyø voïng cuûa toång caùc soá haïng thì baèng toång caùc giaù trò kyø voïng. Do vaäy,
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
( SE
)
xy
xx
) ( ˆ β
E S = = β . Nhö theo Giaû thieát 3.3. Do ñoù, E(Sxy) = βSxx, nghóa laø
vaäy β laø öôùc löôïng khoâng thieân leäch cuûa β. Chöùng minh töông töï cho α^ . Caàn nhaän thaáy raèng vieäc chöùng minh ñoä khoâng thieân leäch phuï thuoäc chuû yeáu vaøo Giaû thieát 3.4. Neáu E(Xtut) ≠ 0, βˆ coù theå bò thieân leäch. BAØI TAÄP 3.3
~ X1). Löu yù raèng β
Söû duïng Phöông trình (3.9) ñeå chöùng minh raèng αˆ laø khoâng thieân leäch. Neâu roõ caùc giaû thuyeát caàn thieát khi chöùng minh. Maëc daàu ñoä khoâng thieân leäch luoân laø moät tính chaát luoân ñöôïc mong muoán, nhöng töï baûn thaân ñoä khoâng thieân leäch khoâng laøm cho thoâng soá öôùc löôïng “toát”, vaø moät öôùc löôïng khoâng thieân leäch khoâng chæ laø tröôøng hôïp caù bieät. ~ Haõy xem xeùt ví duï sau veà moät thoâng soá öôùc löôïng khaùc laø β = (Y2 – Y1)/(X2 –
ñôn giaûn laø ñoä doác cuûa ñöôøng thaúng noái hai ñieåm (X1, Y1)
~ vaø (X2, Y2). Raát deã nhaän thaáy raèng β
u
+
+
( X + βα
)
)
2
1
1
2
=
=
=
+
~ β
β
Y 2 X
Y 1 X
u 2 X
( X + βα X
u X
u 1 X
− −
laø khoâng thieân leäch
− −
2
1
2
1
2
1
− −
~ cuûa caùc thoâng soá öôùc löôïng khoâng thieân leäch nhö treân. Bôûi vì β loaïi boû caùc giaù trò quan saùt töø 3 ñeán n, moät caùch tröïc quan ñaây khoâng theå laø moät thoâng soá öôùc löôïng “toát”. Trong Baøi taäp 3.6, taát caû caùc giaù trò quan saùt ñöôïc söû duïng theå thieát laäp caùc thoâng soá öôùc löôïng khoâng thieân leäch khaùc, nhöng töông töï nhö treân ñaây khoâng phaûi laø laø thoâng soá öôùc löôïng khoâng thieân leäch toát nhaát. Do ñoù, raát caàn coù nhöõng tieâu chuaån boå sung ñeå ñaùnh giaù “ñoä toát” cuûa moät thoâng soá öôùc löôïng.
Nhö ñaõ noùi tröôùc ñaây, caùc giaù trò X laø khoâng ngaãu nhieân vaø E(u2) = E(u1) = 0. ~ Do ñoù, β laø khoâng thieân leäch. Thöïc ra, ta coù theå xaây döïng moät chuoãi voâ haïn
Ramu Ramanathan
15
Thuïc Ñoan/Haøo Thi
Tieâu chuaån thöù hai caàn xem xeùt laø tính nhaát quaùn, ñaây laø moät tính chaát cuûa maãu lôùn ñaõ ñöôïc ñònh nghóa trong Phaàn 2.6 (Ñònh nghóa 2.10). Giaû söû ta choïn ngaãu nhieân moät maãu coù n phaàn töû vaø ñi tìm αˆ vaø βˆ . Sau ñoù choïn moät maãu lôùn hôn vaø öôùc löôïng laïi caùc thoâng soá naøy. Laëp laïi quaù trình naøy nhieàu laàn ñeå coù ñöôïc moät chuoãi nhöõng thoâng soá öôùc löôïng. Tính nhaát quaùn laø tính chaát ñoøi hoûi caùc thoâng soá öôùc löôïng vaãn phuø hôïp khi côõ maãu taêng leân voâ haïn. Öôùc ~ löôïng β ñöôïc trình baøy ôû treân roõ raøng laø khoâng ñaït ñöôïc tính nhaát quaùn bôûi vì khi côõ maãu taêng leân khoâng aûnh höôûng gì ñeán thoâng soá naøy. Tính chaát 3.4 phaùt bieåu caùc ñieàu kieän ñeå moät öôùc löôïng coù tính nhaát quaùn.
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
TÍNH CHAÁT 3.4 (Tính Nhaát Quaùn)
Theo Giaû thieát (3.2), (3.3) vaø (3.4), öôùc löôïng bình phöông toái thieåu coù tính chaát nhaát quaùn. Do ñoù, ñieàu kieän ñeå ñaït ñöôïc tính nhaát quaùn laø E(ut) = 0, Cov(Xt, ut) = 0 vaø Var(Xt) ≠ 0.
(Neáu ñoäc giaû khoâng quan taâm, coù theå boû qua phaàn naøy.)
CHÖÙNG MINH
xu
ˆ = ββ
+
Töø Phöông trình (3.15) vaø (3.10)
(3.17)
S S
/ /
n n
xx
Theo quy luaät soá lôùn (Tính chaát 2.7a), Sxu/n ñoàng quy vôùi kyø voïng cuûa chính noù, ñoù laø Cov(X, u). Töông töï, Sxx/n ñoàng quy vôùi Var(X). Do vaäy daãn tôùi ñieàu, neáu n hoäi tuï ñeán voâ cuøng, β seõ ñoàng quy vôùi β + [Cov(X,u)/Var(X), vaø seõ baèng β neáu Cov(X,u) = 0 – nghóa laø neáu X vaø u khoâng töông quan. Nhö vaäy, βˆ laø öôùc löôïng nhaát quaùn cuûa β.
Maëc duø βˆ laø khoâng thieân leäch vaø nhaát quaùn, vaãn coù nhöõng tieâu chuaån caàn boå sung bôûi ñeå coù theå xaây döïng öôùc löôïng nhaát quaùn vaø khoâng thieân leäch khaùc. Baøi taäp 3.6 laø moät ví duï veà loaïi öôùc löôïng ñoù. Tieâu chuaån söû duïng tieáp theo laø tính hieäu quaû (ñònh nghóa trong Phaàn 2.6). Noùi moät caùch ñôn giaûn, öôùc löôïng khoâng thieân leäch coù tính hieäu quaû hôn neáu öôùc löôïng naøy coù phöông sai nhoû hôn. Ñeå thieát laäp tính hieäu quaû, caàn coù caùc giaû thieát sau veà ut.
GIAÛ THIEÁT 3.5 (Phöông sai cuûa sai soá khoâng ñoåi) Taát caû giaù trò u ñöôïc phaân phoái gioáng nhau vôùi cuøng phöông sai σ2 Var
)
u (
2 σ=
=
, sao cho . Ñieàu naøy ñöôïc goïi laø phöông sai cuûa sai soá khoâng ñoåi
( uE
)
t
2 t (phaân taùn ñeàu).
GIAÛ THIEÁT 3.6 (Ñoäc Laäp Theo Chuoãi) Giaù trò u ñöôïc phaân phoái ñoäc laäp sao cho Cov(ut, us) = E(utus) = 0 ñoái vôùi moïi t ≠ s. Ñaây ñöôïc goïi laø chuoãi ñoäc laäp.
Ramu Ramanathan
16
Thuïc Ñoan/Haøo Thi
Caùc giaû thieát treân ngaàm chæ raèng caùc phaàn dö phaân coù phaân phoái gioáng nhau vaø phaân phoái ñoäc laäp (iid). Töø Hình 1.2 ta thaáy raèng öùng vôùi moät giaù trò
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
X seõ coù moät giaù trò phaân phoái Y ñeå xaùc ñònh phaân phoái coù ñieàu kieän. Sai soá ut laø ñoä leäch töø trung bình coù ñieàu kieän α + βXt. Giaû thieát 3.5 ngaàm ñònh raèng phaân phoái cuûa ut coù cuøng phöông sai (σ2) vôùi phaân phoái cuûa us cho moät quan saùt khaùc s. Hình 3.4a laø moät ví duï veà phöông sai cuûa sai soá thay ñoåi (hoaëc khoâng phaân taùn ñeàu) khi phöông sai thay ñoåi taêng theo giaù trò quan saùt X. Giaû thuyeát 3.5 ñöôïc giaûm nheï trong Chöông 8. Phaàn 3.6 Phuï chöông coù trình baøy moâ taû ba chieàu cuûa giaû thuyeát naøy.
Giaû thieát 3.6 (seõ ñöôïc giaûm nheï trong Chöông 9) ngaàm ñònh raèng laø ut vaø us ñoäc laäp vaø do vaäy khoâng coù moái töông quan. Cuï theå laø, caùc sai soá lieân tieáp nhau khoâng töông quan nhau vaø khoâng taäp trung. Hình 3.4b laø moät ví duï veà töï töông quan khi giaû thuyeát treân bò vi phaïm. Chuù yù raèng khi caùc giaù trò quan saùt keá tieáp nhau taäp trung laïi, thì coù khaû naêng caùc sai soá seõ coù töông quan.
HÌNH 3.4 Ví Duï veà Phöông Sai Cuûa Sai Soá Thay Ñoåi vaø Töï Hoài Quy
Y
X
Y
X
a. Phöông sai cuûa sai soá thay ñoåi
b. Töï hoài quy
TÍNH CHAÁT 3.5 (Hieäu quaû, BLUE vaø Ñònh lyù Gauss-Markov)
Ramu Ramanathan
17
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Theo Giaû thieát 3.2 ñeán 3.6, öôùc löôïng bình phöông toái thieåu thoâng thöôøng (OLS) laø öôùc löôïng tuyeán tính khoâng thieân leäch coù hieäu quaû nhaát trong caùc öôùc löôïng. Vì theá phöông phaùp OLS ñöa ra Öôùc Löôïng Khoâng Thieân leäch Tuyeán Tính Toát Nhaát (BLUE).
Keát quaû naøy (ñöôïc chöùng minh trong Phaàn 3.A.4) ñöôïc goïi laø Ñònh lyù Gauss–Markov, theo lyù thuyeát naøy öôùc löôïng OLS laø BLUE; nghóa laø trong taát caû caùc toå hôïp tuyeán tính khoâng thieân leäch cuûa Y, öôùc löôïng OLS cuûa α vaø β coù phöông sai beù nhaát.
Toùm laïi, aùp duïng phöông phaùp bình phöông toái thieåu (OLS) ñeå öôùc löôïng heä soá hoài quy cuûa moät moâ hình mang laïi moät soá tính chaát mong muoán sau: öôùc löôïng laø (1) khoâng thieân leäch, (2) coù tính nhaát quaùn vaø (3) coù hieäu quaû nhaát. Ñoä khoâng thieân leäch vaø tính nhaát quaùn ñoøi hoûi phaûi keøm theo Giaû thuyeát E(ut) = 0 vaø Cov(Xt, ut) = 0. Yeâu caàu veà tính hieäu quaû vaø BLUE, thì caàn coù theâm giaû thuyeát, Var(ut) = σ2 vaø Cov(ut, us) = 0, vôùi moïi t ≠ s.
3.4 Ñoä Chính Xaùc cuûa Öôùc Löôïng vaø Möùc Ñoä Thích Hôïp cuûa Moâ Hình
351
.52
13875
,0
ˆ =β
351,52
,0
13875351
X
=
+
Söû duïng caùc döõ lieäu trong ví duï veà ñòa oác ta öôùc löôïng ñöôïc thoâng soá nhö sau . Caâu hoûi cô baûn laø caùc öôùc löôïng naøy toát nhö theá ˆ =α vaø
ˆ Yt
vôùi naøo vaø möùc ñoä thích hôïp cuûa haøm hoài quy maãu
döõ lieäu ra sao. Phaàn naøy seõ thaûo luaän phöông phaùp xaùc ñònh thoâng soá ño löôøng ñoä chính xaùc cuûa caùc öôùc löôïng cuõng nhö ñoä phuø hôïp.
Ñoä Chính Xaùc cuûa Caùc Öôùc Löôïng
Töø lyù thuyeát xaùc suaát ta bieát raèng phöông sai cuûa moät bieán ngaãu nhieân ño löôøng söï phaân taùn xung quanh giaù trò trung bình. Phöông sai caøng beù, ôû möùc trung bình, töøng giaù trò rieâng bieät caøng gaàn vôùi giaù trò trung bình. Töông töï, khi ñeà caäp ñeán khoaûng tin caäy, ta bieát raèng phöông sai cuûa bieán ngaãu nhieân caøng nhoû, khoaûng tin caäy cuûa caùc tham soá caøng beù. Nhö vaäy, phöông sai cuûa moät öôùc löôïng laø thoâng soá ñeå chæ ñoä chính xaùc cuûa moät öôùc löôïng. Do ñoù vieäc tính toaùn phöông sai cuûa αˆ vaø βˆ laø luoân caàn thieát.
2
Do αˆ vaø βˆ thuoäc vaøo caùc giaù trò Y, maø Y laïi phuï thuoäc vaøo caùc bieán ngaãu nhieân u1, u2, …, un, neân chuùng cuõng laø bieán ngaãu nhieân vôùi phaân phoái töông öùng. Sau ñaây caùc phöông trình ñöôïc ruùt ra trong Phaàn 3.A.6 ôû phaàn phuï luïc cuûa chöông naøy.
(3.18)
Var
E
=
=
) ( ˆ ββ −
)ˆ( 2 σβ = β &&
2 σ xxS
Ramu Ramanathan
18
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
X
2
2 t
2
(3.19)
Var
E
=
σ
2 )ˆ( = σα ˆ α
[ ( ) ˆ − αα
]
xx
2
ˆ
(3.20)
Cov
E
)ˆ,ˆ( σβα
σ
=
=
−
−
−=
[ )( ( ˆ ββαα
∑= nS ] )
ˆˆ βα
X xxS
βα ˆ ˆ X +=
t
ˆ u
=
trong ñoù Sxx ñöôïc ñònh nghóa theo Phöông trình (3.11) vaø σ2 laø phöông sai cuûa sai soá. Caàn löu yù raèng neáu Sxx taêng, giaù trò phöông sai vaø ñoàng phöông sai (trò tuyeät ñoái) seõ giaûm. Ñieàu naøy cho thaáy söï bieán thieân ôû X caøng cao vaø côõ maãu caøng lôùn thì caøng toát bôûi vì ñieàu ñoù cho chöùng toû ñoä chính cuûa caùc thoâng soá ñöôïc öôùc löôïng.
ˆ Y t
t
t
n
laø moät öôùc löôïng cuûa ut, vaø laø phaàn dö öôùc löôïng.
ˆ u
2
2
Caùc bieåu thöùc treân laø phöông sai cuûa toång theå vaø laø aån soá bôûi vì σ2 laø aån soá. Tuy nhieân, caùc thoâng soá naøy coù theå ñöôïc öôùc löôïng bôûi vì σ2 coù theå ñöôïc ˆ Y laø ñöôøng thaúng öôùc löôïng. öôùc löôïng döïa treân maãu. Löu yù raèng t βα ˆ ˆ X −− Do ñoù, Moät öôùc löôïng deã thaáy cuûa σ2 laø ˆ 2∑ ut / nhöng öôùc löôïng naøy ngaãu nhieân bò thieân leäch. Moät öôùc löôïng khaùc cuûa σ2 ñöôïc cho sau ñaây (xem chöùng minh ôû Phaàn 3.A.7)
(3.21)
s
=
ˆ σ
=
2 t 2
∑ n −
)
( ∑ − x xi
n – 1, ñaõ ñöôïc thaûo luaän trong Phaàn 2.7. n – 1 ñöôïc aùp duïng do Lyù do chia töû soá cho n – 2 thì töông töï nhö tröôøng hôïp chia chi-square cho coù
ñieàu kieän laø baèng 0. Ñeå aùp duïng chia cho n – 2, caàn coù hai ñieàu kieän bôûi Phöông trình (3.4) vaø (3.5). Caên baäc hai cuûa phöông sai öôùc löôïng ñöôïc goïi laø sai soá chuaån cuûa phaàn dö hay sai soá chuaån cuûa hoài quy. Söû duïng öôùc löôïng naøy, ta tính ñöôïc caùc öôùc löôïng cuûa phöông sai vaø ñoàng phöông sai cuûa αˆ vaø βˆ . Caên baäc hai cuûa phöông sai ñöôïc goïi laø sai soá chuaån cuûa heä soá hoài quy vaø kyù hieäu αˆs vaø βˆs . Phöông sai öôùc löôïng vaø ñoàng phöông sai cuûa heä soá hoài quy öôùc löôïng baèng
s
(3.22)
2 ˆ
β =
2 ˆσ xxS
2 t
2 ˆσ
X
(3.23)
2 ˆ α
2σ ˆ
−=
(3.24)
s ˆˆ βα
s ∑= nS xx X xxS
Ramu Ramanathan
19
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
ˆ Y t
t
t YY t
theå tính ñöôïc phaàn dö
Toùm laïi: Tröôùc tieân, caàn tính heä soá hoài quy öôùc löôïng αˆ vaø βˆ baèng caùch aùp duïng Phöông trình (3.9) vaø (3.10). Keát quaû cho cho moái quan heä öôùc löôïng βα ˆ ˆ X += giöõa Y vaø X. sau ñoù tính giaù trò döï baùo cuûa Yt theo . Töø ñoù, ta coù ˆ− . Sau ñoù tính toaùn öôùc löôïng cuûa phöông tuˆ theo sai cuûa ut döïa theo Phöông trình (3.21). Thay keát quaû vaøo Phöông trình (3.18), (3.19) vaø (3.20), ta ñöôïc giaù trò phöông sai vaø ñoàng phöông sai cuûa αˆ vaø βˆ .
Caàn löu yù raèng ñeå coâng thöùc tính phöông sai cuûa phaàn dö s2 ñöôïc cho trong Phöông trình 3.21 coù yù nghóa, caàn coù ñieàu kieän n > 2. Khoâng coù giaû thuyeát naøy, phöông sai ñöôïc öôùc löôïng coù theå khoâng xaùc ñònh ñöôïc hoaëc aâm. Ñieàu kieän toång quaùt hôn ñöôïc phaùt bieåu trong Giaû thuyeát 3.7, vaø baét buoäc phaûi tuaân theo.
GIAÛ THIEÁT 3.7 (n > 2) Soá löôïng quan saùt (n) phaûi lôùn hôn soá löôïng caùc heä soá hoài quy ñöôïc öôùc löôïng (k). Trong tröôøng hôïp hoài quy tuyeán tính ñôn bieán, thì ñieàu kieän n > 2 khoâng coù.
Ví duï 3.2
285
,37
,0
01873
Sau ñaây laø sai soá chuaån trong ví duï veà giaù nhaø, Sai soá chuaån cuûa phaàn dö = s = σˆ = 39,023 Sai soá chuaån cuûa
ˆ = ˆ =
ˆ = αα s ˆ = ββ s
,0
671
−=
Sai soá chuaån cuûa
ˆ = βαβ s ˆˆ
Ñoàng phöông sai giöõa αˆ vaø
Thöïc haønh maùy tính Phaàn 3.1 cuûa Phuï chöông D seõ cho keát quaû töông töï.
Maëc duø coù caùc ñaïi löôïng ño löôøng soá hoïc veà ñoä chính xaùc cuûa caùc öôùc löôïng, töï thaân caùc ño löôøng naøy khoâng söû duïng ñöôïc bôûi vì caùc ño löôøng naøy coù theå lôùn hoaëc nhoû moät caùch tuøy tieän baèng caùch ñôn giaûn laø thay ñoåi ñôn vò ño löôøng (xem theâm ôû Phaàn 3.6). Caùc ño löôøng naøy ñöôïc söû duïng chuû yeáu trong vieäc kieåm ñònh giaû thuyeát, ñeà taøi naøy seõ ñöôïc thaûo luaän chi tieát ôû Phaàn 3.5.
Ñoä Thích Hôïp Toång Quaùt
Ramu Ramanathan
20
Thuïc Ñoan/Haøo Thi
Hình 3.1 cho thaáy roõ raèng khoâng coù ñöôøng thaúng naøo hoaøn toaøn “thích hôïp” vôùi caùc döõ lieäu bôûi vì coù nhieàu giaù trò döï baùo bôûi ñöôøng thaúng caùch xa vôùi giaù trò thöïc teá. Ñeå coù theå ñaùnh giaù moät moái quan heä tuyeán tính moâ taû nhöõng giaù trò quan saùt coù toát hôn moät moái quan heä tuyeán tính khaùc hay khoâng, caàn phaûi coù
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
2
moät ño löôøng toaùn hoïc ñoä thích hôïp. Phaàn naøy seõ phaùt trieån caùc thoâng soá ño löôøng ñoù.
)1
)
2 ˆ Yσ
( Yt
Y n − − Khi thöïc hieän döï baùo veà moät bieán phuï thuoäc Y, neáu ta chæ coù nhöõng thoâng tin veà caùc giaù trò quan saùt cuûa Y coù ñöôïc töø moät soá phaân phoái xaùc suaát, thì coù leõ caùch toát nhaát coù theå laø laø öôùc löôïng giaù trò trung bình Y vaø phöông sai söû . Neáu caàn döï baùo, moät caùch ñôn giaûn, ta coù theå duïng
] (
[ = ∑
)2
TSS
−
=
( n
)1
ˆ Yσ
=
βα ˆ ˆ X +
söû duïng giaù trò trung bình bôûi vì khoâng coøn thoâng tin naøo khaùc. Sai soá khi döï YYt − . Bình phöông giaù trò naøy vaø tính toång bình baùo quan saùt thöù t baèng phöông cho taát caû maãu, ta tính ñöôïc toång phöông sai cuûa Yt so vôùi Y laø ( ∑ − Y Y . Ñaây laø toång bình phöông toaøn phaàn (TSS). Ñoä leäch chuaån cuûa maãu cuûa Y ño löôøng ñoä phaân taùn cuûa Yt xung quanh giaù trò trung bình cuûa Y, noùi caùch khaùc laø ñoä phaân taùn cuûa sai soá khi söû duïng Y laøm bieán döï baùo, vaø ñöôïc cho nhö sau
ˆ Y t
t
ˆ u
=
Giaû söû ta cho raèng Y coù lieân quan ñeán moät bieán X khaùc theo Phöông trình (3.1). Ta coù theå hy voïng raèng bieát tröôùc giaù trò X seõ giuùp döï baùo Y toát hôn laø chæ duøng Y . Cuï theå hôn laø, neáu ta coù caùc öôùc löôïng αˆ vaø βˆ vaø bieát ñöôïc giaù trò cuûa X laø Xt, nhö vaäy öôùc löôïng cuûa Yt seõ laø . Sai soá cuûa öôùc
t
ˆ YY − t
t
(
n
)2
ESS
=
−
. Bình phöông giaù trò sai soá naøy vaø tính toång caùc sai soá löôïng naøy laø
498
,88
. Giaù trò naøy ño löôøng ñoä phaân taùn cuûa sai soá khi söû duïng
ˆ =Yσ
023,39 ˆ =σ vaø
cho toaøn boä maãu, ta coù ñöôïc toång bình phöông sai soá (ESS), hay toång caùc bình phöông phaàn dö, laø ESS = ∑ 2ˆtu . Sai soá chuaån cuûa caùc phaàn dö laø tYˆ ˆ σ Yσˆ ñöôïc cho ôû treân ñeå xem xeùt laøm bieán döï baùo vaø thöôøng ñöôïc so saùnh vôùi möùc ñoä giaûm xuoáng laø bao nhieâu. Bôûi vì ESS caøng nhoû caøng toát, vaø möùc ñoä giaûm xuoáng caøng nhieàu. Trong ví duï ñöa ra, ø, giaûm hôn phaân nöûa so vôùi giaù trò ban ñaàu.
Phöông phaùp naøy khoâng hoaøn toaøn toát laém, tuy nhieân bôûi vì caùc sai soá chuaån raát nhaïy caûm ñoái vôùi ñôn vò ño löôøng Y neân raát caàn coù moät thoâng soá ño löôøng khaùc khoâng nhaïy caûm vôùi ñôn vò ño löôøng. Vaán ñeà naøy seõ ñöôïc ñeà caäp sau ñaây.
HÌNH 3.5 Caùc Thaønh Phaàn cuûa Y
Ramu Ramanathan
21
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Y
(
t YX ,
)t
tY
ˆ Y
βα ˆ ˆ X +
=
tuˆ
Y
Yt −
Y
Yt −ˆ
Y
X
0
X
tX
tYˆ so vôùi Y (laø giaù trò trung bình
Thoâng soá ño löôøng toång bieán thieân cuûa
)2
( ˆ∑ − Y Yt
tYˆ ) cho toaøn maãu laø (RSS). Phaàn 3.A.8 cho thaáy
2
2
Y
Y
cuûa . Ñöôïc goïi laø toång bình phöông hoài quy
ˆ u − = −
(3.25)
)
( Y t
2 t
( ˆ Y t
∑
) ∑ +
Y
)
Y
)
ˆ u
−
=
−
+
∑
( Y t
ˆ( Y t
t
Xβα ˆ ˆ +
. Hình 3.5 Do vaäy, TSS = RSS + ESS. Löu yù raèng
minh hoïa caùc thaønh phaàn treân. Phöông trình (3.25) phaùt bieåu raèng caùc thaønh phaàn cuõng ñöôïc bình phöông. Neáu moái quan heä giöõa X vaø Y laø “chaët cheõ”, caùc ñieåm phaân taùn (Xt, Yt) seõ naèm gaàn ñöôøng thaúng . noùi caùch khaùc ESS seõ caøng nhoû vaø RSS caøng lôùn. Tyû soá
RSS TSS
ESS TSS
−= 1
ñöôïc goïi laø heä soá xaùc ñònh ña bieán vaø kyù hieäu laø R2. Thuaät ngöõ ña bieán khoâng aùp duïng trong hoài quy ñôn bieán bôûi vì chæ coù duy nhaát moät bieán phuï ñoäc laäp X. Tuy nhieân, do bieåu thöùc R2 trong hoài quy ñôn bieán cuõng gioáng nhö trong hoài quy ña bieán neân ôû ñaây chuùng ta duøng cuøng thuaät ngöõ
2
R
0
1
≤ R
2 ≤
ˆ u 1 −= 1 −= =
(3.26)
Y
ESS TSS
RSS TSS
2 t −
)
∑ ( Y t
∑
Ramu Ramanathan
22
Thuïc Ñoan/Haøo Thi
Roõ raøng raèng, R2 naèm giöõa khoaûng töø 0 ñeán 1. R2 khoâng coù thöù nguyeân vì caû töû soá vaø maãu soá ñeàu coù cuøng ñôn vò. Ñieåm quan saùt caøng gaàn ñöôøng thaúng öôùc löôïng, “ñoä thích hôïp” caøng cao, nghóa laø ESS caøng nhoû vaø R2 caøng lôùn. Do vaäy, R2 laø thoâng soá ño löôøng ñoä thích hôïp, R2 caøng cao caøng toát. ESS coøn
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
ñöôïc goïi laø bieán thieân khoâng giaûi thích ñöôïc bôûi vì
tuˆ laø aûnh höôûng cuûa nhöõng bieán khaùc ngoaøi Xt vaø khoâng coù trong moâ hình. RSS laø bieán thieân giaûi thích ñöôïc. Nhö vaäy, TSS, laø toång bieán thieân cuûa Y, coù theå phaân thaønh hai thaønh phaàn: (1) RSS, laø phaàn giaûi thích ñöôïc theo X; vaø (2) ESS, laø phaàn khoâng giaûi thích ñöôïc. Giaù trò R2 nhoû nghóa laø coù nhieàu söï bieán thieân ôû Y khoâng theå giaûi thích ñöôïc baèng X. Ta caàn phaûi theâm vaøo nhöõng bieán khaùc coù aûnh höôûng ñeán Y.
)
(ˆ t rY
ˆ tYY t
Ngoaøi yù nghóa laø moät tyû leä cuûa toång bieán thieân cuûa Y ñöôïc giaûi thích qua moâ hình, R2 coøn coù moät yù nghóa khaùc. Ñoù laø thoâng soá ño löôøng moái töông quan giöõa giaù trò quan saùt Yt vaø giaù trò döï baùo . Caàn xem laïi phaàn trình baøy
2
Cov
2
2
R
=
=
=
veà heä soá töông quan cuûa maãu vaø cuûa toång theå ôû Phaàn 2.3 vaø 3.5. Phaàn 3.A.9 trình baøy
(3.26a)
r ˆ YY
RSS TSS
Var
YY ( t t Var )
Y ( t
)ˆ )ˆ( Y t
Nhö vaäy, bình phöông heä soá töông quan ñôn bieán giöõa giaù trò quan saùt Yt vaø tYˆ baèng phöông trình hoài quy thì seõ cho ra keát quaû baèng vôùi giaù giaù trò döï baùo trò R2 ñöôïc ñònh nghóa trong Phöông trình (3.26a). Keát quaû naøy vaãn ñuùng trong tröôøng hôïp coù nhieàu bieán giaûi thích, mieãn laø trong hoài quy coù moät soá haïng haèng soá.
Coù moät thaéc maéc phoå bieán veà ñoä thích hôïp toång theå, ñoù laø “baèng caùch naøo ñeå xaùc ñònh raèng R2 laø cao hay thaáp?”. Khoâng coù moät quy ñònh chuaån hay nhanh choùng ñeå keát luaän veà R2 nhö theá naøo laø cao hay thaáp. Vôùi chuoãi döõ lieäu theo thôøi gian, keát quaû R2 thöôøng lôùn bôûi vì coù nhieàu bieán theo thôøi gian chòu aûnh höôûng xu höôùng vaø töông quan vôùi nhau raát nhieàu. Do ñoù, giaù trò quan saùt R2 thöôøng lôùn hôn 0.9. R2 beù hôn 0.6 vaø 0.7 ñöôïc xem laø thaáp. Tuy nhieân, ñoái vôùi döõ lieäu cheùo, ñaïi dieän cho daïng cuûa moät yeáu toá thay ñoåi vaøo moät thôøi ñieåm naøo ñoù, thì R2 thöôøng thaáp. Trong nhieàu tröôøng hôïp, R2 baèng 0.6 hoaëc 0.7 thì chöa haún laø xaáu. Ñaây ñôn giaûn chæ laø thoâng soá ño löôøng veà tính ñaày ñuû cuûa moâ hình. Ñieàu quan troïng hôn laø neân ñaùnh giaù moâ hình xem daáu cuûa heä soá hoài quy coù phuø hôïp vôùi caùc lyù thuyeát kinh teá, tröïc giaùc vaø kinh nghieäm cuûa ngöôøi nghieân cöùu hay khoâng.
Ví duï 3.3
Ramu Ramanathan
23
Thuïc Ñoan/Haøo Thi
Trong baøi taäp veà giaù nhaø, TSS, ESS vaø R2 coù caùc giaù trò sau (xem laïi keát quaû ôû Phaàn thöïc haønh maùy tính 3.1):
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
R2 = 0,82052
TSS = 101.815 ESS = 18.274
Nhö vaäy, 82,1% ñoä bieán thieân cuûa giaù nhaø trong maãu ñöôïc giaûi thích bôûi dieän tích söû duïng töông öùng. Trong chöông 4, seõ thaáy raèng theâm vaøo caùc bieán giaûi thích khaùc, nhö soá löôïng phoøng nguû vaø phoøng taém seõ caûi thieän ñoä thích hôïp cuûa moâ hình.
3.5 Kieåm Ñònh Giaû Thuyeát Thoáng Keâ
Nhö ñaõ ñeà luùc ñaàu, kieåm ñònh giaû thuyeát thoáng keâ laø moät trong nhöõng nhieäm vuï chính cuûa nhaø kinh teá löôïng. Trong moâ hình hoài quy (3.1), neáu β baèng 0, giaù trò döï baùo cuûa Y seõ ñoäc laäp vôùi X, nghóa laø X khoâng coù aûnh höôûng ñoái vôùi Y. Do ñoù, caàn coù giaû thuyeát β = 0, vaø ta kyø voïng raèng giaû thuyeát naøy seõ bò baùc boû. Heä soá töông quan (ρ) giöõa hai bieán X vaø Y ño löôøng ñoä töông öùng giöõa hai bieán. Öôùc löôïng maãu cuûa ρ ñöôïc cho trong Phöông trình (2.11). Neáu ρ = 0, caùc bieán khoâng coù töông quan nhau. Do ñoù cuõng caàn kieåm ñònh giaû thuyeát ρ = 0. Phaàn naøy chæ thaûo luaän phöông phaùp kieåm ñònh giaû thuyeát ñoái vôùi α vaø β. Kieåm ñònh giaû thuyeát ñoái vôùi p seõ ñöôïc trình baøy ôû phaàn sau. Caàn löu yù raèng, tröôùc khi tieáp tuïc phaàn tieáp theo, baïn neân xem laïi Phaàn 2.8 veà kieåm ñònh giaû thuyeát vaø Phaàn 2.7 veà caùc loaïi phaân phoái.
Kieåm ñònh giaû thuyeát bao goàm ba böôùc cô baûn sau: (1) thieát laäp hai giaû thuyeát traùi ngöôïc nhau (Giaû thuyeát khoâng vaø Giaû thuyeát ngöôïc laïi), (2) ñöa ra kieåm ñònh thoáng keâ vaø phaân phoái xaùc suaát cho giaû thuyeát khoâng, vaø (3) ñöa ra quy luaät ra quyeát ñònh ñeå baùc boû hay chaáp nhaän giaû thuyeát khoâng. Trong ví duï veà giaù nhaø, Giaû thuyeát khoâng laø Ho : β = 0. Bôûi vì chuùng ta kyø voïng raèng β seõ döông, Giaû thuyeát ngöôïc laïi laø H1: β ≠0. Ñeå thöïc hieän kieåm ñònh naøy, βˆ vaø sai soá chuaån öôùc löôïng s ñöôïc söû duïng ñeå ñöa ra thoáng keâ kieåm ñònh. Ñeå ñöa ra phaân phoái maãu cho α vaø β, maø ñieàu naøy aûnh höôûng giaùn tieáp ñeán caùc soá haïng sai soá ngaãu nhieân u1, u2, …un (xem Phöông trình 3.15), caàn boå sung moät giaû thuyeát veà phaân phoái cuûa ut. GIAÛ THIEÁT 3.8 (Tính Chuaån Taéc cuûa Sai Soá) Moïi giaù trò sai soá ut tuaân theo phaân phoái chuaån N(0, σ2) , nghóa laø maät ñoä coù ñieàu kieän cuûa Y theo X tuaân theo phaân phoái N(α + βX, σ2).
Ramu Ramanathan
24
Thuïc Ñoan/Haøo Thi
Nhö vaäy, caùc soá haïng sai soá u1, u2, …un ñöôïc giaû ñònh laø ñoäc laäp vaø coù phaân phoái chuaån gioáng nhau vôùi giaù trò trung bình baèng khoâng vaø phöông sai baèng σ2. Giaû thieát 3.8 laø giaû thieát caên baûn trong kieåm ñònh giaû thuyeát thoáng keâ. Baûng 3.2 seõ trình baøy toùm taét taát caû caùc giaû thieát ñaõ ñöôïc ñöa ra. Nhöõng soá haïng sai soá thoûa caùc Giaû thieát töø 3.2 ñeán 3.8 thì ñöôïc xem laø sai soá ngaãu nhieân hay sai soá do nhieãu traéng.
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
BAÛNG 3.2 Caùc Giaû Thieát cuûa Moâ Hình Hoài Quy Tuyeán Tính Ñôn Bieán 3.1 Moâ hình hoài quy laø ñöôøng thaúng vôùi aån soá laø caùc heä soá α vaø β; ñoù laø
Yt = α + βXt + ut, vôùi t = 1, 2, 3…, n.
3.2 Taát caû caùc giaù trò quan saùt X khoâng ñöôïc gioáng nhau; phaûi coù ít nhaát moät
giaù trò khaùc bieät.
0.
3.3 Sai soá ut laø bieán ngaãu nhieân vôùi trung bình baèng khoâng; nghóa laø, E(ut) =
2
3.4 Xt ñöôïc cho vaø khoâng ngaãu nhieân, ñieàu naøy ngaàm ñònh raèng khoâng töông
quan vôùi ut; nghóa laø Cov (Xt, ut) = E(Xtut) – E(Xt)E(ut)= 0. ) ( 2 σ=tuE
us).
3.5 ut coù phöông sai khoâng ñoåi vôùi moïi t; nghóa laø Var(ut) = 3.6 ut vaø us coù phaân phoái ñoäc laäp ñoái vôùi moïi t ≠ s, sao cho Cov(ut, us) = E(ut
3.7 Soá löôïng quan saùt (n) phaûi lôùn hôn soá löôïng heä soá hoài quy ñöôïc öôùc löôïng
(ôû ñaây n > 2).
3.8 ut tuaân theo phaân phoái chuaån ut ~ N(0, σ2), nghóa laø öùng vôùi giaù trò Xt cho
tröôùc, Yt ~ N(α + βXt, σ2).
Xaùc Ñònh Trò Thoáng Keâ Kieåm Ñònh Phaàn naøy chöùng minh raèng kieåm ñònh thoáng keâ
=
−
ˆ
0
t c
) βββ ( ˆ s
tuaân theo
phaân phoái Student t, theo giaû thuyeát khoâng, vôùi baäc töï do laø n – 2 (bôûi vì ta ñang öôùc löôïng hai tham soá α vaø β). Löu yù raèng Giaû thuyeát 3.7 raát caàn ñeå chaéc chaén raèng baäc töï do laø döông.
CHÖÙNG MINH
(Ñoäc giaû khoâng quan taâm ñeán nguoàn goác vaán ñeà, coù theå boû qua phaàn naøy).
2
2
Tröôùc heát caàn xem xeùt caùc tính chaát sau TÍNH CHAÁT 3.6
ˆ)2
n
σ
−
)
] 2 σσ
( ∑ 2 ˆ u t n–2.
coù phaân phoái chi-bình phöông vôùi baäc töï do a. αˆ vaø βˆ coù phaân phoái chuaån. [ b. ( =
2ˆσ .
c. αˆ vaø βˆ ñöôïc phaân phoái ñoäc laäp vôùi
Ramu Ramanathan
25
Thuïc Ñoan/Haøo Thi
Tính chaát 3.6a xuaát phaùt töø thöïc teá laø αˆ vaø βˆ laø nhöõng toå hôïp tuyeát tính cuûa ut vaø ut coù phaân phoái chuaån. Ñeå chöùng minh tính chaát b vaø c, neân tham
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
2
ˆ u
2 t
),
khaûo taøi lieäu Hogg vaø Graig (1978, trang 296-298). Taän duïng caùc keát qua ñoù ta ñöôïc
),
~ˆ α N
~
X
2 ( , ˆασα
, ( ˆβσβ
2
2 n −
2
∑ σ
2
~ˆ β N
2 ˆασ vaø
ˆβσ laø phöông sai cuûa αˆ vaø βˆ theo Phöông trình (3.18) vaø (3.19). Baèng caùch chuaån hoùa phaân phoái cuûa thoâng soá öôùc löôïng – nghóa laø tröø cho trung bình vaø chia cho ñoä leäch chuaån) – ta ñöôïc
2
) ˆ σ
( n
)
β
)
~
2 nX
2
−
~
N
),1,0(
~
N
),1,0(
2
2 − σ
( ˆ αα − ˆ ασ ˆ
ˆ
( ˆ β − ˆ βσ
trong ñoù
21
2
)
)
)
t
~
=
÷
=
=
nt
2
−
2
( ˆ ββ − s
( ˆ − ββ ˆ σ
ˆ σ σ
( ˆ ββσ − ˆ σσ
ˆ β
ˆ β
ˆ β
Trong phaàn 2.7, phaân phoái t ñöôïc ñònh nghóa laø tyû soá cuûa soá chuaån chuaån hoùa treân caên baäc hai cuûa moät chi-square ñoäc laäp vôùi noù. Thay vaøo cho β vaø aùp duïng phöông trình (3.18), (3.19) vaø (3.22), ta ñöôïc
ˆ
s
=
ˆ β
ˆ ˆ σσ = σ S
σσσ ˆ β = σ S
xx
xx
trong ñoù
βˆs laø sai soá chuaån öôùc löôïng cuûa βˆ theo Phöông trình (3.22). t ñöôïc trình baøy ôû treân laø trò thoáng keâ kieåm ñònh döïa treân quy luaät ra quyeát ñònh ñöôïc thieát laäp sau naøy. Kieåm ñònh naøy ñöôïc goïi laø kieåm ñònh t. Caùc böôùc kieåm ñònh thoáng keâ phaân ra trong hai tröôøng hôïp kieåm ñònh moät phía vaø kieåm ñònh hai phía ñöôïc trình baøy sau ñaây.
Quy Taéc Ra Quyeát Ñònh Kieåm ñònh t-test moät phía
BÖÔÙC 1 BÖÔÙC 2
H0: β = β0 Kieåm ñònh thoáng keâ laø
−
=
ˆ
0
H1: β ≠β0 t c
) βββ ( ˆ s Theo giaû thuyeát khoâng, kieåm ñònh thoáng keâ coù phaân phoái t vôùi
Ramu Ramanathan
26
Thuïc Ñoan/Haøo Thi
, ñöôïc tính döïa treân maãu.
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
BÖÔÙC 3
BÖÔÙC 4
baäc töï do laø n – 2. Neáu tc tính ñöôïc laø “lôùn”, ta coù theå nghi ngôø raèng β seõ khoâng baèng β0. Ñieàu naøy daãn ñeán böôùc tieáp theo. Trong baûng tra phaân phoái t ôû trang bìa tröôùc cuûa saùch, tra baäc töï do laø n – 2. Vaø choïn möùc yù nghóa (α) vaø xaùc ñònh ñieåm t*n–2(α) sao cho P(t > t*) = α. Baùc boû H0 neáu tc > t*. Neáu giaû thuyeát ngöôïc laïi β < β0 , tieâu chuaån kieåm ñònh ñeå baùc boû H0 laø neáu tc < –t*.
Kieåm ñònh treân ñöôïc minh hoïa baèng hình aûnh qua Hình 3.6 (kyù hieäu α ñöôïc söû duïng ñeå chæ möùc yù nghóa ñeå traùnh nhaàm laãn vôùi α chæ tung ñoä). Neáu tc rôi vaøo dieän tích in ñaäm trong hình veõ (ñöôïc goïi laø vuøng tôùi haïn) nghóa laø tc >t*. Trong tröôøng hôïp ñoù, giaû thuyeát khoâng seõ bò baùc boû vaø keát luaän ñöôïc raèng β lôùn hôn β0 raát nhieàu.
HÌNH 3.6 Kieåm Ñònh Moät Phía vôùi H0: β = β0
H1: β ≠β0
f(tn-2)
Dieän tích a
t*n-2(a)
0 Chaáp nhaän Ho
tn-2 Baùc boû Ho
Ví duï 3.4
tc =
ˆ s ββ ˆ
, laø kieåm ñònh thoáng keâ Trong ví duï veà giaù nhaø, ta coù β0 = 0. Do ñoù,
13875 ,0
ˆ =βs
01873 ,0 ñôn giaûn vaø laø tyû soá giöõa heä soá hoài quy öôùc löôïng treân sai soá chuaån. Tyû soá ˆ =β ñöôïc goïi laø trò thoáng keâ t. Caùc öôùc löôïng laø , vaø theo ví duï 3.2 ta . Do ñoù, trò thoáng keâ t ñöôïc tính seõ laø tc = 0,13875/0,01873 bieát
= 7.41. Baäc töï do baèng n – 2 = 14 – 2 = 12. Cho möùc yù nghóa laø 1%, nghóa laø α = 1%. Tra baûng phaân phoái t, ta ñöôïc t*n–2=2,681 . Do tc > t*, giaû thuyeát H0 bò baùc boû vaø keát luaän ñöôïc raèng β lôùn hôn zero moät caùch ñaùng keå vôùi möùc yù nghóa 1%. Löu yù raèng heä soá naøy vaãn coù yù nghóa trong tröôøng hôïp möùc yù nghóa chæ laø 0,05% bôûi vì t*12(0,0005) = 4,318.
Ramu Ramanathan
27
Thuïc Ñoan/Haøo Thi
Trò thoáng keâ t ñoái vôùi αˆ ñöôïc cho bôûi tc = 52,351/37,285 = 1.404 nhoû hôn t*12(0,0005) = 1.782. Do ñoù khoâng theå baùc boû H0 nhöng thay vaøo ñoù coù theå coù theå keát luaän raèng α khoâng lôùn hôn zero xeùt veà maët thoáng keâ vôùi möùc yù nghóa 5%. Caùc ñieåm αˆ khoâng nghóa ôû hai ñieåm sau. Thöù nhaát, X = 0 thì hoaøn
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
toaøn naêm ngoaøi khoaûng maãu vaø do ñoù öôùc löôïng Yˆ khi X = 0 khoâng ñaùng tin
caäy (xem theâm Phaàn 3.9). Thöù nhì, töø Hình 3.1 coù theå thaáy raèng ñaëc ñieåm hai bieán laø khoâng ñaày ñuû ñeå giaûi thích ñoä bieán thieân giaù cuûa caùc giaù trò quan saùt. Trong chöông 4 seõ cho thaáy αˆ bao haøm caû aûnh höôûng trung bình cuûa bieán bò boû soùt vaø tính phi tuyeán, khi X baèng 0. Caùc aûnh höôûng treân seõ laøm cho αˆ khoâng coù yù nghóa.
Moät Soá Löu YÙùù khi Söû Duïng Kieåm Ñònh t-Test
Maëc duø kieåm ñònh t-test raát höõu ích trong vieäc xaùc ñònh yù nghóa thoáng keâ cuûa caùc heä soá, tuy nhieân raát deã nhaàm laãn giöõa caùc yù nghóa cuûa kieåm ñònh. Ví duï, ôû Ví duï 3.4 kieåm ñònh t-test ñoái vôùi α khoâng theå baùc boû giaû thuyeát khoâng laø α = 0. Nhö vaäy coù phaûi kieåm ñònh naøy “chöùng minh” raèng α = 0 hay khoâng? Caâu traû lôøi laø khoâng. Coù theå chaéc chaén raèng, theo taäp döõ lieäu vaø moâ hình ñöôïc moâ taû, khoâng coù baèng chöùng naøo cho thaáy α > 0. Trong chöông 4, seõ ñeà caäp kieåm ñònh t-test cho nhieàu heä soá hoài quy. Neáu moät trong nhöõng heä soá naøy khoâng coù yù nghóa (nghóa laø, khoâng theå baùc boû giaû thuyeát raèng heä soá baèng 0), ñieàu ñoù khoâng coù nghóa laø bieán töông öùng khoâng coù aûnh höôûng gì ñeán bieán phuï thuoäc hoaëc bieán ñoù khoâng quan troïng. Vaán ñeà naøy seõ ñöôïc thaûo luaän ñaày ñuû trong chöông sau. Trong chöông 5 seõ thaáy raèng khi moâ hình thay ñoåi, möùc yù nghóa cuûa heä soá cuõng thay ñoåi. Do ñoù, caàn thöïc hieän kyõ caùc kieåm ñònh giaû thuyeát ñöa ra vaø khoâng neân voäi vaõ keát luaän maø khoâng xeùt ñeán moâ hình vaø nhöõng phaân tích theâm veà caùc kieåm ñònh chuaån ñoaùn caàn thieát ñeå ñöa ra moät keát luaän yù nghóa (oån ñònh theo ñaëc ñieåm moâ hình).
Phöông Phaùp p-value trong Kieåm Ñònh Giaû thuyeát
p-value = P(t>tc ) = P(sai laàm loaïi I)
Kieåm ñònh t-test coù theå ñöôïc thöïc hieän theo moät phöông phaùp khaùc töông ñöông. Tröôùc tieân tính xaùc suaát ñeå bieán ngaãu nhieân t lôùn hôn trò quan saùt tc, nghóa laø
Ramu Ramanathan
28
Thuïc Ñoan/Haøo Thi
Xaùc suaát naøy (ñöôïc goïi laø p-value) laø phaàn dieän tích beân phaûi tc trong phaân phoái t (xem Hình 3.7) vaø laø xaùc suaát sai laàm loaïi I – nghóa laø xaùc suaát loaïi boû giaû thuyeát H0. Xaùc suaát naøy caøng cao cho thaáy haäu quaû cuûa vieäc loaïi boû sai laàm giaû thuyeát ñuùng H0 caøng nghieâm troïng. p-value beù nghóa laø haäu quaû cuûa
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
vieäc loaïi boû giaû thuyeát ñuùng H0 laø khoâng nghieâm troïng (nghóa laø, xaùc suaát xaûy ra sai laàm loaïi I laø thaáp) vaø do ñoù coù theå yeân taâm khi baùc boû H0. Nhö vaäy, quy luaät ra quyeát ñònh laø khoâng baùc boû H0 neáu p -value quaù lôùn, ví duï: lôùn hôn 0,1, 0,2, 0,3. Noùi caùch khaùc, neáu p-value lôùn hôn möùc yù nghóa α, coù theå keát luaän raèng heä soá hoài quy khoâng lôùn hôn β0 ôû möùc yù nghóa α. Neáu p-value nhoû hôn α, giaû thuyeát H0 bò baùc boû vaø keát luaän ñöôïïc raèng β lôùn hôn β0 moät caùch ñaùng keå.
Ñeå thaáy ñöôïc söï töông ñöông cuûa hai phöông phaùp, löu yù raèng treân Hình 3.7 neáu xaùc suaát P(t>tc ) beù hôn möùc yù nghóa α, thì ñieåm töông öùng laø tc phaûi naèm beân phaûi ñieåm t*n-2(α). Nghóa laø tc rôi vaøo mieàn baùc boû. Töông töï, neáu xaùc suaát P(t>tc ) lôùn hôn möùc yù nghóa α, thì ñieåm töông öùng laø tc phaûi naèm beân traùi ñieåm t*n-2(α) vaø do ñoù rôi vaøo mieàn chaáp nhaän. Sau ñaây laø caùc böôùc boå sung trong phöông phaùp p-value nhö sau:
HÌNH 3.7 Kieåm Ñònh Giaû thuyeát theo Phöông Phaùp p-value
f(tn-2)
Baùc boû Ho neáu p- value< a
t*
0
tc
tn-2
BÖÔÙC 3a Tính xaùc suaát (kyù hieäu p-value) ñeå t lôùn hôn tc , nghóa laø tính
phaàn dieän tích beân phaûi giaù trò tc.
BÖÔÙC 4a Baùc boû H0 vaø keát luaän raèng heä soá coù yù nghóa neáu p-value beù hôn
möùc yù nghóa ñöôïc choïn.
Ramu Ramanathan
29
Thuïc Ñoan/Haøo Thi
Toùm laïi, β ñöôïc xem laø lôùn hôn β0 moät caùch ñaùng keå neáu trò thoáng keâ t lôùn hay p-value laø beù, möùc ñoä nhö theá naøo laø lôùn vaø beù seõ ñöôïc quyeát ñònh bôûi ngöôøi nghieân cöùu. Phöông phaùp phoå bieán trong kieåm ñònh giaû thuyeát laø xaùc ñònh giaù trò moác t*. Tuy nhieân theo höông phaùp tính p-value, laïi caàn tính toaùn phaàn dieän tích moät ñaàu öùng vôùi giaù trò tc cho tröôùc. Ngaøy caøng coù nhieàu phaàn meàm maùy tính tính toaùn saün p-value (chöông trình SHAZAM vaø ESL ñöôïc giôùi thieäu trong saùch naøy) vaø do ñoù phöông phaùp naøy deã öùng duïng deã daøng. Tuy nhieân, caàn caån thaän kieåm tra laïi giaù trò p-value laø duøng cho kieåm moät phía hay kieåm ñònh hai phía.
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Ví duï 3.4a
Ñeå aùp duïng phöông phaùp p-value cho ví duï veà giaù nhaø, ta tính xaùc suaát ñeå t lôùn hôn giaù trò quan saùt β = 7.41. Söû duïng ESL ñeå tính toaùn ta ñöôïc p < 0,0001 (tham khaûo phaàn keát quaû trong phaàn Thöïc haønh maùy tính 3.1). Ñieàu ñoù coù nghóa laø, neáu ta baùc boû giaû thuyeát khoâng, thì cô hoäi ñeå xaûy ra sai laàm loaïi I beù hôn 0,01%, vaø do ñoù hoaøn toaøn yeân taâm khi baùc boû Ho vaø keát luaän ñöôïc raèng β lôùn hôn 0. Ñoái vôùi tham soá α, p-value baèng 0,093, nghóa laø P(t>1,404) = 0,093. Neáu H0: α = 0 bò baùc boû, xaùc suaát ñeå xaûy ra sai laàm loaïi I laø 9,3%, lôùn hôn 5%. Do ñoù, khoâng theå baùc boû H0 ôû möùc yù nghóa 5%, nghóa laø ta coù cuøng keát luaän nhö trong phöông phaùp ñaàu, ñoù laø ôû möùc yù nghóa 5%, α khoâng lôùn hôn zero xeùt veà maët thoáng keâ. Nhö vaäy phöông phaùp p-value coù moät öu ñieåm laø, ta bieát ñöôïc chính xaùc möùc ñoä maø heä soá coù yù nghóa vaø coù theå ñaùnh giaù xem möùc yù nghóa naøy ñuû thaáp hay khoâng ñeå xem xeùt baùc boû H0. Cuoái cuøng, khoâng caàn lo laéng ñoái vôùi caùc giaù trò 0,01, 0,05 vaø 0,1.
Kieåm Ñònh t-test Hai Phía Bao goàm caùc böôùc sau: BÖÔÙC 1
H0: β = β0
H1: β ≠β0
BÖÔÙC 2
−
=
ˆ
0
t c
) βββ ( ˆ s Theo giaû thuyeát khoâng, kieåm ñònh thoáng keâ coù phaân phoái t laø tn-2.
Kieåm ñònh thoáng keâ laø , ñöôïc tính döïa treân maãu.
BÖÔÙC 3
Trong baûng tra phaân phoái t ôû trang bìa tröôùc cuûa saùch, tra baäc töï do laø n – 2 vaø choïn möùc yù nghóa (α) vaø xaùc ñònh ñieåm t*n–2(α) sao cho P(t>t*) = α/2 (phaân nöûa möùc yù nghóa).
BÖÔÙC 3a AÙp duïng phöông phaùp ρ- value, tính giaù trò p
ρ- value = P(t > tc hoaëc t < –tc ) = 2P(t > |tc |)
Ramu Ramanathan
30
Thuïc Ñoan/Haøo Thi
do phaân phoái t ñoái xöùng.
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
BÖÔÙC 4
Baùc boû H0 neáu |tc |> t* vaø keát luaän β khaùc vôùi β0 moät caùch ñaùng keå ôû möùc yù nghóa α.
BÖÔÙC 4a Baùc boû H0 neáu p-value < α, ôû möùc yù nghóa naøy.
Kieåm ñònh treân ñöôïc minh hoïa baèng hình aûnh qua Hình 3.8. Baäc töï do trong tröôøng hôïp naøy baèng n–2. Neáu trò thoáng keâ t (tc ) rôi vaøo vuøng dieän tích ñen, giaû thuyeát khoâng bò baùc boû vaø keát luaän ñöôïc raèng β khaùc vôùi β0. giaù trò t* = 2 ñöôïc söû duïng laø quy luaät ñeå ñaùnh giaù möùc yù nghóa cuûa trò thoáng keâ t ôû möùc 5% (kieåm ñònh hai phía). Bôûi vì t* gaàn baèng 2 vôùi baäc töï do laø 25.
HÌNH 3.8 Kieåm Ñònh Hai Phía vôùi H0: β = β0 H1: β ≠β0
f(tn-2)
Dieän tích a/2
Dieän tích a/2
-t* n-2(a/2)
t*n-2(a/2)
Baùc boû Ho
tn-2 Baùc boû Ho
0 Chaáp nhaän Ho
Ramu Ramanathan
31
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Ví duï 3.5
.1
404
t
ˆ =α
=
41.7 055.3 vaø . Tra baûng giaù trò t, ta coù )005.0(* 12
>
)404.1 ˆ tPα (2 =
Theo caùch tính naøy tc trong ví duï giaù nhaø coù giaù trò nhö caùch tính theo t-test, ˆ =β , ñieàu naøy coù nghóa laø dieän tích cuûa caû 2 phía töông öùng vôùi giaù trò 3.055 laø 0.01. Bôûi ñoái vôùi βˆ thì tc>t* do ñoù ta coù theå loaïi giaû thuyeát H0 vaø keát luaän ñöôïc raèng β khaùc t * vôùi ôû möùc yù nghóa 1%. Ñoái vôùi αˆ thì lôùn hôn giaù trò tc. Do = 179.2)025.0( 12 ñoù ta khoâng theå baùc boû giaû thuyeát H0 (löu yù raèng ta ñang duøng kieåm ñònh giaù trò α ôû möùc yù nghóa 5%). Töø böôùc 3a ta coù theå suy ra ñöôïc giaù trò p-value ñoái = 0.186 (löu yù giaù trò p-value töông öùng vôùi tc trong vôùi tröôøng hôïp kieåm ñònh 2 phía seõ gaáp 2 laàn giaù trò cuûa noù trong tröôøng hôïp kieåm ñònh 1 phía). Do sai laàm loaïi I coù giaù trò 18.6% laø khoâng theå chaáp nhaän ñöôïc neân ta khoâng theå baùc boû giaû thuyeát H0: α = 0. Ñieàu naøy coù nghóa laø α khoâng coù yù nghóa veà thoáng keâ trong khi β laïi coù.
BAØI TAÄP 3.4
Trong ví duï giaù nhaø, haõy kieåm ñònh giaû thuyeát H0: β = 0.1 vaø giaû thuyeát H1: β ≠ 0.1 laàn löôït ôû möùc yù nghóa 0.05 vaø 0.01.
BAØI TAÄP 3.5
Chöùng minh raèng neáu moät heä soá coù yù nghóa ôû möùc 1% thì heä soá naøy cuõng seõ coù yù nghóa ôû möùc cao hôn.
BAØI TAÄP 3.6
Haõy chöùng minh raèng neáu moät heä soá khoâng coù yù nghóa ôû möùc 10% thì heä soá naøy cuõng seõ khoâng coù yù nghóa ôû baát kyø möùc yù nghóa naøo thaáp hôn 10%.
Kieåm Ñònh
2σ
2σ khoâng phoå bieán 2σ goàm caùc
2σ 0
2σ =
Maëc duø thoáng keâ kieåm ñònh möùc yù nghóa phöông sai sai soá nhöng vaãn ñöôïc trình baøy ñaày ñuû trong phaàn naøy. Kieåm ñònh böôùc sau:
BÖÔÙC 1
H1:
H0:
)2 = n ( −
BÖÔÙC 2
2σ 2σ ≠ 0 2 ˆ σ 2 ˆ σ 0 Chi-square vôùi baäc töï do n-2. Neáu Q coù giaù trò “lôùn” ta coù theå 2σ nghi ngôø raèng 0
2σ khoâng baèng
Ramu Ramanathan
32
Thuïc Ñoan/Haøo Thi
Trò kieåm ñònh laø . Sau ñoù tra baûng phaân phoái Qc
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
BÖÔÙC 3
n-2(α) sao cho dieän tích beân phaûi baèng α.
BÖÔÙC 4
n-2(α).
Trong baûng tra phaân phoái Chi-square ôû trang bìa tröôùc cuûa saùch, tra giaù trò cuûa Q* Baùc boû H0 ôû möùc yù nghóa α neáu Qc> Q*
/(
S
)
Nguyeân nhaân toång quaùt laøm cho kieåm ñònh naøy khoâng phoå bieán laø do ngöôøi 2σ söû duïng trong kieåm ñònh khoâng coù thoâng tin sô caáp ban ñaàu veà giaù trò cuûa giaû thuyeát H0. Kieåm Ñònh Ñoä Thích Hôïp
2 r = xy
2 xy
SS xx
yy
Ta coù theå thöïc hieän kieåm ñònh ñoä thích hôïp. Goïi p laø heä soá töông quan toång theå giöõa X vaø Y ñöôïc ñònh nghóa ôû Phöông trình (2.7). Theo phöông trình (2.11), ta thaáy giaù trò öôùc löôïng p2 ñöôïc xaùc ñònh bôûi trong
2
)
ñoù Sxx vaø Sxy ñöôïc ñònh nghóa theo Phöông trình (3.8) vaø (3.9), vaø
2
2
Y t ) TSS S = − = =
(3.27)
yy
∑
∑
( ∑ n
Y t ( YY − t
xy baèng vôùi R2 (ñieàu naøy chæ ñuùng trong tröôøng hôïp hoài qui ñôn bieán maø thoâi). ÔÛ Phaàn kieåm ñònh giaû thuyeát 2.8 trình baøy phöông phaùp kieåm ñònh giaû thuyeát cho raèng X vaø Y khoâng coù moái töông quan. Kieåm ñònh naøy goïi laø kieåm ñònh F (F-test). Kieåm ñònh F- test goàm caùc böôùc sau:
H1: ρxy ≠ 0
ÔÛ Phaàn 3.A.10 ngöôøi ta ñaõ chöùng minh raèng r2
BÖÔÙC 1 BÖÔÙC 2
H0: ρxy = 0 Trò thoáng keâ kieåm ñònh laø Fc = R2(n – 2)/(1 – R2). Fc cuõng coù theå ñöôïc tính theo coâng thöùc sau Fc = RSS(n – 2)/ESS. Theo giaû thuyeát H0, trò thoáng keâ naøy tuaân theo phaân phoái F vôùi 1 baäc töï do ôû töû soá vaø n – 2 baäc töï do ôû maãu soá.
BÖÔÙC 3
Tra baûng F theo 1 baäc töï ôû töû soá vaø n – 2 baäc töï do ôû maãu soá tìm giaù trò F* 1, n – 2 (α) sao cho phaàn dieän tích veà phía phaûi cuûa F* laø α, möùc yù nghóa.
BÖÔÙC 4 Baùc boû giaû thuyeát H0 (taïi möùc yù nghóa α) neáu Fc > F*.
Ramu Ramanathan
33
Thuïc Ñoan/Haøo Thi
Neân löu yù raèng giaû thuyeát H0 ôû treân seõ khoâng hôïp leä khi coù nhieàu giaù trò X. Nhö seõ ñöôïc trình baøy ôû chöông 4, kieåm ñònh F vaãn ñöôïc söû duïng nhöng H0 seõ khaùc.
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Ví duï 3.6
Trong ví duï giaù nhaø, R2 = 0,82052. Fc = 0,82052(14 – 2)/(1 – 0,82052) = 54,86. Theo ví duï 3.5, ESS = 18.274, vaø RSS = TSS – ESS = 83.541. Vì vaäy Fc coøn coù theå ñöôïc tính theo coâng thöùc khaùc nhö ôû böôùc 2: Fc = 83.541 (14 – 2)/18.274 = 54,86. Baäc töï do cuûa töû soá laø 1, cuûa maãu soá laø 12. Vôùi möùc yù 1, 12(0.05) = 4,75. Vì Fc > F* chuùng ta nghóa α = 5%, tra baûng A.4b ta ñöôïc F* baùc boû (taïi möùc yù nghóa 5%) giaû thuyeát H0 cho raèng X vaø Y khoâng töông quan. Thöïc ra, vì Fc > F* 1, 12(0.01) (tra baûng A.4a), giaû thuyeát H0 cuõng bò baùc boû taïi möùc yù nghóa 1%. Nhö vaäy, maëc duø giaù trò R2 khaù nhoû hôn 1, noù cuõng khaùc 0 moät ñaùng keå.
Trình Baøy Caùc Keát Quaû Hoài Quy Caùc keát quaû cuûa phaân tích hoài quy ñöôïc trình baøy theo nhieàu caùch. Theo caùch thoâng thöôøng, ngöôøi ta seõ vieát phöông trình öôùc löôïng keøm vôùi caùc trò thoáng keâ t ôû döôùi moãi heä soá hoài quy nhö sau:
=
+
GIAÙ
,0
SQFT
351,52 (1,404)
13875 (7,41)
821.0
2 =R
023.39=σ
.f.d =
12
=
+
GIAÙ
SQFT
351,52 (37.29)
,0 13875 (0.019)
Neáu nhieàu moâ hình hoài quy ñöôïc öôùc löôïng, vieäc trình baøy keát quaû ôû daïng baûng nhö Baûng 4.2 seõ thuaän tieän hôn.
Vieäc taùch toång caùc bình phöông toaøn phaàn ra thaønh caùc thaønh phaàn thöôøng
ñöôïc toùm taét ôû daïng baûng Phaân Tích Phöông Sai (ANOVA) Baûng 3.3.
Moät caùch khaùc laø ñieàn caùc sai soá chuaån döôùi caùc heä soá hoài quy:
3.6 Thang Ño vaø Ñôn Vò Ño
Giaû söû chuùng ta ñaõ tính GIAÙ theo ñôn vò ñoàng ñoâla thay vì theo ngaøn ñoàng ñoâla. Coät GIAÙ ôû baûng 3.1 seõ chöùa caùc giaù trò nhö 199.900, 228.000, v.v. Nhöõng öôùc löôïng cuûa heä soá hoài quy, caùc sai soá chuaån cuûa chuùng, R2, v.v. seõ bò aûnh höôûng nhö theá naøo bôûi söï thay ñoåi ñôn vò naøy? Caâu hoûi naøy seõ ñöôïc khaûo saùt ôû ñaây vì GIAÙ vaø SQFT ñöôïc tính ôû caùc ñôn vò khaùc nhau. Ñaàu tieân chuùng ta chaïy laïi moâ hình.
Ramu Ramanathan
34
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
GIAÙ = α + βSQFT + u
Goïi GIAÙ* laø giaù tính theo ñoâ la thöôøng. Nhö vaäy GIAÙ* = 1.000 GIAÙ. Nhaân moïi soá haïng trong phöông trình vôùi 1.000 vaø thay GIAÙ* vaøo veá traùi. Chuùng ta coù
GIAÙ* = 1.000α + 1.000βSQFT + 1.000u = GIAÙ* = α* + β*SQFT + u*
Neáu chuùng ta aùp duïng phöông phaùp OLS cho phöông trình naøy vaø cöïc tieåu t)2, chuùng ta seõ tìm ñöôïc caùc giaù trò öôùc löôïng cuûa α* vaø β*. Deã daøng hoùa Σ (u* nhaän thaáy raèng caùc heä soá hoài quy môùi seõ baèng caùc heä soá cuõ nhaân vôùi 1,000. Nhö vaäy, thay ñoåi thang ño cuûa chæ bieán phuï thuoäc trong moâ hình hoài quy laøm cho thang ño cuûa moãi heä soá hoài quy thay ñoåi theo töông öùng. Vì u* = 1,000u, caùc phaàn dö vaø sai soá chuaån cuõng seõ ñöôïc nhaân leân 1.000. Toång caùc bình phöông seõ ñöôïc nhaân theâm 1 trieäu (1.000 bình phöông). Caàn löu yù raèng caùc trò thoáng keâ t, F, vaø R2 seõ khoâng bò aûnh höôûng vì chuùng laø caùc tæ soá trong ñoù yeáu toá thang ño seõ trieät tieâu.
BAÛNG 3.3 Phaân Tích Phöông Sai
Nguoàn
F
Toång bình phöông (SS)
Baäc töï do (d.f.)
Bình phöông trung bình (SS÷d.f.)
83.541
1
Hoài quy (RSS)
2)
= 83.541
=
) Y ˆ( Yt −∑ , 86 54 ( − RSS 2n ESS
N – 2 = 12 1.523
2)
Sai soá (ESS) ∑ 2ˆtu = 18.274 Toång (TSS) Y = 101.815 N – 1 = 13 7.832 Yt −∑ (
Taùc ñoäng cuûa vieäc thay ñoåi thang ño cuûa moät bieán ñoäc laäp seõ ra sao? Giaû söû SQFT ñöôïc tính theo ñôn vò traêm meùt vuoâng thay vì theo meùt vuoâng thoâng thöôøng, nhöng GIAÙ ñöôïc tính theo ñôn vò ngaøn ñoâla nhö tröôùc. Goïi SQFT’ laø bieán tính theo traêm meùt vuoâng. Vaäy SQFT= 100SQFT’. Thay vaøo phöông trình ban ñaàu ta coù:
GIAÙ = α + β100SQFT’ + u
100 ˆ β 'ˆ = β
Ramu Ramanathan
35
Thuïc Ñoan/Haøo Thi
Roõ raøng theo phöông trình naøy, neáu chuùng ta hoài quy GIAÙ theo moät haèng soá vaø SQFT’, heä soá duy nhaát seõ bò aûnh höôûng laø heä soá cuûa SQFT. Neáu β’ laø heä soá cuûa SQFT’, thì . Sai soá chuaån cuûa noù cuõng seõ nhaân vôùi 100. Tuy nhieân, taát caû caùc soá ño khaùc – ESS, giaù trò thoáng keâ t, F, R2 chaúng haïn seõ
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
khoâng bò aûnh höôûng. Toùm laïi, trong moät moâ hình hoài quy tuyeán tính, neáu thang ño cuûa moät bieán ñoäc laäp thay ñoåi caùc heä soá hoài quy cuûa noù vaø caùc sai soá chuaån töông öùng seõ thay ñoåi töông öùng nhöng caùc trò thoáng keâ khaùc seõ khoâng thay ñoåi. Coù lyù do chính ñaùng ñeå thay ñoåi thang ño cuûa caùc giaù trò sao cho caùc soá sau khi thay ñoåi seõ khoâng lôùn cuõng khoâng quaù nhoû vaø töông töï vôùi caùc giaù trò cuûa caùc bieán khaùc. Ñieàu naøy laø vì caùc soá coù giaù trò lôùn seõ laán aùt caùc sai soá vaø caùc soá nhoû seõ gaây ra sai soá laøm troøn, ñaëc bieät laø khi tính giaù trò toång bình phöông, vieäc naøy seõ laøm aûnh höôûng xaáu ñeán ñoä chính xaùc cuûa keát quaû.
Ñeå hieåu moät caùch thöïc teá haäu quaû cuûa vieäc thay ñoåi ñôn vò, haõy Thöïc Haønh
Maùy Tính phaàn 3.2 ôû phuï luïc D. BAØI TAÄP 3.7
Giaû söû chuùng ta ñaët moät bieán môùi X* = SQFT – 1.000 (nghóa laø, X* laø phaàn dieän tích vuoâng treân 1.000) vaø öôùc löôïng moâ hình GIAÙ = a + bX* + v. Giaûi thích baèng caùch naøo baïn coù theå tìm ñöôïc aˆ vaø bˆ töø αˆ vaø βˆ maø khoâng phaûi öôùc löôïng laïi moâ hình môùi.
3.7 ÖÙng duïng: Öôùc Löôïng Ñöôøng Engel Bieåu Dieãn Quan Heä Giöõa Chi Tieâu cho
Chaêm Soùc Söùc Khoûe vaø Thu Nhaäp.
Trong phaàn naøy, chuùng ta seõ trình baøy moät öùng duïng “taäp döôït” vôùi moâ hình hoài quy hai bieán. Döõ lieäu ñöôïc söû duïng laø chuoãi döõ lieäu cheùo cho 50 bang vaø quaän Columbia (n = 51), döõ lieäu ñöôïc thu thaäp töø cuoán Toùm Löôïc Thoáng Keâ Myõ naêm 1995 (Statistical Abstract of the US). Caùc giaù trò cuûa döõ lieäu thöïc coù ôû taäp tin DATA3-2. Caùc bieán laø: EXPHLTH = Chi tieâu toång hôïp (ñôn vò tyû ñoâla) cho chaêm soùc söùc khoûe cuûa bang vaøo naêm 1993, Baûng 153, trang 111, khoaûng töø 0,998- 9,029.
INCOME = Thu nhaäp caù nhaân (ñôn vò tyû ñoâla) cuûa bang vaøo naêm 1993,
Baûng 712, trang 460, khoaûng töø 9,3-64,1.
Ramu Ramanathan
36
Thuïc Ñoan/Haøo Thi
Moâ hình laø ñöôøng Engel tìm ñöôïc ôû ví duï 1.4 vaø ñöôïc aùp duïng vôùi toång chi tieâu cho chaêm soùc söùc khoûe cuûa Myõ laø haøm soá theo toång thu nhaäp caù nhaân. Phaàn ÖÙng Duïng Maùy Tính 3.3 (xem phuï luïc baûng D.1) trình baøy höôùng daãn ñeå tìm ra keát quaû. Baûn chuù thích cuûa baùo caùo in töø maùy tính, söû duïng chöông trình ESL vaø taäp tin PS3-3.ESL, ñöôïc trình baøy ôû baûng 3.4. Phaàn ñöôïc in ñaäm laø nhaäp löôïng cuûa chöông trình vaø caùc phaàn in nghieâng laø caùc nhaän xeùt veà keát quaû. Baïn neân tìm hieåu caùc chuù thích naøy caån thaän vaø söû duïng chöông trình hoài quy baïn coù ñeå chaïy laïi caùc keát quaû naøy (taäp tin PS3-3.SHZ chöùa caùc doøng leänh ñeå söû duïng phaàn meàm SHAZAM). Döôùi ñaây laø moâ hình öôùc löôïng cuøng vôùi trò thoáng keâ maãu t trong ngoaëc ñôn, vaø p-value (giaù trò xaùc suaát p) trong ngoaëc vuoâng:
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
=
+
EXPHLTH ,0 141652 INCOME
σˆ = 2,547
176496 ,0 (0.378) [0.707] d.f. = 49 (49.272) [<0.0001] F = 2.428 R2 = 0,98
Moâ hình raát thích hôïp vôùi soá lieäu vì 98% söï bieán ñoåi cuûa chi tieâu cho chaêm soùc söùc khoûe ñöôïc giaûi thích bôûi bieán thu nhaäp. Nhö ñaõ giaûi thích ôû Baûng 3.3, soá haïng haèng soá khoâng coù yù nghóa veà maët thoáng keâ vaø phuø hôïp vôùi tieâu chuaån lyù thuyeát ñeà ra trong ví duï 1.4, chæ ra raèng α = 0. Ñeå bieát theâm chi tieát, xem caùc chuù thích ôû Baûng 3.4
3.8 Khoaûng Tin Caäy
Nhö ñaõ ñöôïc chæ ra ôû Phaàn 2.9, moät caùch ñeå xem xeùt tröïc tieáp ñeán vieäc öôùc löôïng α vaø β trong ñieàu kieän khoâng chaéc chaén laø xaùc ñònh khoaûng tin caäy. Nhö vaäy, ví duï, thay vì noùi βˆ = 0,139 chuùng ta coù theå noùi raèng vôùi möùc xaùc
^ β seõ naèm trong khoaûng töø 0,09 ñeán 0,17. Töø keát quaû caùc giaù
suaát cho tröôùc, trò thoáng keâ kieåm ñònh ôû phaàn 3.5 ta coù:
~
nt
2
−
nt
2
−
ˆ ββ − sβ
n-2(0,025) laø ñieåm name treân phaân phoái t vôùi n-2 baäc töï do sao cho
~ vaø ˆ − αα sα
Ramu Ramanathan
37
Thuïc Ñoan/Haøo Thi
Ñaët t* P(t>t*) = 0,025. Ñieàu naøy töông ñöông vôùi P(- t*≤ t ≤ t*) = 0,95. Nhö vaäy,
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
(1) exphlth (2) income
BAÛNG 3.4 Baùo Caùo töø Maùy Tính Keøm Theo Chuù Giaûi cho Phaàn 3.7 Caùc leänh ESL ñöôïc in ñaäm vaø caùc nhaän xeùt ñöôïc in nghieâng. Danh saùch bieán (0) Haèng soá
(Ñoà thò cuûa möùc chi tieâu theo thu nhaäp cho thaáy coù söï quan heä chaët cheõ giöõa hai
bieán)
income
| o | | + | | | o | + | o | o | o | o o + o | o o | o | o ooo | oooo o + ooo | +−−−−−−−−−+−−−−−−−−−+−−−−−−−−−+−−−−−−−−−+−−−−−−−−−+−−−−−−−−−+ 9,3 683,5
exphlth 94,178 78,648 52,764 26,881 0,998
Öôùc löôïng OLS vôùi 51 quan saùt 1-51 Bieán phuï thuoäc EXPHLTH
Ramu Ramanathan
38
Thuïc Ñoan/Haøo Thi
Bieán (0) haèng (1) income Heä soá 0,176496 0,141652 Sai soá chuaån 0,467509 0,002875 T-stat 0,377525 49,271792 2 Prob(t >T) 0,707414 <0,0001***
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Giaù trò öôùc löôïng cuûa heä soá cuûa bieán thu nhaäp laø βˆ = 0,141652 vaø öôùc löôïng cuûa soá haïng haèng soá laø αˆ = 0,176496. Trò thoáng keâ t (heä soá chia cho sai soá chuaån) cuûa bieán thu nhaäp laø 49,271792, ñaây laø giaù trò raát yù nghóa. 2Prob(t >T) laø vuøng dieän tích ôû hai ñaàu phaân phoái t chaën bôûi giaù trò kieåm ñònh t vaø laø giaù trò p-value hoaëc xaùc suaát sai laàm loaïi I (ñoái vôùi kieåm ñònh 2 phía). Neáu p-value nhoû (trong tröôøng hôïp naøy, nhoû hôn 0,0001), chuùng ta “an toaøn” khi baùc boû giaû thuyeát Ho raèng β = 0, vaø keát luaän raèng heä soá cuûa bieán thu nhaäp laø khaùc 0 ñaùng keå. Giaù trò p-value cuûa soá haïng haèng soá baèng 0,707414 gôïi yù raèng neáu chuùng ta baùc boû giaû thuyeát Ho cho raèng α = 0, chuùng ta coù theå phaïm phaûi sai laàm loaïi I trong 70,7 % soá laàn. Vì möùc sai laàm naøy quaù cao, chuùng ta khoâng theå baùc boû giaû thuyeát Ho. Nhö vaäy chuùng ta keát luaän raèng soá haïng haèng soá khoâng khaùc 0 ñaùng keå. Löu yù raèng trong ví duï 1.4, vieäc suy dieãn lyù thuyeát ra ñöôøng Engel aùm chæ raèng khoâng coù soá haïng haèng soá. Soá haïng haèng soá khoâng coù yù nghóa laø phuø hôïp vôùi keát quaû theo lyù thuyeát. Xu höôùng chi tieâu caän bieân cho vieäc chaêm soùc söùc khoûe laáy töø thu nhaäp laø 0,141652; nghóa laø, vôùi moãi khoaûn taêng thu nhaäp 100 ñoâla, chuùng ta coù theå kyø voïng caùc caù nhaân seõ chi trung bình 14,17 ñoâla cho chaêm soùc söùc khoûe.
Giaù trò R2 (R-square) chæ ra raèng 98% söï bieán ñoåi cuûa chi tieâu ñöôïc giaûi thích bôûi bieán thu nhaäp. Söï khaùc nhau giöõa giaù trò R2 Hieäu chænh vaø Khoâng hieäu chænh seõ ñöôïc giaûi thích ôû chöông 4 cuøng vôùi caùc giaù trò thoáng keâ maãu ñeå choïn moâ hình.
Giaù trò thoáng keâ maãu Durbin-Watson vaø heä soá töông quan chuoãi baäc nhaát seõ ñöôïc giaûi thích ôû chöông 9, nhaèm giaûi quyeát söï vi phaïm giaû thieát 3.6 cho raèng caùc soá haïng sai soá cuûa hai quan saùt laø khoâng töông quan. Giaù trò trung bình cuûa bieán phuï thuoäc laø Y vaø S.D. laø ñoä leäch chuaån cuûa Sy
15,068863 S.D. cuûa bieán phuï thuoäc 17,926636
317,898611 2,547102
0,980 Sai soá chuaån cuûa phaàn dö R- hieäu chænh 0,980
Giaù trò trung bình cuûa bieán phuï thuoäc Toång bình phöông sai soá (ESS) R- bình phöôngkhoâng hieäu chænh Trò thoáng keâ F 2427,709468 <0,0001
p-value = Prob(F>2427.709) Heä soá töï töông quan baäc nhaát
2,209485 -0,121
Trò thoáng keâ Durbin- Watson
Giaù trò thoáng keâ maãu ñeå choïn moâ hình
Ramu Ramanathan
39
Thuïc Ñoan/Haøo Thi
SGMASQ 6,487727 AIC 6,741876 FPE 6,742147
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
SHIBATA 6,722193
HQ GCV 6,939901 6,752532 SCHWARZ RICE 7,272471 6,7638
(giaù trò “thích hôïp” = exphlth quan saùt tröø phaàn dö)
(In giaù trò chi tieâu thöïc vaø döï baùo, giaù trò phaàn dö.
(löu caùc öôùc löôïng phaàn dö trong maùy vaøo ut.)
?genr ut=uhat Generated var. no. 3 (ut) ?genr =exphlth-ut Generated var. no. 4 (yhat) ?print –o exphlth yhat ut; Daáu hieäu –o chæ in ra ôû daïng baûng)
exphlth yhat ut
Ramu Ramanathan
40
Thuïc Ñoan/Haøo Thi
Obs 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 0,998 1,499 4,285 1,573 2,021 2,26 1,953 2,103 3,428 2,277 3,452 3,485 3,433 3,747 4,4 3,878 5,197 4,118 6,111 6,903 6,187 7,341 7,999 8,041 12,216 10,066 9,029 10,384 10,635 1,493862 1,763001 2,598749 2,131297 1,720505 2,343775 1,989644 2,244618 3,179523 2,910384 3,731965 4,057766 3,476992 4,652705 4,666871 3,916114 4,341071 4,426062 5,672601 7,301601 5,686766 7,485749 8,533975 7,967367 13,250993 11,027054 8,84561 9,256401 10,276297 -0,49586172 -0,26400087 1,686251 -0,55829655 0,30049479 -0,08377483 -0,03664435 -0,1416183 0,24847729 -0,63338356 -0,27996523 -0,57276526 -0,0439923 -0,90570543 -0,26687065 -0,03811407 0,85592937 -0,30806194 0,43839884 -0,39850129 0,50023362 -0,14474913 -0,53497529 0,07353344 -1,034993 -0,96105374 0,1833899 1,127599 0,35870284
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
*
30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 12,06 13,014 14,194 15,154 14,502 16,203 15,949 15,129 16,401 23,421 6,682 20,104 18,241 25,741 27,136 33,456 34,747 41,521 44,811 49,816 67,033 94,178 10,318793 10,276297 13,619289 16,96228 14,32755 13,477637 14,68168 16,395672 15,701576 20,985202 20,036133 19,002072 18,56295 30,093438 27,756177 31,042507 37,516012 36,439456 40,320726 49,0465 64,004971 96,995765 1,741207 2,737703 0,57471128 -1,80828 0,17445035 2,725363 1,26732 -1,256672 0,69942416 2,435798 -13,354133 1,101928 -0,32194997 -4,352438 -0,62017675 2,413493 -2,769012 5,081544 4,490274 0,7694999 3,028029 -2,817765
* st
* st
* ≤−
≤
=
=
−
≤
α
α
t P tP ( ) 95.0 ˆ( α ˆ +≤ αα ˆ − αα s α
β
ˆ ± β st * vaø
Töø ñaây coù theå ruùt ra raèng khoaûng tin caäy 95% cuûa α vaø β laàn löôït laø αα st * ˆ ±
Ví duï 3.7
52,351 ± (2,179x37,285) = (-28,893; 133,595) 0,13875 ± (2,179x0,018373) = (0,099; 0,179)
Trong ví duï veà giaù nhaø, sai soá chuaån cuûa αˆ vaø βˆ laø αˆs = 37,285 vaø βˆs = 0,18373. Ñoàng thôøi, töø baûng t, ta coù t* 12(0,025) = 2,179. Do ñoù, khoaûng tin caäy 95% laø Ñoái vôùi α: Ñoái vôùi β:
Ramu Ramanathan
41
Thuïc Ñoan/Haøo Thi
Löu yù raèng caùc khoaûng tin caäy naøy laø töông ñoái roäng. Ñaây laø daáu hieäu cho thaáy moâ hình hoài quy tuyeán tính thích hôïp raát keùm vôùi taäp döõ lieäu. Moät moâ hình hoài quy thích hôïp seõ cho khoaûng tin caäy heïp hôn.
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
BAØI TAÄP 3.8
Xaùc ñònh khoaûng tin caäy cuûa α vaø β trong Phaàn ÖÙng Duïng 3.7
3.9 Döï Baùo
13875
351
.52
ˆ Y
.0
X
+
=
ˆ ˆ X βα+=
ˆ Y 0
0
Nhö ñaõ ñeà caäp tröôùc ñaây, moät trong caùc öùng duïng phoå bieán cuûa moâ hình hoài quy laø ñeå döï baùo (chuû ñeà naøy seõ ñöôïc thaûo luaän chi tieát hôn ôû chöông 11). Trong ví duï giaù nhaø, chuùng ta coù theå ñaët caâu hoûi giaù baùn döï baùo cuûa moät ngoâi nhaø coù dieän tích 2,000 meùt vuoâng seõ laø bao nhieâu. Moâ hình hoài quy öôùc löôïng . Nhö vaäy, khi X = 2,000, giaù trò döï baùo cuûaY laø laø 52,351 + (2,000x0,13875) = 329,851. Vì giaù ñöôïc tính theo ñôn vò ngaøn ñoâla, giaù trò döï baùo naøy cuõng coù ñôn vò ngaøn ñoâla. Vì vaäy, theo moâ hình, giaù trung bình öôùc löôïng cuûa moät caên hoä dieän tích 2,000 meùt vuoâng laø 329.851 ñoâla. Moät caùch toång quaùt, deã daøng nhaän thaáy neáu X coù giaù trò X0 thì giaù trò döï baùo cuûa Y0 seõ laø . Giaù trò trung bình coù ñieàu kieän cuûa bieán döï ñoaùn Y
cho tröôùc X = X0 laø
=
=
+
=
=
0
0
0
ˆ( XYE X ) E XYE ( X ) )ˆ( α )ˆ( X += βαβ EX 0
ˆY laø giaù trò döï baùo coù ñieàu kieän khoâng thieân leäch cuûa giaù baùn trung
Nhö vaäy
0 bình taïi X0. Khoaûng Tin Caäy cho Giaù Trò Döï Baùo Trung Bình Vì α vaø β ñöôïc öôùc löôïng coù sai soá, giaù trò döï baùo
0
ˆY cuõng chòu sai soá. Ñeå xeùt ñeán yeáu toá naøy, chuùng ta tính sai soá chuaån vaø khoaûng tin caäy cho giaù trò döï baùo trung bình. Döôùi ñaây laø öôùc löôïng cuûa phöông sai cuûa giaù trò döï baùo (xem chöùng minh ôû Phaàn 3.A.11 )
2
X
X
)
(
0
s
=
+
2 ˆ σ
(3.28)
2 Y 0
1 n
− S
xx
* st
]
−
+
ˆ[ Y 0
ˆ, * Yst ˆ 0 Y 0
ˆ Y 0
Ramu Ramanathan
42
Thuïc Ñoan/Haøo Thi
Khoaûng tin caäy cuûa giaù trò döï baùo trung bình laø
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
0
trong ñoù t* laø giaù trò ngöôõng cuûa phaân phoái t. Löu yù raèng khi X0 caøng leäch xa ˆYs caøng lôùn vaø khoaûng tin caäy töông öùng caøng roäng. giaù trò trung bình X , thì Ñieàu naøy coù nghóa raèng neáu döï baùo ñöôïc thöïc hieän quaù xa khoûi phaïm vi cuûa maãu, ñoä tin caäy cuûa döï baùo seõ giaûm ñi. Neáu X0 = X , khoaûng tin caäy seõ heïp nhaát. Hình 3.9 cho yù nieäm veà “daûi tin caäy” vôùi caùc giaù trò X0. HÌNH 3.9 Daûi Khoaûng Tin Caäy cuûa Caùc Giaù Trò Döï Baùo
Khoaûng Tin Caäy cho Döï Baùo Ñieåm Phöông sai maãu trình baøy ôû phaàn tröôùc duøng ñeå döï baùo giaù trò trung bình. Beân caïnh ñoù chuùng ta cuõng muoán tìm phöông sai cuûa sai soá döï baùo cho caùc giaù trò thöïc Y0 töông öùng vôùi X0. Coâng thöùc döôùi ñaây ñöôïc laáy töø Phuï luïc 3.A.12:
2
(
)
X
X
0
Var
=
=
2 s ˆ u
)ˆ( u 0
0
2 s ˆ Y 0
− S
1 ++ n
xx
2 1ˆ σ
>
ˆ u
Y
=
−
0
0
(3.29)
ˆ Y laø sai soá cuûa döï baùo ñieåm. Khoaûng tin caäy ñöôïc tính 0 ˆYs . Khi côõ maãu lôùn, soá haïng thöù hai vaø thöù ba ôû treân seõ
0
Ramu Ramanathan
43
Thuïc Ñoan/Haøo Thi
thay vì theo trong ñoù 0ˆus
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
σˆ2
0ˆus moät giaù trò gaàn baèng σˆ . Ngoaøi ra, t* cuõng gaàn baèng 2 trong tröôøng hôïp möùc yù nghóa 95%. Nhö vaäy, khoaûng tin caäy cuûa maãu coù kích thöôùc lôùn laø
ˆ 0 ±Y
khoâng ñaùng keå so vôùi
Ví duï 3.8
2
2 ˆ
ˆ0us
0Ys = 111,555 vaø
Trong ví duï giaù nhaø, chuùng ta coù =1634,353 vaø khoaûng
tin caäy töông öùng khi X0 = 2.000 seõ laø (307, 353) vaø (242, 418). Khoaûng tin caäy vôùi côõ maãu lôùn laø (252,408). (Xem phaàn Thöïc Haønh Maùy Tính 3.4 ñeå chaïy laïi keát quaû naøy).
t laø giaù trò döï baùo cuûa bieán phuï thuoäc cho quan saùt t, vaø Yt laø giaù trò
Chuùng ta neân choïn loaïi khoaûng tin caäy naøo trong soá hai loaïi treân? Vì quan taâm chính laø sai soá döï baùo ñoái vôùi giaù trò thöïc Y0, phöông trình (3.29) thöôøng ñöôïc söû duïng. Löu yù raèng khoaûng tin caäy theo phöông trình naøy roäng hôn nhieàu khoaûng tin caäy döïa treân phöông trình (3.28) So Saùnh Caùc Giaù Trò Döï Baùo Caùc nhaø phaân tích kinh teá vaø kinh doanh thöôøng söû duïng nhieàu hôn moät moâ hình ñeå döï baùo. Moät soá ño thöôøng duøng ñeå so saùnh naêng löïc döï baùo cuûa caùc moâ hình khaùc nhau laø sai soá bình phöông trung bình (hoaëc ñoâi khi ngöôøi ta söû duïng caên baäc hai cuûa noù, vaø ñöôïc goïi laø caên baäc hai sai soá bình phöông trung bình). Goïi Yf
2
thöïc. Sai soá bình phöông trung bình ñöôïc tính nhö sau:
∑
=
f − YY t t − 2n
( ) MSE RMSE = MSE
Neáu hai moâ hình ñöôïc söû duïng ñeå döï baùo Y, moâ hình naøo coù MSE nhoû hôn
seõ ñöôïc ñaùnh giaù laø moâ hình toát hôn cho muïc ñích döï baùo.
f − YY t
t
=
Moät soá ño höõu ích khaùc laø sai soá phaàn traêm tuyeät ñoái trung bình (MAPE)
∑
100 MAPE 1 n Y t
44
Ramu Ramanathan
Thuïc Ñoan/Haøo Thi
Soá ño naøy chæ coù yù nghóa neáu taát caû caùc giaù trò Y ñeàu döông (xem Phaàn ÖÙng Duïng 3.11). Moät caùch khaùc, chuùng ta coù theå tính sai soá phaàn traêm bình phöông trung bình (MSPE) hoaëc caên cuûa noù
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
2
t
=
RMSPE =
MSPE
MSPE
∑
1 n
f − YY t Y t
100
Moät phöông phaùp khaùc ñeå ñaùnh giaù moâ hình vaø naêng löïc döï baùo cuûa noù laø thöïc hieän döï baùo haäu maãu. Theo phöông phaùp naøy, ngöôøi phaân tích seõ khoâng söû duïng moät soá quan saùt cuoái cuøng (chaúng haïn, 10% soá quan saùt cuoái cuøng) trong vieäc öôùc löôïng moâ hình, nhöng seõ söû duïng caùc öôùc löôïng thoâng soá töø taäp quan saùt ñaàu tieân ñeå döï baùo Yt cho phaàn maãu ñeå daønh. Sau ñoù chuùng ta coù theå tính MSE vaø MAPE cho giai ñoaïn haäu maãu. Moâ hình naøo coù caùc giaù trò ño löôøng naøy thaáp hôn seõ toát hôn cho muïc ñích döï baùo.
3.10 Tính Nhaân Quaû trong Moâ Hình Hoài Quy
Ramu Ramanathan
45
Thuïc Ñoan/Haøo Thi
Khi ñònh moâ hình ôû daïng Y = α + βX + u, chuùng ta ngaàm giaû ñònh raèng X gaây ra Y. Maëc duø R2 ño ñoä thích hôïp, noù khoâng theå ñöôïc söû duïng ñeå xaùc ñònh tính nhaân quaû. Noùi caùch khaùc, vieäc X vaø Y töông quan chaët vôùi nhau khoâng coù nghóa raèng söï thay ñoåi X daãn ñeán söï thay ñoåi Y hay ngöôïc laïi. Ví duï, heä soá töông quan giöõa soá löôïng kaênguru cuûa UÙc vaø toång daân soá nöôùc naøy coù theå laø raát cao. Phaûi chaêng ñieàu naøy coù nghóa raèng söï thay ñoåi moät bieán seõ laøm cho bieán kia thay ñoåi? Roõ raøng laø khoâng, vì ôû ñaây chuùng ta coù moät tröôøng hôïp töông quan giaû taïo. Neáu chuùng ta hoài quy moät trong caùc bieán vôùi bieán coøn laïi, chuùng ta seõ coù söï hoài qui giaû taïo. Laáy moät ví duï khaùc thöïc teá hôn, giaû söû chuùng ta hoài quy soá löôïng vuï troäm trong moät thaønh phoá vôùi soá haïng haèng soá vaø soá nhaân vieân caûnh saùt (X) vaø sau ñoù quan saùt thaáy heä soá goùc öôùc löôïng coù giaù trò döông, coù nghóa raèng coù töông quan thuaän giöõa X vaø Y. Phaûi chaêng ñieàu naøy coù nghóa raèng vieäc taêng soá löôïng caûnh saùt seõ laøm taêng soá vuï troäm, do ñoù ngaàm keùo theo phaûi coù chính saùch giaûm löïc löôïng caûnh saùt? Roõ raøng keát luaän naøy laø khoâng theå chaáp nhaän ñöôïc. Ñieàu xaûy ra coù theå laø moái quan heä nhaân quaû laø ngöôïc laïi, coù nghóa laø thaønh phoá neân thueâ theâm caûnh saùt vì soá vuï troäm taêng leân, vaø nhö vaäy vieäc hoài quy X theo Y laø hôïp lyù hôn. Tuy nhieân, trong thöïc teá, hai bieán seõ ñöôïc xaùc ñònh keát hôïp vaø do ñoù chuùng ta neân ñònh roõ hai phöông trình, moät vôùi Y theo X vaø caùc bieán khaùc vaø phöông trình coøn laïi vôùi X theo Y vaø caùc bieán khaùc. Vieäc xaùc ñònh ñoàng thôøi caùc bieán seõ ñöôïc trình baøy chi tieát ôû chöông 13. Nhö seõ thaáy ôû chöông naøy caùc öôùc löôïng thu ñöôïc baèng caùch boû qua tính ñoàng thôøi seõ bò sai leäch vaø khoâng nhaát quaùn. Cuõng coù theå laø söï töông quan cao quan saùt ñöôïc giöõa X vaø Y coù theå hoaøn toaøn laø do caùc bieán khaùc vaø khoâng bieán naøo trong soá chuùng coù theå tröïc tieáp gaây ra caùc bieán coøn laïi. Nhöõng ví duï naøy nhaán maïnh taàm quan troïng cuûa vieäc caân nhaéc kyõ löôõng baûn chaát cô cheá haønh vi tieàm aån laø gì, töùc laø, quaù trình phaùt döõ lieäu laø gì (DGP), vaø laäp moâ hình moät caùch phuø hôïp. Lyù thuyeát kinh teá, kieán thöùc cuûa nhaø phaân tích veà caùc haønh vò tieàm aån, kinh nghieäm quaù khöù,
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
v.v. phaûi gôïi yù moâ hình neân phaûi ñöôïc xaùc ñònh nhö theá naøo. Tuy nhieân, coù theå kieåm ñònh phöông höôùng cuûa söï nhaân quaû moät caùch roõ raøng (chi tieát seõ trình baøy ôû chöông 10). Ñoäc giaû quan taâm ñeán vaán ñeà naøy coù theå tham khaûo baøi vieát cuûa Granger (1969) vaø Sims (1972).
Ñeå minh hoïa taàm quan troïng cuûa vieäc xaùc ñònh chính xaùc söï nhaân quaû, giaû
söû chuùng ta ñaûo ngöôïc vò trí cuûa X vaø Y vaø öôùc löôïng moâ hình:
(3.1’)
Xt = α* + β*Yt + vt
t
'
X
X
−
' + βα
+
Y t
t
' v t
t
* α * β
1 * β
v * β
−=
+
=
Lieäu chuùng ta coù theå tìm ñöôïc ñöôøng thaúng gioáng nhö tröôùc khoâng? Caâu traû lôøi, noùi chung, laø khoâng. Vì thuû tuïc bình phöông nhoû nhaát ñöôïc aùp duïng cho phöông trình (3.1) seõ cöïc tieåu hoùa toång bình phöông cuûa caùc ñoä leäch ñöùng töø ñöôøng thaúng (xem hình 3.10). Traùi laïi, ñöôøng thaúng nghòch cöïc tieåu hoùa toång bình phöông cuûa caùc ñoä leäch ngang vt. Tìm Yt theo Xt, Phöông trình (3.1’) coù theå ñöôïc vieát laïi nhö sau:
Ramu Ramanathan
46
Thuïc Ñoan/Haøo Thi
Vieäc cöïc tieåu hoùa ∑ 2ˆtu , laøm töông töï nhö vôùi phöông trình (3.1), vaø cöïc tieåu hoùa ∑ 2ˆtv seõ thöôøng cho ra caùc keát quaû khaùc nhau. Cuï theå hôn, giaù trò öôùc löôïng cuûa β’ seõ khaùc vôùi giaù trò β töø phöông trình (3.1). HÌNH 3.10 Cöïc Tieåu Hoùa Toång Bình Phöông theo Truïc Tung vaø Truïc Hoaønh
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Ví duï 3.9
=
+
Quan heä öôùc löôïng khi ∑ 2ˆtu ñöôïc cöïc tieåu hoùa laø (xem Phaàn Thöïc Haønh Maùy Tính 3.5) GIAÙ , 351 52 , 13875 0 SQFT
=
+
SQFT
, 385 33
, 913666 5
GIAÙ
Khi quan heä nhaân quaû ñöôïc ñaûo ngöôïc vaø ∑ 2ˆtv ñöôïc cöïc tieåu hoùa, chuùng ta coù
−=
+
−=
+
Nghòch ñaûo quan heä öôùc löôïng thöù hai vaø bieåu dieãn GIAÙ laø haøm cuûa SQFT, ta coù
GIAÙ SQFT 645.5 169.0 SQFT 385.33 913666 .5 1 913666 .5
Löu yù raèng daáu cuûa soá haïng haèng soá bò nghòch ñaûo vaø ñoä doác laø hoaøn toaøn khaùc.
S
xy
*
=
=
ˆ β
'
S
1 ˆ β
yy
Nhö vaäy döôùi ñieàu kieän gì thì hai ñöôøng öôùc löôïng seõ nhö nhau? Ñeå traû lôøi caâu hoûi naøy, ñaàu tieân aùp duïng OLS cho phöông trình (3.1’); nghóa laø cöïc tieåu hoùa∑ 2ˆtv . Hoaùn ñoåi X vaø Y trong phöông trình 3.10, ta coù:
yy S
xy
S / ˆ ' =β . Öôùc löôïng bình phöông nhoû nhaát laøm cöïc tieåu
'β baèng β , ñieàu kieän laø
xx
ˆ =β . Ñeå S / xy S
Vaø do ñoù ∑ 2ˆtu laø
xy
yy
=
=
2 xy SS xx
yy
xx
xy
Nhöng veá traùi cuûa phöông trình thöù hai laø rxy
2, bình phöông cuûa heä soá hoài quy ñôn giöõa X vaø Y (ñònh nghóa ôû phöông trình 2.11). Nhö vaäy, ñieàu kieän caàn laø X vaø Y phaûi töông quan hoaøn haûo. Tính chaát 2.4d noùi raèng neáu toàn taïi söï töông quan hoaøn haûo giöõa hai bieán, thì phaûi toàn taïi moät quan heä tuyeán tính
Ramu Ramanathan
47
Thuïc Ñoan/Haøo Thi
S S S 1 hoaëc S S
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
chính xaùc giöõa chuùng. Vì vaäy, söï thích hôïp giöõa X vaø Y phaûi hoaøn haûo thì chuùng ta môùi nhaän ñöôïc cuøng moät ñöôøng hoài quy cho duø chuùng ta aùp duïng OLS cho phöông trình (3.1) hay (3.1’). Nhìn chung, söï töông quan giöõa X vaø Y seõ khoâng hoaøn haûo, chính vì vaäy chuùng ta seõ khoâng nhaän ñöôïc cuøng moät ñöôøng thaúng hoài quy. Ñieàu naøy nhaán maïnh taàm quan troïng cuûa vieäc xaùc ñònh ñuùng höôùng quan heä nhaân quaû thay vì vieäc choïn thieáu suy xeùt bieán X vaø Y.
Nhö ñaõ ñöôïc minh hoïa tröôùc ñaây trong ví duï veà toäi phaïm, quan heä nhaân quaû coù theå theo caû hai chieàu, tình huoáng naøy ñöôïc goïi laø phaûn hoài. Quan heä giöõa giaù baùn vaø löôïng baùn cuõng laø ví duï cuûa hieän töôïng naøy. Vì giaù vaø löôïng baùn ñöôïc xaùc ñònh cuøng luùc bôûi quan heä töông taùc giöõa cung vaø caàu, cho neân caùi naøy coù theå aûnh höôûng caùi kia. Töông töï, hieän töôïng phaûn hoài cuõng ñöôïc tìm thaáy trong quan heä giöõa thu nhaäp toång hôïp vaø tieâu duøng hay ñaàu tö. Nhöõng tình huoáng naøy seõ ñöôïc trình baøy ôû chuû ñeà moâ hình hoài quy heä phöông trình ôû chöông 13.
3.11 Öùng Duïng: Quan Heä giöõa Baèng Saùng Cheá vaø Chi Phí cho Hoaït Ñoäng
Nghieân Cöùu vaø Phaùt Trieån (R&D) Phaàn naøy seõ trình baøy moät ví duï “dieãn taäp” khasùc veà phaân tích hoài quy. Döõ lieäu duøng trong ví duï naøy ôû taäp tin DATA3.3, maø seõ ñeà caäp ñeán caùc bieán sau: PATENTS = Soá öùng duïng baèng saùng cheá ñöôïc ghi nhaän, ñôn vò ngaøn, giao
ñoäng töø 84,5 - 189,4
R&D =
Chi phí cho nghieân cöùu vaø phaùt trieån, ñôn vò tæ ñoâla 1992, ñöôïc xaùc ñònh baèng tæ soá giöõa chi phí theo ñoâla hieän haønh vaø chæ soá giaûm phaùt toång saûn phaåm quoác noäi goäp (GDP), giao ñoäng töø 57,94 ñeán 166,7.
Döõ lieäu theo naêm laáy trong voøng 34 naêm töø 1960 ñeán 1993 cho toaøn boä
nöôùc Myõ. Nguoàn ñöôïc trình baøy ôû phuï luïc D.
Neáu moät quoác gia chi nhieàu hôn cho hoaït ñoäng nghieân cöùu vaø phaùt trieån, chuùng ta coù theå kyø voïng raèng quoác gia naøy seõ ñaït ñöôïc nhieàu caûi tieán ñöôïc baûo veä thoâng qua luaät baèng saùng cheá hôn. Do ñoù, chuùng ta kyø voïng toàn taïi moät quan heä döông giöõa soá löôïng baèng saùng cheá ñöôïc ban boá vaø chi tieâu cho R&D. Maëc duø hieäu quaû cuûa hoaït ñoäng nghieân cöùu vaø phaùt trieån seõ treã vaøi naêm sau khi döï aùn ñöôïc baét ñaàu, ñeå ñôn giaûn hoùa chuùng ta boû qua hieän töôïng naøy. ÔÛ nhöõng chöông sau chuùng ta seõ khaûo saùt hieäu öùng treã cuûa caùc bieán ñoäc laäp vaø seõ quay laïi ví duï naøy.
Ramu Ramanathan
48
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Moâ hình hoài quy tuyeán tính öôùc löôïng ñöôïc trình baøy döôùi ñaây keøm vôùi caùc trò thoáng keâ maãu t trong ngoaëc ñôn (Phaàn Thöïc Haønh Maùy Tính 3.6 höôùng daãn caùch laäp laïi keát quaû cuûa phaàn naøy vaø Baûng 3.5 trình baøy keát quaû.)
&
=
+
SAÙNGCHEÁ
, 0 DR792 (13,97)
d.f. = 32 ˆ =σ
, 34 571 (5,44) R2 = 0,859 Fc (1,32) =195,055
Ñeå kieåm ñònh moâ hình veà söï yù nghóa toång theå, chuùng ta söû duïng trò thoáng keâ F, coù giaù trò baèng 195,055. Theo giaû thuyeát H0 thì soá baèng saùng cheá vaø chi phí cho R&D laø khoâng töông quan, Fc tuaân theo phaân phoái F vôùi baäc töï do ôû töû soá laø 1 vaø baäc töï do ôû maãu soá laø 32 (= 34 – 2). Töø baûng A.4a (cuõng ôû trong bìa sau) chuùng ta coù nhaän xeùt raèng giaù trò ngöôõng F(1,32) ôû möùc yù nghóa 1% naèm giöõa 7,31 vaø 7,56. Vì Fc cao hôn nhieàu so vôùi giaù trò naøy, chuùng ta keát luaän raèng soá baèng saùng cheá vaø chi phí cho R&D laø töông quan ñaùng keå. Keát luaän naøy ñöôïc cuõng coá theâm thoâng qua giaù trò thoáng keâ maãu t. Kieåm ñònh hai ñaàu vôùi möùc yù nghóa 1%, baûng t trong bìa tröôùc cuûa quyeån saùch (hay Baûng A.2) cho thaáy giaù trò ngöôõng vôùi 32 baäc töï do naèm giöõa 2,704 vaø 2,75. Vì giaù trò quan saùt tc cao hôn nhöõng giaù trò naøy nhieàu chuùng ta keát luaän raèng caû soá haïng tung ñoä goác vaø ñoä doác coù giaù trò khaùc 0 moät caùch ñaùng keå. Soá ño ñoä thích hôïp R2 cho bieát moâ hình giaûi thích ñöôïc 85,9% söï bieán ñoåi cuûa bieán phuï thuoäc. Maëc duø ñaây döôøng nhö laø moät söï thích hôïp toát, tuy nhieân chuùng ta thaáy töø hình 3.11 raèng moâ hình khoâng hoaøn toaøn theå hieän söï bieán ñoåi thöïc teá veà soá baèng saùch cheá. Ñöôøng thaúng hoài quy laø ñöôøng lieàn vaø noù khoâng ñaïi dieän ñaày ñuû baûn chaát ñöôøng cong cuûa döõ lieäu quan saùt. Chính vì ñieàu naøy moâ hình seõ döï baùo raát keùm soá löôïng baèng saùng cheá taïi nhieàu naêm.
=
t
t
Ñieåm naøy ñöôïc neâu ra roõ hôn ôû Baûng 3.5, baûng naøy coù nhieàu trò thoáng keâ )tYˆ , coät naêm laø giaù trò höõu ích khaùc. Coät thöù tö laø giaù trò trung bình öôùc löôïng( phaàn dö ñöôïc tính baèng giaù trò quan saùt tröø ñi giaù trò trung bình öôùc löôïng ( ˆ u vaø coät cuoái cuøng laø sai soá phaàn traêm tuyeät ñoái (APE), ñöôïc tính baèng 100
. Giaù trò döï baùo trình baøy ôû baûng 3.5 ñöôïc laøm troøn ñeán 1 chöõ
)t ˆ YY − t Yu /ˆ
t
172,11
soá thaäp phaân. Vì döõ lieäu goác veà soá baèng saùng cheá chæ coù moät soá thaäp phaân, neân vieäc coá gaéng coù ñöôïc caùc giaù trò döï baùo coù ñoä chính xaùc ñeán hôn moät soá thaäp phaân laø khoâng coù yù nghóa. HÌNH 3.11 Soá Baèng Saùng Cheá Theo Chi Phí cho R&D cuûa Nöôùc Myõ
Soá baèng saùng cheá (ngaøn)
Ramu Ramanathan
49
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Chí phí R&D (tyû)
BAÛNG 3.5 Baùo Caùo Maùy Tính coù Chuù Thích cho Phaàn Öùng Duïng cuûa
Phaàn 3.11.
(2) R&D
(1) Naêm
(3) PATENTS (SAÙNG CHEÁ)
Caùc leänh ELS ñöôïc in ñaäm vaø caùc lôøi nhaän xeùt ñöôïc in nghieâng Danh saùch caùc bieán (0) Haèng soá Thôøi ñoaïn: 1, quan saùt lôùn nhaát: 34, phaïm vi quan saùt: suoát 1960-1993, hieän haønh 1960-1993 ?ols PATENTS const R&D; (Öôùc löôïng moâ hình theo OLS)
Öôùc löôïng theo OLS söû duïng 34 quan saùt töø 1960-1993 Bieán phuï thuoäc – PATENTS
T stat 5,437521 13,966211
Sai soá chuaån 6,357873 0,056704
Heä soá 34,571064 0,791935
2Prob(t > |T|) < 0,0001*** < 0,0001***
Bieán (0) Haèng (2) R&D
119,238235
29,305827
3994,300257
11,172371
0,859
0,855
Giaù trò trung bình cuûa bieán phuï thuoäc Toång bình phöông sai soá (ESS) R-bình phöông khoâng hieäu chænh Trò thoáng keâ F
195,055061
<0,0001
Trò Durbin-Watson
0,233951
0,945
S.D. cuûa bieán phuï thuoäc Sai soá chuaån cuûa phaàn dö R- bình phöông hieäu chænh p-value = Prob(F>2427.709) Heä soá töï töông quan baäc nhaát
Ramu Ramanathan
50
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Caùc giaù trò thoáng keâ ñeå choïn moâ hình
FPE SHIBATA
132,164347 131,300527
132,146377 144,560215 133,143342
AIC SCHWARZ RICE
124,821883 136,255226 132,623251
SGMASQ HQ GCV
(löu caùc öôùc löôïng phaàn dö )
(tính giaù trò “gaén”)
(laøm troøn ñeán moät soá thaäp
(in caùc giaù trò ôû daïng
?genr ut=uhat ?genr temp = PATENTS -ut genr fitted = int (0.5+ (10*temp))/10 phaân) ?genr error = PATENTS – fitted (tính sai soá döï baùo) ?genr abspcerr = int (0.5 + (1000*abs(error)/PATENTS))/100 (tính sai soá % tuyeät ñoái vaø laøm troøn ñeán hai chöõ soá thaäp phaân) ?print –o R&D PATENTS fitted error abspcerr; baûng)
R&D
Patens
OBS
Fited
1960 1961 1962 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992
57,94 60,59 64,44 70,66 76,83 80,00 84,82 86,84 88,81 88,28 85,29 83,18 85,07 86,72 85,45 83,41 87,44 90,11 94,50 99,28 103,64 108,77 113,96 121,72 133,33 144,78 148,39 150,90 154,36 157,19 161,86 164,54 166,70
84,5 S6,2 90,4 91,1 93,2 100,4 93,5 93,0 98,7 104,4 109,4 111,1 105,3 109,6 107,4 108,0 110,0 109,0 109,3 108,9 113,0 114,5 118,4 112,4 120,6 127,1 133,0 139,8 151,9 166,3 176,7 178,4 187,2
80,5 82,6 85,6 90,5 95,4 97,9 101,7 103,3 104,9 104,5 102,1 100,4 101,9 103,2 102,2 100,6 103,8 105,9 109,4 113,2 116,6 120,7 124,8 131,0 140,2 149,2 152,1 154,1 156,8 159,1 162,8 164,9 166,6
Error ABSPCER R 4,73 6,35 5,31 0,66 2,36 2,49 8,77 11,08 6,28 0,10 6,67 9,63 3,23 5,84 4,84 6,85 5,64 2,84 0,09 3,95 3,19 5,41 5,41 16,55 -16,25 17,39 14,36 10,23 3,23 4,33 7,87 7,57 11,00
4,0 5,6 4,8 0,6 -2,2 2,5 -8,2 -10,3 -6,2 -0,1 7,3 10,7 3,4 6,4 5,2 7,4 6,2 3,1 -0,1 -4,3 -3,5 -6,2 -6,4 -18,5 -19,6 -22,1 -19,1 -14,3 -4,9 7,2 13,9 13,5 20,6
Ramu Ramanathan
51
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
1993
165,20
189,4
155,4
24,0
12,67
Nhieàu giaù trò APE lôùn hôn 5%, vaø trong moät soá naêm chuùng vöôït qua 10%, ñaây laø tæ leä khaù lôùn. Chuùng ta cuõng quan saùt thaáy raèng caùc ñieåm phaân taùn co cuïm laïi vôùi nhau trong caùc naêm töø 1966-1977, chæ ra raèng moät yeáu toá naøo ñoù khaùc hôn laø chi phí R&D gaây ra söï thay ñoåi veà soá baèng saùng cheá. Do ñoù, quan saùt kyõ hôn caùc keát quaû chæ cho thaáy söï xaùc ñònh sai moâ hình. Trong chöông 6, chuùng ta seõ duøng taäp döõ lieäu naøy ñeå öôùc löôïng moâ hình ñöôøng cong vaø seõ xem xeùt xem lieäu vieäc xaùc ñònh naøy coù theå hieän toát hôn caùc bieán ñoåi quan saùt ñöôïc veà soá baèng saùng cheá khoâng.
TOÙM TAÉT Maëc duø moâ hình hoài quy tuyeán tính ñôn hai bieán ñöôïc söû duïng trong chöông naøy, nhöng haàu heát caùc khía caïnh cô baûn cuûa vieäc tieán haønh phaân tích thöïc nghieäm ñaõ ñöôïc ñeà caäp. Thaät höõu ích khi toùm taét laïi caùc keát quaû ñaõ ñöôïc thaûo luaän töø ñaàu ñeán giôø.
Moät moâ hình hoài quy tuyeán tính ñôn laø Yt = α + βXt + ut (t = 1, 2, …, n). Xt vaø Yt laø quan saùt thöù t laàn löôït cuûa bieán ñoäc laäp vaø bieán phuï thuoäc, α vaø β laø caùc thoâng soá cuûa toång theå khoâng bieát seõ ñöôïc öôùc löôïng töø döõ lieäu cuûa X vaø Y, ut soá haïng sai soá khoâng quan saùt ñöôïc, ñaây laø caùc bieán ngaãu nhieân vôùi caùc tính chaát ñöôïc ñeà caäp döôùi ñaây, n laø toång soá quan saùt. Ñoä doác (β) ñöôïc dieãn dòch laø aûnh höôûng caän bieân cuûa söï taêng moät ñôn vò giaù trò Xt leân Yt , α + βXt laø trò trung bình coù ñieàu kieän cuûa Y cho tröôùc X = Xt.
Thuû tuïc bình phöông nhoû nhaát thoâng thöôøng (OLS) cöïc tieåu hoùa toång bình vaø tính toaùn caùc öôùc löôïng (kyù hieäu αˆ vaøβˆ ø) cuûa soá haïng phöông sai soá ∑ 2 tuˆ tung ñoä goác α vaø ñoä doác β. Yeâu caàu duy nhaát ñeå thöïc hieän vieäc öôùc löôïng caùc thoâng soá theo OLS laø n coù giaù trò nhoû nhaát baèng 2 vaø ít nhaát moät trong nhöõng giaù trò cuûa X laø khaùc nhau – nghóa laø, khoâng phaûi taát caû caùc giaù trò cuûa X laø nhö nhau.
= 0
vaø ∑
Neáu ut laø bieán ngaãu nhieân coù giaù trò trung bình baèng 0, vaø Xt cho tröôùc vaø khoâng ngaãu nhieân, thì E(ut) = 0 vaø E(Xtut) = 0. Caùc phöông trình chuaån laø ∑ = 0 ˆt ˆtu tuX . Lôøi giaûi cuûa caùc phöông trình naøy cho keát quaû laø caùc öôùc löôïng theo OLS cuûa α vaø β.
Döôùi caùc giaû ñònh vöøa neâu ra, caùc öôùc löôïng theo OLS laø khoâng thieân leäch vaø nhaát quaùn. Söï nhaát quaùn ñöôïc giöõ nguyeân ngay caû neáu Xt laø ngaãu nhieân, mieãn laø Cov(X, u) = 0 vaø 0 < Var(X) < ∝ - nghóa laø, mieãn laø X vaø u khoâng töông quan vaø X khoâng laø haèng soá.
Ramu Ramanathan
52
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
βˆ
Neáu caùc giaù trò u tuaân theo phaân phoái ñoäc laäp vaø töông töï nhau (iid) vôùi moät phöông sai xaùc ñònh, αˆ vaø βˆ cuõng seõ laø caùc öôùc löôïng khoâng thieân leäch tuyeán tính toát nhaát (BLUE); töùc laø, trong soá taát caû toå hôïp tuyeán tính khoâng thieân leäch cuûa caùc giaù trò cuûa Y, αˆ vaø βˆ coù phöông sai nhoû nhaát. Keát quaû naøy ñöôïc goïi laø ñònh lyù Gauss-Markov vaø coù nghóa raèng, ngoaøi tính chaát khoâng thieân leäch vaø nhaát quaùn, caùc öôùc löôïng theo OLS cuõng laø caùc öôùc löôïng hieäu quaû nhaát. Neáu caùc giaù trò cuûa u tuaân theo phaân phoái chuaån ñoäc laäp vaø töông töï nhau N(0, σ2), caùc öôùc löôïng theo OLS cuõng laø caùc öôùc löôïng thích hôïp nhaát (MLE).
ˆ u
Töø αˆ vaø βˆ , giaù trò döï baùo cuûa Yt (kyù hieäu laø −
βα ˆ ˆ X +=
=
baèng
, vaø phaàn dö ñöôïc öôùc löôïng baèng
Y t
t
ˆ Y t
t
tYˆ ) thu ñöôïc ˆ Y . Sai soá chuaån t
cuûa caùc phaàn dö laø moät öôùc löôïng cuûa ñoä leäch chuaån σ vaø ñöôïc tính theo coâng
)2
−
] 2/1
2 (ˆ nut
[ ˆ ∑ = σ
. Töø caùc keát quaû naøy, ta coù theå suy ra sai soá chuaån thöùc cuûa αˆ vaø βˆ ( αˆs vaø βˆs ). Caùc sai soá chuaån caøng nhoû, ñoä chính xaùc cuûa caùc öôùc löôïng cuûa caùc thoâng soá caøng lôùn. Söï bieán ñoåi cuûa X caøng lôùn caøng toát vì ñieàu naøy coù khuynh höôùng caûi thieän ñoä chính xaùc cuûa caùc öôùc löôïng rieâng leû.
Caùc böôùc tieán haønh kieåm ñònh ñoái thuyeát moät ñaàu veà β ñöôïc tieán haønh nhö
sau:
BÖÔÙC 1 BÖÔÙC 2
ˆ
t c
H0: β = β0 H1: β > β0 ) βββ ( ˆ , trong ñoù βˆs laø sai s 0 / − = Trò thoáng keâ kieåm ñònh laø soá chuaån öôùc löôïng cuûa βˆ . Theo giaû thuyeát H0, giaù trò naøy tuaân theo phaân phoái t vôùi n – 2 baäc töï do.
BÖÔÙC 3
Tra baûng t vôùi giaù trò öùng vôùi n – 2 baäc töï do vaø moät möùc yù nghóa cho tröôùc (chaúng haïn α), vaø tìm ñieåm t* n-2(α) sao cho P(t> t*) = α.
BÖÔÙC 4
Baùc boû H0 taïi möùc yù nghóa α neáu tc > t*. Neáu giaû thuyeát ngöôïc laïi H1 laø β < β0, H0 seõ bò baùc boû neáu tc < - t*.
Kieåm ñònh coù theå ñöôïc thöïc hieän theo moät caùch töông ñöông. Caùc böôùc 3 vaø 4 ñöôïc ñieàu chænh nhö sau: BÖÔÙC 3a Tính xaùc suaát (kyù hieäu laø p-value) sao cho t > |tc|. BÖÔÙC 4a Baùc boû H0 vaø keát luaän laø heä soá coù yù nghóa neáu p-value nhoû hôn
moät möùc yù nghóa naøo ñoù (α).
Ramu Ramanathan
53
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Caùc böôùc kieåm ñònh giaû thuyeát ngöôïc laïi H1 coù tính hai phía ñöôïc thöïc hieän
=
nhö sau: BÖÔÙC 1 BÖÔÙC 2
H0: β = β0 H1: β ≠ β0 Trò thoáng keâ kieåm ñònh laø
. Theo giaû thuyeát H0,
ˆ
tc
) βββ ( ˆ s 0 / −
BÖÔÙC 3
giaù trò tuaân theo phaân phoái t vôùi n – 2 baäc töï do. Tra baûng t vôùi giaù trò öùng vôùi n – 2 baäc töï do vaø moät möùc yù nghóa cho tröôùc (chaúng haïn α), vaø tìm ñieåm t* n-2(α/2) sao cho P(t> t*) = α/2 (moät nöõa cuûa möùc yù nghóa). Baùc boû H0 taïi möùc yù nghóa α neáu |tc| > t*.
BÖÔÙC 4
Caùc böôùc hieäu chænh ñeå thöïc hieän kieåm ñònh theo phöông phaùp p-value nhö
sau: BÖÔÙC 3a Tính p-value = 2P(t > |tc|). BÖÔÙC 4a Baùc boû H0 neáu p-value nhoû hôn moät möùc yù nghóa naøo ñoù (α).
Trò thoáng keâ ño löôøng ñoä thích hôïp cuûa moät moâ hình laø R2 = 1- (ESS/TSS),
trong ñoù
2
ESS
vaø
. R2 coù giaù trò töø 0 ñeán 1. Giaù trò naøy caøng
2ˆtu
− t YY
TSS ∑ =
∑=
cao ñoä thích hôïp caøng toát. R2 mang hai yù nghóa: (1) noù laø tyû leä cuûa toång phöông sai cuûa Y maø moâ hình giaûi thích, vaø (2) noù laø bình phöông cuûa heä soá )tYˆ töông quan giöõa giaù trò quan saùt (Yt) cuûa bieán phuï thuoäc vaø giaù trò döï baùo ( .
H1: ρxy ≠ 0
Kieåm ñònh veà ñoä thích hôïp toång theå cuûa moâ hình coù theå ñöôïc thöïc hieän baèng caùch söû duïng giaù trò R2. Caùc böôùc ñöôïc tieán haønh nhö sau (ρxy laø heä soá töông quan cuûa toång theå cuûa hai bieán X vaø Y): BÖÔÙC 1 BÖÔÙC 2
BÖÔÙC 3
H0: ρxy = 0 Trò thoáng keâ kieåm ñònh laø Fc = R2(n – 2)/(1 – R2). Theo giaû thuyeát H0, trò thoáng keâ naøy tuaân theo phaân phoái F vôùi 1 baäc töï do ôû töû soá vaø n – 2 baäc töï do ôû maãu soá. Tra baûng F theo töû soá 1 baäc töï do vaø maãu soá n – 2 baäc töï do vaø moät möùc yù nghóa cho tröôùc (chaúng haïn α) tìm gí trò F* sao cho: P(F>F*) = α. Baùc boû giaû thuyeát H0 (taïi möùc yù nghóa α) neáu Fc > F*.
BÖÔÙC 4 Khoaûng tin caäy 95% cuûa β ñöôïc xaùc ñònh nhö sau:
Ramu Ramanathan
54
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
* st
* st
−
+
ˆ
ˆ
( ˆ β
)β
ˆ, β β
Y
ˆ ˆ X βα+
=
. Phöông
Döï baùo coù ñieàu kieän cuûa Y, cho tröôùc X baèng X0, laø
0
sai cuûa noù (pheùp ño ñoä tin caäy cuûa döï baùo) tæ leä thuaän vôùi khoaûng caùch cuûa X0 so vôùi giaù trò trung bình X . Nhö vaäy, X0 caøng xa khoûi giaù trò trung bình cuûa X, giaù trò döï baùo caøng keùm tin caäy.
Thay ñoåi thang ño cuûa bieán phuï thuoäc daãn ñeán thay ñoåi töông öùng thang ño cuûa moãi heä soá hoài quy. Tuy nhieân, caùc giaù trò R2 vaø trò thoáng keâ t seõ khoâng ñoåi. Neáu thang ño cuûa moät bieán ñoäc laäp thay ñoåi, heä soá hoài quy cuûa noù vaø caùc heä sai soá chuaån töông öùng bò thay ñoåi cuøng thang ño, tuy nhieân taát caû caùc trò thoáng keâ khaùc khoâng thay ñoåi.
Vieäc xaùc ñònh chính xaùc quan heä nhaân quaû laø heát söùc quan troïng trong moâ hình hoài quy. Giaû thieát chuaån laø X gaây ra Y. Tuy nhieân, neáu X vaø Y ñöôïc traùo ñoåi, vaø moâ hình ñöôïc öôùc löôïng baèng Xt = α* + β*Yt + vt, ñöôøng thaúng hoài quy noùi chung seõ khaùc vôùi ñöôøng ñöôïc xaùc ñònh töø moâ hình Yt = α + βXt + ut. THUAÄT NGÖÕ
Analysis of variance (ANOVA) Best linear unbiased estimator (BLUE) Coefficient of multiple determination Conditional mean of Y given X
Critical region Data-generating process (DGP) Engel curve Error sum of square (ESS) Estimated residual Explained variation Feedback Fitted straight line F-test Gauss-Markov theorem Goodness of fit Heteroscedasticity Homoscedasticity
Phaân tích phöông sai Öôùc löôïng khoâng thieân leäch tuyeán tính toát nhaát Heä soá xaùc ñònh boäi Giaù trò trung bình ñieàu kieän cuûa Y bieát tröôùc X Vuøng ngöôõng (vuøng tôùi haïn) Quaù trình phaùt döõ lieäu Ñöôøng cong Engel Toång bình phöông sai soá Phaàn dö öôùc löôïng Söï bieán ñoåi giaûi thích ñöôïc Phaûn hoài Ñöôøng thaúng thích hôïp Kieåm ñònh F Ñònh lyù Gauss-Markov Ñoä khôùp Phöông sai cuûa sai soá thay ñoåi Ñoàng phöông sai sai soá (tính chaát phöông sai cuûa sai soá khoâng thay
Ramu Ramanathan
55
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Joinly determined Linear estimator Marginal effect of X on Y Mean absolute percent error (MAPE)
Mean squared error (MSE) Mean squared percentage error (MSPE) Method of least square
Nonlinear regression model Normal equation Ordinary least squares (OLS) Population parameters Population regression function Population regression line Population variance Postsample forecast p-value Regression coefficients Regression sum of squares (RSS) Residual Root mean squared error
ñoåi) Ñöôïc xaùc ñònh cuøng luùc Öôùc löôïng tuyeán tính Hieäu öùng caän bieân cuûa X leân Y Sai soá phaàn traêm tuyeät ñoái trung bình Sai soá bình phöông trung bình Sai soá phaàn traêm bình phöông trung bình Phöông phaùp bình phöông toái thieåu Moâ hình hoài quy phi tuyeán Phöông trình chuaån Bình phöông toái thieåu thöôøng Tham soá cuûa toång theå Haøm hoài quy cuûa toång theå Ñöôøng hoài quy cuûa toång theå Phöông sai cuûa toång theå Döï baùo haäu maãu Giaù trò p Caùc heä soá hoài quy Toång bình phöông hoài quy Phaàn dö Caên baäc hai cuûa sai soá bình phöông trung bình Öôùc löôïng cuûa maãu Ñöôøng hoài quy cuûa maãu Haøm hoài quy cuûa maãu Bieåu ñoà phaân taùn cuûa maãu Töông quan chuoãi Ñoäc laäp chuoãi Khaùc 0 moät caùch ñaùng keå Lôùn hôn 0 moät caùch ñaùng keå Moâ hình hoài quy tuyeán tính ñôn Töông quan giaû taïo Hoài quy giaû taïo Sai soá chuaån cuûa heä soá hoài quy
Sample estimate Sample regression line Sample regression function Sample scatter diagram Serial correlation Serial independence Significanly different from zero Significanly greater from zero Simple linear regression model Spurious correlation Spurious regression Standard error of a regression coefficient Standard error of the regression Standard error of the residuals Statistically insignificant
Sai soá chuaån cuûa hoài quy Sai soá chuaån cuûa phaàn dö Khoâng coù yù nghóaveà thoáng keâ
Ramu Ramanathan
56
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Statistically not greater than zero Statistically significant Sum of squares of the residuals (ESS)
Khoâng lôùn hôn 0 veà maët thoáng keâ Coù yù nghóa veà thoáng keâ Toång bình phöông cuûa caùc phaàn dö Toång bình phöông toaøn phaàn Phöông sai toång Trò thoáng keâ t Kieåm ñònh t Bieán ñoåi khoâng giaûi thích ñöôïc Sai soá thay ñoåi ngaãu nhieân Sai soá do nhieãu traéng
t], ta coù phöông sai cuûa sai soá thay ñoåi (phaân taùn khoâng nhö nhau). Hình 3.A.2
Y
f(u)
α + βX
0
Total sum of squares (TSS) Total variance t-statistic t-test Unexplained variation Well-behaved errors White-noise errors 3.A PHUÏ LUÏC Chöùng Minh Caùc Phöông Trình (cid:0) 3.A.1 Bieåu dieãn 3 chieàu cuûa moâ hình tuyeán tính ñôn Hình 3.A.1 bieåu dieãn baèng ñoà thò caùc giaû thieát lieät keâ trong baûng 3.2 cho tröôøng hôïp cuûa moâ hình hoài qui 2 bieán ñôn. Truïc X vaø Y ñaïi dieän cho caùc giaù trò cuûa caùc bieán X vaø Y. Truïc Z laø haøm maät ñoä xaùc suaát f(u) cuûa sai soá ngaãu nhieân u. Ñöôøng thaúng α + βX laø trung bình coù ñieàu kieän cuûa Y vôùi X cho tröôùc, ñöôïc giaû söû laø tuyeán tính. Caùc phaân phoái thoáng keâ ñöôïc veõ xung quanh ñöôøng trung bình cho 3 giaù trò X1, X2 vaø X3 laø caùc phaân phoái coù ñieàu kieän töông öùng. Nhö ñaõ ñeà caäp trong baøi, giaû thieát raèng Var(ut) = σ2 ñöôïc goïi laø phöông sai cuûa sai soá khoâng ñoåi, coù nghóa “phaân taùn nhö nhau”. Hình 3.A.1 moâ taû tính baát bieán cuûa phöông sai cuûa sai soá cho taát caû caùc quan saùt. Neáu caùc phöông sai naøy khoâng baát bieán maø thay ñoåi theo t [nhö vaäy, Var(ut) ≠σ2 minh hoïa tröôøng hôïp phöông sai cuûa sai soá thay ñoåi trong ñoù phöông sai taêng khi X taêng. Tröôøng hôïp naøy ñöôïc xem xeùt chi tieát hôn trong chöông 8. Hình 3.A.1 Bieåu dieãn ñoà thò cuûa Moâ hình Hoài Qui Tuyeán Tính Ñôn
X1
X2
X
X3
Ramu Ramanathan
57
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Hình 3.A.2 Minh hoïa veà phöông sai cuûa sai soá khoâng ñoåi
Y
f(u)
α + βX
0
X1
X2
X
X3
(cid:0) 3.A.2 Caùc Keát Quaû Cuûa Pheùp Tính Toång
Caùc tính chaát 3.1 vaø 3.2 ñöôïc chöùng minh ôû ñaây
TÍNH CHAÁT 3.1
2
2
=
=
=
−
S
− XX
− XnX
X
X
) 2
xx
t
2 t
2 t
t
∑
∑
∑
( ∑
1 n
( ) ( )
Chöùng minh
2
2
2
(
)
(
)
(
)
=
+
−
=
−
+
− XX
X
XX2
X
]
t
[ XXX2X t
2 t
2 t
t
∑
∑
∑
∑
∑
2
. Do ñoù, bieåu thöùc trôû thaønh
. Hôn nöõa
t
2 t
XnXt =∑
2
∑ + (+
. Keát hôïp soá haïng thöù hai vaø ba trong bieåu thöùc ta ñöôïc phaàn thöù
( ) XnXX2X
2 t
Nhö tröôùc ñaây, X nhö nhau vôùi moãi giaù trò t. Do vaäy, bieåu thöùc treân = ∑ − −∑ nhaát cuûa tính chaát. Ta bieát raèng X = (∑Xt)/n. Thay vaøo, ta coù phaàn thöù hai cuûa tính chaát.
) XnXnX2X
TÍNH CHAÁT 3.2
(
)
=
−
−
=
=
−
S
)( YYXX
− YXnYX
]nY ) /
xy
t
t
t
t
YX t t
t
t
∑
∑
∑
[ ( )( ∑∑ X
Chöùng minh
Ramu Ramanathan
58
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
(
)
(
)
−
−
+
=
t
t
t
t
t
=
+
t − YXnYXXYYX
∑
t
t
t
t
)( − YYXX ∑
∑
+
−
t
− YXXYYXYX ∑ − YXnYnXXnYYX t − YXnYX t
t
∑ − = ∑ = ∑
X
Y
vaø
, ta coù ñaúng thöùc thöù hai.
Thay
/) nX t
/) nY t
(∑=
∑= (
Ramu Ramanathan
59
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
(cid:0) 3.A.3 Chöùng Minh Caùc Phöông Trình Chuaån Baèng Pheùp Bình Phöông Nhoû Nhaát Trong phaàn naøy ta aùp duïng phöông phaùp bình phöông nhoû nhaát, ñöôïc trình baøy trong phaàn 3.2 vaø chöùng minh caùc phöông trình chuaån (3.4) vaø (3.5). Tieâu chuaån bình phöông nhoû nhaát laø choïn giaù trò cuûa αˆ vaø βˆ laøm toái thieåu toång bình phöông sai soá:
= nt
= nt
2
2 t
t
∑∑ ˆ = u
= 1t
= 1t
) ( )ˆ,ˆ( =βα ˆ ˆ β−α− ESS X Y t
)
2 t
∑
Ñeå toái thieåu ESS vôùi αˆ vaø βˆ , ta cho ñaïo haøm rieâng (xem phaàn 2.A.3 veà ñaïo haøm rieâng) ∂ESS/∂ αˆ vaø ∂ESS/∂ βˆ baèng 0 vaø giaûi phöông trình naøy. Ta coù
=
=
=
2
2
X
0
t
Y t
t
∑
∑
∑
∂ ESS ˆ α∂
( ˆ ∂ u ˆ α∂
∂
∑
( )( ˆ ˆ β−α− ˆ u2 ) =− 1 ) =− 1 (ˆ u t ˆ ∂ u t ˆ α∂
=
=
=
−
=
=
2
2
X
X
0
t
t
Y t
t
t
∑
∑
∑
∂ ESS ˆ β∂
2 )ˆ( u t ˆ β∂
Töø ñoù ta thu ñöôïc caùc phöông trình sau:
) ( )( ) ˆ ˆ β−α− ˆ u2 (ˆ − Xu t ˆ ∂ u t ˆ β∂
=
X
0
Y t
t
( )
=
Y t
0XX t
t
∑ ∑
Laáy toång töøng soá haïng vaø löu yù raèng αˆ vaø βˆ khoâng phuï thuoäc vaøo t vaø laø thöøa soá chung coù theå ñöa ra ngoaøi caùc toång, ta ñöôïc
nY
X
t
t
( ) ˆ ˆ β−α− ˆ ˆ β−α−
ˆ β+
∑ X
X
t
t
XY t
2 t
∑ ∑
∑
Phöông trình ñaàu tieân töông ñöông vôùi phöông trình (3.4) vaø phöông trình thöù 2 töông ñöông vôùi phöông trình (3.5).
(cid:0) 3.A.4 Öôùc Löôïng Khoâng Thieân Leäch Tuyeán Tính Toát Nhaát (Blue) Vaø Ñònh Lyù Gauss-Markov
Töø lyù thuyeát thoáng keâ ta bieát raèng moät trong nhöõng tính chaát mong muoán cho moät öôùc löôïng laø öôùc löôïng tuyeán tính khoâng thieân leäch phöông sai nhoû nhaát (xem ñònh nghóa 2.8). Noùi caùch khaùc, giöõa caùc toå hôïp tuyeán tính cuûa bieán phuï thuoäc khoâng thieân leäch, ta choïn moät bieán coù
Ramu Ramanathan
60
Thuïc Ñoan/Haøo Thi
ˆ α= ˆ ˆ β+α= ∑
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
phöông sai nhoû nhaát. Ñaây laø öôùc löôïng khoâng thieân leäch toát nhaát (BLUE). Trong phaàn naøy ta chöùng minh ñònh lyù Gauss-Markov, ñònh lyù naøy cho raèng öôùc löôïng OLS ruùt ra trong phaàn 3.2 cuõng coù tính chaát BLUE.
Ñaàu tieân löu yù raèng öôùc soá OLS βˆ thöïc söï coù theå ñöôïc bieåu dieãn nhö laø toå hôïp tuyeán
tính cuûa Yt. Ñeå thaáy ñieàu naøy, ta vieát laïi phöông trình (3.12) döôùi ñaây.
t
Y t
=
−
S
(3.12)
xy
YX t t
∑
( )( ∑∑ X n
)
X
, keát quaû naøy coù theå ñöôïc bieåu dieãn nhö
Löu yù
t /∑= nX
(
)
=
− YXYX t
t
t
− YXX t
t
∑
∑
∑
Vì βˆ = Sxy/Sxx töø phöông trình (3.10), ta coù
=
ω
=βˆ
Y t
Y t
t
∑
∑
− XX t S
xx
phuï thuoäc vaøo Xt. Baây giôø xem
Ñaây laø toå hôïp tuyeán tính cuûa Yt vôùi troïng soá
=ω t
xx
− XX t S ~
, vôùi at coù tính khoâng ngaãu
tYa t
toång toå hôïp tuyeán tính cuûa caùc giaù trò cuûa Y coù daïng ∑=β ~ nhieân. Öôùc löôïng khoâng thieân leäch toát nhaát (BLUE) coù 2 tính chaát: (1) β
khoâng thieân leäch vaø
) laø nhoû nhaát.
~ (2) Var( β Chöùng minh Goïi dt = at - ωt laø hieäu soá caùc troïng soá (löu yù raèng dt chæ phuï thuoäc caùc bieán X vaø do ñoù ñöôïc xem laø khoâng ngaãu nhieân). Vaäy at = ωt + dt. Tieáp theo laø
+ω t
Yd t t
∑
∑
( ~ =β ˆ +β= ) Yd t t
+β=
+β=
β+α
E
d
X
t
t
t
∑
∑
∑
α+β=
t
Xd t
t
∑ β+
d ∑
Ramu Ramanathan
61
Thuïc Ñoan/Haøo Thi
) ( ) ) ( ) ~ +β=β ( YdE t t ( YEd t
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Ñeå coù tính khoâng thieân leäch, keát quaû naøy phaûi baèng β, ñieàu naøy xaûy ra khi vaø chæ khi
vaø
∑dt = 0
∑dtXt = 0
ñöôïc xaùc ñònh bôûi Var[∑(ωt + dt)Yt]. Töø tính chaát 2.A.5c, phöông
~ Phöông sai cuûa öôùc löôïng β sai cuûa toång caùc bieán ngaãu nhieân ñoäc laäp laø toång caùc phöông sai (tính ñoäc laäp ñöôïc ñaûm baûo bôûi giaû thieát 3.6). Hôn nöõa, do giaû thieát 3.5 veà phöông sai cuûa sai soá khoâng ñoåi, ut vaø do vaäy Yt coù phöông sai khoâng ñoåi σ2. Töø ñoù ta coù
) = σ2∑(ωt + dt)2 = σ2∑ω2
t + σ2∑d2
t + 2σ2∑ωtdt
=
, do caùc ñieàu kieän veà tính khoâng thieân leäch
0
d
Soá haïng thöù 3 baèng 0 vì ∑ωtdt =
t
xx
−∑ XX t S
, soá haïng ñaàu ñoäc laäp vôùi caùc bieán choïn dt. Bôûi vì soá haïng thöù hai laø toång caùc bình
~ Var( β
∑dt = 0 vaø ∑dtXt = 0 laøm cho moãi soá haïng trong toång baèng 0. Trong bieåu thöùc treân veà phöông ~ sai cuûa β phöông, chæ coù caùch duy nhaát ñeå toái thieåu soá haïng naøy laø choïn moãi giaù trò cuûa ds baèng 0. Ñieàu naøy laøm cho at = ωt vaø do ñoù , nhö vaäy yù noùi raèng öôùc löôïng OLS thöïc söï laø BLUE vaø do vaäy seõ coù hieäu quaû nhaát. Ñieàu naøy xaùc minh ñònh lyù Gauss-Markov. Cuõng neân löu yù raèng vieäc chöùng minh ñònh lyù caàn ñeán caùc giaû thieát 3.5 vaø 3.6 veà phöông sai cuûa sai soá khoâng ñoåi vaø tính ñoäc laäp theo chuoãi. Neáu moät trong 2 giaû thieát naøy bò vi phaïm, thì phöông phaùp OLS khoâng cho öôùc löôïng hieäu quaû. (cid:0) 3.A.5 Öôùc Löôïng Thích Hôïp Nhaát Lyù do cuûa phöông phaùp öôùc löôïng thích hôïp nhaát ñöôïc dieãn taû chi tieát trong phaàn 2.A.4. Baïn ñoïc coù theå xem phaàn ñoù tröôùc khi baét ñaàu phaàn naøy. Trong phaàn ñoù, phöông phaùp naøy ñaõ ñöôïc aùp duïng cho tröôøng hôïp öôùc löôïng giaù trò trung bình vaø phöông sai cuûa moät phaân phoái chuaån. ÔÛ ñaây ta aùp duïng kyõ thuaät töông töï vaøo baøi toaùn hoài qui. Bôûi vì nguyeân lyù thích hôïp nhaát ñoøi hoûi kieán thöùc veà caùc phaân phoái trong baøi toaùn, neân ta caàn giaû thieát 3.7. Caùc böôùc ñeå xaùc ñònh moät öôùc löôïng thích hôïp nhaát raát deã hieåu. Tröôùc tieân, laäp haøm thích hôïp lieân keát haøm maät ñoä cuûa caùc quan saùt vôùi caùc thoâng soá chöa bieát. Ñeå cöïc ñaïi haøm naøy, laáy vi phaân rieâng phaàn logarit cuûa haøm thích hôïp cho moãi thoâng soá chöa bieát vaø cho baèng 0. Keá ñeán giaûi caùc ñieàu kieän baäc nhaát ñeå tìm caùc öôùc löôïng thích hôïp nhaát. Haøm maät ñoä cuûa u ñöôïc xaùc ñònh theo [xem phöông trình (2.4)]
2
1
σ
2 2 /
=
)( uf
− ue
(
)
σ
π
2
Bôûi vì caùc quan saùt laø ñoäc laäp nhau, haøm thích hôïp cuûa u1, u2,…,un laø
Ramu Ramanathan
62
Thuïc Ñoan/Haøo Thi
~ ˆ β≡β
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
L(α, β, σ2) = f(u1) f(u2) f(u3).. .. ..f(un)
2
1
/(
)
∑−
σ
u
2
2 t
=
e
n
(
)
σ
π
2
2
2 1
(
)
/(
)
∑−
β−α−
σ
X
2
Y t
t
=
e
n
(
)
σ
π
2
Thöïc hieän cöïc ñaïi hoùa logarit cuûa haøm thích hôïp thì deã hôn, giaù trò cöïc ñaïi seõ baèng vôùi giaù trò lôùn nhaát L bôûi vì loga coù tính chaát taêng ñeàu; nghóa laø neáu a>b, thì ln(a)>ln(b).
2
)
(
X
t
Y t
ln
ln
ln(
−=
−σ
) −π
L
n
n
2
2
∑
β−α− σ 2
−=
−σ
n
n
2
SSE 22 σ
Trong ñoù SSE = ∑ (Yt - α -βXt)2. α vaø β chæ xuaát hieän trong soá haïng SSE. Do ñoù, lnL lôùn nhaát baèng vôùi SEE nhoû nhaát (bôûi vì coù daáu aâm tröôùc SSE). Nhöng SSE nhoû nhaát nghóa laø caùc öôùc löôïng bình phöông nhoû nhaát. Do ñoù, caùc öôùc löôïng bình phöông nhoû nhaát cuõng laø MLE vôùi ñieàu kieän caùc sai soá cuûa u tuaân theo phaân phoái N(0,σ2). Bôûi vì caùc öôùc löôïng thích hôïp nhaát laø ñoàng nhaát vaø hieäu quaû moät caùch tieäm caän, neân caùc öôùc löôïng OLS cuõng vaäy. Ñeå coù MLE cuûa σ2, laáy vi phaân rieâng phaàn lnL theo σ vaø cho baèng 0. Ta coù
ln ln( ) −π
∂
−=
+
0
3 =
)
n σ
SEE σ
2
(ln L σ∂
. Nhö ñaõ phaùt bieåu tröôùc ñoù, giaù trò naøy khoâng thieân leäch. Moät öôùc löôïng
nu2 /ˆ t
2σˆ ñaõ xaùc ñònh trong
cho n-2 vaø duøng
)
uXX
t
t
Giaûi phöông trình tìm σ2 ta ñöôïc σ2 = SSE/n. Nhöng SSE phuï thuoäc vaøo α vaø β. Tuy nhieân, ta coù theå duøng caùc öôùc löôïng cuûa chuùng αˆ vaø βˆ . Do ñoù ta thu ñöôïc MLE cuûa phöông sai cuûa ut baèng vôùi khoâng thieân leäch coù theå tìm ñöôïc baèng caùch chia ∑ 2 tuˆ phöông trình (3.21). Ñieàu kieän khoâng thieân leäch ñöôïc chöùng minh trong phuïc luïc phaàn 3.A.7. (cid:0) 3.A.6 Tìm Caùc Phöông Sai Cuûa Caùc Öôùc Löôïng Töø phöông trình (3.10), ta coù βˆ = Sxy/Sxx. Vì X laø khoâng ngaãu nhieân theo giaû thieát 3.4, Sxx cuõng khoâng ngaãu nhieân vaø do ñoù Var(βˆ ) = Var(Sxy)S2 xx. Töø phöông trình (3.15), Sxy = βSxx + Sxu vaø −∑ ( do ñoù Var(Sxy) = Var(Sxu). Töø phöông trình (3.16) ta löu yù raèng Sxu = . Tính chaát
2.A.5c cho thaáy phöông sai cuûa toång caùc bieán ngaãu nhieân laø toång cuûa caùc phöông sai vôùi ñieàu
Ramu Ramanathan
63
Thuïc Ñoan/Haøo Thi
~ ∑=σ
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
kieän ñoàng phöông sai (covariance) caùc soá haïng baèng 0. Theo giaû thieát 3.6, ut vaø us laø khoâng töông quan vôùi moïi t≠s vaø ñoàng phöông sai baèng 0. Do ñoù,
2
(
)
(
)
)
(
)
(
)
=
=
=
SVar
Var
− uXX
Var
− uXX
− XX
uVar
[ (
]
]
xu
t
t
t
t
t
t
[ ∑
∑
∑
2
2
(
)
σ
XX
=σ2Sxx. Töø ñoù seõ coù
t
Vôùi giaû thieát 3.5, Var(ut) = σ2. Do ñoù, Var(Sxu) = ∑ −
2
2
σ
)
xx
xy
=
=
)ˆ( =β
Var
σ S
S
xx
( SVar 2 xx
2
)
− 2n
laø moät öôùc löôïng khoâng thieân leäch cuûa σ2.
S 2 S xx Vaäy ta ñaõ chöùng minh phöông trình (3.18). Thuû tuïc ñeå chöùng minh caùc phöông trình (3.19) vaø (3.20) cuõng töông töï vaø seõ laø baøi taäp cho baïn ñoïc. (cid:0) 3.A.7 Öôùc Löôïng Khoâng Thieân Leäch Cuûa Phöông Sai Cuûa Soá Haïng Sai Soá Theo phöông trình (3.21), s2 =
u2 /()ˆ t
=σ ∑ ˆ (
Ñieàu naøy ñöôïc chöùng minh nhö sau.
(
)
ˆ ˆ β−α−
−
ˆ β−
ˆ β−
Y
X
X
ˆ = Yu t
t
= YX t
t
t
+β+α=
Y
uX
Duøng phöông trình (3.9) cho αˆ . Vì Yt ñöôïc xaùc ñònh bôûi phöông trình (3.1), vôùi u baèng ∑ut/n. Do ñoù, nhoùm taát caû caùc soá haïng β ta coù,
ˆ
(
)
)
β+α=
ˆ u
uX
ˆ β−β+ X
X
t
+ uX t
t
t
)(
)
(
)
( +β+α− ˆ( β−β−
=
−
− XX
u
u
t
t
Toång bình phöông cuûa
tuˆ ñöôïc xaùc ñònh theo
2
2
2
−
u
− XX
2
− uXX
u
2 t
t
t
t
t
∑
2
)
)
− u ˆ( β−β+
∑ ˆ( β−β−
2
S
S
∑ ∑ ˆ = u = S
uu
xx
xu
Duøng kyù hieäu töông töï nhö trong phöông trình (3.11) vaø (3.16). Töø phöông trình (3.15), Sxu= Sxy - βSxx = Sxx( βˆ - β). Thay keát quaû naøy vaøo phöông trình treân vaø keát hôïp caùc soá haïng thöù hai vaø ba ta coù
2
( ) ( ) ) ( )( ) ˆ( ) β−β+ ˆ( β−β−
S
S
2 t
uu
xx
=∑ ˆ u
Ramu Ramanathan
64
Thuïc Ñoan/Haøo Thi
) ˆ( β−β−
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
Ñeå tính giaù trò kyø voïng cuûa toång bình phöông cuûa sai soá, ta caàn E(Suu)vaø E[( βˆ - β)2]. Töø tính chaát 2.11b ta löu yù raèng
E(Suu) = (n-1)Var(u) = (n-1)σ2. Hôn nöõa,
2
E[( βˆ - β)2] = Var( βˆ ) =
σ S
xx
Töø phöông trình (3.18). Ñaët taát caû caùc keát quaû, ta coù
2
2
2 ])
−
=
2 =σ−σ−
E
2n
2 t
ES xx
uu
( ) =∑ ˆ u
Chia cho n-2 ta coù keát quaû mong muoán
2
2 )ˆ( =σ
σ=
E
E
∑ ˆ 2 u t − 2n
2σˆ laø öôùc löôïng khoâng thieân leäch cuûa σ2.
Vaäy, (cid:0) 3.A.8 Chöùng Minh Phöông Trình 3.25 Giaù trò toång bình phöông ñöôïc vieát laïi nhö sau:
2
2
) ( ( ˆ[( β−β ) σ− ( SE ) 1n
=
−
−
t
t
t
t
2
2
ˆ ˆ ( )
+
+
−
−
∑
=
2
∑ − YY t
t
t
t
t
∑
=
−
, hai soá haïng ñaàu tieân laø hai soá haïng coù trong phöông trình (3.25). Baây giôø taát
t
ˆ ˆ YYu t t
ˆ ( ) ( ) ˆ)( YYYY t ) − YY ∑ ( + YYYY ˆ( )ˆ ∑ − YY
−
−
=
=
.
t
t
t
t
∑
Vôùi caû ñieàu ta caàn laø phaûi chöùng minh raèng ∑
ˆ ( ) ) ˆ(ˆ − 0YYu ˆ)( YYYY t
−
=
− YX
t
t
t
t
t
t
∑
∑
∑
∑
∑− ˆ ˆ uYXu
) ) ˆ β+α ˆ α= ˆ β+ ˆ u ˆ(ˆ YYu t ˆ(ˆ u t
=
=
. Töø phöông trình (3.5),
X
0
t
Y t
t
∑
( ) ˆ ˆ β−α− ˆ u
=
=
, vaäy keát quaû ñöôïc chöùng minh.
0
t
Y t
t
∑
( ˆ Xu t ) XX t
Töø phöông trình chuaån ñaàu tieân (3.4), ∑ ˆ ∑ ˆ β−α− (cid:0) 3.A.9 Chöùng Minh Phöông Trình 3.26a Ñeå chöùng minh phöông trình (3.26a), tröôùc tieân ta tìm ñoàng phöông sai maãu (kyù hieäu bôûi Cov) giöõa Yt vaø
tYˆ . Töø phöông trình (2.10),
Ramu Ramanathan
65
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
−
−
Cov (Yt,
tYˆ ) =
t
t
∑
1 − 1n
YX =β+α ˆ ˆ
. Vaäy,
Löu yù raèng trung bình cuûa
tYˆ cuõng laø Y bôûi vì
(
)ˆ
)
)
−
−
+
=
−
+
ˆ( ˆ uYY
ˆ( − YY
= YYYY t
t
t
t
t
t
Do ñoù,
2
)
)
−
∑
∑
+
Cov (Yt,
tYˆ ) =
ˆ(ˆ YYu t t − 1n
ˆ( − YY t − 1n
tYˆ baèng
Phaàn tröôùc ñaõ cho thaáy soá haïng thöù nhaát baèng 0. Do vaäy, ñoàng phöông sai cuûa Yt vaø vôùi soá haïng thöù hai, laø RSS/(n-1);
Cov (Yt,
tYˆ ) =
RSS − 1n
Ta cuõng coù
2
)
=
vaø
Var(
Var(Yt) =
tYˆ ) =
TSS − 1n
ˆ( −∑ YY t − 1n
RSS − 1n
Töø phöông trình (2.7) ta nhôù laïi bình phöông cuûa heä soá töông quan ñôn giöõa Yt vaø
tYˆ ñöôïc xaùc
2
)ˆ,
(
Cov
YY t t
=2 YYr ˆ
)ˆ(
)
( YVar t
YVar t
ñònh bôûi Thay theá ñoàng phöông sai vaø phöông sai töø bieåu thöùc vöøa ruùt ra vaø boû n-1, ta coù
2 RSS
2
=
=
R
2 =ˆ r YY
TSS
RSS
RSS TSS
tYˆ ñöôïc döï baùo bôûi moâ
( ) ˆ)( YYYY
xy = R2 Cho Moâ Hình Hoài Qui Ñôn
Vaäy, bình phöông cuûa töông quan ñôn giöõa giaù trò quan saùt Yt vaø giaù trò hình hoài qui laø nhö nhau vaø laø R2 ñöôïc ñònh nghóa trong phöông trình (3.26). (cid:0) 3.A.10 Chöùng Minh Raèng r2
Ramu Ramanathan
66
Thuïc Ñoan/Haøo Thi
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
2
xy = S2 xy/(SxxSyy). Syy ˆ ˆ β+α= X
. Vì
vaø
ˆ Y t
t
)
) ˆ( t YY
ˆ − YY
− XX
Trong phaàn naøy ta seõ chöùng minh raèng trong tröôøng hôïp moâ hình hoài qui ñôn, R2 cuõng baèng vôùi bình phöông cuûa töông quan ñôn giöõa X vaø Y. Töø phöông trình (2.11), r2 baèng vôùi toång bình phöông TSS. Hôn nöõa, RSS = ∑ − (ˆ β= Y
. Do ñoù,
, ta coù
X
t
t
2
2
2
2
ˆ ˆ β+α=
=
=
=
− XX
t
t
xx
∑
∑
) βˆ ( ) βˆ RSS ˆ( − YY S
S
xy
(
)
ˆ β=
S
S
RSS
xx
xy
S
xx
ˆ β=
Töø phöông trình (3.10), βˆ = Sxy/Sxx. Thay keát quaû naøy cho moät soá haïng βˆ ôû treân, ta thu ñöôïc
ˆ β
Thay theá Sxy töø keát quaû naøy vaø löu yù raèng Syy = TSS, ta coù
2
=
=
=
=
2 r xy
xy TSS
2 xy SS xx
yy
xx
S S R SS xy xy TSS S
2
)ˆ(
ˆ(
)]
=
ˆ[ − XYEYE
Keát quaû ñaõ ñöôïc chöùng minh. (cid:0) 3.A.11 Chöùng Minh Phöông Trình 3.28
0
0
2
2
]
ˆ[(
)]
ˆ ˆ[ β+α=
=
ˆ( β−β
E
X
E
X
0
0
=
)ˆ( +α
β−α− X 0 )ˆ( +β
) +α−α )ˆ,ˆ( βα
Var
X
Var
X2
Cov
2 0
0
YVar 0
Trong pheùp bieán ñoåi treân, ta ñaõ duøng tính chaát 2.4a. Thay töø phöông trình (3.18), (3.19) vaø (3.20), ta ñöôïc
2 t
2
)ˆ(
+
−
2 o
xx
xx
xx
σ= ∑ nS
2
2
+
X 2 X YVar 0 XX 0 S 1 S
Xn
S
Xn
2 t
2 t
xx
X −∑
=∑ X
. Löu yù raèng vaø thay vaøo bieåu thöùc phöông sai, ta Vôùi Sxx =
2
2
(
)
+
X
X
2
2
0
)ˆ(
σ=
+
+
YVar 0
1 n
2 − XX2XX 0 S
1 n
− 0 S
xx
xx
σ=
Ramu Ramanathan
67
Thuïc Ñoan/Haøo Thi
ñöôïc
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
=
ˆ ˆ β+α=
ˆ u
X
0
0
ˆ Y 0
0
ˆ − YY 0 bình. Do ñoù ta coù
)
)ˆ(
)ˆ,
=
+
−
uVar
2
)ˆ( 0
( YVar 0
YVar 0
( YYCov 0
0
0Yˆ ) ñöôïc xaùc ñònh bôûi phöông trình (3.26).
0Yˆ ) = 0, bôûi vì u0 khoâng töông quan vôùi caùc soá dö khaùc vaø do ñoù khoâng
Keát quaû naøy laø phöông trình (3.28). 3.A.12 Chöùng minh phöông trình 3.29 Goïi laø giaù trò döï baùo cuûa trung laø sai soá taïi ñieåm döï baùo cuûa Y0, vôùi
2
(
)
X
X
2
σ=
uVar
)ˆ( 0
1 ++ n
− 0 S
xx
1
Ramu Ramanathan
68
Thuïc Ñoan/Haøo Thi
Vì Y0 = α + βX0 + u0, Var(Y0) = σ2. Maët khaùc, Var( Cuoái cuøng, Cov(Y0, töông quan vôùi αˆ vaø βˆ . Vaäy ta coù
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
BAØI TAÄP
Caâu Hoûi Lyù Thuyeát 3.1. Taát caû caùc toång ôû boán bieåu thöùc döôùi ñaây ñöôïc duøng ñeå tính caùc döõ lieäu maãu, khoâng phaûi daønh cho taäp hôïp hoaøn chænh. Haõy chæ ra nhöõng bieåu thöùc sai vaø ñuùng. Giaûi thích taïi sao nhöõng bieåu thöùc ñoù ñuùng hay sai.
0
=
0
=
a.
t
n ˆ1 t u ∑ = t n ˆ uX ∑ = t 1 t
0
=
b.
t
n t u ∑ = 1
0
=
c.
uX t
t
n ∑ = t 1
0
d.
)ˆ( =tuE
3.2. Coù söï khaùc bieät gì giöõa soá haïng sai soá vaø phaàn dö? Haõy giaûi thích söï khaùc . Vaø giaûi bieät giöõa ut vaø E(ut). Sau ñoù, haõy chöùng minh raèng
+
t
t
u X += βα thích giaù trò kyø voïng nghóa laø gì cuõng nhö neâu ra nhöõng giaû thieát caàn thieát ñeå chöùng minh bieåu thöùc treân. 3.3. Cho moâ hình tuyeán tính ñôn bieán , haõy chöùng minh raèng Y t
)
=
döôùi nhöõng giaû thieát nhaát ñònh, phöông phaùp öôùc löôïng OLS cho keát quaû caùc öôùc löôïng khoâng cheäch? Coù nghóa laø caàn phaûi chöùng minh raèng
YE ( t
^ YE ( t
. Haõy neâu ra nhöõng giaû thieát caàn thieát cho vieäc chöùng minh
ñoù.
3.4. Neâu ra nhöõng giaû thieát caàn thieát cho moãi phaùt bieåu sau. Ñoàng thôøi giaûi
thích lyù do taïi sao nhöõng giaû thieát naøy caàn cho phaùt bieåu ñoù. e. Ñeå öôùc löôïng α vaø β baèng phöông phaùp OLS f. Ñeå chöùng minh raèng caùc öôùc löôïng cuûa caùc thoâng soá theo phöông
phaùp OLS laø khoâng cheäch vaø nhaát quaùn.
g. Ñeå chöùng minh raèng caùc öôùc löôïng theo phöông phaùp OLS laø hieäu
quaû
h. Ñeå thöïc hieän kieåm ñònh t vaø F
3.5. Nhöõng caâu hoûi sau laø ñuùng hay sai? Neáu nhöõng caâu hoûi naøy chæ ñuùng moät phaàn, baïn haõy chæ ra phaàn ñuùng ñoù. Giaûi thích lyù do taïi sao nhöõng caâu (phaàn) ñoù ñuùng. i. Caùc öôùc löôïng heä soá goùc theo phöông phaùp OLS seõ chính xaùc hôn
neáu caùc giaù trò X gaàn vôùi trò trung bình maãu cuûa chuùng hôn. j. Neáu Xt vaø ut töông quan, caùc öôùc löôïng vaãn seõ khoâng cheäch. k. Caùc öôùc löôïng khoâng theå laø öôùc löôïng khoâng leäch tuyeán tính toát nhaát
Ramu Ramanathan
69
Thuïc Ñoan/Haøo Thi
(BLUE) tröø phi taát caû caùc giaù trò ut tuaân theo phaân boá chuaån.
Nhaäp moân kinh teá löôïng vôùi caùc öùng duïng Chöông 3: Moâ hình hoài quy tuyeán tính ñôn
Chöông trình Giaûng daïy Kinh teá Fulbright Nieân khoùa 2003-2004
Phöông phaùp phaân tích Baøi ñoïc
l. Neáu caùc soá haïng sai soá khoâng tuaân theo phaân phoái chuaån thì caùc kieåm
ñònh t vaø F khoâng theå ñöôïc thöïc hieän.
m. Neáu phöông sai cuûa ut lôùn thì caùc khoaûng tin caäy cuûa caùc öôùc löôïng
seõ lôùn (roäng) hôn.
n. Neáu phöông sai cuûa X lôùn thì caùc khoaûng tin caäy öôùc löôïng seõ heïp
hôn.
o. Khi trò soá p-value lôùn thì heä soá seõ khaùc 0 moät caùch ñaùng keå. p. Neáu baïn choïn moät möùc yù nghóa cao hôn thì heä soá hoài qui seõ coù khaû
Ramu Ramanathan
70
Thuïc Ñoan/Haøo Thi
naêng coù yù nghóa hôn.