Final experimental results

 

The following experiments were performed by using a Final Algorithm that combines three approximation methods ,mentioned above, and provides minimum run time with minimum LOD error rate.

Input files

Num loci

Num people

Complexity

Complexity threshold

Exact LOD score/likelihood

Approx LOD score/likelihood

Exact run time (sec)

Approx run time (sec)

LOD error rate

datafile27_25_9.dat pedfile27_25_9.dat

27

25

9

6

-0.551546/ 1159.696562

-0.550202/ 564.502678

26.046

3.812

0.24%

datafile40_70_9.dat pedfile40_70_9.dat

40

70

9

6

-14.013844/ 9272.252685

-14.013844/ 5092.260239

98.875

35.4

0%

datafile76_20_7.dat pedfile76_20_7.dat

76

20

7

6

-11.597888/ 7049.623917

-11.597888/ 3497.377749

14.703

5.422

0%

datafile30_20_8.dat pedfile30_20_8.dat

30

20

8

6

0.009217/ 685.974775

0.009321/ 218.039512

103.437

14.64

1.12%

datafile23_80_8.dat pedfile23_80_8.dat

23

80

8

7

-50.655849/ 8746.564028

-50.693937/ 4392.057199

127

36.4

0.075%

datafile20_22_9.dat pedfile20_22_9.dat

20

22

9

7

-0.675122/ 431.563191

-0.673222/ 257.030192

2317.469

62.484

0.28%

datafile20_21_9.dat pedfile20_21_9.dat

20

21

9

7

0.029953/ 633.82178

0.029262/ 148.016724

5392.064

3.79

2.3%

datafile60_20_8.dat pedfile60_20_8.dat

60

20

8

7

0.001860/ 1484.125933 0.001882/ 341.179421

135.218

40.046 1.18%

datafile63_20_11.dat pedfile63_20_11.dat

63

20

11

7

-0.610024/ 4476.642004 -0.610024/ 2125.180642

1049.812

44.375

0%

datafile30_150_10.dat pedfile30_150_10.dat

30 150 10 7 -28.22995/ 17687.62947 -28.228052/ 5408.491191 400.485 280.2 0.006%

datafile80_20_18.dat pedfile80_20_18.dat

80 20 18 7 -1.788475/ 5161.509564 -1.788475/ 1283.198118 2564.031 219.234 0%

datafile50_70_10.dat pedfile50_70_10.dat

50 70 10 7 -8.22144/ 11316.609554 -8.22144/ 5875.223639 752.797 53.4 0%

datafile23_80_9.dat pedfile23_80_9.dat

23 80 9 7 -24.952202/ 7661.63884 -24.951648/ 3887.824927 26.8 10 0.002%

datafile56_90_10.dat pedfile56_90_10.dat

56 90 10 8 -6.233115/ 16609.07515 -6.233115/ 8757.641992 602.266 200.969 0%

 

    LOD error rate is calculated as:    |(Exact LOD score - Approx LOD score) / Exact LOD score| * 100%

 

The following experiments  show that the maximum LOD score when approximate program is activated is accepted in the same point as when exact program is activated.

Exact max LOD score    - the maximum LOD score that exact program gives as a result.

Approx max LOD score - the maximum LOD score that approximate program gives as a result.

Exact max theta    - the recombination value ( distance) between iterated locus and its left neighbor on the map, that gives the maximum LOD score in exact program activation.

Approx max theta - the recombination value ( distance) between iterated locus and its left neighbor on the map, that gives the maximum LOD score in approximate program activation.

In all this input files there are three evaluations in the interval.

 

Approx max theta

Exact max theta

Approx max LOD score

Exact max LOD score

Num People Num Loci Input files
0.046 0.046

-8.068408

-8.068408

20 76 datafile76_20_7_3.dat pedfile76_20_7_3.dat
0.046 0.046

-37.005818

-36.946805

80 23 datafile23_80_8_3.dat pedfile23_80_8_3.dat

0.1

0.1

1.067207

1.050892

25 23 datafile23_25_7_3.dat pedfile23_25_7_3.dat
0.1 0.1

-6.233115

-6.233115

90 56 datafile56_90_10_3.dat pedfile56_90_10_3.dat

 

Back to main page