# Nhập liệu
likert=file.choose()
datalikert=read.csv(likert,header=TRUE)
save(datalikert,file="datalikert.rda")
attach(datalikert)
is.data.frame(datalikert)
## [1] TRUE
### CÂU 1
# Thống kê các biến
library(psych)
ABCDEFtb=data.frame(Atb,Btb,Ctb,Dtb,Etb,Ftb)
describe(ABCDEFtb)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## Atb 1 198 4.06 0.72 4.00 4.08 1.19 2.20 5 2.80 -0.15 -1.14 0.05
## Btb 2 198 3.95 0.83 4.00 4.04 1.11 1.25 5 3.75 -0.89 0.37 0.06
## Ctb 3 198 3.68 0.87 3.62 3.71 0.93 1.50 5 3.50 -0.11 -0.61 0.06
## Dtb 4 198 3.80 0.78 4.00 3.85 0.49 1.00 5 4.00 -0.68 0.71 0.06
## Etb 5 198 3.95 0.81 4.00 4.04 0.74 1.25 5 3.75 -0.97 1.02 0.06
## Ftb 6 198 3.85 0.74 3.90 3.84 1.04 2.00 5 3.00 -0.04 -1.12 0.05
# Tính Mode biến A1
library(DescTools)
##
## Attaching package: 'DescTools'
## The following objects are masked from 'package:psych':
##
## AUC, ICC, SD
Mode(A1)
## [1] 5
## attr(,"freq")
## [1] 94
table(A1)
## A1
## 3 4 5
## 65 39 94
### NHẬN XÉT CÂU 1: Biến có giá trị trung bình cao nhất là Atb = 4.06
# Biến có độ phân tán lớn nhất là biến Btb
# Mode của biến A1 là 5 với 94 quan sát
### CÂU 2
# Mã hóa biến Giới tính
GTINHMH=GTINH
GTINHMH=replace(GTINHMH,GTINH==1,"Nam")
GTINHMH=replace(GTINHMH,GTINH==2,"NỮ")
GTcode=data.frame(GTINH,GTINHMH)
GTcode
## GTINH GTINHMH
## 1 2 NỮ
## 2 2 NỮ
## 3 2 NỮ
## 4 1 Nam
## 5 2 NỮ
## 6 1 Nam
## 7 1 Nam
## 8 1 Nam
## 9 2 NỮ
## 10 1 Nam
## 11 1 Nam
## 12 1 Nam
## 13 2 NỮ
## 14 1 Nam
## 15 2 NỮ
## 16 1 Nam
## 17 1 Nam
## 18 1 Nam
## 19 2 NỮ
## 20 1 Nam
## 21 1 Nam
## 22 1 Nam
## 23 1 Nam
## 24 1 Nam
## 25 1 Nam
## 26 1 Nam
## 27 1 Nam
## 28 1 Nam
## 29 1 Nam
## 30 2 NỮ
## 31 2 NỮ
## 32 2 NỮ
## 33 1 Nam
## 34 2 NỮ
## 35 2 NỮ
## 36 2 NỮ
## 37 1 Nam
## 38 2 NỮ
## 39 1 Nam
## 40 1 Nam
## 41 1 Nam
## 42 2 NỮ
## 43 1 Nam
## 44 2 NỮ
## 45 1 Nam
## 46 1 Nam
## 47 2 NỮ
## 48 1 Nam
## 49 2 NỮ
## 50 2 NỮ
## 51 1 Nam
## 52 2 NỮ
## 53 2 NỮ
## 54 1 Nam
## 55 2 NỮ
## 56 2 NỮ
## 57 1 Nam
## 58 2 NỮ
## 59 2 NỮ
## 60 2 NỮ
## 61 2 NỮ
## 62 2 NỮ
## 63 2 NỮ
## 64 2 NỮ
## 65 2 NỮ
## 66 1 Nam
## 67 1 Nam
## 68 2 NỮ
## 69 2 NỮ
## 70 2 NỮ
## 71 2 NỮ
## 72 2 NỮ
## 73 2 NỮ
## 74 1 Nam
## 75 2 NỮ
## 76 2 NỮ
## 77 2 NỮ
## 78 2 NỮ
## 79 1 Nam
## 80 2 NỮ
## 81 2 NỮ
## 82 2 NỮ
## 83 1 Nam
## 84 2 NỮ
## 85 1 Nam
## 86 1 Nam
## 87 1 Nam
## 88 2 NỮ
## 89 1 Nam
## 90 2 NỮ
## 91 1 Nam
## 92 1 Nam
## 93 2 NỮ
## 94 1 Nam
## 95 2 NỮ
## 96 2 NỮ
## 97 1 Nam
## 98 2 NỮ
## 99 2 NỮ
## 100 1 Nam
## 101 2 NỮ
## 102 2 NỮ
## 103 1 Nam
## 104 1 Nam
## 105 1 Nam
## 106 1 Nam
## 107 1 Nam
## 108 1 Nam
## 109 2 NỮ
## 110 1 Nam
## 111 2 NỮ
## 112 1 Nam
## 113 2 NỮ
## 114 2 NỮ
## 115 1 Nam
## 116 2 NỮ
## 117 1 Nam
## 118 1 Nam
## 119 1 Nam
## 120 1 Nam
## 121 1 Nam
## 122 1 Nam
## 123 1 Nam
## 124 1 Nam
## 125 1 Nam
## 126 1 Nam
## 127 1 Nam
## 128 1 Nam
## 129 1 Nam
## 130 1 Nam
## 131 2 NỮ
## 132 2 NỮ
## 133 1 Nam
## 134 1 Nam
## 135 2 NỮ
## 136 2 NỮ
## 137 1 Nam
## 138 2 NỮ
## 139 1 Nam
## 140 1 Nam
## 141 2 NỮ
## 142 1 Nam
## 143 2 NỮ
## 144 1 Nam
## 145 1 Nam
## 146 1 Nam
## 147 1 Nam
## 148 1 Nam
## 149 1 Nam
## 150 1 Nam
## 151 1 Nam
## 152 1 Nam
## 153 1 Nam
## 154 1 Nam
## 155 1 Nam
## 156 1 Nam
## 157 1 Nam
## 158 1 Nam
## 159 1 Nam
## 160 1 Nam
## 161 1 Nam
## 162 1 Nam
## 163 1 Nam
## 164 2 NỮ
## 165 2 NỮ
## 166 1 Nam
## 167 2 NỮ
## 168 1 Nam
## 169 1 Nam
## 170 1 Nam
## 171 1 Nam
## 172 2 NỮ
## 173 1 Nam
## 174 2 NỮ
## 175 1 Nam
## 176 2 NỮ
## 177 1 Nam
## 178 1 Nam
## 179 1 Nam
## 180 2 NỮ
## 181 1 Nam
## 182 1 Nam
## 183 2 NỮ
## 184 2 NỮ
## 185 1 Nam
## 186 1 Nam
## 187 1 Nam
## 188 2 NỮ
## 189 1 Nam
## 190 1 Nam
## 191 1 Nam
## 192 1 Nam
## 193 1 Nam
## 194 1 Nam
## 195 1 Nam
## 196 1 Nam
## 197 1 Nam
## 198 1 Nam
# Bảng tần số Giới tính
table(GTINHMH)
## GTINHMH
## Nam NỮ
## 121 77
table(GTINHMH)/sum(table(GTINHMH))
## GTINHMH
## Nam NỮ
## 0.6111111 0.3888889
### NHẬN XÉT CÂU 2: Trong tổng thể 198 người; có 121 nam (61,1%) và 77 nữ (38,9%)
### CÂU 3
# Mã hóa biến Độ tuổi
DTUOIMH=DTUOI
DTUOIMH=replace(DTUOIMH,DTUOI==1,"< 20 tuổi")
DTUOIMH=replace(DTUOIMH,DTUOI==2,"Từ 20 - 30 tuổi")
DTUOIMH=replace(DTUOIMH,DTUOI==3,"Từ 31 - 50 tuổi")
DTUOIMH=replace(DTUOIMH,DTUOI==4,"Trên 50 tuổi")
DTcode=data.frame(DTUOI,DTUOIMH)
DTcode
## DTUOI DTUOIMH
## 1 1 < 20 tuổi
## 2 2 Từ 20 - 30 tuổi
## 3 1 < 20 tuổi
## 4 2 Từ 20 - 30 tuổi
## 5 1 < 20 tuổi
## 6 2 Từ 20 - 30 tuổi
## 7 1 < 20 tuổi
## 8 1 < 20 tuổi
## 9 2 Từ 20 - 30 tuổi
## 10 1 < 20 tuổi
## 11 2 Từ 20 - 30 tuổi
## 12 1 < 20 tuổi
## 13 2 Từ 20 - 30 tuổi
## 14 1 < 20 tuổi
## 15 4 Trên 50 tuổi
## 16 3 Từ 31 - 50 tuổi
## 17 2 Từ 20 - 30 tuổi
## 18 1 < 20 tuổi
## 19 3 Từ 31 - 50 tuổi
## 20 4 Trên 50 tuổi
## 21 2 Từ 20 - 30 tuổi
## 22 3 Từ 31 - 50 tuổi
## 23 3 Từ 31 - 50 tuổi
## 24 4 Trên 50 tuổi
## 25 3 Từ 31 - 50 tuổi
## 26 2 Từ 20 - 30 tuổi
## 27 3 Từ 31 - 50 tuổi
## 28 1 < 20 tuổi
## 29 3 Từ 31 - 50 tuổi
## 30 3 Từ 31 - 50 tuổi
## 31 2 Từ 20 - 30 tuổi
## 32 3 Từ 31 - 50 tuổi
## 33 3 Từ 31 - 50 tuổi
## 34 4 Trên 50 tuổi
## 35 3 Từ 31 - 50 tuổi
## 36 2 Từ 20 - 30 tuổi
## 37 3 Từ 31 - 50 tuổi
## 38 1 < 20 tuổi
## 39 3 Từ 31 - 50 tuổi
## 40 3 Từ 31 - 50 tuổi
## 41 2 Từ 20 - 30 tuổi
## 42 1 < 20 tuổi
## 43 3 Từ 31 - 50 tuổi
## 44 1 < 20 tuổi
## 45 4 Trên 50 tuổi
## 46 2 Từ 20 - 30 tuổi
## 47 3 Từ 31 - 50 tuổi
## 48 1 < 20 tuổi
## 49 3 Từ 31 - 50 tuổi
## 50 4 Trên 50 tuổi
## 51 1 < 20 tuổi
## 52 2 Từ 20 - 30 tuổi
## 53 1 < 20 tuổi
## 54 3 Từ 31 - 50 tuổi
## 55 4 Trên 50 tuổi
## 56 4 Trên 50 tuổi
## 57 4 Trên 50 tuổi
## 58 3 Từ 31 - 50 tuổi
## 59 3 Từ 31 - 50 tuổi
## 60 4 Trên 50 tuổi
## 61 2 Từ 20 - 30 tuổi
## 62 4 Trên 50 tuổi
## 63 3 Từ 31 - 50 tuổi
## 64 3 Từ 31 - 50 tuổi
## 65 3 Từ 31 - 50 tuổi
## 66 3 Từ 31 - 50 tuổi
## 67 3 Từ 31 - 50 tuổi
## 68 4 Trên 50 tuổi
## 69 3 Từ 31 - 50 tuổi
## 70 4 Trên 50 tuổi
## 71 2 Từ 20 - 30 tuổi
## 72 3 Từ 31 - 50 tuổi
## 73 3 Từ 31 - 50 tuổi
## 74 4 Trên 50 tuổi
## 75 3 Từ 31 - 50 tuổi
## 76 2 Từ 20 - 30 tuổi
## 77 4 Trên 50 tuổi
## 78 4 Trên 50 tuổi
## 79 4 Trên 50 tuổi
## 80 3 Từ 31 - 50 tuổi
## 81 1 < 20 tuổi
## 82 2 Từ 20 - 30 tuổi
## 83 3 Từ 31 - 50 tuổi
## 84 3 Từ 31 - 50 tuổi
## 85 2 Từ 20 - 30 tuổi
## 86 3 Từ 31 - 50 tuổi
## 87 2 Từ 20 - 30 tuổi
## 88 1 < 20 tuổi
## 89 1 < 20 tuổi
## 90 3 Từ 31 - 50 tuổi
## 91 3 Từ 31 - 50 tuổi
## 92 1 < 20 tuổi
## 93 2 Từ 20 - 30 tuổi
## 94 3 Từ 31 - 50 tuổi
## 95 1 < 20 tuổi
## 96 3 Từ 31 - 50 tuổi
## 97 2 Từ 20 - 30 tuổi
## 98 4 Trên 50 tuổi
## 99 4 Trên 50 tuổi
## 100 2 Từ 20 - 30 tuổi
## 101 2 Từ 20 - 30 tuổi
## 102 2 Từ 20 - 30 tuổi
## 103 1 < 20 tuổi
## 104 2 Từ 20 - 30 tuổi
## 105 2 Từ 20 - 30 tuổi
## 106 1 < 20 tuổi
## 107 2 Từ 20 - 30 tuổi
## 108 2 Từ 20 - 30 tuổi
## 109 2 Từ 20 - 30 tuổi
## 110 4 Trên 50 tuổi
## 111 3 Từ 31 - 50 tuổi
## 112 2 Từ 20 - 30 tuổi
## 113 3 Từ 31 - 50 tuổi
## 114 3 Từ 31 - 50 tuổi
## 115 4 Trên 50 tuổi
## 116 3 Từ 31 - 50 tuổi
## 117 2 Từ 20 - 30 tuổi
## 118 3 Từ 31 - 50 tuổi
## 119 4 Trên 50 tuổi
## 120 3 Từ 31 - 50 tuổi
## 121 4 Trên 50 tuổi
## 122 3 Từ 31 - 50 tuổi
## 123 2 Từ 20 - 30 tuổi
## 124 3 Từ 31 - 50 tuổi
## 125 4 Trên 50 tuổi
## 126 2 Từ 20 - 30 tuổi
## 127 4 Trên 50 tuổi
## 128 3 Từ 31 - 50 tuổi
## 129 2 Từ 20 - 30 tuổi
## 130 4 Trên 50 tuổi
## 131 3 Từ 31 - 50 tuổi
## 132 3 Từ 31 - 50 tuổi
## 133 2 Từ 20 - 30 tuổi
## 134 3 Từ 31 - 50 tuổi
## 135 3 Từ 31 - 50 tuổi
## 136 4 Trên 50 tuổi
## 137 2 Từ 20 - 30 tuổi
## 138 3 Từ 31 - 50 tuổi
## 139 3 Từ 31 - 50 tuổi
## 140 4 Trên 50 tuổi
## 141 3 Từ 31 - 50 tuổi
## 142 4 Trên 50 tuổi
## 143 3 Từ 31 - 50 tuổi
## 144 2 Từ 20 - 30 tuổi
## 145 3 Từ 31 - 50 tuổi
## 146 3 Từ 31 - 50 tuổi
## 147 3 Từ 31 - 50 tuổi
## 148 4 Trên 50 tuổi
## 149 4 Trên 50 tuổi
## 150 3 Từ 31 - 50 tuổi
## 151 2 Từ 20 - 30 tuổi
## 152 1 < 20 tuổi
## 153 3 Từ 31 - 50 tuổi
## 154 1 < 20 tuổi
## 155 3 Từ 31 - 50 tuổi
## 156 1 < 20 tuổi
## 157 2 Từ 20 - 30 tuổi
## 158 3 Từ 31 - 50 tuổi
## 159 1 < 20 tuổi
## 160 2 Từ 20 - 30 tuổi
## 161 3 Từ 31 - 50 tuổi
## 162 3 Từ 31 - 50 tuổi
## 163 1 < 20 tuổi
## 164 3 Từ 31 - 50 tuổi
## 165 2 Từ 20 - 30 tuổi
## 166 2 Từ 20 - 30 tuổi
## 167 3 Từ 31 - 50 tuổi
## 168 3 Từ 31 - 50 tuổi
## 169 2 Từ 20 - 30 tuổi
## 170 4 Trên 50 tuổi
## 171 2 Từ 20 - 30 tuổi
## 172 4 Trên 50 tuổi
## 173 4 Trên 50 tuổi
## 174 1 < 20 tuổi
## 175 1 < 20 tuổi
## 176 3 Từ 31 - 50 tuổi
## 177 3 Từ 31 - 50 tuổi
## 178 1 < 20 tuổi
## 179 3 Từ 31 - 50 tuổi
## 180 2 Từ 20 - 30 tuổi
## 181 2 Từ 20 - 30 tuổi
## 182 2 Từ 20 - 30 tuổi
## 183 4 Trên 50 tuổi
## 184 4 Trên 50 tuổi
## 185 3 Từ 31 - 50 tuổi
## 186 1 < 20 tuổi
## 187 3 Từ 31 - 50 tuổi
## 188 1 < 20 tuổi
## 189 2 Từ 20 - 30 tuổi
## 190 3 Từ 31 - 50 tuổi
## 191 1 < 20 tuổi
## 192 1 < 20 tuổi
## 193 2 Từ 20 - 30 tuổi
## 194 1 < 20 tuổi
## 195 3 Từ 31 - 50 tuổi
## 196 2 Từ 20 - 30 tuổi
## 197 4 Trên 50 tuổi
## 198 4 Trên 50 tuổi
# Mã hóa biến Thu nhập
TNHAPMH=TNHAP
TNHAPMH=replace(TNHAPMH,TNHAP==1,"< 10 triệu")
TNHAPMH=replace(TNHAPMH,TNHAP==2,"Từ 10 - 20 triệu")
TNHAPMH=replace(TNHAPMH,TNHAP==3,"Trên 20 triệu")
TNcode=data.frame(DTUOI,TNHAPMH)
TNcode
## DTUOI TNHAPMH
## 1 1 Trên 20 triệu
## 2 2 Trên 20 triệu
## 3 1 Từ 10 - 20 triệu
## 4 2 Từ 10 - 20 triệu
## 5 1 < 10 triệu
## 6 2 Từ 10 - 20 triệu
## 7 1 Trên 20 triệu
## 8 1 Từ 10 - 20 triệu
## 9 2 Trên 20 triệu
## 10 1 Trên 20 triệu
## 11 2 Từ 10 - 20 triệu
## 12 1 Từ 10 - 20 triệu
## 13 2 Từ 10 - 20 triệu
## 14 1 Từ 10 - 20 triệu
## 15 4 Trên 20 triệu
## 16 3 < 10 triệu
## 17 2 Trên 20 triệu
## 18 1 < 10 triệu
## 19 3 Trên 20 triệu
## 20 4 Trên 20 triệu
## 21 2 Trên 20 triệu
## 22 3 Trên 20 triệu
## 23 3 Từ 10 - 20 triệu
## 24 4 Từ 10 - 20 triệu
## 25 3 < 10 triệu
## 26 2 Từ 10 - 20 triệu
## 27 3 Trên 20 triệu
## 28 1 Từ 10 - 20 triệu
## 29 3 Trên 20 triệu
## 30 3 Trên 20 triệu
## 31 2 Từ 10 - 20 triệu
## 32 3 Từ 10 - 20 triệu
## 33 3 Từ 10 - 20 triệu
## 34 4 < 10 triệu
## 35 3 Từ 10 - 20 triệu
## 36 2 Trên 20 triệu
## 37 3 Trên 20 triệu
## 38 1 Từ 10 - 20 triệu
## 39 3 Trên 20 triệu
## 40 3 Từ 10 - 20 triệu
## 41 2 Từ 10 - 20 triệu
## 42 1 Từ 10 - 20 triệu
## 43 3 Trên 20 triệu
## 44 1 Trên 20 triệu
## 45 4 Trên 20 triệu
## 46 2 Trên 20 triệu
## 47 3 Trên 20 triệu
## 48 1 Từ 10 - 20 triệu
## 49 3 Từ 10 - 20 triệu
## 50 4 Từ 10 - 20 triệu
## 51 1 Trên 20 triệu
## 52 2 Trên 20 triệu
## 53 1 Từ 10 - 20 triệu
## 54 3 Từ 10 - 20 triệu
## 55 4 Trên 20 triệu
## 56 4 Từ 10 - 20 triệu
## 57 4 Trên 20 triệu
## 58 3 Trên 20 triệu
## 59 3 Từ 10 - 20 triệu
## 60 4 Từ 10 - 20 triệu
## 61 2 Trên 20 triệu
## 62 4 Trên 20 triệu
## 63 3 Từ 10 - 20 triệu
## 64 3 < 10 triệu
## 65 3 Trên 20 triệu
## 66 3 Từ 10 - 20 triệu
## 67 3 Từ 10 - 20 triệu
## 68 4 Từ 10 - 20 triệu
## 69 3 Trên 20 triệu
## 70 4 Trên 20 triệu
## 71 2 Từ 10 - 20 triệu
## 72 3 Từ 10 - 20 triệu
## 73 3 Từ 10 - 20 triệu
## 74 4 Từ 10 - 20 triệu
## 75 3 Trên 20 triệu
## 76 2 Từ 10 - 20 triệu
## 77 4 Từ 10 - 20 triệu
## 78 4 Từ 10 - 20 triệu
## 79 4 Từ 10 - 20 triệu
## 80 3 Trên 20 triệu
## 81 1 Từ 10 - 20 triệu
## 82 2 Từ 10 - 20 triệu
## 83 3 Từ 10 - 20 triệu
## 84 3 Trên 20 triệu
## 85 2 Từ 10 - 20 triệu
## 86 3 Trên 20 triệu
## 87 2 Từ 10 - 20 triệu
## 88 1 Từ 10 - 20 triệu
## 89 1 Trên 20 triệu
## 90 3 Từ 10 - 20 triệu
## 91 3 Từ 10 - 20 triệu
## 92 1 Trên 20 triệu
## 93 2 Từ 10 - 20 triệu
## 94 3 Từ 10 - 20 triệu
## 95 1 Trên 20 triệu
## 96 3 Từ 10 - 20 triệu
## 97 2 Từ 10 - 20 triệu
## 98 4 Trên 20 triệu
## 99 4 Trên 20 triệu
## 100 2 Từ 10 - 20 triệu
## 101 2 Trên 20 triệu
## 102 2 Từ 10 - 20 triệu
## 103 1 Trên 20 triệu
## 104 2 Trên 20 triệu
## 105 2 Từ 10 - 20 triệu
## 106 1 Trên 20 triệu
## 107 2 Từ 10 - 20 triệu
## 108 2 Trên 20 triệu
## 109 2 Từ 10 - 20 triệu
## 110 4 Từ 10 - 20 triệu
## 111 3 < 10 triệu
## 112 2 Trên 20 triệu
## 113 3 Từ 10 - 20 triệu
## 114 3 Từ 10 - 20 triệu
## 115 4 Trên 20 triệu
## 116 3 Trên 20 triệu
## 117 2 Từ 10 - 20 triệu
## 118 3 Từ 10 - 20 triệu
## 119 4 Trên 20 triệu
## 120 3 Trên 20 triệu
## 121 4 Trên 20 triệu
## 122 3 Từ 10 - 20 triệu
## 123 2 Trên 20 triệu
## 124 3 Từ 10 - 20 triệu
## 125 4 Trên 20 triệu
## 126 2 Từ 10 - 20 triệu
## 127 4 Từ 10 - 20 triệu
## 128 3 Trên 20 triệu
## 129 2 Từ 10 - 20 triệu
## 130 4 Trên 20 triệu
## 131 3 Từ 10 - 20 triệu
## 132 3 Trên 20 triệu
## 133 2 Trên 20 triệu
## 134 3 Từ 10 - 20 triệu
## 135 3 Trên 20 triệu
## 136 4 Từ 10 - 20 triệu
## 137 2 < 10 triệu
## 138 3 Trên 20 triệu
## 139 3 Từ 10 - 20 triệu
## 140 4 Từ 10 - 20 triệu
## 141 3 Trên 20 triệu
## 142 4 Từ 10 - 20 triệu
## 143 3 Từ 10 - 20 triệu
## 144 2 Trên 20 triệu
## 145 3 < 10 triệu
## 146 3 < 10 triệu
## 147 3 Từ 10 - 20 triệu
## 148 4 Từ 10 - 20 triệu
## 149 4 Từ 10 - 20 triệu
## 150 3 Từ 10 - 20 triệu
## 151 2 Từ 10 - 20 triệu
## 152 1 Trên 20 triệu
## 153 3 Từ 10 - 20 triệu
## 154 1 Trên 20 triệu
## 155 3 Từ 10 - 20 triệu
## 156 1 Từ 10 - 20 triệu
## 157 2 Trên 20 triệu
## 158 3 Từ 10 - 20 triệu
## 159 1 Từ 10 - 20 triệu
## 160 2 Trên 20 triệu
## 161 3 Từ 10 - 20 triệu
## 162 3 < 10 triệu
## 163 1 Trên 20 triệu
## 164 3 Từ 10 - 20 triệu
## 165 2 Từ 10 - 20 triệu
## 166 2 Trên 20 triệu
## 167 3 Từ 10 - 20 triệu
## 168 3 < 10 triệu
## 169 2 Trên 20 triệu
## 170 4 Từ 10 - 20 triệu
## 171 2 < 10 triệu
## 172 4 Từ 10 - 20 triệu
## 173 4 Từ 10 - 20 triệu
## 174 1 Từ 10 - 20 triệu
## 175 1 Từ 10 - 20 triệu
## 176 3 Từ 10 - 20 triệu
## 177 3 < 10 triệu
## 178 1 Từ 10 - 20 triệu
## 179 3 < 10 triệu
## 180 2 Trên 20 triệu
## 181 2 Từ 10 - 20 triệu
## 182 2 Từ 10 - 20 triệu
## 183 4 Từ 10 - 20 triệu
## 184 4 Từ 10 - 20 triệu
## 185 3 Trên 20 triệu
## 186 1 Từ 10 - 20 triệu
## 187 3 Từ 10 - 20 triệu
## 188 1 Trên 20 triệu
## 189 2 Từ 10 - 20 triệu
## 190 3 < 10 triệu
## 191 1 Trên 20 triệu
## 192 1 Từ 10 - 20 triệu
## 193 2 Từ 10 - 20 triệu
## 194 1 Từ 10 - 20 triệu
## 195 3 Trên 20 triệu
## 196 2 Từ 10 - 20 triệu
## 197 4 Từ 10 - 20 triệu
## 198 4 Từ 10 - 20 triệu
# Xác định nhóm tuổi và thu nhập có số lượng lớn nhất
library(gmodels)
## Registered S3 method overwritten by 'gdata':
## method from
## reorder.factor DescTools
CrossTable(DTUOIMH,TNHAPMH)
##
##
## Cell Contents
## |-------------------------|
## | N |
## | Chi-square contribution |
## | N / Row Total |
## | N / Col Total |
## | N / Table Total |
## |-------------------------|
##
##
## Total Observations in Table: 198
##
##
## | TNHAPMH
## DTUOIMH | < 10 triệu | Trên 20 triệu | Từ 10 - 20 triệu | Row Total |
## ----------------|------------------|------------------|------------------|------------------|
## < 20 tuổi | 2 | 15 | 19 | 36 |
## | 0.284 | 0.178 | 0.021 | |
## | 0.056 | 0.417 | 0.528 | 0.182 |
## | 0.125 | 0.203 | 0.176 | |
## | 0.010 | 0.076 | 0.096 | |
## ----------------|------------------|------------------|------------------|------------------|
## Trên 50 tuổi | 1 | 14 | 23 | 38 |
## | 1.396 | 0.003 | 0.249 | |
## | 0.026 | 0.368 | 0.605 | 0.192 |
## | 0.062 | 0.189 | 0.213 | |
## | 0.005 | 0.071 | 0.116 | |
## ----------------|------------------|------------------|------------------|------------------|
## Từ 20 - 30 tuổi | 2 | 20 | 29 | 51 |
## | 1.092 | 0.046 | 0.050 | |
## | 0.039 | 0.392 | 0.569 | 0.258 |
## | 0.125 | 0.270 | 0.269 | |
## | 0.010 | 0.101 | 0.146 | |
## ----------------|------------------|------------------|------------------|------------------|
## Từ 31 - 50 tuổi | 11 | 25 | 37 | 73 |
## | 4.411 | 0.191 | 0.199 | |
## | 0.151 | 0.342 | 0.507 | 0.369 |
## | 0.688 | 0.338 | 0.343 | |
## | 0.056 | 0.126 | 0.187 | |
## ----------------|------------------|------------------|------------------|------------------|
## Column Total | 16 | 74 | 108 | 198 |
## | 0.081 | 0.374 | 0.545 | |
## ----------------|------------------|------------------|------------------|------------------|
##
##
table(DTUOIMH,TNHAPMH)
## TNHAPMH
## DTUOIMH < 10 triệu Trên 20 triệu Từ 10 - 20 triệu
## < 20 tuổi 2 15 19
## Trên 50 tuổi 1 14 23
## Từ 20 - 30 tuổi 2 20 29
## Từ 31 - 50 tuổi 11 25 37
table(DTUOIMH,TNHAPMH)/sum(table(DTUOIMH,TNHAPMH))
## TNHAPMH
## DTUOIMH < 10 triệu Trên 20 triệu Từ 10 - 20 triệu
## < 20 tuổi 0.010101010 0.075757576 0.095959596
## Trên 50 tuổi 0.005050505 0.070707071 0.116161616
## Từ 20 - 30 tuổi 0.010101010 0.101010101 0.146464646
## Từ 31 - 50 tuổi 0.055555556 0.126262626 0.186868687
### NHẬN XÉT CÂU 3: Nhóm tuổi và nhóm thu nhập có số lượng quan sát lớn nhất
# là từ 31 - 50 tuổi với thu nhập từ 10 - 20 triệu; có 37 trường hợp/đối tượng (18,7%)
### CÂU 4
library(gmodels)
CrossTable(GTINHMH,DTUOIMH)
##
##
## Cell Contents
## |-------------------------|
## | N |
## | Chi-square contribution |
## | N / Row Total |
## | N / Col Total |
## | N / Table Total |
## |-------------------------|
##
##
## Total Observations in Table: 198
##
##
## | DTUOIMH
## GTINHMH | < 20 tuổi | Trên 50 tuổi | Từ 20 - 30 tuổi | Từ 31 - 50 tuổi | Row Total |
## -------------|-----------------|-----------------|-----------------|-----------------|-----------------|
## Nam | 24 | 21 | 35 | 41 | 121 |
## | 0.182 | 0.213 | 0.471 | 0.292 | |
## | 0.198 | 0.174 | 0.289 | 0.339 | 0.611 |
## | 0.667 | 0.553 | 0.686 | 0.562 | |
## | 0.121 | 0.106 | 0.177 | 0.207 | |
## -------------|-----------------|-----------------|-----------------|-----------------|-----------------|
## NỮ | 12 | 17 | 16 | 32 | 77 |
## | 0.286 | 0.334 | 0.741 | 0.459 | |
## | 0.156 | 0.221 | 0.208 | 0.416 | 0.389 |
## | 0.333 | 0.447 | 0.314 | 0.438 | |
## | 0.061 | 0.086 | 0.081 | 0.162 | |
## -------------|-----------------|-----------------|-----------------|-----------------|-----------------|
## Column Total | 36 | 38 | 51 | 73 | 198 |
## | 0.182 | 0.192 | 0.258 | 0.369 | |
## -------------|-----------------|-----------------|-----------------|-----------------|-----------------|
##
##
### NHẬN XÉT CÂU 4:Trong tất cả các nhóm tuổi thì giới tính nam đều chiếm tỷ trọng cao hơn
#mã hoá dữ liệu giới tính
GTINHMH=GTINH
GTINHMH=replace(GTINHMH,GTINH==1,"Nam")
GTINHMH=replace(GTINHMH,GTINH==2,"Nu")
GTcode=data.frame(GTINH,GTINHMH)
GTcode
## GTINH GTINHMH
## 1 2 Nu
## 2 2 Nu
## 3 2 Nu
## 4 1 Nam
## 5 2 Nu
## 6 1 Nam
## 7 1 Nam
## 8 1 Nam
## 9 2 Nu
## 10 1 Nam
## 11 1 Nam
## 12 1 Nam
## 13 2 Nu
## 14 1 Nam
## 15 2 Nu
## 16 1 Nam
## 17 1 Nam
## 18 1 Nam
## 19 2 Nu
## 20 1 Nam
## 21 1 Nam
## 22 1 Nam
## 23 1 Nam
## 24 1 Nam
## 25 1 Nam
## 26 1 Nam
## 27 1 Nam
## 28 1 Nam
## 29 1 Nam
## 30 2 Nu
## 31 2 Nu
## 32 2 Nu
## 33 1 Nam
## 34 2 Nu
## 35 2 Nu
## 36 2 Nu
## 37 1 Nam
## 38 2 Nu
## 39 1 Nam
## 40 1 Nam
## 41 1 Nam
## 42 2 Nu
## 43 1 Nam
## 44 2 Nu
## 45 1 Nam
## 46 1 Nam
## 47 2 Nu
## 48 1 Nam
## 49 2 Nu
## 50 2 Nu
## 51 1 Nam
## 52 2 Nu
## 53 2 Nu
## 54 1 Nam
## 55 2 Nu
## 56 2 Nu
## 57 1 Nam
## 58 2 Nu
## 59 2 Nu
## 60 2 Nu
## 61 2 Nu
## 62 2 Nu
## 63 2 Nu
## 64 2 Nu
## 65 2 Nu
## 66 1 Nam
## 67 1 Nam
## 68 2 Nu
## 69 2 Nu
## 70 2 Nu
## 71 2 Nu
## 72 2 Nu
## 73 2 Nu
## 74 1 Nam
## 75 2 Nu
## 76 2 Nu
## 77 2 Nu
## 78 2 Nu
## 79 1 Nam
## 80 2 Nu
## 81 2 Nu
## 82 2 Nu
## 83 1 Nam
## 84 2 Nu
## 85 1 Nam
## 86 1 Nam
## 87 1 Nam
## 88 2 Nu
## 89 1 Nam
## 90 2 Nu
## 91 1 Nam
## 92 1 Nam
## 93 2 Nu
## 94 1 Nam
## 95 2 Nu
## 96 2 Nu
## 97 1 Nam
## 98 2 Nu
## 99 2 Nu
## 100 1 Nam
## 101 2 Nu
## 102 2 Nu
## 103 1 Nam
## 104 1 Nam
## 105 1 Nam
## 106 1 Nam
## 107 1 Nam
## 108 1 Nam
## 109 2 Nu
## 110 1 Nam
## 111 2 Nu
## 112 1 Nam
## 113 2 Nu
## 114 2 Nu
## 115 1 Nam
## 116 2 Nu
## 117 1 Nam
## 118 1 Nam
## 119 1 Nam
## 120 1 Nam
## 121 1 Nam
## 122 1 Nam
## 123 1 Nam
## 124 1 Nam
## 125 1 Nam
## 126 1 Nam
## 127 1 Nam
## 128 1 Nam
## 129 1 Nam
## 130 1 Nam
## 131 2 Nu
## 132 2 Nu
## 133 1 Nam
## 134 1 Nam
## 135 2 Nu
## 136 2 Nu
## 137 1 Nam
## 138 2 Nu
## 139 1 Nam
## 140 1 Nam
## 141 2 Nu
## 142 1 Nam
## 143 2 Nu
## 144 1 Nam
## 145 1 Nam
## 146 1 Nam
## 147 1 Nam
## 148 1 Nam
## 149 1 Nam
## 150 1 Nam
## 151 1 Nam
## 152 1 Nam
## 153 1 Nam
## 154 1 Nam
## 155 1 Nam
## 156 1 Nam
## 157 1 Nam
## 158 1 Nam
## 159 1 Nam
## 160 1 Nam
## 161 1 Nam
## 162 1 Nam
## 163 1 Nam
## 164 2 Nu
## 165 2 Nu
## 166 1 Nam
## 167 2 Nu
## 168 1 Nam
## 169 1 Nam
## 170 1 Nam
## 171 1 Nam
## 172 2 Nu
## 173 1 Nam
## 174 2 Nu
## 175 1 Nam
## 176 2 Nu
## 177 1 Nam
## 178 1 Nam
## 179 1 Nam
## 180 2 Nu
## 181 1 Nam
## 182 1 Nam
## 183 2 Nu
## 184 2 Nu
## 185 1 Nam
## 186 1 Nam
## 187 1 Nam
## 188 2 Nu
## 189 1 Nam
## 190 1 Nam
## 191 1 Nam
## 192 1 Nam
## 193 1 Nam
## 194 1 Nam
## 195 1 Nam
## 196 1 Nam
## 197 1 Nam
## 198 1 Nam
### Câu 5
#biểu đồ tròn Pie
table(GTINHMH)
## GTINHMH
## Nam Nu
## 121 77
pie(table(GTINHMH),col=rainbow(2),main="Biểu đồ tròn của giới tính")
pie(table(GTINHMH),col=c("green","orange"),main="Biểu đồ tròn của giới tính")
#mã hoá dữ liệu độ tuổi
DTUOIMH=DTUOI
DTUOIMH=replace(DTUOIMH,DTUOI==1,"< 20 tuoi")
DTUOIMH=replace(DTUOIMH,DTUOI==2,"Tu 20 - 30 tuoi")
DTUOIMH=replace(DTUOIMH,DTUOI==3,"Tu 31 - 50 tuoi")
DTUOIMH=replace(DTUOIMH,DTUOI==4,"Tren 50 tuoi")
DTcode=data.frame(DTUOI,DTUOIMH)
DTcode
## DTUOI DTUOIMH
## 1 1 < 20 tuoi
## 2 2 Tu 20 - 30 tuoi
## 3 1 < 20 tuoi
## 4 2 Tu 20 - 30 tuoi
## 5 1 < 20 tuoi
## 6 2 Tu 20 - 30 tuoi
## 7 1 < 20 tuoi
## 8 1 < 20 tuoi
## 9 2 Tu 20 - 30 tuoi
## 10 1 < 20 tuoi
## 11 2 Tu 20 - 30 tuoi
## 12 1 < 20 tuoi
## 13 2 Tu 20 - 30 tuoi
## 14 1 < 20 tuoi
## 15 4 Tren 50 tuoi
## 16 3 Tu 31 - 50 tuoi
## 17 2 Tu 20 - 30 tuoi
## 18 1 < 20 tuoi
## 19 3 Tu 31 - 50 tuoi
## 20 4 Tren 50 tuoi
## 21 2 Tu 20 - 30 tuoi
## 22 3 Tu 31 - 50 tuoi
## 23 3 Tu 31 - 50 tuoi
## 24 4 Tren 50 tuoi
## 25 3 Tu 31 - 50 tuoi
## 26 2 Tu 20 - 30 tuoi
## 27 3 Tu 31 - 50 tuoi
## 28 1 < 20 tuoi
## 29 3 Tu 31 - 50 tuoi
## 30 3 Tu 31 - 50 tuoi
## 31 2 Tu 20 - 30 tuoi
## 32 3 Tu 31 - 50 tuoi
## 33 3 Tu 31 - 50 tuoi
## 34 4 Tren 50 tuoi
## 35 3 Tu 31 - 50 tuoi
## 36 2 Tu 20 - 30 tuoi
## 37 3 Tu 31 - 50 tuoi
## 38 1 < 20 tuoi
## 39 3 Tu 31 - 50 tuoi
## 40 3 Tu 31 - 50 tuoi
## 41 2 Tu 20 - 30 tuoi
## 42 1 < 20 tuoi
## 43 3 Tu 31 - 50 tuoi
## 44 1 < 20 tuoi
## 45 4 Tren 50 tuoi
## 46 2 Tu 20 - 30 tuoi
## 47 3 Tu 31 - 50 tuoi
## 48 1 < 20 tuoi
## 49 3 Tu 31 - 50 tuoi
## 50 4 Tren 50 tuoi
## 51 1 < 20 tuoi
## 52 2 Tu 20 - 30 tuoi
## 53 1 < 20 tuoi
## 54 3 Tu 31 - 50 tuoi
## 55 4 Tren 50 tuoi
## 56 4 Tren 50 tuoi
## 57 4 Tren 50 tuoi
## 58 3 Tu 31 - 50 tuoi
## 59 3 Tu 31 - 50 tuoi
## 60 4 Tren 50 tuoi
## 61 2 Tu 20 - 30 tuoi
## 62 4 Tren 50 tuoi
## 63 3 Tu 31 - 50 tuoi
## 64 3 Tu 31 - 50 tuoi
## 65 3 Tu 31 - 50 tuoi
## 66 3 Tu 31 - 50 tuoi
## 67 3 Tu 31 - 50 tuoi
## 68 4 Tren 50 tuoi
## 69 3 Tu 31 - 50 tuoi
## 70 4 Tren 50 tuoi
## 71 2 Tu 20 - 30 tuoi
## 72 3 Tu 31 - 50 tuoi
## 73 3 Tu 31 - 50 tuoi
## 74 4 Tren 50 tuoi
## 75 3 Tu 31 - 50 tuoi
## 76 2 Tu 20 - 30 tuoi
## 77 4 Tren 50 tuoi
## 78 4 Tren 50 tuoi
## 79 4 Tren 50 tuoi
## 80 3 Tu 31 - 50 tuoi
## 81 1 < 20 tuoi
## 82 2 Tu 20 - 30 tuoi
## 83 3 Tu 31 - 50 tuoi
## 84 3 Tu 31 - 50 tuoi
## 85 2 Tu 20 - 30 tuoi
## 86 3 Tu 31 - 50 tuoi
## 87 2 Tu 20 - 30 tuoi
## 88 1 < 20 tuoi
## 89 1 < 20 tuoi
## 90 3 Tu 31 - 50 tuoi
## 91 3 Tu 31 - 50 tuoi
## 92 1 < 20 tuoi
## 93 2 Tu 20 - 30 tuoi
## 94 3 Tu 31 - 50 tuoi
## 95 1 < 20 tuoi
## 96 3 Tu 31 - 50 tuoi
## 97 2 Tu 20 - 30 tuoi
## 98 4 Tren 50 tuoi
## 99 4 Tren 50 tuoi
## 100 2 Tu 20 - 30 tuoi
## 101 2 Tu 20 - 30 tuoi
## 102 2 Tu 20 - 30 tuoi
## 103 1 < 20 tuoi
## 104 2 Tu 20 - 30 tuoi
## 105 2 Tu 20 - 30 tuoi
## 106 1 < 20 tuoi
## 107 2 Tu 20 - 30 tuoi
## 108 2 Tu 20 - 30 tuoi
## 109 2 Tu 20 - 30 tuoi
## 110 4 Tren 50 tuoi
## 111 3 Tu 31 - 50 tuoi
## 112 2 Tu 20 - 30 tuoi
## 113 3 Tu 31 - 50 tuoi
## 114 3 Tu 31 - 50 tuoi
## 115 4 Tren 50 tuoi
## 116 3 Tu 31 - 50 tuoi
## 117 2 Tu 20 - 30 tuoi
## 118 3 Tu 31 - 50 tuoi
## 119 4 Tren 50 tuoi
## 120 3 Tu 31 - 50 tuoi
## 121 4 Tren 50 tuoi
## 122 3 Tu 31 - 50 tuoi
## 123 2 Tu 20 - 30 tuoi
## 124 3 Tu 31 - 50 tuoi
## 125 4 Tren 50 tuoi
## 126 2 Tu 20 - 30 tuoi
## 127 4 Tren 50 tuoi
## 128 3 Tu 31 - 50 tuoi
## 129 2 Tu 20 - 30 tuoi
## 130 4 Tren 50 tuoi
## 131 3 Tu 31 - 50 tuoi
## 132 3 Tu 31 - 50 tuoi
## 133 2 Tu 20 - 30 tuoi
## 134 3 Tu 31 - 50 tuoi
## 135 3 Tu 31 - 50 tuoi
## 136 4 Tren 50 tuoi
## 137 2 Tu 20 - 30 tuoi
## 138 3 Tu 31 - 50 tuoi
## 139 3 Tu 31 - 50 tuoi
## 140 4 Tren 50 tuoi
## 141 3 Tu 31 - 50 tuoi
## 142 4 Tren 50 tuoi
## 143 3 Tu 31 - 50 tuoi
## 144 2 Tu 20 - 30 tuoi
## 145 3 Tu 31 - 50 tuoi
## 146 3 Tu 31 - 50 tuoi
## 147 3 Tu 31 - 50 tuoi
## 148 4 Tren 50 tuoi
## 149 4 Tren 50 tuoi
## 150 3 Tu 31 - 50 tuoi
## 151 2 Tu 20 - 30 tuoi
## 152 1 < 20 tuoi
## 153 3 Tu 31 - 50 tuoi
## 154 1 < 20 tuoi
## 155 3 Tu 31 - 50 tuoi
## 156 1 < 20 tuoi
## 157 2 Tu 20 - 30 tuoi
## 158 3 Tu 31 - 50 tuoi
## 159 1 < 20 tuoi
## 160 2 Tu 20 - 30 tuoi
## 161 3 Tu 31 - 50 tuoi
## 162 3 Tu 31 - 50 tuoi
## 163 1 < 20 tuoi
## 164 3 Tu 31 - 50 tuoi
## 165 2 Tu 20 - 30 tuoi
## 166 2 Tu 20 - 30 tuoi
## 167 3 Tu 31 - 50 tuoi
## 168 3 Tu 31 - 50 tuoi
## 169 2 Tu 20 - 30 tuoi
## 170 4 Tren 50 tuoi
## 171 2 Tu 20 - 30 tuoi
## 172 4 Tren 50 tuoi
## 173 4 Tren 50 tuoi
## 174 1 < 20 tuoi
## 175 1 < 20 tuoi
## 176 3 Tu 31 - 50 tuoi
## 177 3 Tu 31 - 50 tuoi
## 178 1 < 20 tuoi
## 179 3 Tu 31 - 50 tuoi
## 180 2 Tu 20 - 30 tuoi
## 181 2 Tu 20 - 30 tuoi
## 182 2 Tu 20 - 30 tuoi
## 183 4 Tren 50 tuoi
## 184 4 Tren 50 tuoi
## 185 3 Tu 31 - 50 tuoi
## 186 1 < 20 tuoi
## 187 3 Tu 31 - 50 tuoi
## 188 1 < 20 tuoi
## 189 2 Tu 20 - 30 tuoi
## 190 3 Tu 31 - 50 tuoi
## 191 1 < 20 tuoi
## 192 1 < 20 tuoi
## 193 2 Tu 20 - 30 tuoi
## 194 1 < 20 tuoi
## 195 3 Tu 31 - 50 tuoi
## 196 2 Tu 20 - 30 tuoi
## 197 4 Tren 50 tuoi
## 198 4 Tren 50 tuoi
### Câu 6
#biểu đồ cột Ba
DTUOI.freq<-table(DTUOIMH)
DTUOI.freq
## DTUOIMH
## < 20 tuoi Tren 50 tuoi Tu 20 - 30 tuoi Tu 31 - 50 tuoi
## 36 38 51 73
barplot(DTUOI.freq,xlab="Nhóm tuổi",ylab="Số lượng",main="Biểu đồ cột của các nhóm tuổi",col=c("red","green","blue","purple"))
# thống kê các biến
library(psych)
ABCDEFtb=data.frame(Atb,Btb,Ctb,Dtb,Etb,Ftb)
describe(ABCDEFtb)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## Atb 1 198 4.06 0.72 4.00 4.08 1.19 2.20 5 2.80 -0.15 -1.14 0.05
## Btb 2 198 3.95 0.83 4.00 4.04 1.11 1.25 5 3.75 -0.89 0.37 0.06
## Ctb 3 198 3.68 0.87 3.62 3.71 0.93 1.50 5 3.50 -0.11 -0.61 0.06
## Dtb 4 198 3.80 0.78 4.00 3.85 0.49 1.00 5 4.00 -0.68 0.71 0.06
## Etb 5 198 3.95 0.81 4.00 4.04 0.74 1.25 5 3.75 -0.97 1.02 0.06
## Ftb 6 198 3.85 0.74 3.90 3.84 1.04 2.00 5 3.00 -0.04 -1.12 0.05
### Câu 7
#Biểu đồ Histogram
hist(Ftb)
hist(Ftb,col="green",xlab="Giá trị Ftb",border="yellow",ylab="Tần số Ftb",main="Đồ thị Histogram Ftb",xlim=c(2,5),ylim=c(0,1))
lines(density(na.omit(Ftb)),col="purple",lwd=2)
#Nhận xét
# Về hình dạng: Biến Ftb cguwa có hình dạng của hình chuông
# Sự phân bố của biến Ftb chưa đạt phân phối chuẩn
### Câu 8
#Biểu đồ Boxplot
boxplot(Ftb,ylab="Giá trị Ftb",col="lightgreen")
boxplot.stats(Ftb)
## $stats
## [1] 2.0 3.0 3.9 4.4 5.0
##
## $n
## [1] 198
##
## $conf
## [1] 3.7428 4.0572
##
## $out
## numeric(0)
#Nhận xét: Không có giá trị ngoại lai
### Câu 9
#Biểu đồ Scatter Plot giữa Atb và Ftb
plot(Ftb~Atb,xlab="Giá trị Atb",ylab="Giá trị Ftb",col="purple",pch=19)
plot(Atb,Ftb,xlab="Giá trị Atb",ylab="Giá trị Ftb",col="blue",pch=21)
abline(lm(Ftb~Atb))
#Nhận xét: Atb và Ftb cùng chiều, Atb càng tăng thì Ftb cũng tăng
#Mức độ tuyến tính của Atb và Ftb thấp
### Câu 10
library(psych)
fa=data.frame(A1,A2,A3,A4,A5)
alpha(fa)
##
## Reliability analysis
## Call: alpha(x = fa)
##
## raw_alpha std.alpha G6(smc) average_r S/N ase mean sd median_r
## 0.93 0.94 0.95 0.74 14 0.0083 4.1 0.72 0.78
##
## 95% confidence boundaries
## lower alpha upper
## Feldt 0.91 0.93 0.95
## Duhachek 0.92 0.93 0.95
##
## Reliability if an item is dropped:
## raw_alpha std.alpha G6(smc) average_r S/N alpha se var.r med.r
## A1 0.96 0.96 0.97 0.86 25.2 0.0049 0.0057 0.87
## A2 0.90 0.91 0.91 0.71 9.7 0.0125 0.0270 0.67
## A3 0.90 0.91 0.91 0.71 9.7 0.0125 0.0270 0.67
## A4 0.91 0.92 0.92 0.73 11.1 0.0111 0.0361 0.72
## A5 0.90 0.91 0.91 0.70 9.5 0.0121 0.0277 0.67
##
## Item statistics
## n raw.r std.r r.cor r.drop mean sd
## A1 198 0.74 0.73 0.60 0.60 4.1 0.89
## A2 198 0.94 0.94 0.95 0.90 4.0 0.83
## A3 198 0.94 0.94 0.95 0.90 4.0 0.83
## A4 198 0.89 0.90 0.89 0.84 4.1 0.74
## A5 198 0.94 0.94 0.95 0.90 4.0 0.73
##
## Non missing response frequency for each item
## 2 3 4 5 miss
## A1 0.00 0.33 0.20 0.47 0
## A2 0.02 0.30 0.35 0.33 0
## A3 0.02 0.30 0.35 0.33 0
## A4 0.01 0.21 0.45 0.33 0
## A5 0.01 0.24 0.48 0.28 0
#Nhận xét:
# Giá trị Cronbanh's Alpha = 0.93 > 0.6; như vậy nhân tố A đảm bảo độ tin cậy
# Hệ số tương quan biến tổng các biến/thang đo đều >0.3 (Đạt yêu cầu)
# => Nên không bất kì loại biến/thang đo nào
### Câu 11
library(psych)
abcde=data.frame(A1,A2,A3,A4,A5,B1,B2,B3,B4,C1,C2,C3,C4,D1,D2,D3,E1,E2,E3,E4)
KMO(abcde)
## Kaiser-Meyer-Olkin factor adequacy
## Call: KMO(r = abcde)
## Overall MSA = 0.67
## MSA for each item =
## A1 A2 A3 A4 A5 B1 B2 B3 B4 C1 C2 C3 C4 D1 D2 D3
## 0.77 0.70 0.61 0.75 0.68 0.78 0.81 0.53 0.56 0.58 0.53 0.86 0.84 0.72 0.74 0.88
## E1 E2 E3 E4
## 0.60 0.82 0.74 0.61
cortest.bartlett(abcde)
## R was not square, finding R from data
## $chisq
## [1] 5048.334
##
## $p.value
## [1] 0
##
## $df
## [1] 190
#Nhận xét:
# Giá trị KMO = 0.67 > 0.5; đảm bảo yêu cầu phân tích EFA
# Giá trị Sig. (P-Value) của Barlett's = 0,000 < 5%;
# Do đó kết quả EFA có thể sử dụng để phân tích
###Câu 12
fitabcde1<-principal(abcde,cor=TRUE,nfactors=19,rotate="none")
fitabcde1
## Principal Components Analysis
## Call: principal(r = abcde, nfactors = 19, rotate = "none", cor = TRUE)
## Standardized loadings (pattern matrix) based upon correlation matrix
## PC1 PC2 PC3 PC4 PC5 PC6 PC7 PC8 PC9 PC10 PC11 PC12
## A1 0.45 -0.26 0.51 0.05 0.21 -0.20 -0.25 -0.52 0.11 -0.04 0.16 0.03
## A2 0.64 -0.46 0.53 0.01 -0.06 0.09 0.09 0.07 0.02 -0.08 -0.13 -0.16
## A3 0.63 -0.47 0.53 0.03 -0.05 0.10 0.10 0.07 0.03 -0.10 -0.12 -0.16
## A4 0.56 -0.56 0.46 0.00 0.03 0.03 0.04 0.08 -0.07 0.13 0.10 0.27
## A5 0.65 -0.53 0.44 0.02 0.00 -0.06 0.11 0.16 -0.06 0.02 0.02 0.13
## B1 0.54 0.16 -0.43 0.61 0.06 0.24 -0.07 -0.02 0.08 0.00 -0.05 0.07
## B2 0.45 0.27 -0.27 0.44 0.13 -0.52 0.18 0.13 -0.24 0.03 0.18 -0.09
## B3 0.61 0.19 -0.46 0.59 0.07 0.04 0.01 0.00 0.00 -0.05 0.00 0.02
## B4 0.57 0.15 -0.41 0.61 0.08 0.24 -0.03 -0.08 0.04 0.00 -0.02 -0.03
## C1 0.51 0.59 0.03 -0.50 0.21 0.16 -0.09 0.13 -0.11 0.04 0.11 -0.11
## C2 0.44 0.52 0.09 -0.52 0.27 0.28 -0.13 0.06 -0.07 0.14 0.20 -0.01
## C3 0.38 0.64 -0.02 -0.27 0.23 -0.22 0.15 0.07 0.30 0.15 -0.27 0.22
## C4 0.48 0.57 0.08 -0.35 0.26 -0.17 0.09 -0.09 0.07 -0.26 -0.09 -0.13
## D1 0.40 0.44 -0.01 -0.13 -0.61 -0.08 -0.19 0.19 0.18 -0.29 0.20 0.12
## D2 0.53 0.33 0.08 -0.02 -0.59 -0.13 -0.11 -0.12 0.03 0.39 -0.07 -0.17
## D3 0.17 0.67 0.17 -0.15 -0.32 0.14 0.26 -0.30 -0.36 -0.11 -0.12 0.17
## E1 -0.29 0.49 0.66 0.43 0.08 -0.02 -0.15 0.08 -0.06 0.00 -0.04 0.02
## E2 -0.34 0.40 0.45 0.29 -0.05 0.16 0.50 -0.08 0.25 0.07 0.28 -0.06
## E3 -0.27 0.50 0.68 0.40 0.03 0.00 -0.12 0.10 -0.01 -0.04 -0.10 -0.05
## E4 -0.30 0.48 0.66 0.42 0.08 -0.01 -0.17 0.07 -0.05 0.03 -0.02 0.05
## PC13 PC14 PC15 PC16 PC17 PC18 PC19 h2 u2 com
## A1 -0.11 0.01 -0.01 0.00 0.01 -0.01 0.00 1 2.1e-06 5.2
## A2 -0.10 -0.05 0.03 -0.07 0.06 0.06 -0.01 1 4.4e-04 3.4
## A3 -0.10 0.04 0.05 0.09 -0.06 -0.02 0.02 1 1.1e-03 3.5
## A4 0.20 -0.06 0.07 0.06 0.05 0.00 -0.02 1 3.3e-06 4.1
## A5 0.04 0.06 -0.14 -0.07 -0.05 -0.04 0.01 1 1.4e-04 3.3
## B1 0.04 0.19 0.04 0.00 0.05 0.01 0.01 1 1.6e-06 3.8
## B2 -0.08 0.02 0.07 -0.01 0.00 0.00 0.00 1 1.3e-06 5.8
## B3 -0.02 -0.07 -0.14 0.04 -0.01 0.05 -0.03 1 6.1e-04 3.3
## B4 0.02 -0.13 0.06 -0.03 -0.05 -0.05 0.01 1 4.5e-04 3.5
## C1 -0.07 -0.04 -0.07 0.05 0.09 -0.03 0.04 1 4.8e-05 4.1
## C2 -0.04 0.04 0.05 -0.04 -0.07 0.03 -0.04 1 5.6e-05 5.0
## C3 -0.14 -0.02 0.02 0.00 0.00 0.00 0.00 1 2.1e-06 4.7
## C4 0.33 0.01 0.00 -0.01 -0.01 0.00 0.00 1 2.9e-07 5.1
## D1 -0.05 -0.01 0.02 0.00 0.00 0.00 0.00 1 2.8e-06 4.6
## D2 0.11 0.01 -0.02 0.00 0.00 0.00 0.00 1 9.2e-07 4.1
## D3 -0.07 0.01 0.01 0.00 0.00 0.00 0.00 1 4.9e-07 4.1
## E1 0.02 0.01 -0.01 0.07 -0.05 0.04 0.02 1 2.1e-03 3.4
## E2 0.01 0.00 -0.02 0.00 0.01 0.00 0.00 1 3.3e-07 6.2
## E3 -0.03 0.02 -0.01 0.00 0.04 -0.06 -0.07 1 1.0e-06 3.2
## E4 0.05 -0.05 0.02 -0.07 0.01 0.01 0.04 1 2.1e-03 3.5
##
## PC1 PC2 PC3 PC4 PC5 PC6 PC7 PC8 PC9 PC10 PC11
## SS loadings 4.59 4.23 3.45 2.65 1.17 0.71 0.62 0.54 0.43 0.41 0.39
## Proportion Var 0.23 0.21 0.17 0.13 0.06 0.04 0.03 0.03 0.02 0.02 0.02
## Cumulative Var 0.23 0.44 0.61 0.75 0.80 0.84 0.87 0.90 0.92 0.94 0.96
## Proportion Explained 0.23 0.21 0.17 0.13 0.06 0.04 0.03 0.03 0.02 0.02 0.02
## Cumulative Proportion 0.23 0.44 0.61 0.75 0.80 0.84 0.87 0.90 0.92 0.94 0.96
## PC12 PC13 PC14 PC15 PC16 PC17 PC18 PC19
## SS loadings 0.31 0.24 0.08 0.07 0.04 0.03 0.02 0.01
## Proportion Var 0.02 0.01 0.00 0.00 0.00 0.00 0.00 0.00
## Cumulative Var 0.98 0.99 0.99 0.99 1.00 1.00 1.00 1.00
## Proportion Explained 0.02 0.01 0.00 0.00 0.00 0.00 0.00 0.00
## Cumulative Proportion 0.98 0.99 0.99 0.99 1.00 1.00 1.00 1.00
##
## Mean item complexity = 4.2
## Test of the hypothesis that 19 components are sufficient.
##
## The root mean square of the residuals (RMSR) is 0
## with the empirical chi square 0 with prob < NA
##
## Fit based upon off diagonal values = 1
fitabcde2<-principal(abcde,cor=TRUE,nfactors=5,rotate="none")
fitabcde2
## Principal Components Analysis
## Call: principal(r = abcde, nfactors = 5, rotate = "none", cor = TRUE)
## Standardized loadings (pattern matrix) based upon correlation matrix
## PC1 PC2 PC3 PC4 PC5 h2 u2 com
## A1 0.45 -0.26 0.51 0.05 0.21 0.57 0.430 2.9
## A2 0.64 -0.46 0.53 0.01 -0.06 0.90 0.099 2.8
## A3 0.63 -0.47 0.53 0.03 -0.05 0.90 0.103 2.8
## A4 0.56 -0.56 0.46 0.00 0.03 0.83 0.169 2.9
## A5 0.65 -0.53 0.44 0.02 0.00 0.90 0.095 2.7
## B1 0.54 0.16 -0.43 0.61 0.06 0.88 0.121 3.0
## B2 0.45 0.27 -0.27 0.44 0.13 0.56 0.437 3.5
## B3 0.61 0.19 -0.46 0.59 0.07 0.96 0.036 3.1
## B4 0.57 0.15 -0.41 0.61 0.08 0.90 0.099 2.9
## C1 0.51 0.59 0.03 -0.50 0.21 0.89 0.110 3.2
## C2 0.44 0.52 0.09 -0.52 0.27 0.82 0.176 3.5
## C3 0.38 0.64 -0.02 -0.27 0.23 0.68 0.323 2.4
## C4 0.48 0.57 0.08 -0.35 0.26 0.75 0.251 3.1
## D1 0.40 0.44 -0.01 -0.13 -0.61 0.74 0.256 2.7
## D2 0.53 0.33 0.08 -0.02 -0.59 0.76 0.245 2.6
## D3 0.17 0.67 0.17 -0.15 -0.32 0.63 0.368 1.9
## E1 -0.29 0.49 0.66 0.43 0.08 0.95 0.047 3.1
## E2 -0.34 0.40 0.45 0.29 -0.05 0.56 0.440 3.7
## E3 -0.27 0.50 0.68 0.40 0.03 0.95 0.048 2.9
## E4 -0.30 0.48 0.66 0.42 0.08 0.94 0.055 3.1
##
## PC1 PC2 PC3 PC4 PC5
## SS loadings 4.59 4.23 3.45 2.65 1.17
## Proportion Var 0.23 0.21 0.17 0.13 0.06
## Cumulative Var 0.23 0.44 0.61 0.75 0.80
## Proportion Explained 0.29 0.26 0.21 0.16 0.07
## Cumulative Proportion 0.29 0.55 0.76 0.93 1.00
##
## Mean item complexity = 2.9
## Test of the hypothesis that 5 components are sufficient.
##
## The root mean square of the residuals (RMSR) is 0.03
## with the empirical chi square 71.75 with prob < 0.99
##
## Fit based upon off diagonal values = 0.99
fitabcde3<-principal(abcde,cor=TRUE,nfactors=5,rotate="none")
fitabcde3
## Principal Components Analysis
## Call: principal(r = abcde, nfactors = 5, rotate = "none", cor = TRUE)
## Standardized loadings (pattern matrix) based upon correlation matrix
## PC1 PC2 PC3 PC4 PC5 h2 u2 com
## A1 0.45 -0.26 0.51 0.05 0.21 0.57 0.430 2.9
## A2 0.64 -0.46 0.53 0.01 -0.06 0.90 0.099 2.8
## A3 0.63 -0.47 0.53 0.03 -0.05 0.90 0.103 2.8
## A4 0.56 -0.56 0.46 0.00 0.03 0.83 0.169 2.9
## A5 0.65 -0.53 0.44 0.02 0.00 0.90 0.095 2.7
## B1 0.54 0.16 -0.43 0.61 0.06 0.88 0.121 3.0
## B2 0.45 0.27 -0.27 0.44 0.13 0.56 0.437 3.5
## B3 0.61 0.19 -0.46 0.59 0.07 0.96 0.036 3.1
## B4 0.57 0.15 -0.41 0.61 0.08 0.90 0.099 2.9
## C1 0.51 0.59 0.03 -0.50 0.21 0.89 0.110 3.2
## C2 0.44 0.52 0.09 -0.52 0.27 0.82 0.176 3.5
## C3 0.38 0.64 -0.02 -0.27 0.23 0.68 0.323 2.4
## C4 0.48 0.57 0.08 -0.35 0.26 0.75 0.251 3.1
## D1 0.40 0.44 -0.01 -0.13 -0.61 0.74 0.256 2.7
## D2 0.53 0.33 0.08 -0.02 -0.59 0.76 0.245 2.6
## D3 0.17 0.67 0.17 -0.15 -0.32 0.63 0.368 1.9
## E1 -0.29 0.49 0.66 0.43 0.08 0.95 0.047 3.1
## E2 -0.34 0.40 0.45 0.29 -0.05 0.56 0.440 3.7
## E3 -0.27 0.50 0.68 0.40 0.03 0.95 0.048 2.9
## E4 -0.30 0.48 0.66 0.42 0.08 0.94 0.055 3.1
##
## PC1 PC2 PC3 PC4 PC5
## SS loadings 4.59 4.23 3.45 2.65 1.17
## Proportion Var 0.23 0.21 0.17 0.13 0.06
## Cumulative Var 0.23 0.44 0.61 0.75 0.80
## Proportion Explained 0.29 0.26 0.21 0.16 0.07
## Cumulative Proportion 0.29 0.55 0.76 0.93 1.00
##
## Mean item complexity = 2.9
## Test of the hypothesis that 5 components are sufficient.
##
## The root mean square of the residuals (RMSR) is 0.03
## with the empirical chi square 71.75 with prob < 0.99
##
## Fit based upon off diagonal values = 0.99
#Nhận xét:
# Số nhân tố được trích : 5 nhóm nhân tố
# Giá trị Eigenvalue = 1.83 > 1 (đạt yêu cầu)
# Giá trị tổng phương sai trích = 0.80 = 80% > 50% (đạt yêu cầu)
# Cấu trúc nhân tố có phù hợp với dữ liệu ban đầu vì các biến/thang đo sắp xếp theo đúng nhóm nhân tố
###Câu13
library(psych)
vars=cbind(Atb,Btb,Ctb,Dtb,Etb,Ftb)
pairs.panels(vars)
library(Hmisc)
##
## Attaching package: 'Hmisc'
## The following objects are masked from 'package:DescTools':
##
## %nin%, Label, Mean, Quantile
## The following object is masked from 'package:psych':
##
## describe
## The following objects are masked from 'package:base':
##
## format.pval, units
result=rcorr(as.matrix(vars),type="pearson")
result$r
## Atb Btb Ctb Dtb Etb Ftb
## Atb 1.00000000 0.06221638 0.03185639 0.04014488 -0.08600054 0.2222299
## Btb 0.06221638 1.00000000 0.16097635 0.22801725 -0.10776531 0.2511952
## Ctb 0.03185639 0.16097635 1.00000000 0.50335736 0.02399959 0.5915879
## Dtb 0.04014488 0.22801725 0.50335736 1.00000000 0.13044247 0.5733869
## Etb -0.08600054 -0.10776531 0.02399959 0.13044247 1.00000000 0.2096923
## Ftb 0.22222985 0.25119517 0.59158791 0.57338686 0.20969230 1.0000000
result$P
## Atb Btb Ctb Dtb Etb Ftb
## Atb NA 0.3838835582 6.559338e-01 5.744306e-01 0.228316785 0.0016512737
## Btb 0.383883558 NA 2.347636e-02 1.234270e-03 0.130738634 0.0003574484
## Ctb 0.655933822 0.0234763641 NA 4.085621e-14 0.737161764 0.0000000000
## Dtb 0.574430580 0.0012342702 4.085621e-14 NA 0.066996016 0.0000000000
## Etb 0.228316785 0.1307386340 7.371618e-01 6.699602e-02 NA 0.0030269715
## Ftb 0.001651274 0.0003574484 0.000000e+00 0.000000e+00 0.003026971 NA
#Nhận xét:
# Biến Ftb có mối tương quan dương với Atb, Btb, Ctb ,Dtb ,Etb vì hệ số tương quan pearson của các mối quan hệ này >0
# Mối tương quan giữa Ftb và các biến Atb, Btb, Ctb ,Dtb, Etb có ý nghĩa thống kê vì giá trị Sig, (P-Value của các mối quan hệ này < 5%)