exact_pvalue function

P-value for the two exponential comparison in Han et al.(2012)

This function computes the exact p-value from the likelihood ratio test

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

Returns

a2: Beta distribution quantile computed using bisec.R pval: p-value

Examples

exact_pvalue(1, 302.04, 2, 25, 1)