R语言基础绘图教程——第5章:直方图和柱状图

#读入数据
data = read.table("barplot.txt",header=T)
#绘制条形图，仔细喊下面没一行代码都生成一个图，看他们的差别会知道参数是干嘛的。
barplot(data[,2])
barplot(data[,2],names.arg = data[,1])
barplot(data[,2],names.arg = data[,1],main="条形图",xlab="分组",ylab="统计量")
barplot(data[,2],names.arg = data[,1],main="条形图",xlab="分组",ylab="统计量",col="blue")
barplot(data[,2],names.arg = data[,1],main="条形图",xlab="分组",ylab="统计量",col=c("grey","red","blue","orange","green")

barplot(data[,2])

# 添加横轴分组名称
barplot(data[,2],names.arg = data[,1])

# 添加横纵坐标的名称和图片名称
barplot(data[,2],names.arg = data[,1],main="条形图",xlab="分组",ylab="统计量")

# 更改条形颜色
barplot(data[,2],names.arg = data[,1],main="条形图",xlab="分组",ylab="统计量",col="blue")

为不同的组设置不同的颜色
barplot(data[,2],names.arg = data[,1],main="条形图",xlab="分组",ylab="统计量",col=c("grey","red","blue","orange","green"))

#转换数据
data2 = t(data[,c(2,3)])
#绘制柱状图
barplot(as.matrix(data2))

#调整图形
barplot(as.matrix(data2),
        names.arg = data[,1],main="条形图",xlab="分组",ylab="统计量",
        col=c("blue","red"),
        legend=c("Low Dose","High Dose"),
        ylim=c(0,50))
box() #边框

###------水平柱状图
par(las=2)#调整水平轴数字方向
barplot(as.matrix(data2),
        names.arg = data[,1],main="条形图",xlab="统计量",
        col=c("blue","red"),
        legend=c("Low Dose","High Dose"),
        xlim=c(0,50),
        horiz=TRUE,
        space=0.5)

box() #边框

#非堆积柱状图
barplot(as.matrix(data2),
        names.arg = data[,1],main="条形图",xlab="分组",ylab="统计量",
        col=c("blue","red"),
        legend=c("Low Dose","High Dose"),
        ylim=c(0,30),
        beside=TRUE)
box() #边框

ggplot(data,aes())+ geom_bar()

library(ggplot2)
csub <- read.table("barplot_using_ggplot2.txt", sep="\t", header=T)
head(csub)
ggplot(csub,aes(x=Year,y=Anomaly10y))+
  geom_bar(stat="identity",position="identity")

#positive and negative values with different color
csub$positive <- csub$Anomaly10y >= 0
ggplot(csub,aes(x=Year,y=Anomaly10y,fill=positive))+
  geom_bar(stat="identity",position="identity")

#change more params
ggplot(csub, aes(Year,Anomaly10y, fill=positive))+
  geom_bar(stat="identity", position="identity", colour="black", size=0.25) +
  scale_fill_manual(values=c("#CCEEFF", "#FFDDDD"), guide=F)

library(ggplot2)
#readin data
data = read.table("histogram1.txt", sep="\t", header=T)
#draw histogram 
ggplot(data, aes(waiting)) +
  geom_histogram()
#set binwidth 5
ggplot(data, aes(waiting)) +
  geom_histogram(binwidth=5, fill="white", color="black")
#how to set the binwidth to make the histogram have 15 groups?
binsize <- diff(range(data$waiting))/15
ggplot(data, aes(waiting)) +
  geom_histogram(binwidth=binsize, fill="white", color="black")

#set the value of "origin"
ggplot(data, aes(waiting)) +
  geom_histogram(binwidth=8, fill="white", color="black", origin=30)
ggplot(data, aes(waiting)) +
  geom_histogram(binwidth=8, fill="white", color="black", origin=40)

library(ggplot2)

#readin data
data = read.table("histogram2.txt", sep="\t", header=T)

#draw histogram with 
ggplot(data, aes(bwt)) +
  geom_histogram()

#draw histogram with in different subset
data$smoke = factor(data$smoke)
ggplot(data, aes(bwt, fill=smoke)) +
  geom_histogram(position="identity", alpha=0.4)