#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))

```