火山图(2)

解决了这些问题,火山图表达的意思也就一目了然了,横着看差异,竖着看(统计学)意义。

三条虚线将其分成了六个区,两条竖线是差异倍数线,多取1.2/1.5/2,绝对值越大,差异表达越显著;一条横线是显著性水平线,多取0.05/0.01,如图所示:
①基因或蛋白表达显著,差异性明显下调
②基因或蛋白表达显著,差异较小或无差异
③基因或蛋白表达显著,差异性明显上调
④基因或蛋白表达不显著,差异性明显下调
⑤基因或蛋白表达不显著,差异较小或无差异
⑥基因或蛋白表达不显著,差异性明显上调
其实在我们实际看图的时候,图的左上角和右上角是差异基因集中的地方,这才是我们关注的重点。
#医学生[超话]##医学sci##生物信息学##R语言绘图##生信分析##火山图#
原文节选自:医学僧的科研日记(id:zzudoctor)

火山图(1)

火山图(Volcano plot)是散点图的一种,可以直观展示两组样品的的差异程度及其统计学显著性,快速直观地识别变化幅度较大且具有统计学意义的数据点。常应用于转录组研究,也能应用于基因组,蛋白质组,代谢组等。

如图,这样的火山图在表达什么意思呢?

我们先把它简单分解一下,横坐标为logFC,纵坐标为-log10(adj P value),三条虚线将其分成了六部分。
在理解每一部分代表什么意思之前,我们需要先解决几个问题:

1.什么是logFC:FC全拼为fold change,翻译成中文就是差异倍数,简单来说就是基因在一组样品中的表达值的均值除以其在另一组样品中的表达值的均值。所以火山图只适合展示两组样品之间的比较。那么这样就会出现一个问题,两个数相除,大于1代表上调,小于1代表下调,最大值可以无穷大,最小值确为0,这样在图里展现出来就非常的不美观,下调的数据点会挤在很小的一个空间里。所以这时候就需要对数据做一些转换,也就是我们的横坐标logFC,上调的数据转换后Log (fold change)都大于等于1,下调的数据转换后Log (fold change)都小于等于-1。这样无论是展示还是描述是不是都方便了许多。

2.什么是-log10(adj P value):我们都知道P value是用来衡量检验得到的结果是否具有统计学差异,常用P-value<0.05作为统计检验显著的常规标准。adjustedP-value是校正后的P值,这涉及到一个统计学问题,数据较多时经过多次检验,会引入假阳性结果,需要做一个多重假设检验校正。这个校正怎么做呢?统计学家提出了许多计算方法,如Bonferroni correction ,holm, hochberg, hommel, BH, BY, fdr等。BH是我们比较常用的一个校正方法,获得的值是假阳性率 FDR (false discovery rate)。FDR筛选时就可以不用遵循0.05这个标准了。我们可以设置FDR<0.05表示我们容许数据中存在至多5%假阳性率;FDR<0.1表示我们对假阳性率的容忍度至多是10%。由于FDR值是0-1之间,数值越小表示显著性越大。而经过-Log 10 (adjusted P-value)转换后表示数值越大越显著,更利于我们展示结果。此外,adjusted P-value, q value, fdr一般代表相同的含义,都是多重假设检验校正后的P-value,可能的区别就在于校正算法的不同。

解决了这些问题,火山图表达的意思也就一目了然了,横着看差异,竖着看(统计学)意义。#医学生[超话]##医学sci##生物信息学##R语言绘图##火山图#

原文节选自:医学僧的科研日记(id:zzudoctor)

小提琴图

小提琴图 (Violin Plot)是用来展示多组数据的分布状态以及概率密度。这种图表结合了箱形图和密度图的特征,主要用来显示数据的分布形状。跟箱形图类似,但是在密度层面展示更好。在数据量非常大不方便一个一个展示的时候小提琴图特别适用。

如此优雅的名字和形状,经常出现在高分文章中。

那么小提琴图该怎么看呢?前面说什么来着,它就是核密度图和箱线图的结合体,那不就简单多了!

诸位请看:小提琴图的内部是箱线图(有的图中位数会用白点表示,但归根结底都是箱线图的变化);外部包裹的就是核密度图,某区域图形面积越大,某个值附近分布的概率越大。

通过箱线图,可以查看有关数据的基本分布信息,例如中位数,平均值,四分位数,以及最大值和最小值,但不会显示数据在整个范围内的分布。如果数据的分布有多个峰值(也就是数据分布极其不均匀),那么箱线图就无法展现这一信息,这时候小提琴图的优势就展现出来了!

那么大家以后看到这样的图是不是也明白了作者想要侧重表达什么了呢!#生物信息学##R语言绘图##生信分析##小提琴图# 节选自:医学僧的科研日记(ID:zzudoctor)


