findoptimal {AlgebraicHaploPackage} R Documentation

## 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.

### 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

Jan wolfertz

### 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)
```

[Package AlgebraicHaploPackage version 1.2 Index]