條形圖和餅圖：可視化基于頻率的分類變量
- 條形圖
barplot()
mtcars[1:5,]
cyl.freq <- table(mtcars$cyl)
cyl.freq
barplot(cyl.freq)
- 根據(jù)其他分類變量進(jìn)一步分割每一條頻率
cyl.freq.matrix <- table(mtcars$am,mtcars$cyl)
cyl.freq.matrix
barplot(cyl.freq.matrix,beside=TRUE,horiz=TRUE,las=1,
main="Performance car counts\nby transmission and cylinders",
names.arg=c("V4","V6","V8"),legend.text=c("auto","manual"),
args.legend=list(x="bottomright"))
- ggplot2 繪制條形圖
library("ggplot2")
qplot(factor(mtcars$cyl),geom="bar")
qplot(factor(mtcars$cyl),geom="blank",fill=factor(mtcars$am),xlab="",
ylab="",main="Performance car counts\nby transmission and cylinders") +
geom_bar(position="dodge") + scale_x_discrete(labels=c("V4","V6","V8"))+
scale_y_continuous(breaks=seq(0,12,2)) +
scale_fill_grey(name="Trans.",labels=c("auto","manual")) +
theme_bw() + coord_flip()
- 餅圖
pie()
pie(table(mtcars$cyl),labels=c("V4","V6","V8"),
col=c("white","gray","black"),main="Performance cars by cylinders")
直方圖：可視化連續(xù)變量的分布
hist()
mtcars$hp
hist(mtcars$hp)
# 自主設(shè)置間距：
hist(mtcars$hp,breaks=seq(0,400,25),col="gray",main="Horsepower",xlab="HP")
abline(v=c(mean(mtcars$hp),median(mtcars$hp)),lty=c(2,3),lwd=2)
legend("topright",legend=c("mean HP","median HP"),lty=c(2,3),lwd=2)
- ggplot2繪制直方圖
library("ggplot2")
qplot(mtcars$hp)
qplot(mtcars$hp,geom="blank",main="Horsepower",xlab="HP") +
geom_histogram(color="black",fill="white",breaks=seq(0,400,25),closed="right") +
geom_vline(mapping=aes(xintercept=c(mean(mtcars$hp),median(mtcars$hp)),
linetype=factor(c("mean","median"))),show.legend=TRUE) +
scale_linetype_manual(values=c(2,3)) + labs(linetype="")
箱線圖：五位數(shù)概括法的可視化表示
- 獨(dú)立箱線圖
boxplot()， 低于下四分位數(shù)或高于上四分位數(shù)1.5陪四分位差（IQR）的觀測(cè)值定義為異常點(diǎn)猜丹，boxplot的rang參數(shù)可以控制此分類懂缕，默認(rèn)range=1.5
hist(quakes$mag)
boxplot(quakes$mag)
- 并列箱線圖
比較不同組的五數(shù)概況分布~
stations.fac <- cut(quakes$stations,breaks=c(0,50,100,150)) #定義分組因子
stations.fac[1:5]
boxplot(quakes$mag~stations.fac,xlab="# stations detected",ylab="Magnitude",col="gray")
- ggplot2繪制箱線圖
library("ggplot2")
qplot(stations.fac,quakes$mag,geom="boxplot",
xlab="# stations detected",ylab="Magnitude")
散點(diǎn)圖：識(shí)別兩個(gè)不同連續(xù)型數(shù)值變量的觀測(cè)值之間的關(guān)系
- 單一散點(diǎn)圖
plot(),points()
iris[1:5,]
plot(iris[,4],iris[,3],type="n",xlab="Petal Width (cm)",
ylab="Petal Length (cm)") # type="n"生成正確維數(shù)的繪圖區(qū)域
points(iris[iris$Species=="setosa",4],
iris[iris$Species=="setosa",3],pch=19,col="black") # 繪制物種為setosa的點(diǎn)圖
points(iris[iris$Species=="virginica",4],
iris[iris$Species=="virginica",3],pch=19,col="gray") # 繪制物種為virginica的點(diǎn)圖
points(iris[iris$Species=="versicolor",4],
iris[iris$Species=="versicolor",3],pch=1,col="black") # 繪制物種為versicolor的點(diǎn)圖
legend("topleft",legend=c("setosa","virginica","versicolor"),
col=c("black","gray","black"),pch=c(19,19,1)) # 設(shè)置圖例
- 上面散點(diǎn)圖的簡(jiǎn)化版：先設(shè)置向量指定每個(gè)觀測(cè)個(gè)體所需點(diǎn)的字符和顏色
iris_pch <- rep(19,nrow(iris))
iris_pch[iris$Species=="versicolor"] <- 1
iris_col <- rep("black",nrow(iris))
iris_col[iris$Species=="virginica"] <- "gray"
plot(iris[,4],iris[,3],col=iris_col,pch=iris_pch,
xlab="Petal Width (cm)",ylab="Petal Length (cm)")
- 散點(diǎn)圖矩陣
pairs()
pairs(iris[,1:4],pch=iris_pch,col=iris_col,cex=0.75)
- ggplot 來繪制散點(diǎn)圖
library("ggplot2")
qplot(iris[,4],iris[,3],xlab="Petal width",ylab="Petal length",
shape=iris$Species) +
scale_shape_manual(values=4:6) + labs(shape="Species")
- 使用GGally包的
ggpairs()繪制散點(diǎn)圖矩陣
library("GGally")
ggpairs(iris,mapping=aes(col=Species),axisLabels="internal")