发布     👍 0 举报 写留言 🖊   
✋热门推荐
  • 看了知乎上的影评和解析感觉自己看了个假的电影“孤儿院是什么”“孤儿院就是关没有小孩的大人的”“那如果没有我,你也会被关进孤儿院吗”“是的”“你刚才吹的是什么?
  • ①13.9!查名一猫睫毛膏女防水纤长卷翘浓$zPGkYwKmPZf$ ②18!李宁瑜伽球防爆加厚普拉提球$K0VrYwKHeVG$ ③19.9!洁柔抽纸Face
  • さすが和辻家、ソファーが豪華原文镜像♥️ 朋友圈里可爱的句子 2.0 ver. ♥️1.请你一定要对号入座 我说的所有星光月亮 清泉小溪 蘸着糖的奶油和蜂蜜 这
  • 看一看外围市场表现,美股三大指数小幅收跌,A50期指小幅反弹,从情绪上来讲,前面两天跌的差不多了,市场量也比较小,今天大概率是能收红小反弹的,但是不要幻想太多,
  • 当听到流星陨落的消息,落下的眼泪在为你惋惜,亲爱的你,或者我不认识的你们,多一点与世界抗争的勇气,于千万人之中,我在默默支持你!#晚安公主#太多正在与抑郁症纠缠
  • 【玛4:4】“你们当记念我仆人摩西的律法,就是我在何烈山为以色列众人所吩咐他的律例典章。【玛4:3】你们必践踏恶人,在我所定的日子,他们必如灰尘在你们脚掌之下。
  • 微博创作的这个特点也决定了博文内容的良莠不齐,优质内容与无病呻吟并存、真实报道与无端指责同台,加大了平台方的管理难度。当然,如果你有独到的见解,也可以跨领域进行
  • 人就是要有短暂的幸福感满足感才行今天给妈妈打电话 正好被问起有没有 于是面临贫困的我突然多了二百块 嘻嘻 不用在余额宝提钱辣突然发现我们的互联网+老师挺通情达理
  • #子平果老占星[超话]#为什么说算命不是迷信#星座运势12月##紫薇星盘# 八字算命主要根据人的出生时间来推断命的格局与人生运势,这种预测方法是中华国学智慧的结
  • #今日早盘观点#短期来看,经济增长尚难言见底、通胀预期存在上行风险、股票市场供给增加、中美签署第一阶段协议时间尚不确定、当前盈利兑现却并不充分等因素影响下,预计
  • 尊重降温穿上秋裤,想知道双十一剁手的你们为自己买秋裤了吗?前来交作业[嘻嘻],今天是用一件卫衣抵御成都10℃天气的我。
  • 奇怪我为什么会喜欢千玺,把一个弟弟当偶像,就因为他长得帅吗,慢慢地我不再去对呛身边那些对他恶语相向的人,因为我觉得知道我喜欢千玺还在我面前这么说是一件很不礼貌的
  • 结果我等到现在[微笑] 我是还没长大是吧[微笑] 还是小朋友是吧[微笑] 那q弹的肚肚是啥玩意[微笑]:妈我觉得我长大了。我问爸爸“点解我有肚煲嘅?
  • 同时,本团支部为每一位参与活动的同学准备了抽奖环节,他们乐在其中,并祝我们的团日活动完满举行。大家提前准备好活动所需,主动的去负责一个环节,认认真真的向经过的同
  • 【投稿】洞娘您好 之前借平台求了一下本地群没有蹲到 自己建了一个 主要就想找近一点的娃娘大家一起玩 互相安利 偶尔可以约个外拍约dp之类的 介绍一下自己 我坐标
  • 真正的朋友不只是锦上添花,更多是雪中送炭,不要拒绝真诚的话,更不要拒绝一颗真诚的心广州市深圳市长沙市上海市成都市天津市杭州市武汉市青岛市济南市北京市无[赞啊]痕
  • 知难而上,兵荒马乱成了日常,但,”因为是凭借我自己意愿选择的生活,所以这些就是我正常的生活,和窗外的树木,清晨的咖啡一样正常”(大概记忆)。勇敢地尝试每一个可能
  • 上次看莫扎特是在飞机上,当戏说历史的闹剧看的,觉得戏里的莫扎特太美国太不古典,这次感觉还好,电影主题甚得吾心:平庸者无论用尽什么方法也战胜不了天才,最终只能沦为
  • “四有”是国家对公民的基本要求,也是提高整个中华民族的思想道德素质和科学文化素质的基本内容。因此,要实现社会主义现代化,就要培养一代有理想、有道德、有文化、有纪
  • 生日快乐我的大震瓜,希望你健康幸福,未开多多指教,想爆你丑照还是留点面子吧[吃瓜],我还想跟你一起摆地摊[笑cry]各种不正经的想法...似兄弟的情谊[笑cry