binaries/src/clustalw/src/tree/dayhoff.h

   1 /**
   2  * Author: Mark Larkin
   3  *
   4  * Copyright (c) 2007 Des Higgins, Julie Thompson and Toby Gibson.
   5  */
   6 #ifndef DAYHOFFPAMS2_H
   7 #define DAYHOFFPAMS2_H
   8 /* Mark tidy up Nov 2005 */
   9 /* DAYHOFF.H
  10
  11    Table of estimated PAMS (actual no. of substitutions per 100 residues)
  12    for a range of observed amino acid distances from 75.0% (the first entry
  13    in the array), in 0.1% increments, up to 93.0%.
  14
  15    These values are used to correct for multiple hits in protein alignments.
  16    The values below are for observed distances above 74.9%.  For values above
  17    93%, an arbitrary value of 1000 PAMS (1000% substitution) is used.
  18
  19    These values are derived from a Dayhoff model (1978) of amino acid
  20    substitution and assume average amino acid composition and that amino
  21    acids replace each other at the same rate as in the original Dayhoff model.
  22
  23    Up to 75% observed distance, use Kimura's emprical formula to derive
  24    the correction.  For 75% or greater, use this table.  Kimura's formula
  25    is accurate up to about 75% and fails completely above 85%.
  26 */
  27 namespace clustalw
  28 {
  29
  30 int dayhoff_pams[]={
  31   195,   /* 75.0% observed d; 195 PAMs estimated = 195% estimated d */
  32   196,   /* 75.1% observed d; 196 PAMs estimated */
  33                   197,    198,    199,    200,    200,    201,    202,  203,
  34   204,    205,    206,    207,    208,    209,    209,    210,    211,  212,
  35   213,    214,    215,    216,    217,    218,    219,    220,    221,  222,
  36   223,    224,    226,    227,    228,    229,    230,    231,    232,  233,
  37   234,    236,    237,    238,    239,    240,    241,    243,    244,  245,
  38   246,    248,    249,    250,    /* 250 PAMs = 80.3% observed d */
  39                                   252,    253,    254,    255,    257,  258,
  40   260,    261,    262,    264,    265,    267,    268,    270,    271,  273,
  41   274,    276,    277,    279,    281,    282,    284,    285,    287,  289,
  42   291,    292,    294,    296,    298,    299,    301,    303,    305,  307,
  43   309,    311,    313,    315,    317,    319,    321,    323,    325,  328,
  44   330,    332,    335,    337,    339,    342,    344,    347,    349,  352,
  45   354,    357,    360,    362,    365,    368,    371,    374,    377,  380,
  46   383,    386,    389,    393,    396,    399,    403,    407,    410,  414,
  47   418,    422,    426,    430,    434,    438,    442,    447,    451,  456,
  48   461,    466,    471,    476,    482,    487,    493,    498,    504,  511,
  49   517,    524,    531,    538,    545,    553,    560,    569,    577,  586,
  50   595,    605,    615,    626,    637,    649,    661,    675,    688,  703,
  51   719,    736,    754,    775,    796,    819,    845,    874,    907,  945,
  52          /* 92.9% observed; 945 PAMs */
  53   988    /* 93.0% observed; 988 PAMs */
  54 };
  55 }
  56 #endif
  57