Phasing trios

compute all kinds of p values for exact coreTDT

1 2 3 | ```
pvalue_calculator(y1, y2, n1, n2, theta1 = 0.25, theta2 = 0.5)
loglr_comp(x1, x2, n1, n2, theta1 = 0.25, theta2 = 0.5)
loglr_comp_2side(x1, x2, n1, n2, theta1 = 0.25, theta2 = 0.5)
``` |

`y1` |
integer, N112 |

`y2` |
integer, N122 |

`n1` |
integer, N11 |

`n2` |
integer, N12 |

`theta1` |
float, probaliblity of N112 given N11 |

`theta2` |
float, probaliblity of N122 given N12 |

`x1` |
integer, N112 |

`x2` |
integer, N122 |

`pvalue_pr` |
pvalues computed from probabaility model |

`pvalue_lr` |
pvalues from likelihood ratio test with restricted alternative hypothesis |

`pvalue_lr2` |
pvalues from likelihood ratio test |

Questions? Problems? Suggestions? Tweet to @rdrrHQ 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.