R语言学习笔记（十七）：高级绘图ggplot2

#ggplot2

library(ggplot2)
#散点图
ggplot(data=mtcars,aes(x=wt,y=mpg))+geom_point()+labs(title="Automation Data",x="Weight",y="Miles Per Gallon")


#散点图+回归曲线
ggplot(data=mtcars,aes(x=wt,y=mpg))+geom_point(pch=17,color="blue",size=2)+geom_smooth(method="lm",color="red",linetype=2)+labs(title="Automobile Data",x="Weight",y="Miles Per Gallon")



#分组---出错
mtcars$am<-factor(mtcars$am,levels=c(0,1),labels=c("Automation","Manual"))
mtcars$vs<-factor(mtcars$vs,levels=c(0,1),labels=c("V-Engine","Straight Engine"))
mtcars$cyl<-factor(mtcars$cyl)

ggplot(data=mtcars,aes=c(x=hp,y=mpg,shape=cyl,color=cyl))+geom_point(size=3)+facet_grid(am~vs)+labs(title="Automobile Data by Engine Type",x="Horsepower",y="Miles Per Gallon")


#直方图
data(singer,package="lattice")
ggplot(singer,aes(x=height))+geom_histogram()


#箱线图
ggplot(singer,aes(x=voice.part,y=height))+geom_boxplot()

#复杂的叠加图：箱线+散点+地毯
data(Salaries,package="car")
library(ggplot2)
ggplot(Salaries,aes(x=rank,y=salary))+geom_boxplot(fill="cornflowerblue",color="black",notch = TRUE)+geom_point(position="jitter",color="blue",alpha=.5)+geom_rug(side="1",color="black")


#图形组合：小提琴+箱线
data(singer,package="lattice")
ggplot(singer,aes(x=voice.part,y=height))+geom_violin(fill="lightblue")+geom_boxplot(fill="lightgreen",width=.2)



#分组 - 密度图
data(Salaries,package="car")
library(ggplot2)
ggplot(data=Salaries,aes(x=salary,fill=rank))+geom_density(alpha=.3)



#分组 - 散点图
ggplot(Salaries,aes(x=yrs.since.phd,y=salary,color=rank,shape=sex))+geom_point()


#分组 -直方图
ggplot(Salaries,aes(x=rank,fill=sex))+geom_bar(position="stack")+labs(title='position="stack"')

ggplot(Salaries,aes(x=rank,fill=sex))+geom_bar(position="dodge")+labs(title='position="dodge"')

ggplot(Salaries,aes(x=rank,fill=sex))+geom_bar(position="fill")+labs(title='position="fill"')

#分组 - 刻画
data(singer,package="lattice")
library(ggplot2)
ggplot(data=singer,aes(x=height))+geom_histogram()+facet_wrap(~voice.part,nrow=4)

ggplot(Salaries,aes(x=yrs.since.phd,y=salary,color=rank,shape=rank))+geom_point()+facet_grid(.~sex)

#分组 - 刻画 -核密度
ggplot(data=singer,aes(x=height,fill=voice.part))+geom_density()+facet_grid(voice.part~.)

#添加光滑曲线
data(Salaries,package="car")
library(ggplot2)
ggplot(data=Salaries,aes(x=yrs.since.phd,y=salary))+geom_smooth()+geom_point()

ggplot(data=Salaries,aes(x=yrs.since.phd,y=salary,linetype=sex,shape=sex,color=sex))+geom_smooth(method=lm,formula=y~poly(x,2),se=FALSE,size=1)+geom_point(size=2)


#坐标轴
data(Salaries,package="car")
library(ggplot2)
ggplot(data=Salaries,aes(x=rank,y=salary,fill=sex))+geom_boxplot()+scale_x_discrete(breaks=c("AsstProf","AssocProf","Prof"),labels=c("Assistant
Professor","Assistant
Professor","Full
Professor"))+scale_y_continuous(breaks=c(50000,100000,150000,200000),labels=c("$50K","$100K","$150K","$200K"))+labs(title="Faculty Salary by Rank and Sex",x="",y="")

ggplot(data=Salaries,aes(x=rank,y=salary,fill=sex))+geom_boxplot()+scale_x_discrete(breaks=c("AsstProf","AssocProf","Prof"),labels=c("Assistant
Professor","Assistant
Professor","Full
Professor"))+scale_y_continuous(breaks=c(50000,100000,150000,200000),labels=c("$50K","$100K","$150K","$200K"))+labs(title="Faculty Salary by Rank and Sex",x="",y="")+theme(legend.position = c(.1,.8))


#标尺
ggplot(mtcars,aes(x=wt,y=mpg,size=disp))+geom_point(shape=21,color="black",fill="cornsilk")+labs(x="Weight",y="Miles Per Gallon",title="Bubble Chart",size="Engine
Displacement")

data(Salaries,package="car")
ggplot(data=Salaries,aes(x=yrs.since.phd,y=salary,color=rank))+scale_color_manual(values=c("orange","olivedrab","navy"))+geom_point(size=2)


#主题
data(Salaries,package="car")
library(ggplot2)
mytheme<-theme(plot.title = element_text(face="bold.italic",size="14",color="brown"),
               axis.title = element_text(face="bold.italic",size=10,color="brown"),
               axis.text = element_text(face="bold",size=9,color="darkblue"),
               panel.background = element_rect(fill="white",color="darkblue"),
               panel.grid.major.y = element_line(color="grey",linetype=1),
               panel.grid.minor.y = element_line(color="grey",linetype=2),
               panel.grid.minor.x = element_blank(),
               legend.position = "top")

ggplot(Salaries,aes(x=rank,y=salary,fill=sex))+geom_boxplot()+labs(title="Salary by Rank and Sex",x="Rank",y="Salary")+mytheme


#多重图
data(Salaries,package="car")
library(ggplot2)
p1<-ggplot(data=Salaries,aes(x=rank))+geom_bar()
p2<-ggplot(data=Salaries,aes(x=sex))+geom_bar()
p3<-ggplot(data=Salaries,aes(x=yrs.since.phd,y=salary))+geom_point()

library(gridExtra)
grid.arrange(p1,p2,p3,ncol=3)


保存图片
ggsave(file="d:\mygraph.pdf")