datasets1 <- generateAnovaDatasets(2, 4, 4, 4)
datasets2 <- generateAnovaDatasets(2, 4, 4, 4)
datasets <- bindDatasets(datasets1, datasets2)
data <- datasets$generateData(1)
stopifnot(data$numberOfDatasets == 4,
all.equal(datasets2$generateData(1), datasets$generateData(3),
check.attributes = FALSE),
all.equal(datasets2$sphericalRandomEffects(1), datasets$sphericalRandomEffects(3)),
all.equal(datasets2$createXMatrix(data), datasets$createXMatrix(data)),
all.equal(datasets2$createZMatrix(data), datasets$createZMatrix(data)))
preparedDataset <- prepareMixedEffectDataset(Reaction ~ Days + (Days|Subject), sleepstudy)
datasets1 <- generateMixedEffectDatasets(2, preparedDataset)
datasets2 <- generateMixedEffectDatasets(2, preparedDataset)
datasets <- bindDatasets(datasets1, datasets2)
data <- datasets$generateData(1)
stopifnot(data$numberOfDatasets == 4,
all.equal(datasets2$generateData(1), datasets$generateData(3),
check.attributes = FALSE),
all.equal(datasets2$sphericalRandomEffects(1), datasets$sphericalRandomEffects(3)),
all.equal(datasets2$createXMatrix(data), datasets$createXMatrix(data)),
all.equal(datasets2$createZMatrix(data), datasets$createZMatrix(data)))
Run the code above in your browser using DataLab