完成品 以下のようなVolcano plotをPythonで書いていきます。 Script 以下をコピペで上記のグラフが出てきます。 import pandas as pd url = "https://raw.githubusercontent.com/vappiah/bioinfoscripts/main/testvolcano.csv" c = pd . read_csv ( url ) def volcanoplot ( de , foldchange , pvalue , featurename , label = 10 ) : import numpy as np import pandas as pd import matplotlib import matplotlib . pyplot as plt import seaborn as sns import joblib from adjustText import adjust_text topfeature = [ ] de = de . sort_values ( pvalue ) de [ "logpvalue" ] = np . log10 ( de [ pvalue ] ) * - 1 conditionlist = [ ( de [ foldchange ] > 1 ) & ( de [ pvalue ] < 0.1 ) , ( de [ foldchange ] < - 1 ) & ( de [ pvalue ] < 0.1 ) , ] choicelist = [ 'UP' , 'DOWN' ] de [ 'diffexpressed' ] = np . select ( conditionlist , choicelist , default = 'NaN' ) ...