import sys import os f1 = open(sys.argv[1],'r') f2 = open(sys.argv[2],'r') fout = open(sys.argv[3],'w') dic1 = {} for line1 in f1: tmp1 = line1.strip().split("\t") Gene_ID = tmp1[0] anno = tmp1[1:] dic1[Gene_ID] = anno dic2 = {} for line2 in f2: tmp2 = line2.strip().split("\t") Gene_ID = tmp2[0] gene_name = tmp2[1] dic2[Gene_ID] = gene_name for key in dic2: if key in dic1: fout.write(key + "\t" + gene_name + "\t" + anno + "\n") else: fout.write(key + "\t" + gene_name + "\t" + "NA" + "\n" ) f1.close() f2.close() fout.close(),修改fout处,实现将dic1的键值部分添加到dic2的键和值后面
时间: 2024-01-02 19:04:07 浏览: 17
可以将最后的 `fout.write` 的语句改为:
```
if key in dic1:
fout.write(key + "\t" + dic2[key] + "\t" + "\t".join(dic1[key]) + "\n")
else:
fout.write(key + "\t" + dic2[key] + "\t" + "NA" + "\n" )
```
这样就可以将 `dic1` 的键值部分添加到 `dic2` 的键和值后面了。
相关问题
import sys import os f1 = open(sys.argv[1],'r') f2 = open(sys.argv[2],'r') fout = open(sys.argv[3],'w') dic1 = {} for line1 in f1: tmp1 = line1.strip().split("\t") Gene_ID = tmp1[0] anno = tmp1[1:] dic1[Gene_ID] = anno dic2 = {} for line2 in f2: tmp2 = line2.strip().split("\t") Gene_ID = tmp2[0] gene_name = tmp2[1] dic2[Gene_ID] = gene_name for key in dic2: if key in dic1: fout.write(key + "\t" + gene_name + "\t" + anno + "\n") else: fout.write(key + "\t" + gene_name + "\t" + "NA" + "\n" ) f1.close() f2.close() fout.close(),修改fout处,实现将dic1的anno添加到dic2的键和值后一列
要将 `dic1` 的 `anno` 添加到 `dic2` 的键和值后一列,可以将 `fout.write()` 函数中的第三个参数改为 `gene_name + "\t" + dic1[key][0] + "\t" + dic1[key][1] + "\t" + dic1[key][2] + "\n"`,假设 `dic1` 中每个键对应的值有三个元素。修改后的代码如下:
```
import sys
import os
f1 = open(sys.argv[1],'r')
f2 = open(sys.argv[2],'r')
fout = open(sys.argv[3],'w')
dic1 = {}
for line1 in f1:
tmp1 = line1.strip().split("\t")
Gene_ID = tmp1[0]
anno = tmp1[1:]
dic1[Gene_ID] = anno
dic2 = {}
for line2 in f2:
tmp2 = line2.strip().split("\t")
Gene_ID = tmp2[0]
gene_name = tmp2[1]
dic2[Gene_ID] = gene_name
for key in dic2:
if key in dic1:
fout.write(key + "\t" + gene_name + "\t" + dic1[key][0] + "\t" + dic1[key][1] + "\t" + dic1[key][2] + "\n")
else:
fout.write(key + "\t" + gene_name + "\t" + "NA" + "\n" )
f1.close()
f2.close()
fout.close()
```
这样修改后,输出文件的格式为 `Gene_ID gene_name anno[0] anno[1] anno[2]`,其中 `anno[0]`、`anno[1]` 和 `anno[2]` 分别对应 `dic1` 中每个键对应的值的三个元素。
import sys import os f1 = open(sys.argv[1],'r') f2 = open(sys.argv[2],'r') f3 = open(sys.argv[3],'r') fout = open(sys.argv[4],'w') dic1 = {} for line1 in f1: tmp1 = line1.strip().split() Gene_ID = tmp1[0] # anno = tmp1[1:] # dic1[Gene_ID] = anno dic2 = {} for line2 in f2: tmp2 = line2.strip().split("\t") Gene_ID = tmp2[0] gene_name = tmp2[6] dic2[Gene_ID] = gene_name dic3 = {} for line3 in f3: tmp3 = line3.strip().split("\t") new_name = tmp3[0] raw_name = tmp3[1] dic3[new_name] = raw_name for key in dic2: if key in dic1: fout.write(key + "\t" + dic2[key] + "\n") else: fout.write(key + "\t" + dic2[key] + "\t" + "NA" + "\n" ) f1.close() f2.close() fout.close()修改代码,其中f1中只有一列数据,需要根据f1文件的内容去f2文件检索找到对应的第7列信息,由于第7列每行的字符串由逗号分割,需去掉逗号并进行换行成为新的一列,再根据新的一列信息去和f3第1列匹配,最后打印出f3中与匹配结果对应的第二列信息
import sys
import os
f1 = open(sys.argv[1],'r')
f2 = open(sys.argv[2],'r')
f3 = open(sys.argv[3],'r')
fout = open(sys.argv[4],'w')
dic2 = {}
for line2 in f2:
tmp2 = line2.strip().split("\t")
Gene_ID = tmp2[0]
gene_info = tmp2[6].replace(",", "\n")
dic2[Gene_ID] = gene_info
dic3 = {}
for line3 in f3:
tmp3 = line3.strip().split("\t")
new_name = tmp3[0]
raw_name = tmp3[1]
dic3[new_name] = raw_name
for line1 in f1:
Gene_ID = line1.strip()
if Gene_ID in dic2:
gene_info = dic2[Gene_ID].strip()
if gene_info in dic3:
fout.write(Gene_ID + "\t" + gene_info + "\t" + dic3[gene_info] + "\n")
else:
fout.write(Gene_ID + "\t" + gene_info + "\t" + "NA" + "\n")
else:
fout.write(Gene_ID + "\t" + "NA" + "\t" + "NA" + "\n")
f1.close()
f2.close()
fout.close()
相关推荐
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![7z](https://img-home.csdnimg.cn/images/20210720083312.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)