# 读取猪只档案
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)
ped文件批量修改个体号
最后编辑于 :
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 例如:我需要在当前目录下所有的css文件路径前追加域名 我想的是用sed替换去实现,鲍哥的思路是用for循环 1....
- for f in * ; do mv -- "$f" "PRE_$f" ; done 将‘PRE_’修改即可 fo...
- package test2; import java.io.BufferedReader; import java...