-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathhmm_pi_a_b.cpp
46 lines (43 loc) · 2.37 KB
/
hmm_pi_a_b.cpp
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
double pre_pi[3] = {1,9.56820074913016E-77,1.81858419938713E-41};
double pre_a[3][3] = {{0.938655874201866,0.128135097379301,0.195148174983384},{0.0347142240780744,0.869569485812218,3.22171776979643E-25},{0.0266299017200585,0.00229541680848012,0.804851825016615}};
double pre_b[8][3] = {
{2.05436603274678E-15,0.0125713759559856,5.05122127258106E-13},
{2.52998394807497E-22,0.114285235965998,2.22910357695258E-27},
{2.27719498162092E-07,0.873138340658857,2.53255883998954E-25},
{0.999999772280499,5.04741915858667E-06,0.00287017823161201},
{5.84207610274566E-21,3.58200708685093E-37,0.98198902187119},
{7.82574783508609E-51,1.45012517083579E-169,0.0151407998966923},
{0,0,0},
{0,0,0}};
double fall_pi[3] = {1,2.02885E-62,1.4382E-69};
double fall_a[3][3]={{0.847339186,0.094980361,0.117705798},{0.075112028,0.883376935,0.029549676},{0.077548786,0.021642704,0.852744525}};
double fall_b[8][3]={
{1.16527E-35,3.66807E-31,0.015534493},
{1.76139E-20,0.000531981,0.166125668},
{5.75741E-07,6.65886E-05,0.81833829},
{0.999946355,0.000510303,1.5493E-06},
{5.3069E-05,0.465112339,3.67054E-12},
{1.87443E-14,0.273749465,1.1375E-11},
{1.3259E-10,0.221482264,3.84995E-25},
{1.53108E-27,0.038547061,1.69949E-30}};
double run_pi[3] = {5.23E-193,1,8.79E-132};
double run_a[3][3] = {{0.911397594,0.004520484,0.117142582},{0.016535984,0.928022794,0.105564851},{0.072066422,0.067456722,0.777292568}};
double run_b[8][3] = {
{0.376659084,2.89177E-36,7.32707E-22},
{0.623340916,0.0000000000177971,0.000132496},
{4.75775E-11,2.94127E-19,0.579184439},
{1.64936E-22,5.24812E-05,0.420680016},
{3.64058E-23,0.329060907,3.04952E-06},
{2.42714E-31,0.497695775,2.88081E-15},
{5.99981E-26,0.173190837,6.12921E-17},
{0,0,0}};
double walk_pi[3] = {0.999999970460279,1.37693917802767E-19,2.95397206340818E-08};
double walk_a[3][3] = {{0.970315824155936,0.136526478038046,0.191485954538908},{0.0261673499961692,0.779338299831794,0.0212671498421256},{0.00351682584789484,0.0841352221301586,0.787246895618965}};
double walk_b[8][3] = {{0,0,0},
{0,0,0},
{3.36988852249168E-11,7.05490068038612E-17,0.497929004197568},
{0.9999999999658,0.00403935434381675,0.502070854220998},
{5.00984836242739E-13,0.995960645656183,1.41581432624145E-07},
{0,0,0},
{0,0,0},
{0,0,0}};