重测序及SNP分析

FST值的数据可视化

2022-06-16  本文已影响0人  爱吃海椒的妹妹

一、准备工作;

用excel创建如下图的excel数据,然后保存成CSV格式


image.png

创建工作目录

getwd()

把准备好的文件放在这个工作目录下

一、读入数据;

 dt1<- read.csv("test.csv",sep=",", header = T)

read.csv和read.csv2与read.table相同,但默认值除外。
它们用于读取“逗号分隔值”文件(’.csv’)或(read.csv2)在使用逗号作为小数点和分号作为字段分隔符的国家/地区中使用的变体。
类似地,read.delim和read.delim2用于读取默认分隔符为TAB的文件。 请注意,这些变体中的header =
TRUE和fill = TRUE,并且注释字符已禁用。
————————————————
版权声明:本文为CSDN博主「慎终君」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/wmm131333/article/details/100691190

A.read.delim和read.table的区别为:

读取数据速度不同、要求不同、空串不同。

一、读取数据速度不同

1、read.delim:read.delim的读取数据速度比read.table的读取数据速度更快。

2、read.table:read.table的读取数据速度比read.delim的读取数据速度更慢。

二、要求不同

1、read.delim:read.delim不要求所有列都对等,会按最大列,或指定的列数填充。

2、read.table:read.table严格要求所有列都对等。
三、空串不同

1、read.delim:read.delim空串时,可指定相应的字符或数字填充进行读取。

2、read.table:read.table不允许有空串,必须有相应的字符或数字才可以进行读取。

B.sep和header的含义见

初学者R语言:外部数据导入 - 简书 (jianshu.com)

geom_point()可以设置散点的大小,形状,颜色等,如geom_point(size=6,shape=1)具体见R语言绘图基础篇-散点图(geom_point) - 知乎 (zhihu.com)

二、加载ggplot2包;

library(ggplot2)

三、建立颜色与数据的关系,这里让点按照数据大小显示不同的颜色;

p1<-ggplot(dt1, aes(No, Fst,color= factor(Gene)))+geom_point()

四、设置x轴和y轴范围

设置x轴范围

p1<-p1+scale_x_continuous(limits = c(0, 114),breaks = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104[105],105,106,107,108,109,110,111,112,113,114))

1:114指的是1,2,3,4......114

设置y轴范围

p1<-p1+scale_y_continuous(limits = c(-0.065377,0.910141),breaks = c(0,0.25,0.50,0.75,1),labels = c("0.00","0.25","0.50","0.75","1.00"))

五、设置图例、坐标轴、图表的标题;

p1<-p1+labs(y="Fst",x="SNP")

六、自定义图表主题,对图表做精细调整;

top.mar=0.2

right.mar=0.2

bottom.mar=0.2

left.mar=0.2

mytheme<-theme(panel.grid.major =element_blank(),

              panel.grid.minor = element_blank(),

              panel.background = element_blank(),

              panel.border = element_blank(),

              axis.line.y = element_line(color = "black"),

              axis.line.x = element_line(color = "white"),

              axis.title.y = element_text(size = rel(1.2)),

              axis.title.x = element_text(size = rel(1.2),color = "black"),

              axis.text.y = element_text(size=rel(1.2),color="black"),

              axis.text.x = element_text(size=rel(1.2),color="black"),

              plot.margin=unit(x=c(top.mar,right.mar,bottom.mar,left.mar),units="inches"))



mytheme<-theme(panel.grid.major =element_blank(),

                panel.grid.minor = element_blank(),

                panel.background = element_blank(),

                panel.border = element_blank(),

                axis.line.y = element_line(color = "black"),

                axis.line.x = element_line(color = "black"),

                #axis.title.x = element_text(size = rel(1.2)),

                axis.title.y = element_text(size = rel(1.2)),

                axis.text.y = element_text(size=rel(1.2),color="black"),

                #axis.text.x = element_text(size=rel(1.2),color="black"),

                plot.margin=unit(x=c(top.mar,right.mar,bottom.mar,left.mar),units="inches"))

p1<-p1+mytheme
ggsave(file="4.pdf",width = 15,height = 4)

作者:EwanH
链接:https://www.jianshu.com/p/bb0beec0ed63
来源:简书
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

作者:EwanH
链接:https://www.jianshu.com/p/7b6701ae593b
来源:简书
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

上一篇下一篇

猜你喜欢

热点阅读