lme_permute2 {crctStepdown} | R Documentation |
Returns the test statistic from a specified null hypothesis and model under a single new permutation
lme_permute2(
fitlist,
data,
null_par = rep(0, length(fitlist)),
cl_var = "cl",
rand_func = NULL
)
fitlist |
A list of glm model objects fitted under the null hypotheses |
data |
A data frame containing the data used to fit the models in fitlist |
null_par |
A vector of the same length as fitlist specifying the value(s) of the treatment effect parameter(s) under the null hypotheses |
cl_var |
String specifying the name of the column identifying the clusters/cluster-time |
rand_func |
String of 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 |
A vector of the length of fitlist with the test statistics for each model and null hypothesis
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)
nullfitlist <- list()
for(i in 1:length(fitlist)){
nullfitlist[[i]] <- est_null_model(fitlist[[i]],
data,
tr_var = "treat",
null_par = 0)
}
out <- lme_permute2(nullfitlist,
data=data,
cl_var = "cl")