setwd("D:/Dropbox/Nghien cuu dinh luong/Dao tao/Bai giang/Khoa hoc R")
data<-read.csv("Data.csv")
data$gioi[data$SEX==1]<-"Nam"
data$gioi[data$SEX==0]<-"Nu"
data$levl[data$Hvan==1]<-"THPT"
data$levl[data$Hvan==2]<-"DH"
data$levl[data$Hvan==3]<-"SDH"
hoc<-table(data$levl,data$gioi)
data[1:5,]
## stt SEX Hvan Income AGE VAY NUMP MARITAL.ST STYPE.H TRANO PP
## 1 1 1 1 6 29 1 2 1 2 0 3.8
## 2 2 0 2 5 48 0 2 1 2 1 3.2
## 3 3 1 1 5 32 1 2 1 2 0 3.6
## 4 4 1 2 4 35 1 2 1 3 1 3.0
## 5 5 0 3 11 58 0 4 1 3 1 3.0
## Name LV gioi levl
## 1 Nguy?n Xuân Phú 9 Nam THPT
## 2 Hoàng Trung Hi?u 10 Nu DH
## 3 ?inh Lê Th?ng 4 Nam THPT
## 4 Nguy?n Tr?ng Phú 8 Nam DH
## 5 Tr?n ??i D??ng 5 Nu SDH
Nguồn: nghiencuudinhluong.com
library(fBasics)
## Loading required package: timeDate
## Loading required package: timeSeries
##
## Rmetrics Package fBasics
## Analysing Markets and calculating Basic Statistics
## Copyright (C) 2005-2014 Rmetrics Association Zurich
## Educational Software for Financial Engineering and Computational Science
## Rmetrics is free software and comes with ABSOLUTELY NO WARRANTY.
## https://www.rmetrics.org --- Mail to: info@rmetrics.org
bangmau <- qualiPalette(3, "Set3")
barplot(hoc,col=bangmau,beside=T,horiz=T,las=1)
legend("right",legend=rownames(hoc),box.lty=1,cex=0.5,fill=bangmau)
# Vẽ với ggplot2
library(ggplot2)
ggplot(data,aes(gioi,fill=levl))+geom_bar(position="dodge")