conf_int_search {crctStepdown} R Documentation

## Randomisation test confidence interval limit search

### Description

A multi-variate Robbins-Monroe search process to estimate the upper of lower limits for a confidence set for parameters from a list of fitted model objects.

### Usage

conf_int_search(
fitlist,
data,
actual_tr,
start,
nsteps = 1000,
alpha = 0.05,
plots = TRUE,
cl_var = "cl",
rand_func = NULL,
verbose = TRUE,
type = "rw"
)


### Arguments

 fitlist A list of p fitted mer class model objects data A data frame used to fit the models in fitlist actual_tr A vector of length p with the original point estimates of the treatment effects start A vector of length p with a set of starting values for the search. The uncorrected confidence set limits (tr_eff +/- 2*SE) usually provide a good starting point. nsteps Number of steps for the search process. alpha Numeric value. The process searches for the 100(1-2*alpha) confidence intervals plots Logical value indicating whether to plot the search process. Default to TRUE cl_var String indicating the name of the column in data with the IDs of the clusters rand_func The name of a function that re-randomises the clusters. The function should produce a data frame that identifies the clusters in the treatment group under the new randomisation scheme. The data frame can either have a single column with name cl_var or two columns of cl_var and t identifying the cluster ID and time period a cluster joins the treatment group. If NULL then clusters are randomised in a 1:1 ratio to treatment and control verbose Logical indicating whether to provide verbose output showing progress and estimates type Method of correction: options are "rw" = Romano-Wolf randomisation test based stepdown, "h" = Holm standard stepdown, "b" = Bonferroni, or "none" for no correction.

### Details

A version of the search process proposed by Garthwaite (1996) adapted for multiple limits. Given a set of estimates of the upper or lower, the process calculates the test statistics for the two-sided null hypotheses that the treatment effects equal these values, and then conducts a single iteration randomisation test of the same null hypotheses. The estimates are then updated based on whether the actual test statistic is higher or lower than the randomisation test statistic it would be compared to under the resampling stepdown approach of Romano & Wolf (2005). At the limits of the confidence set all values should be rejected in a two-sided hypothesis test with a family-wise error rate of alpha, which provides a probabilistic basis for the search process. See Watson (2021) for more details.

### Value

A vector of length p with the estimates of the limits

### Examples

out <- twoarm_sim()
data <- out[[1]]
fit1 <- lme4::glmer(y1 ~ treat + (1|cl) ,
data=data,
family="poisson")

fit2 <- lme4::glmer(y2 ~ treat + (1|cl),
data=data,
family="poisson")
fitlist <- list(fit1,fit2)
tr_eff <- rep(NA,2)
for(i in 1:2){
res <- summary(fitlist[[i]])
tr_eff[i] <- res\$coefficients["treat",'Estimate']
}
conf_int_search(fitlist,
data = data,
actual_tr=tr_eff,
start=tr_eff+1,
nsteps=100,
alpha=0.025,
plots = FALSE,
cl_var = "cl",
verbose = FALSE)


[Package crctStepdown version 0.2.1 Index]