Chose the most likely solution of the three potential onces.

Share:

Description

Starting with a 3x3 matrix, three potential haplotypes 2x2 matrices will be calculted and evaluated. The most likly one is choosen. A dicrete solution values are not enforced. This reduces increases the right prediction on real data.

Usage

1
findoptimal(A, B, C, D, mmorg, exact = 1e-05)

Arguments

A

Firts entry of the 2x2 matrix.

B

Second numeric entry of teh 2x2 matrix.

C

Third numeric entry of the 2x2 matrix.

D

Last numeric entry of the 2x2 matrix.

mmorg

3x3 matrix of the pair of original snip pairs.

exact

exact is a parameter when data are assumed to be equal. Acutually 1e-5 is taken. Should not be larger than the invers of four times the number of people.

Details

It chose the 2x2 model of haplotypes with the smallest prediction error.

Value

AA

Coefficient of x^3

BB

Coefficient of x^2

CC

Coefficient of x

DD

Coefficient of the intercept

Author(s)

Jan wolfertz

Examples

1
2
3
4
5
dd2=matrix(c(4,0,0,0,30,0,0,0,23),ncol=3,byrow=TRUE)
A=c(38+0i,2+12.1655i,2-12.1655i)
B=c(0+0i,36-12.1655i,36+12.1655i)
C=c(0+0i,36-12.1655i,36+12.1655i)
D=c(76+0i,40+12.1655i,40-12.1655i)