品种特异位点，利用样本vcf参照打分

m0_51268845

于 2021-11-05 16:32:37 发布

阅读量262

点赞数

文章标签：生物信息学

本文链接：https://blog.csdn.net/m0_51268845/article/details/121165893

版权

import pandas as pd
df = pd.read_csv("labuladuo_loci.txt")
quanshe = pd.read_table("北京犬舍.vcf")

loci_list = quanshe['ID'].values.tolist()
vareity_list = quanshe.columns.tolist()
score = {}
for variety in variety_list:
   score[variety] = 0
quanshe.set_index(["ID"], inplace=True)
df.set_index(["ID"], inplace=True)
vareity_list.pop(2)
for locus in loci_list:
   for variety in vareity_list:
       if(quanshe.loc[locus,variety]==df.loc[locus,"ID"]):
           score[variety] = score[variety] + 1