setwd("E:/bioinformatics/fanzhimin")

tg <- ToothGrowth

library(ggplot2)

summarySE <- function(data=NULL, measurevar, groupvars=NULL, na.rm=FALSE,

+                       conf.interval=.95, .drop=TRUE) {

+     library(plyr)

+ length2 <- function (x, na.rm=FALSE) {

+     if (na.rm) sum(!is.na(x))

+     else       length(x)

+ }

+ datac <- ddply(data, groupvars, .drop=.drop,

+                .fun = function(xx, col) {

+                    c(N    = length2(xx[[col]], na.rm=na.rm),

+                      mean = mean   (xx[[col]], na.rm=na.rm),

+                      sd   = sd     (xx[[col]], na.rm=na.rm)

+                    )

+                },

+                measurevar

+ )

+ 

+ datac <- plyr::rename(datac, c("mean" = measurevar))

+ datac$se <- datac$sd / sqrt(datac$N) 

+ ciMult <- qt(conf.interval/2 + .5, datac$N-1)

+ datac$ci <- datac$se * ciMult

+ 

+ return(datac)

+ }

tgc <- summarySE(tg, measurevar="len", groupvars=c("supp","dose"))

ggplot(tgc, aes(x=dose, y=len, colour=supp))+geom_point()+geom_line()+ geom_errorbar(aes(ymin=len-se, ymax=len+se), width=.1)+ggtitle("R2-36")

library(ggplot2)

D<-ggplot(midwest,aes(x=area,y=poptotal))+geom_point(aes(col=state,size=popdensity))+geom_smooth(method='loess')+xlim(c(0,0.1))+ylim(c(0,500000))+labs(title="R2-36")

mytheme<- theme(panel.grid = element_blank(), panel.background=element_blank(),axis.text=element_blank(),axis.ticks = element_blank())+theme(legend.position = "bottom")

D+mytheme

install.packages("RColorBrewer")

library(RColorBrewer)

 data1 <- read.csv("task3.csv", header = TRUE)

data1 <- data1[order(data1$PTS),]

row.names(data1) <- data1$Name

 data1 <- data1[,2:20]

 data2 <- data.matrix(data1)

heatmap(data2, Rowv=NA, Colv=NA, col = brewer.pal(9, "Blues"), scale="column", margins=c(5,10))