bisec function

Bisection algorithm in Beta distribution

Running bisection algorithm to search for a2, the minimizer of (log((a2)^dea1*(1-a2)^dea2-delta))^2

bisec(delta, dea1, dea2, upbd, lowbd)

Arguments

• delta: Test statistic in Han et al. (2012), delta = (ttot1/(ttot1+ttot2))^dea1*(ttot2/(ttot1+ttot2))^dea2;
• dea1: first parameter in Beta distribution (number of events from the first arm)
• dea2: second parameter in Beta distribution (number of events from the second arm)
• upbd: upper bound of a2
• lowbd: lower bound of a2

Returns

a2

Examples

bisec(-74.4824, 33, 98, 1, 0.252)