Nothing
test.best <- function() {
data(williams)
expect_true(all(getDesign(searchCrossOverDesign(s=10, v=5, p=5, start.designs=list(williams5t), n=c(200,10)))==williams5t))
}
test.best()
test.evaluation <- function() {
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)))
general.carryover(design1, model=1)
general.carryover(design2, model=1)
general.carryover(design3, model=1)
design.efficiency(design1)
design.efficiency(design2)
design.efficiency(design3)
model <- 1
v <- 2
H <- Crossover:::linkMatrix(model=1, v=2)
for (design in list(design1, design2, design3)) {
rcD1 <- Crossover:::rcd_R(design, v, model=model) # R-Code
rcD2 <- Crossover:::rcd(design, v, model=model) # C-Code
all(rcD1==rcD2)
X1 <- rcD1
#X1 <- Crossover:::getRCDesignMatrix(rcD1, v*v+v) # R-Code
A1 <- Crossover:::infMatrix_R(X1, v, model=model) # R-Code
A1b <- Crossover:::infMatrix_R(X1, v, model=model, method=2) # R-Code
all(A1==A1b)
Csub <- contrMat(n = rep(1, v), type = "Tukey")
C <- as.matrix(cbind(Csub, matrix(0, dim(Csub)[1], v)))
general.carryover(design, model=model)
diag(C %*% ginv(t(H) %*% A1 %*% H) %*% t(C))
diag(ginv(t(H) %*% A1 %*% H) %*% t(C) %*% C)
A2 <- Crossover:::infMatrix(X1, v, model=model) # C-Code
A2
rcD1
X1
A1
t(H)%*%A1%*%H
}
general.carryover(design1, model=model)
}
test.evaluation()
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.