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
| void set_trinucl_index(void) {
// The 64 possible trinucleotides -- Look at the commented function "print_table_of_all_trinucl" at the bottom of this file
char const *trinuc_arr[] = {
"aaa", "aac", "aag", "aat", "aca",
"acc", "acg", "act", "aga", "agc",
"agg", "agt", "ata", "atc", "atg",
"att", "caa", "cac", "cag", "cat",
"cca", "ccc", "ccg", "cct", "cga",
"cgc", "cgg", "cgt", "cta", "ctc",
"ctg", "ctt", "gaa", "gac", "gag",
"gat", "gca", "gcc", "gcg", "gct",
"gga", "ggc", "ggg", "ggt", "gta",
"gtc", "gtg", "gtt", "taa", "tac",
"tag", "tat", "tca", "tcc", "tcg",
"tct", "tga", "tgc", "tgg", "tgt",
"tta", "ttc", "ttg", "ttt"
};
// Initialisation of the dictionnary that will give for a trinucleotide its corresponding index in RFDD-Base
hcreate(64);
// Element of the dictionnary of trinucleotides
ENTRY trinuc_dic_e;
int i;
for (i=0; i < 64; i++) {
trinuc_dic_e.key = trinuc_arr[i];
trinuc_dic_e.data = (void*) i+1;
hsearch(trinuc_dic_e, ENTER);
}
} |
Partager