火山简图画法,画简单的火山

　　1、导入数据

　　importpandasaspd #数据分析

　　importnumpyasnp #科学计算

　　进口海船#统计可视化

　　#读取数据

　　df=pd.read_csv( ./dataset_volcano.txt ，sep=t )

　　结果=pd .数据帧()

　　结果[x]=df[logFC]

　　结果[y]=df[P。值]

　　结果[-log10(pvalue)]=-df[ P . Value ].申请(NP。日志10)http://www . Sina.com/2、设置阈值

　　#设置pvalue和差异倍数的阈值

　　截止值=0.0000001

　　3、设置分组

　　#分组为向上、正常、向下

　　结果。loc[(结果。xcut _ off _ logFC)(结果。ycut _ off _ pvalue)， group]=up

　　结果。loc[(结果。x-cut _ off _ logFC)(结果。ycut _ off _ pvalue)， group]=down

　　结果。loc[(结果。x=-cut _ off _ logFC)(结果。x=cut _ off _ logFC)(结果。y=cut _ off _ pvalue)， group ]= normal 4、绘制散点图

　　#绘制散点图

　　ax=sns.scatterplot(x=x ，y=-log10(pvalue)，

　　hue=group ，

　　hue_order=(下,正常,上),

　　调色板=(#377EB8 ，灰色， #E41A1C )，

　　=0.5,

　　s=15http:///

　　#确定坐标轴显示范围

　　xmin=-6

　　xmax=10

　　ymin=7

　　ymax=13

　　ax.spines[右]。set_visible(False)#去掉右边框

　　ax.spines[top].set_visible(False)#去掉上边框

　　ax.vlines(-cut_off_logFC，ymin，ymax，color=dimgrey ，linestyle=dashed ，linewidth=1)#画竖直线

　　ax.vlines(cut_off_logFC，ymin，ymax，color=dimgrey ，linestyle=dashed ，linewidth=1)#画竖直线

　　斧头。hlines(-NP。log10(cut _ off _ pvalue)，xmin，xmax，color=dimgrey ，linestyle=dashed ，linewidth=1)#画竖水平线

　　ax.set_xticks(range(xmin，xmax，4))#设置x轴刻度

　　ax.set_yticks(range(ymin，ymax，2))#设置英语字母表的第25个字母轴刻度

　　ax.set_ylabel(-log10(pvalue)，fontweight=bold)#设置英语字母表的第25个字母轴标签

　　ax.set_xlabel(log2(foldchange)，fontweight=bold)#设置x轴标签以上就是大蟒中画火山图的方法，希望能对大家有所帮助，更多知识尽在盛行信息技术软件开发工作室106 .5、设置散点图

郑重声明：本文由网友发布，不代表盛行IT的观点，版权归原作者所有，仅为传播更多信息之目的，如有侵权请联系，我们将第一时间修改或删除，多谢。