aboutsummaryrefslogtreecommitdiff
path: root/wqflask/utility/helper_functions.py
blob: cf16879f97a44f44ae73f8c940d0bc391f9430be (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
from __future__ import absolute_import, print_function, division

from base.trait import GeneralTrait
from base import data_set
from base.species import TheSpecies

from wqflask import user_manager

from flask import Flask, g

import logging
logger = logging.getLogger(__name__ )


def get_species_dataset_trait(self, start_vars):
    #assert type(read_genotype) == type(bool()), "Expecting boolean value for read_genotype"
    self.dataset = data_set.create_dataset(start_vars['dataset'])
    logger.debug("After creating dataset")
    self.species = TheSpecies(dataset=self.dataset)
    logger.debug("After creating species")
    self.this_trait = GeneralTrait(dataset=self.dataset,
                                   name=start_vars['trait_id'],
                                   cellid=None,
                                   get_qtl_info=True)
    logger.debug("After creating trait")

    #if read_genotype:
    #self.dataset.group.read_genotype_file()
    #self.genotype = self.dataset.group.genotype


def get_trait_db_obs(self, trait_db_list):
    if isinstance(trait_db_list, basestring):
        trait_db_list = trait_db_list.split(",")

    self.trait_list = []
    for trait in trait_db_list:
        data, _separator, hmac = trait.rpartition(':')
        data = data.strip()
        assert hmac==user_manager.actual_hmac_creation(data), "Data tampering?"
        trait_name, dataset_name = data.split(":")
        dataset_ob = data_set.create_dataset(dataset_name)
        trait_ob = GeneralTrait(dataset=dataset_ob,
                               name=trait_name,
                               cellid=None)
        self.trait_list.append((trait_ob, dataset_ob))

def get_species_groups():

    species_query = "SELECT SpeciesId, MenuName FROM Species"
    species_ids_and_names = g.db.execute(species_query).fetchall()

    species_and_groups = []
    for species_id, species_name in species_ids_and_names:
        this_species_groups = {}
        this_species_groups['species'] = species_name
        groups_query = "SELECT InbredSetName FROM InbredSet WHERE SpeciesId = %s" % (species_id)
        groups = [group[0] for group in g.db.execute(groups_query).fetchall()]

        this_species_groups['groups'] = groups
        species_and_groups.append(this_species_groups)

    return species_and_groups