AlgebraicHaploPackage (version 1.2)

findoptimal: Chose the most likely solution of the three potential onces.

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

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.

Value

AA

Coefficient of x^3

BB

Coefficient of x^2

CC

Coefficient of x

DD

Coefficient of the intercept

Details

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

Examples

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)