exact_pvalue: P-value for the two exponential comparison in Han et al.(2012)
Description
This function computes the exact p-value from the likelihood ratio test
Usage
exact_pvalue(ttot1,ttot2,dea1,dea2,mono)
Arguments
ttot1
total time on test 1
ttot2
total time on test 2
dea1
number of death 1
dea2
number of death 2
mono
0: 2-sided hypothesis: H0: lam1 is equal to lam2; H1: lam1 is not equal to lam2
1: 1-sided hypothesis: H0: lam1 is greater than or equal to lam2; H1: lam1 is less than lam2
2: 1-sided hypothesis: H0: lam1 is less than or equal to lam2; H1: lam1 is greater than lam2
Value
a2: Beta distribution quantile computed using bisec.R
pval: p-value