This function computes the boundary of the decision region in a manner that can be employed in the field, as a table, for example. See section 4.2 of the reference below.
computeBoundary(b1, b0, p, glrTables = NULL, tol=1e-7)The acceptance boundary value (corresponds to the boundary \(b_1\) in the appendix of reference)
The rejection boundary value (corresponds to the boundary \(b_0\) in the appendix of reference)
The vector of probabilities, \((p_0, p_1)\) with \(p_0 < p_1\).
A previously computed set of likelihood functions, to speed up computation for the same hypothesis testing problem. Otherwise, it is computed ab initio, resulting in a longer running time.
A numerical tolerance, defaults to 1e-7
The upper boundary that indicates rejection of the null hypothesis
The upper boundary that indicates acceptance of the null hypothesis
The estimated \(\alpha\) and \(\beta\) values corresponding to the two boundaries
This essentially computes the probabilities of hitting the boundaries using a recursion.
Mei-Chiung Shih, Tze Leung Lai, Joseph F. Heyse, and Jie Chen. Sequential Generalized Likelihood Ratio Tests for Vaccine Safety Evaluation (Statistics in Medicine, Volume 29, issue 26, p.2698-2708, 2010.)
Please also consult the website http://med.stanford.edu/biostatistics/ClinicalTrialMethodology/ for further developments.
See Also glrSearch
# NOT RUN {
computeBoundary(b1=2.8, b0=3.3, p=c(.5, .75))
# }
Run the code above in your browser using DataLab