||
废话少说,直接进入正题。
输入格式:
至少包括两列,一列是分组信息,一列是表型数据。缺失值使用NA表示。
保存为txt文件(如,genotype_phnotype_input.txt)。
代码如下,请根据自己的实际情况修改
#加载如下R包
library(ggplot2)
library(Hmisc)
library(dplyr)
library(export)
# data包含基因型和表型,如,genotype1\tgenotype2\tphnotype1\tphnotype2\phnotype3 ...
rawdata <- read.table('/Users/yaya/Desktop/scripts/genotype_phnotype_input.txt',sep='\t',header=T)
# 选择需要画图的标记和表型,如WMC1000,PH_02LY
data_column <- select(rawdata, WMC1000, PH_02LY)
# 去除含有NA数据的行
# data <- filter(data_column, WMC1000 != 'NA', PH_02LY != 'NA')
# 选择基因型为 0 和 1,表型不缺失的行
data <- filter(data_column, WMC1000 == '0' | WMC1000 == '1', PH_02LY != 'NA')
head(data)
cbPalette=c("forestgreen", "darkred")
p <- ggplot(data, aes(x=as.factor(WMC1000), y=PH_02LY)) + geom_boxplot() + geom_jitter(aes(color =as.factor(WMC1000)), width=0.1, height=0, show.legend=FALSE) + scale_color_manual(values = cbPalette,guide=FALSE)+ theme(panel.background = element_rect(fill = 'white'),axis.line = element_line(size=1.0, colour = "black"),axis.text.x = element_text(size = 12, family="Times New Roman", color = "black"), axis.text.y = element_text(size = 12, family="Times New Roman", colour = "black"),axis.title.x=element_blank(),axis.title.y= element_text(size=12, family="Times New Roman", color="black",vjust=1))+ ylab("PH_02LY") + xlab("WMC1000") + scale_x_discrete(breaks=c("0","1"),labels=c("CS type","Non-CS type"))
# 设置输出文件
graph2ppt(x=p, file="/Users/yaya/Desktop/scripts/WMC1000_MTKW.pptx", paper=A4)
# 设置箱线图颜色 geom_boxplot(aes(fill =as.factor(WMC1000))) + scale_color_manual(values = cbPalette,guide=FALSE)
出图效果
编辑好之后,可以导出为pdf等格式,可以按照要投稿的杂志格式组图和修改。
把这次图片的ppt格式和pdf格式放到下面的链接里了,可以看看是R包”export"否值得拥有。
https://pan.baidu.com/s/1W4FHyvSYdIqRspKPwr5vaA
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-1 11:12
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社