面试十大禁忌
使⽤clusterprofile做聚类分析library(clusterProfiler )
#
gene_symbol
EXOSC10
ARHGEF10L
VWA5B1
SRRM1
PTAFR
CSMD2
SH3GLB1
GBP6野猪肉可以吃吗
ZNF326
表白的歌词
AKNAD1
STRIP1
GOLPH3L
太阳的影子
...............
a <- read.table("",colClass = "character")
b <- a[,1]
eg = bitr(b, fromType="SYMBOL", toType="ENTREZID", annoDb="db")
head(eg)
write.table(eg, file="")
idType("db")
gene <- eg[,2]
head(gene)
ggo <- groupGO(gene = gene,
organism = "human",
罗汉果的ont = "BP",
level = 7,波普尔
readable = TRUE)
head(summary(ggo))
write.table(as.matrix(ggo@result), file="")
#导出ggo的结果
ego_cc <- enrichGO(gene = gene,
organism = "human",
ont = "CC",
pAdjustMethod = "BH",
minGSSize = 1,
pvalueCutoff = 1,
qvalueCutoff = 1,
readable = TRUE)
head(summary(ego))
write.table(as.matrix(ego_cc@result), file="")
ego_BP <- enrichGO(gene = gene,
organism = "human",
ont = "BP",
pAdjustMethod = "BH",
minGSSize = 1,
pvalueCutoff = 1,
qvalueCutoff = 1,
readable = TRUE)
write.table(as.matrix(ego_BP@result), file="")
ego_MF <- enrichGO(gene = gene,
双冠龙图片organism = "human",
ont = "MF",
pAdjustMethod = "BH",
minGSSize = 1,
pvalueCutoff = 1,
qvalueCutoff = 1,
readable = TRUE)
write.table(as.matrix( ego_MF@result), file="") kk <- enrichKEGG(gene = gene,
organism="human",
pvalueCutoff = 1,
qvalueCutoff = 1,
minGSSize = 1,
readable = TRUE,
u_internal_data =FALSE)
write.table(as.matrix(kk@result), file="")
# barplot(ggo, drop=TRUE, showCategory=24)
barplot(ego_cc, showCategory=24)
dotplot(ego_cc)
barplot(ego_BP, showCategory=24)
dotplot(ego_BP)
barplot(ego_MF, showCategory=24)
dotplot(ego_MF)
barplot(kk, drop=TRUE, showCategory=24)
dev.off()
>生态工程原理
