library(wooldridge)
ls(package:wooldridge)
## Warning in ls(package:wooldridge): 'package:wooldridge' converted to character
## string
## [1] "admnrev" "affairs" "airfare" "alcohol"
## [5] "apple" "approval" "athlet1" "athlet2"
## [9] "attend" "audit" "barium" "beauty"
## [13] "benefits" "beveridge" "big9salary" "bwght"
## [17] "bwght2" "campus" "card" "catholic"
## [21] "cement" "census2000" "ceosal1" "ceosal2"
## [25] "charity" "consump" "corn" "countymurders"
## [29] "cps78_85" "cps91" "crime1" "crime2"
## [33] "crime3" "crime4" "discrim" "driving"
## [37] "earns" "econmath" "elem94_95" "engin"
## [41] "expendshares" "ezanders" "ezunem" "fair"
## [45] "fertil1" "fertil2" "fertil3" "fish"
## [49] "fringe" "gpa1" "gpa2" "gpa3"
## [53] "happiness" "hprice1" "hprice2" "hprice3"
## [57] "hseinv" "htv" "infmrt" "injury"
## [61] "intdef" "intqrt" "inven" "jtrain"
## [65] "jtrain2" "jtrain3" "jtrain98" "k401k"
## [69] "k401ksubs" "kielmc" "labsup" "lawsch85"
## [73] "loanapp" "lowbrth" "mathpnl" "meap00_01"
## [77] "meap01" "meap93" "meapsingle" "minwage"
## [81] "mlb1" "mroz" "murder" "nbasal"
## [85] "ncaa_rpi" "nyse" "okun" "openness"
## [89] "pension" "phillips" "pntsprd" "prison"
## [93] "prminwge" "rdchem" "rdtelec" "recid"
## [97] "rental" "return" "saving" "school93_98"
## [101] "sleep75" "slp75_81" "smoke" "traffic1"
## [105] "traffic2" "twoyear" "volat" "vote1"
## [109] "vote2" "voucher" "wage1" "wage2"
## [113] "wagepan" "wageprc" "wine"
data(sleep75)
View(sleep75)
# câu 1:
#(i) Khi các yếu tố khác được cố định, có bằng chứng cho thấy nam giới ngủ nhiều hơn nữ giới không?
sleep75$age2 <- sleep75$age^2
model <- lm(sleep ~ totwrk + educ + age + age2 + male, data = sleep75)
summary(model)
##
## Call:
## lm(formula = sleep ~ totwrk + educ + age + age2 + male, data = sleep75)
##
## Residuals:
## Min 1Q Median 3Q Max
## -2378.00 -243.29 6.74 259.24 1350.19
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 3840.83197 235.10870 16.336 <2e-16 ***
## totwrk -0.16342 0.01813 -9.013 <2e-16 ***
## educ -11.71332 5.86689 -1.997 0.0463 *
## age -8.69668 11.20746 -0.776 0.4380
## age2 0.12844 0.13390 0.959 0.3378
## male 87.75243 34.32616 2.556 0.0108 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 417.7 on 700 degrees of freedom
## Multiple R-squared: 0.1228, Adjusted R-squared: 0.1165
## F-statistic: 19.59 on 5 and 700 DF, p-value: < 2.2e-16
print("TRA LOI (i) Nam vs nữ = +87.75 (p = 0.011)Nam ngủ nhiều hơn nữ 87.75 phút/tuần")
## [1] "TRA LOI (i) Nam vs nữ = +87.75 (p = 0.011)Nam ngủ nhiều hơn nữ 87.75 phút/tuần"
#-------------------------
#(ii) Có sự đánh đổi nào giữa thời gian làm việc và thời gian ngủ không? Ước lượng mức đánh đổi đó là bao nhiêu?
print("TRA LOI (ii) Làm việc - ngủ = −0.163 (p < 0.01) làm thêm 1 giờ ngủ giảm 9.01 phút.")
## [1] "TRA LOI (ii) Làm việc - ngủ = −0.163 (p < 0.01) làm thêm 1 giờ ngủ giảm 9.01 phút."
#-------------------------
#(iii) Cần chạy thêm mô hình hồi quy nào khác để kiểm định giả thuyết không (H₀) rằng, khi giữ yếu tố khác cố định, tuổi (age) không có ảnh hưởng đến thời gian ngủ?
print("TRA LOI (iii) Theo model tren: (iii) Tuổi - giac ngu: p > 0.05 -> Không có ý nghĩa thống kê. Model cho thấy tuổi không a/h đến giác ngủ.")
## [1] "TRA LOI (iii) Theo model tren: (iii) Tuổi - giac ngu: p > 0.05 -> Không có ý nghĩa thống kê. Model cho thấy tuổi không a/h đến giác ngủ."
# Giả thiết H0: Giả sử mô hình rút gọn và model đầy đủ là như nhau
model_reduced <- lm(sleep ~ totwrk + educ + male, data = sleep75)
model <- lm(sleep ~ totwrk + educ + age + age2 + male, data = sleep75)
# Kiểm định ANOVA để so sánh 2 mô hình:
anova(model_reduced, model)
## Analysis of Variance Table
##
## Model 1: sleep ~ totwrk + educ + male
## Model 2: sleep ~ totwrk + educ + age + age2 + male
## Res.Df RSS Df Sum of Sq F Pr(>F)
## 1 702 122631662
## 2 700 122147777 2 483885 1.3865 0.2506
print("Tra loi (iii) p-value = 0.2506 > 0.05 -> Không bác bỏ H0 -> Tuổi không a/h đến giấc ngủ -> age vào mô hình không cải thiện đáng kể khả năng giải thích. Tuổi không có ảnh hưởng đáng kể đến thời gian ngủ khi giữ các yếu tố khác cố định.")
## [1] "Tra loi (iii) p-value = 0.2506 > 0.05 -> Không bác bỏ H0 -> Tuổi không a/h đến giấc ngủ -> age vào mô hình không cải thiện đáng kể khả năng giải thích. Tuổi không có ảnh hưởng đáng kể đến thời gian ngủ khi giữ các yếu tố khác cố định."