ped文件批量修改个体号

# 读取猪只档案
a=pd.read_excel(r"E:\新希望六和\文件\测序\20230811-205id-淄博\PLINK_150823_0352\秀杰猪种猪档案2023-10-12+09_20_0820231012.xlsx")
# 读取ped文件
snp_map_file = pd.read_table(open(r"E:\新希望六和\文件\测序\20230811-205id-淄博\PLINK_150823_0352\889647.ped"),delim_whitespace=True,encoding="gb18030",header=None,low_memory=False)
a2=pd.merge(snp_map_file[1],a[["耳牌号","个体号"]],left_on= 1,right_on = "耳牌号",how="left")
fh_914out=pd.concat([a2["个体号"],snp_map_file.iloc[:,2:]], axis=1).reset_index()
fh_914out["index"]=fh_914out["index"]+1
fh_914out.to_csv(r"E:\新希望六和\文件\测序\20230811-205id-淄博\淄博205id修改后.ped",sep='\t',index = None,header=None)
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容