搜索
查看: 1793|回复: 0

[mRNA-seq] 转录组入门-关于plotMA函数的疑问

[复制链接]

2

主题

6

帖子

151

积分

注册会员

Rank: 2

积分
151
发表于 2017-12-14 11:24:14 | 显示全部楼层 |阅读模式
本帖最后由 ferfresh 于 2017-12-14 11:24 编辑

在参考hoptop大神的转录组入门教程中,运行plotMA函数的代码,出来的结果的截然不同。
大神的代码:
[Python] 纯文本查看 复制代码
plotMA(res, ylim = c(-5,5))
topGene <- rownames(res)[which.min(res$padj)]
with(res[topGene, ], {
  points(baseMean, log2FoldChange, col="dodgerblue", cex=2, lwd=2)
  text(baseMean, log2FoldChange, topGene, pos=2, col="dodgerblue")
})

示例结果:

我的代码:
[Python] 纯文本查看 复制代码
plotMA(res, ylim = c(-5,5))
topGene <- rownames(res)[which.min(res$padj)]
with(res[topGene, ], {
  points(baseMean, log2FoldChange, col="dodgerblue", cex=2, lwd=2)
  text(baseMean, log2FoldChange, topGene, pos=2, col="dodgerblue")
})

我的结果:

图中的坐标系完全不一样。不知道哪里出问题。请路过的群友赐教!





上一篇:【py】rosalind1_Counting DNA Nucleotides
下一篇:求助TMB的制作方法和流程
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|手机版|小黑屋|生信技能树 ( 粤ICP备15016384号  

GMT+8, 2019-9-17 07:16 , Processed in 0.029747 second(s), 30 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.