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

Description Usage Arguments Details Value Author(s) Examples

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)


Search within the AlgebraicHaploPackage package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.