как добавить ключи и значения из вложенных словарей в списки - PullRequest
0 голосов
/ 02 мая 2019

Я строю три вложенных словаря для анализа моих больших данных.Я пытаюсь произвольно анализировать значения внутри них, чтобы составить точечный график, поэтому я создаю список, чтобы добавить к ним свои данные, а затем создаю график рассеяния с помощью matplotlib.Моя проблема в том, что я получаю ошибку, когда пытаюсь добавить!Ошибка типа: не подлежащий обработке тип: 'список'.поэтому я запутался, чтобы изменить структуру своих словарей или есть возможность справиться с этим из того, что я создал.

моя структура словарей выглядит следующим образом:

data_geo1:

'ENSG00000268358': {'Sample_19-leish_023_v2': 0, 'Sample_4-leish_012_v3': 0, 'Sample_25-leish027_v2': 0, 'Sample_6-leish_015_v3': 0, 'Sample_23-leish026_v2': 1, 'Sample_20-leish_023_v3': 0, 'Sample_18-leish_022_v3': 0, 'Sample_10-leish_017_v3': 0, 'Sample_13-leish_019_v2': 0, 'Sample_1-Leish_011_v2': 0, 'Sample_11-leish_018_v2': 0, 'Sample_3-leish_012_v2': 0, 'Sample_2-leish_011_v3': 0, 'Sample_29-leish032_v2': 0, 'Sample_8-leish_016_v3': 0, 'Sample_28-leish028_v3': 0, 'Sample_27-leish028_v2': 1, 'Sample_26-leish027_v3': 0, 'Sample_12-leish_018_v3': 0, 'Sample_5-leish_015_v2': 0, 'Sample_16-leish_021_v3': 0, 'Sample_21-leish_024_v2': 0, 'Sample_9-leish_017_v2': 0, 'Sample_24-leish026_v3': 1, 'Sample_22-leish_024_v3': 0, 'Sample_14-leish_019_v3': 0, 'Sample_30-leish032_v3': 0, 'Sample_7-leish_016_v2': 0, 'Sample_15-leish_021_v2': 0, 'Sample_17-leish_022_v2': 1}

data_ali:

{'ENSG00000268358': {'Sample_19-leish_023_v2': 0, 'Sample_16-leish_021_v3': 2, 'Sample_20': 0, 'Sample_24-leish026_v3': 1, 'Sample_6-leish_015_v3': 0, 'Sample_12-leish_018_v3': 0, 'Sample_22-leish_024_v3': 0, 'Sample_23-leish026_v2': 2, 'Sample_25-leish027_v2': 0, 'Sample_18-leish_022_v3': 1, 'Sample_14': 0, 'Sample_2-leish_011_v3': 0, 'Sample_13-leish_019_v2': 0, 'Sample_1-Leish_011_v2': 0, 'Sample_11-leish_018_v2': 0, 'Sample_20-leish_023_v3': 0, 'Sample_3-leish_012_v2': 0, 'Sample_10-leish_017_v3': 1, 'Sample_7': 0, 'Sample_29-leish032_v2': 1, 'Sample_8-leish_016_v3': 0, 'Sample_6': 0, 'Sample_7-leish_016_v2': 0, 'Sample_9': 0, 'Sample_8': 0, 'Sample_27-leish028_v2': 0, 'Sample_26-leish027_v3': 0, 'Sample_5': 1, 'Sample_4': 0, 'Sample_3': 0, 'Sample_19': 0, 'Sample_1': 0, 'Sample_2': 0, 'Sample_9-leish_017_v2': 0, 'Sample_5-leish_015_v2': 0, 'Sample_4-leish_012_v3': 0, 'Sample_21-leish_024_v2': 0, 'Sample_18': 0, 'Sample_13': 0, 'Sample_12': 0, 'Sample_11': 0, 'Sample_10': 1, 'Sample_17': 0, 'Sample_16': 0, 'Sample_15': 1, 'Sample_14-leish_019_v3': 0, 'Sample_30-leish032_v3': 0, 'Sample_28-leish028_v3': 1, 'Sample_15-leish_021_v2': 0, 'Sample_17-leish_022_v2': 0}

здесь вся моя структура кода с начала,как вы видите в последних строках, я пытался создать список и добавить свои значения в список, но мне это не удалось.

import os
import numpy as np
import matplotlib.pyplot as plt

path = "/home/ali/Desktop/data/"
root = "/home/ali/Desktop/SAMPLES/"



data_geo1={}
with open(path+"GSE98212_H_DE_genes_count.txt","rt") as fin: #data for sample 1-30
    h = fin.readline()
    sample1 = h.split()
    sample_names = [s.strip('"') for s in sample1[1:31]]
    for l in fin.readlines():
        l = l.strip().split()
        if l:
            gene1= l[0].strip('"')
            data_geo1[gene1] = {}
            for i, x in enumerate(l[1:31]):
                data_geo1[gene1][sample_names[i]] = int(x)

#print(data_geo1)

data_geo2={}
with open (path+"GSE98212_L_DE_genes_count.txt","rt") as fin:
    h= fin.readline()
    sample2=h.split()
    sample_names=sample2[1:21]
    for l in fin.readlines():
        l = l.strip().split()
        if l:
            gene2= l[0].strip()
            data_geo2[gene2]={}
            for i,x in enumerate (l[1:21]):
                data_geo2[gene2][sample_names[i]]= int(x)
#print(data_geo2)

data_ali={}

for sample_name in os.listdir(root):
    with open(os.path.join(root, sample_name, "counts.txt"), "r") as fin:
        for line in fin.readlines():
            gene, reads = line.split()
            reads = int(reads)
            if gene.startswith('ENSG'):
                data_ali.setdefault(gene, {})[sample_name] = reads
        gene = l[0].strip()
#print(data_ali)

list_samples= data_ali[gene].keys()
#print(list_samples)
for sample in list_samples:
    reads_data_ali = []
    for gene in data_ali.keys():
        reads_data_ali.append(data_ali[gene][sample_name])

Я ожидаю, что результат будет:

[[0, 0], [0, 2], [11, 12], [4, 4], [18, 17], [2, 2], [381, 383], [1019, 1020], [198, 194], [66, 65], [2223, 2230], [30, 30], [0, 0], [33, 34], [0, 0], [411, 409], [804, 803], [11829, 7286], [137, 139], [277, 278], [3475, 3482], [5, 5], [2, 1], [70, 70], [48, 48], [234, 232], [121, 120], [928, 925], [220, 159], [165, 165], [702, 700], [1645, 1643], [79, 78], [1064, 1067], [971, 972], [0, 0]]

Ответы [ 2 ]

0 голосов
/ 02 мая 2019

Что касается вашего описания, я предполагаю, что ваш код создания словарей data_ali и data_geo1 получил правильные выходные данные, поэтому проблема может заключаться в последнем коде создания списка.Я нахожу два вопроса: 1 for gene in data_ali.keys():, в следующем цикле, reads_data_geo1.append(data_geo1[gene1][sample_names]), здесь это [gene1] 2 for sample in list_samples:, поэтому, возможно, вам следует использовать reads_data_ali.append(data_ali[gene][sample]), вы можете изменить имя этих переменных и посмотреть, работает ли он.

0 голосов
/ 02 мая 2019

Вы можете попытаться избежать ключевой ошибки, проверив, существует ли ключ в вашем словаре перед .append (...). Попробуйте взглянуть на словарь .get () метод. Это хорошо для предотвращения ошибок такого типа.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...