折线图-ggplot2 

http://blog.163.com/yugao1986@126/blog/static/6922850820131161531421/

http://blog.sina.com.cn/s/blog_7cffd1400101f50v.html

《R Graphics Cookbook-By Winston Chang》

#======================

折线图

library(ggplot2)  #作图

library(gcookbook)  #案例数据

library(plyr)  #ddply()、desc()等函数

# In addition to the variables mapped to the x- and y-axes, map another

# (discrete) variable to colour or linetype Load plyr so we can use ddply()

# to create the example data set Summarize the ToothGrowth data

tg <- ddply(ToothGrowth, c("supp", "dose"), summarise, length = mean(len))

# Map supp to colour

ggplot(tg, aes(x = dose, y = length, colour = supp)) + geom_line()

 

==============

# save as :

# figure_TSS_5k_methy_level_Male_Adipose.R.pdf

library(ggplot2)  #作图

setwd("E:\\_data\\")

group<-rep("xxxx",200)

window_index<-rep(1,200)

methy<-rep(0.0001,200)

inputfile_name=paste0("level.average")  

mydata <- read.table(inputfile_name, header=FALSE)

for(ii in 1:100)

{

  group[ii]<-"MAL"

  window_index[ii]<-ii

  methy[ii]<-mydata[1,ii]

}

inputfile_name=paste0(".level.average")  

mydata <- read.table(inputfile_name, header=FALSE)

for(ii in 1:100)

{

  group[ii+100]<-"MAC"

  window_index[ii+100]<-ii

  methy[ii+100]<-mydata[1,ii]

}

   

mydata<-data.frame(group=group,window_index=window_index,methy=methy)

ggplot(mydata, aes(x = window_index, y =methy , colour = group)) + geom_line()

==============

dev.new(width=5, height=4)  #窗口大小尺寸

==============