R代码很容得到,关键是数据不会准备,遇到问题看不懂
R代码真的很容易得到,例如GOplot包的代码,百度一下就得到了:http://wencke.github.io/, 没有错这是这么简单,可以直接复制+粘贴
library(GOplot)
# Load the dataset
data(EC)
# Get a glimpse of the data format of the results of the functional analysis...
head(EC$david)
# ...and of the data frame of selected genes
head(EC$genelist)
# Generate the plotting object
circ <- circle_dat(EC$david, EC$genelist)
chord <- chord_dat(circ, EC$genes, EC$process)
# Generate the matrix with a list of selected genes
chord <- chord_dat(data = circ, genes = EC$genes)
# Generate the matrix with selected processes
chord <- chord_dat(data = circ, process = EC$process)
# Create the plot
chord <- chord_dat(data = circ, genes = EC$genes, process = EC$process)
GOChord(chord, space = 0.02, gene.order = 'logFC', gene.space = 0.25, gene.size = 5)
。。。
看到官网给的代码好简单呀,自己直接复制到R,结果也出来了,而且很漂亮:
这是直接用R包自带的数据来做的,肯定非常简单的。现在问题来了,自己怎么准备数据?
这也是很多临床医生比较头痛的问题,当然对R比较熟悉的会自己准备,但是很多初学者搞数据挖掘根本自己看不懂,也不会准备,特别是R报错的时候,更加头痛了,绝大多数据R报错也是数据出现了一些小问题,只是看不懂而已。
现在有非常多论坛或者公众号都经常会分享一些代码,例如(GOplot,WGCNA等等)然而绝多数都是用R包自带的数据,然而对很多人来说并没有卵用,因为需要什么数据都不知道,怎么读入也不知道。很多老师讲课也很喜欢用R自带数据来讲,因为这样比较方便,但是学员就惨了,因为数据不知道怎么样准备与读入,后面的分析再怎么样整都没有卵用。
● 火山图的高级玩法
赞 (0)