R/cf_reh_exact.R

In EQUIVNONINF: Testing for Equivalence and Noninferiority

cf_reh_exact <- function(X1,X2,X3,alpha,SW,TOL,ITMAX)         {


cat("\n"," X1=",X1,"   X2=",X2,"   X3=", X3,"   alpha=",alpha,
"\n"," SW=",SW, " TOL=",TOL," ITMAX=",ITMAX)

ConfAsy <- asympt_confbs(X1,X2,X3,alpha)                   

thet0 <-  4*ConfAsy[1,1]^2
exactres_l <-  exact_confb_l(X1,X2,X3,alpha, SW,TOL,ITMAX, thet0)
C_l_exact <- exactres_l[1,1]
prob_ge <- exactres_l[1,2]
IT <- exactres_l[1,3]

thet0 <- 4*ConfAsy[1,3]^2
exactres_r <- exact_confb_r(X1,X2,X3,alpha,SW,TOL,ITMAX, thet0)
C_r_exact <- exactres_r[1,1]
prob_le <- exactres_r[1,2]
IT <- exactres_r[1,3]  
cat("\n","\n"," C_l_exact", C_l_exact,"  C_r_exact", C_r_exact)  }