Nothing
checkResult <- function (result, refDesign, model) {
allDesigns <- result@misc$designs
df <- data.frame(run=numeric(0), step=numeric(0), crit=numeric(0))
for (i in 1:length(allDesigns)) {
designs <- allDesigns[[i]]
for (j in 1:length(designs)) {
design <- designs[[j]]
crit <- sum(general.carryover(design, model=model)$Var.trt.pair)/2
df <- rbind(df, data.frame(run=i, step=j, crit=crit))
}
}
plot <- ggplot(df, aes(x=step, y=crit)) + geom_point() + facet_wrap( ~ run)
print(plot)
return(df)
}
test.2v.designs <- function() {
if (!"extended" %in% strsplit(Sys.getenv("CROSSOVER_UNIT_TESTS"),",")[[1]]) {
cat("Skipping design tests for v=2.\n")
return()
}
design1 <- t(rbind(c(1,1,2,2),
c(2,2,1,1),
c(1,1,2,2),
c(2,2,1,1),
c(1,2,2,1),
c(2,1,1,2)))
design2 <- t(rbind(c(1,1,2,1),
c(2,2,1,2),
c(1,1,2,1),
c(2,2,1,2),
c(1,2,2,1),
c(2,1,1,2)))
design3 <- t(rbind(c(1,1,2,2),
c(2,2,1,1),
c(1,1,2,2),
c(2,2,1,1),
c(1,1,2,2),
c(2,2,1,1)))
design4 <- t(rbind(c(1,1,2,2),
c(2,2,1,1),
c(1,1,2,2),
c(2,2,1,1),
c(1,1,2,2),
c(2,2,1,1))) # Why is there no design 4 in the book?
design5 <- t(rbind(c(1,2,2,2),
c(2,1,1,1),
c(1,2,2,2),
c(2,1,1,1),
c(1,1,2,2),
c(2,2,1,1)))
design6 <- t(rbind(c(1,2,1,2),
c(2,1,2,1),
c(1,1,2,1),
c(2,2,1,2),
c(1,2,1,2),
c(2,1,2,1)))
design7 <- t(rbind(c(1,1,2,1),
c(2,2,1,2),
c(1,1,1,2),
c(2,2,2,1),
c(1,2,1,1),
c(2,1,2,2)))
design8 <- t(rbind(c(1,2,2,2),
c(2,1,1,1),
c(1,1,2,2),
c(2,2,1,1),
c(1,2,1,2),
c(2,1,2,1)))
designs <- list(design1, design2, design3, design4, design5, design6, design7, design8)
models <- c("Standard additive model",
"Self-adjacency model",
"Proportionality model",
"Placebo model",
"No carry-over into self model",
"Treatment decay model",
"Full set of interactions",
"Second-order carry-over effects")
s <- 6; p <- 4; v <- 2
results <- list()
for (i in c(1,2,3,4,5,6)) {
model <- models[i]
cat("======= ", model, " =======","\n")
result <- searchCrossOverDesign(s=s, p=p, v=v, model=model, start.designs=list(designs[[i]]))
print(result)
print(plot(result))
print(general.carryover(designs[[i]], model=i))
checkResult(result, designs[[i]], model=i)
}
estimable(design3, v, model=3)
searchCrossOverDesign(s=s, p=p, v=v, model=3)
searchCrossOverDesign(s=s, p=p, v=v, model=3, start.designs=list(design3))
}
test.2v.designs()
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.