# 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 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 5
## Biểu đồ tròn GTINH
pie(table(GTINHMH),col=c("green","pink"),sub="Biểu đồ tròn của Giới tính")
### CÂU 6
# 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
## Biểu đồ cột DTUOI
DTUOI.freq<-table(DTUOIMH)
DTUOI.freq
## DTUOIMH
## < 20 tuổi Từ 20 - 30 tuổi Từ 31 - 50 tuổi Trên 50 tuổi
## 36 51 73 38
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","orange"))
### CÂu 7
#Mã hóa biến Ftb
MHFtrungbinh=Ftb
MHFtrungbinh=replace(MHFtrungbinh,Ftb<=2,1)
MHFtrungbinh=replace(MHFtrungbinh,Ftb>2&Ftb<=4,2)
MHFtrungbinh=replace(MHFtrungbinh,Ftb>4,3)
FTB=data.frame(Ftb,MHFtrungbinh)
FTB
## Ftb MHFtrungbinh
## 1 4.2 3
## 2 3.6 2
## 3 4.2 3
## 4 4.6 3
## 5 4.0 2
## 6 4.0 2
## 7 3.6 2
## 8 3.8 2
## 9 4.2 3
## 10 3.8 2
## 11 3.0 2
## 12 3.8 2
## 13 3.2 2
## 14 3.6 2
## 15 3.6 2
## 16 3.6 2
## 17 4.2 3
## 18 3.6 2
## 19 4.0 2
## 20 3.6 2
## 21 2.0 1
## 22 3.2 2
## 23 2.4 2
## 24 3.0 2
## 25 3.0 2
## 26 3.0 2
## 27 4.0 2
## 28 3.4 2
## 29 3.8 2
## 30 3.2 2
## 31 4.2 3
## 32 3.0 2
## 33 4.0 2
## 34 4.8 3
## 35 3.6 2
## 36 3.4 2
## 37 3.0 2
## 38 3.4 2
## 39 4.0 2
## 40 3.2 2
## 41 3.6 2
## 42 3.0 2
## 43 3.0 2
## 44 4.6 3
## 45 4.6 3
## 46 4.6 3
## 47 4.6 3
## 48 4.6 3
## 49 4.6 3
## 50 4.4 3
## 51 4.2 3
## 52 4.6 3
## 53 4.6 3
## 54 5.0 3
## 55 4.4 3
## 56 3.0 2
## 57 3.0 2
## 58 5.0 3
## 59 3.2 2
## 60 3.0 2
## 61 5.0 3
## 62 3.6 2
## 63 5.0 3
## 64 3.0 2
## 65 4.8 3
## 66 4.6 3
## 67 3.2 2
## 68 5.0 3
## 69 3.8 2
## 70 5.0 3
## 71 5.0 3
## 72 5.0 3
## 73 3.8 2
## 74 3.0 2
## 75 3.0 2
## 76 3.0 2
## 77 4.0 2
## 78 4.4 3
## 79 4.4 3
## 80 4.4 3
## 81 4.2 3
## 82 4.6 3
## 83 4.6 3
## 84 5.0 3
## 85 4.4 3
## 86 3.0 2
## 87 3.0 2
## 88 5.0 3
## 89 3.2 2
## 90 3.0 2
## 91 5.0 3
## 92 3.6 2
## 93 5.0 3
## 94 3.0 2
## 95 4.8 3
## 96 4.6 3
## 97 3.2 2
## 98 5.0 3
## 99 3.8 2
## 100 3.0 2
## 101 2.6 2
## 102 3.0 2
## 103 3.0 2
## 104 3.0 2
## 105 5.0 3
## 106 3.2 2
## 107 4.2 3
## 108 4.2 3
## 109 4.0 2
## 110 3.0 2
## 111 4.4 3
## 112 3.0 2
## 113 4.2 3
## 114 3.0 2
## 115 2.6 2
## 116 4.4 3
## 117 4.6 3
## 118 4.6 3
## 119 3.0 2
## 120 3.0 2
## 121 3.0 2
## 122 2.8 2
## 123 4.2 3
## 124 4.6 3
## 125 3.6 2
## 126 4.4 3
## 127 2.8 2
## 128 4.2 3
## 129 4.0 2
## 130 2.8 2
## 131 3.6 2
## 132 3.0 2
## 133 3.8 2
## 134 3.6 2
## 135 4.4 3
## 136 4.4 3
## 137 3.2 2
## 138 4.2 3
## 139 4.0 2
## 140 4.0 2
## 141 3.6 2
## 142 3.4 2
## 143 4.0 2
## 144 3.0 2
## 145 4.0 2
## 146 4.2 3
## 147 4.2 3
## 148 4.4 3
## 149 3.8 2
## 150 4.4 3
## 151 4.0 2
## 152 4.0 2
## 153 3.4 2
## 154 4.4 3
## 155 4.0 2
## 156 5.0 3
## 157 2.6 2
## 158 2.2 2
## 159 4.2 3
## 160 4.6 3
## 161 4.6 3
## 162 5.0 3
## 163 4.4 3
## 164 3.0 2
## 165 3.0 2
## 166 5.0 3
## 167 3.2 2
## 168 3.0 2
## 169 5.0 3
## 170 3.6 2
## 171 5.0 3
## 172 3.0 2
## 173 4.8 3
## 174 4.6 3
## 175 3.2 2
## 176 5.0 3
## 177 3.8 2
## 178 5.0 3
## 179 5.0 3
## 180 5.0 3
## 181 3.8 2
## 182 3.0 2
## 183 3.0 2
## 184 3.0 2
## 185 4.0 2
## 186 4.4 3
## 187 4.4 3
## 188 3.8 2
## 189 3.6 2
## 190 4.6 3
## 191 3.0 2
## 192 4.4 3
## 193 3.0 2
## 194 3.0 2
## 195 3.0 2
## 196 3.4 2
## 197 3.4 2
## 198 4.4 3
## Biểu đồ histogram của biến Ftb
hist(Ftb,col="lightblue",xlab="Giá trị Ftb",border="lightyellow",ylab="Tần số Ftb",main="Đồ thị Histogram Ftb",xlim=c(2,5),ylim=c(0,1))
lines(density(na.omit(Ftb)),col="lightgreen",lwd=2)
# Nhận xét: Ftb không phân bố chuẩn vì đường này không có hình dạng hình chuông
### CÂU 8
## Biểu đồ Boxplot của biến Ftb
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 đồ Boxplot giữa Atb và Ftb
plot(Ftb~Atb,xlab="Giá trị Atb",ylab="Giá trị Ftb",col="lightpink",pch=19)
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
# Nhận xét: Mức độ tuyến tính của Atb và Ftb thấp