# no interaction effect expected
refmean <- 1
treatgroups <- 4
timepoints <- 5
treateff <- 1.5
timeeff <- 0.85
factors_levels_names <- list(treatment=letters[1:treatgroups], time=1:timepoints)
effects_treat_time <- calculate_mean_matrix(refmean = refmean,
nlfA = treatgroups, nlfB = timepoints,
fAeffect = treateff, fBeffect = timeeff,
label_list = factors_levels_names)
effsize(effects_treat_time)
# we add cell specific interaction effect keeping design and main effect coefficients
cellswithinteraction <- matrix(c(rep(2,3), 3:5), 3,2)
#second level of factor A interacts with 3rd, 4th and 5th level of factor B
effects_treat_time_interact <- calculate_mean_matrix(refmean = refmean,
nlfA = treatgroups, nlfB = timepoints,
fAeffect = treateff, fBeffect = timeeff,
label_list = factors_levels_names,
groupswinteraction = cellswithinteraction,
interact=1.3)
effsize(effects_treat_time_interact)
Run the code above in your browser using DataLab