“…4. Simulation from CTBU (α 1 , α 2 , α 3 , θ 1 , θ 2 ) distribution Koshti and Kamalja [19] developed a Matlab function for simulating random pairs of observations from the CT BU(α 1 , α 2 , α 3 , θ 1 , θ 2 ) distribution when α 1 is set to 0. Following their approach, we have developed a R function, rctbu(α 1 , α 2 , α 3 , θ 1 , θ 2 , n), to generate random samples of size n from the CT BU(α 1 , α 2 , α 3 , θ 1 , θ 2 ) distribution.…”