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