# Nhap lieu
likert=file.choose()
datalikert=read.csv(likert,header=TRUE)
save(datalikert,file="datalikert.rda")
attach(datalikert)
is.data.frame(datalikert)
## [1] TRUE
#cau 1
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
# Nhận xét :
#biến có giá trị trung bình lớn nhất là Atb 4.06
#biến có độ phân giải lớn nhất là Ctb
#tính Mode (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
#nhận xét : Mode của A1 là 5 và 94 lần xuất hiện
#cau 1
library(psych)
ABCDEFtb=data.frame(Atb,Btb,Ctb,Dtb,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
## Dtb.1 5 198 3.80 0.78 4.00 3.85 0.49 1.00 5 4.00 -0.68 0.71 0.06
## Etb 6 198 3.95 0.81 4.00 4.04 0.74 1.25 5 3.75 -0.97 1.02 0.06
## Ftb 7 198 3.85 0.74 3.90 3.84 1.04 2.00 5 3.00 -0.04 -1.12 0.05
#Nhận xét :
#biến có giá trị trung bình cao nhất là Atb 4.06
#biến có độ phân giải lớn nhất là Ctb
#tính Mode (AI)
library(DescTools)
Mode (A1)
## [1] 5
## attr(,"freq")
## [1] 94
#bài 5
#Mã hóa giới tính
GTINHMH=GTINH
GTINHMH=replace(GTINH,GTINH==1,"Nam")
GTINHMH=replace(GTINH,GTINH==2,"Nữ")
GTcode=data.frame(GTINH,GTINHMH)
GTcode
## GTINH GTINHMH
## 1 2 Nữ
## 2 2 Nữ
## 3 2 Nữ
## 4 1 1
## 5 2 Nữ
## 6 1 1
## 7 1 1
## 8 1 1
## 9 2 Nữ
## 10 1 1
## 11 1 1
## 12 1 1
## 13 2 Nữ
## 14 1 1
## 15 2 Nữ
## 16 1 1
## 17 1 1
## 18 1 1
## 19 2 Nữ
## 20 1 1
## 21 1 1
## 22 1 1
## 23 1 1
## 24 1 1
## 25 1 1
## 26 1 1
## 27 1 1
## 28 1 1
## 29 1 1
## 30 2 Nữ
## 31 2 Nữ
## 32 2 Nữ
## 33 1 1
## 34 2 Nữ
## 35 2 Nữ
## 36 2 Nữ
## 37 1 1
## 38 2 Nữ
## 39 1 1
## 40 1 1
## 41 1 1
## 42 2 Nữ
## 43 1 1
## 44 2 Nữ
## 45 1 1
## 46 1 1
## 47 2 Nữ
## 48 1 1
## 49 2 Nữ
## 50 2 Nữ
## 51 1 1
## 52 2 Nữ
## 53 2 Nữ
## 54 1 1
## 55 2 Nữ
## 56 2 Nữ
## 57 1 1
## 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 1
## 67 1 1
## 68 2 Nữ
## 69 2 Nữ
## 70 2 Nữ
## 71 2 Nữ
## 72 2 Nữ
## 73 2 Nữ
## 74 1 1
## 75 2 Nữ
## 76 2 Nữ
## 77 2 Nữ
## 78 2 Nữ
## 79 1 1
## 80 2 Nữ
## 81 2 Nữ
## 82 2 Nữ
## 83 1 1
## 84 2 Nữ
## 85 1 1
## 86 1 1
## 87 1 1
## 88 2 Nữ
## 89 1 1
## 90 2 Nữ
## 91 1 1
## 92 1 1
## 93 2 Nữ
## 94 1 1
## 95 2 Nữ
## 96 2 Nữ
## 97 1 1
## 98 2 Nữ
## 99 2 Nữ
## 100 1 1
## 101 2 Nữ
## 102 2 Nữ
## 103 1 1
## 104 1 1
## 105 1 1
## 106 1 1
## 107 1 1
## 108 1 1
## 109 2 Nữ
## 110 1 1
## 111 2 Nữ
## 112 1 1
## 113 2 Nữ
## 114 2 Nữ
## 115 1 1
## 116 2 Nữ
## 117 1 1
## 118 1 1
## 119 1 1
## 120 1 1
## 121 1 1
## 122 1 1
## 123 1 1
## 124 1 1
## 125 1 1
## 126 1 1
## 127 1 1
## 128 1 1
## 129 1 1
## 130 1 1
## 131 2 Nữ
## 132 2 Nữ
## 133 1 1
## 134 1 1
## 135 2 Nữ
## 136 2 Nữ
## 137 1 1
## 138 2 Nữ
## 139 1 1
## 140 1 1
## 141 2 Nữ
## 142 1 1
## 143 2 Nữ
## 144 1 1
## 145 1 1
## 146 1 1
## 147 1 1
## 148 1 1
## 149 1 1
## 150 1 1
## 151 1 1
## 152 1 1
## 153 1 1
## 154 1 1
## 155 1 1
## 156 1 1
## 157 1 1
## 158 1 1
## 159 1 1
## 160 1 1
## 161 1 1
## 162 1 1
## 163 1 1
## 164 2 Nữ
## 165 2 Nữ
## 166 1 1
## 167 2 Nữ
## 168 1 1
## 169 1 1
## 170 1 1
## 171 1 1
## 172 2 Nữ
## 173 1 1
## 174 2 Nữ
## 175 1 1
## 176 2 Nữ
## 177 1 1
## 178 1 1
## 179 1 1
## 180 2 Nữ
## 181 1 1
## 182 1 1
## 183 2 Nữ
## 184 2 Nữ
## 185 1 1
## 186 1 1
## 187 1 1
## 188 2 Nữ
## 189 1 1
## 190 1 1
## 191 1 1
## 192 1 1
## 193 1 1
## 194 1 1
## 195 1 1
## 196 1 1
## 197 1 1
## 198 1 1
pie(table(GTINHMH),col=c("lightpink","lightyellow"),main="Đồ thị tỷ lệ giới tính")
#bài 3
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")
DTMH=data.frame(DTUOI,DTUOIMH)
DTMH
## 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
#bài 6
DTUOI.freq<-table(DTUOIMH)
DTUOI.freq
## DTUOIMH
## < 20 tuổi trên 50 tuổi từ 20-30 tuổi từ 31-50 tuổi
## 36 38 51 73
barplot(DTUOI.freq,xalb="Nhóm tuổi",ylab="số lượng",col=c("green","blue","red","yellow"),main="Biểu đồ cột của các nhóm tuổi")
## Warning in plot.window(xlim, ylim, log = log, ...): "xalb" is not a graphical
## parameter
## Warning in axis(if (horiz) 2 else 1, at = at.l, labels = names.arg, lty =
## axis.lty, : "xalb" is not a graphical parameter
## Warning in title(main = main, sub = sub, xlab = xlab, ylab = ylab, ...): "xalb"
## is not a graphical parameter
## Warning in axis(if (horiz) 1 else 2, cex.axis = cex.axis, ...): "xalb" is not a
## graphical parameter
#bài 7
hist(Ftb,col="blue",xlab="giá trị Ftb",border="green",ylab="Tần số Ftb",main="Đồ thị Histogram Ftb",xlim=c(2,5),ylim=c(0,1))
lines(density(na.omit(Ftb)),col="green",lwd=2)
# nhận xét :phân bố không đuề
#bài 8
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: dữ liệu không có biến ngoại lai
# bài 9
plot(Ftb~Atb,xlab="Giá trị Atb",ylab="Giá trị Ftb",col="green",pch=21)
abline(lm(Ftb~Atb))
