簡體   English   中英

Python:如何將列與文本條目以及彼此進行比較

[英]Python: How to compare columns with text entries and with each other

我是一個剛開始編寫腳本的 Python 新手。 我有一個包含三個樣本的數據框,對於每個樣本,我都有一個肽序列列表,例如:

d = {'Sample 1': ['QSFLEVSYYPMAGYIKEDSIM', 'MLPIQTRIAS', 'AAVACTVLRCLAAEQQTSRSVDEAY'], 'Sample 2': ['QSFLEVSYYPTEIRQMGM', 'AEAARLVLAARIKGDAM', 'AAVACTVLRCLAAEQQTSRSVDEAY'], 'Sample 3': ['AAIGVHLGCTSACVAVYAADRGMKL', 'QSFLEVSYYPTEIRQMGM', 'AAVACTVLRCLAAEQQTSRSVDEAY']}
df = pd.DataFrame(data=d)

我想要做的是找出三個樣品中哪些肽是相互的,並繪制維恩圖

我用兩個樣本做了這個,但后來我找不到太多幫助,所以我“手動”做了維恩繪圖

import matplotlib_venn as vplt
from matplotlib import pyplot as plt
from matplotlib_venn import venn2, venn2_circles
import pandas as pd
import numpy as np

df_txt = pd.read_csv('all_Peptides_sortFN.txt', delimiter='\t')
df_txt.count()
df_txt['Peptides'] = df_txt["FN"] + df_txt["FN_18"]
AllPeps = df_txt.groupby('Peptides').count()
AllPeps.drop_duplicates(subset=None, keep='first', inplace=False)
AllPeps = df_txt.groupby('Peptides').count()
AllPeps.drop_duplicates(subset=None, keep='first', inplace=False)
AllPeps['result'] = np.where(AllPeps['FN'] == AllPeps['FN_18'], '1', '0')
AllPeps.head()
AllPeps.result.value_counts()
plt.figure(figsize=(10,10))
v = venn2(subsets=(3411, 1737, 1150), set_labels = ('Filgrastim', 'Filgrastim 18O'))
plt.title('Comparison of peptides in just digested and 18O-labeled Filgrastim(new batch)')
plt.savefig('VennFN.png')
plt.show()

我感謝每一個建議:)

你試過這個嗎?

from matplotlib_venn import venn3
venn3([set(v) for v in d.values()], set_labels=d.keys())

暫無
暫無

聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM