# nhập liệu
datafile1=read.csv("/Users/theanh/Downloads/Likert.csv",header=TRUE)
save(datafile1,file="datafile1.rda")
attach(datafile1)
is.data.frame(datafile1)
## [1] TRUE
#CÂU 14
mreg=lm(Ftb~Atb+Btb+Ctb)
summary(mreg)
##
## Call:
## lm(formula = Ftb ~ Atb + Btb + Ctb)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.73037 -0.44846 0.09795 0.39376 1.34503
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.76022 0.32504 2.339 0.020361 *
## Atb 0.20139 0.05691 3.539 0.000504 ***
## Btb 0.13162 0.04943 2.663 0.008399 **
## Ctb 0.47633 0.04730 10.070 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5705 on 194 degrees of freedom
## Multiple R-squared: 0.4128, Adjusted R-squared: 0.4038
## F-statistic: 45.47 on 3 and 194 DF, p-value: < 2.2e-16
### NHẬN XÉT CÂU 14
# Phương trình hồi quy:
#Ftb = 0.76022 + 0.20139*Atb + 0.13162*Btb + 0.47633*Ctb + ei
# Mô hình có ý nghĩa thống kê với F = 45,47 và p < 0,001.
# R² = 0,4128; R² hiệu chỉnh = 0,4038, mô hình giải thích 40,38% sự biến thiên của Ftb.
# Atb, Btb và Ctb đều tác động cùng chiều và có ý nghĩa thống kê với Ftb.
# Ctb có tác động mạnh nhất với hệ số 0,47633.
#CÂU 15
mreg=lm(Ftb~Atb+Btb+Ctb)
summary(mreg)
##
## Call:
## lm(formula = Ftb ~ Atb + Btb + Ctb)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.73037 -0.44846 0.09795 0.39376 1.34503
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.76022 0.32504 2.339 0.020361 *
## Atb 0.20139 0.05691 3.539 0.000504 ***
## Btb 0.13162 0.04943 2.663 0.008399 **
## Ctb 0.47633 0.04730 10.070 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5705 on 194 degrees of freedom
## Multiple R-squared: 0.4128, Adjusted R-squared: 0.4038
## F-statistic: 45.47 on 3 and 194 DF, p-value: < 2.2e-16
library(car)
## Loading required package: carData
vif(mreg)
## Atb Btb Ctb
## 1.004380 1.030053 1.027108
require(lmtest)
## Loading required package: lmtest
## Loading required package: zoo
##
## Attaching package: 'zoo'
## The following objects are masked from 'package:base':
##
## as.Date, as.Date.numeric
dwtest(mreg)
##
## Durbin-Watson test
##
## data: mreg
## DW = 1.7709, p-value = 0.04459
## alternative hypothesis: true autocorrelation is greater than 0
#NHẬN XÉT:Mô hình phù hợp và có ý nghĩa thống kê; Atb, Btb, Ctb đều tác động cùng chiều đến Ftb. Adjusted R² = 40,38%
#VIF ≈ 1 nên không có đa cộng tuyến.
#DW = 1,7709 và p = 0,04459 cho thấy có tự tương quan dương.
#CÂU 16
mreg=lm(Ftb~Atb+Btb+Ctb)
summary(mreg)
##
## Call:
## lm(formula = Ftb ~ Atb + Btb + Ctb)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.73037 -0.44846 0.09795 0.39376 1.34503
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.76022 0.32504 2.339 0.020361 *
## Atb 0.20139 0.05691 3.539 0.000504 ***
## Btb 0.13162 0.04943 2.663 0.008399 **
## Ctb 0.47633 0.04730 10.070 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5705 on 194 degrees of freedom
## Multiple R-squared: 0.4128, Adjusted R-squared: 0.4038
## F-statistic: 45.47 on 3 and 194 DF, p-value: < 2.2e-16
op<-par(mfrow=c(2,2))
plot(mreg)
#NHẬN XÉT:
#Đồ thị 1: Phần dư tập trung quanh đường y = 0 → giả định trung bình phần dư bằng 0 được chấp nhận.
#Đồ thị 2: Các phần dư nằm gần đường chuẩn → phần dư tuân theo phân phối chuẩn.
#Đồ thị 3: Phần dư chuẩn hóa phân bố tương đối đồng đều → phương sai phần dư ổn định.
#Đồ thị 4: Các điểm phần dư phân bố phù hợp → giả định của mô hình được chấp nhận.