R2-29 第十次作业综合作业

CISO 2018-01-10 22:54:02 阅读: 1087

# Figure 1

library(ggplot2)

library(Rmisc)

dat <- ToothGrowth

dat2 = summarySE(dat,measurevar="len", groupvars=c("supp","dose"))

pd = position_dodge(0.1)

p=ggplot(dat2, aes(x=dose, y=len, colour=supp, group=supp))+geom_line(position=pd)

p+geom_errorbar(aes(ymin=len-se, ymax=len+se), colour="black", width=.1, position=pd)+geom_point(position=pd,size=3, shape=21, fill="white")+xlab("Dose (mg)") + ylab("Tooth length")+labs(color="Supplement type")+ggtitle("R2-29 CC")

# Figure 2

options(scipen = 999)

data("midwest",package="ggplot2")

p=ggplot(midwest,aes(x=area,y=poptotal))+geom_point(aes(col=state,size=popdensity))+geom_smooth(method="loess",se=T,col="red")+xlim(c(0,0.1))+ylim(c(0,500000))+xlab("Area")+ylab("Population")

p+ggtitle("R2-29 CC")

# Figure 3

library(reshape2)

library(plyr)

library(scales)

dat<- read.csv(file.choose())

dat2 <- melt(dat, id.vars=c("Name"))

dat2 <- ddply(dat2, .(variable), transform, label = rescale(value))

ggplot(dat2, aes(x=variable,y=Name)) + geom_tile(aes(fill=label))+scale_fill_gradient(low="yellow",high ="blue")+xlab("Type")+ ggtitle("R2-29 CC")