#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
#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
#biểu đồ tròn Pie(câu 5)
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
# biểu đồ cột Bar(câu6)
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
#Biểu đồ Histogram(câu7)
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
#Biểu đồ Boxplot(câu8)
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)
#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))
```