Nothing
EEw2s2 <- function(des='',randomize=FALSE) {
if (des=='') {
cat(" ", "\n")
cat("Catalog of D-efficient Estimation Equivalent RS","\n")
cat(" Designs for (2 wp factors and 2 sp factors) ","\n")
cat(" ", "\n")
cat(" Jones and Goos, JQT(2012) pp. 363-374","\n")
cat(" ", "\n")
cat(format("Design Name",width=11),format("whole plots",width=11),format("sub-plots/whole plot",width=21),"\n")
cat("----------------------------------------","\n")
cat(format("EE18R9WP", width=11),format(" 9",width=11),format(" 2",width=21),"\n")
cat(format("EE20R10WP", width=11),format(" 10",width=11),format(" 2",width=21),"\n")
cat(format("EE21R7WP", width=11),format(" 7",width=11),format(" 3",width=21),"\n")
cat(format("EE22R11WP", width=11),format(" 11",width=11),format(" 2",width=21),"\n")
cat(format("EE24R8WP", width=11),format(" 8",width=11),format(" 3",width=21),"\n")
cat(format("EE27R9WP", width=11),format(" 9",width=11),format(" 3",width=21),"\n")
cat(format("EE28R7WP", width=11),format(" 7",width=11),format(" 4",width=21),"\n")
cat(format("EE30R10WP", width=11),format(" 10",width=11),format(" 3",width=21),"\n")
cat(format("EE32R8WP", width=11),format(" 8",width=11),format(" 4",width=21),"\n")
cat(format("EE33R11WP", width=11),format(" 11",width=11),format(" 3",width=21),"\n")
cat(format("EE35R7WP", width=11),format(" 7",width=11),format(" 5",width=21),"\n")
cat(format("EE36R9WP", width=11),format(" 9",width=11),format(" 4",width=21),"\n")
cat(format("EE40R8WP", width=11),format(" 8",width=11),format(" 5",width=21),"\n")
cat(format("EE40R10WP", width=11),format(" 10",width=11),format(" 4",width=21),"\n")
cat(format("EE42R7WP", width=11),format(" 7",width=11),format(" 6",width=21),"\n")
cat(format("EE44R11WP", width=11),format(" 11",width=11),format(" 4",width=21),"\n")
cat(format("EE45R9WP", width=11),format(" 9",width=11),format(" 5",width=21),"\n")
cat(format("EE48R8WP", width=11),format(" 8",width=11),format(" 6",width=21),"\n")
cat(format("EE50R10WP", width=11),format(" 10",width=11),format(" 5",width=21),"\n")
cat(format("EE54R9WP", width=11),format(" 9",width=11),format(" 6",width=21),"\n")
cat(format("EE55R11WP", width=11),format(" 11",width=11),format(" 5",width=21),"\n")
cat(format("EE60R10WP", width=11),format(" 10",width=11),format(" 6",width=21),"\n")
cat(format("EE66R11WP", width=11),format(" 11",width=11),format(" 6",width=21),"\n")
cat(" ","\n")
cat("==> to retrieve a design type EEw2s2('EE21R7WP') etc.","\n")
} else if (des=='EE16R8WP') {v <- c(145, 433, 202, 490, 643, 67, 52, 628, 589, 13, 454, 166, 103, 391, 320, 176)
Full <-expand.grid(WP=c(1:8),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:16)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:8),each=2))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:16)
#this randomizes subplots
s<-2
w<-8
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE18R9WP') { v <- c(181, 505, 164, 488, 3, 651, 544, 220, 608, 284, 24, 672, 79, 727, 62, 710, 504, 342)
Full <-expand.grid(WP=c(1:9),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:18)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:9),each=2))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:18)
#this randomizes subplots
s<-2
w<-9
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE20R10WP') { v <- c(1, 721, 622, 262, 23, 743, 64, 784, 185, 545, 566, 206, 287, 647, 88, 808, 609, 249, 590, 410)
Full <-expand.grid(WP=c(1:10),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:20)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:10),each=2))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:20)
#this randomizes subplots
s<-2
w<-10
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE21R7WP') {v <- c(456, 15, 141, 65, 506, 380, 143, 458, 17, 123, 438, 564, 404, 26, 341, 237, 489, 174, 308, 560, 56)
Full <-expand.grid(WP=c(1:7),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:21)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:7),each=3))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:21)
#this randomizes subplots
s<-3
w<-7
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE22R11WP') {v <- c(595, 199, 2, 794, 267, 663, 818, 26, 71, 863, 886, 94, 436, 634, 690, 294, 625, 229, 219, 417, 341, 737)
Full <-expand.grid(WP=c(1:11),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:22)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:11),each=2))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:22)
#this randomizes subplots
s<-2
w<-11
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE24R8WP') {v <- c(49, 193, 553, 578, 434, 74, 643, 139, 499, 436, 580, 76, 525, 237, 165, 198, 54, 558, 15, 519, 375, 480, 408, 48)
Full <-expand.grid(WP=c(1:8),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:24)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:8),each=3))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:24)
#this randomizes subplots
s<-3
w<-8
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE27R9WP') {WP<-c(1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7, 7, 8, 8, 8, 9, 9, 9)
w1<-c(1.0000, 1.0000, 1.0000, 0.0000, 0.0000, 0.0000, 0.7805, 0.7805, 0.7805, 1.0000, 1.0000, 1.0000, -1.0000, -1.0000, -1.0000, -1.0000, -1.0000, -1.0000,
-1.0000, -1.0000, -1.0000, 0.0000, 0.0000, 0.0000, 1.0000, 1.0000, 1.0000)
w2<-c(-1.0000, -1.0000, -1.0000, 1.0000, 1.0000, 1.0000, 0.0000, 0.0000, 0.0000, -1.0000, -1.0000, -1.0000, -1.0000, -1.0000, -1.0000, -1.0000, -1.0000, -1.0000,
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 0.4752, 0.4752, 0.4752)
s1<-c(-1, 1, 0, 1, 0, -1, -1, 0, 1, 0, 1, -1, 0, -1, 1, -1, 0, 1, -1, 1, 0, 0, -1, 1, -1, 0, 1)
s2<-c(-1.000, 0.000, 1.000, 1.000, -1.000, 1.000, -1.000, 0.000, 0.000, 1.000, 0.000, -1.000, 1.000, 0.000, -1.000, 1.000, -1.000, 0.000, -1.000, -1.000,
1.000, -1.000, 1.000, 1.000, -1.000, 0.827, 0.000)
EE<-data.frame(WP,w1,w2,s1,s2)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:9),each=3))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:27)
#this randomizes subplots
s<-3
w<-9
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE28R7WP') {WP<-c(1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 6, 7, 7, 7, 7)
w1<-c(1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1, 0, 0, 0, 0)
w2<-c(1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, -1, -1, -1, -1, 1, 1, 1, 1)
s1<-c(1.0000, -1.0000, 0.0000, 1.0000, 1.0000, -1.0000, 1.0000, -1.0000, 1.0000, -1.0000, -1.0000, 1.0000, -1.0000, -1.0000, 1.0000, 0.0000, 1.0000, -1.0000,
1.0000, 0.0000, -1.0000, -1.0000, 1.0000, 0.0000, -1.0000, -1.0000, 0.0555, 1.0000)
s2<-c(-1, -1, 1, 0, -1, -1, 1, 1, 1, 1, -1, -1, 1, -1, 0, -1, -1, 1, 1, 0, -1, 0, -1, 1, 0, 1, -1, 1)
EE<-data.frame(WP,w1,w2,s1,s2)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:7),each=4))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:28)
#this randomizes subplots
s<-4
w<-7
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE30R10WP') {WP<-c(1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7, 7, 8, 8, 8, 9, 9, 9, 10, 10, 10)
w1<-c(0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0, -1, -1, -1, 1, 1, 1, -1, -1, -1, 0, 0, 0, -1, -1, -1, 1, 1, 1)
w2<-c(0.0000, 0.0000, 0.0000, -1.0000, -1.0000, -1.0000, 1.0000, 1.0000, 1.0000, 0.0000, 0.0000, 0.0000, 1.0000,
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, -1.0000, -1.0000, -1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
1.0000, -0.2242, -0.2242, -0.2242)
s1<-c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, -1, -1, -1, -1, -1, -1, 1)
s2<-c(0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, 1, 1, -1, 0, -1, 1, -1, -1, 1, -1)
EE<-data.frame(WP,w1,w2,s1,s2)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:10),each=3))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:30)
#this randomizes subplots
s<-3
w<-10
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE32R8WP') {v <- c(433, 145, 1, 577, 194, 626, 482, 50, 499, 211, 643, 67, 596, 164, 452, 20, 101, 605, 317, 245, 374, 518, 302, 14, 335, 119,
623, 263, 64, 568, 352, 424)
Full <-expand.grid(WP=c(1:8),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:32)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:8),each=4))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:32)
#this randomizes subplots
s<-4
w<-8
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE30R10WP') {WP<-c(1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7, 7, 8, 8, 8, 9, 9, 9, 10, 10, 10)
w1<-c(0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0, -1, -1, -1, 1, 1, 1, -1, -1, -1, 0, 0, 0, -1, -1, -1, 1, 1, 1)
w2<-c(0.0000, 0.0000, 0.0000, -1.0000, -1.0000, -1.0000, 1.0000, 1.0000, 1.0000, 0.0000, 0.0000, 0.0000, 1.0000,
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, -1.0000, -1.0000, -1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
1.0000, -0.2242, -0.2242, -0.2242)
s1<-c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, -1, -1, -1, -1, -1, -1, 1)
s2<-c(0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, 1, 1, -1, 0, -1, 1, -1, -1, 1, -1)
EE<-data.frame(WP,w1,w2,s1,s2)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:10),each=3))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:30)
#this randomizes subplots
s<-3
w<-10
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE33R11WP') {WP<-c(1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7, 7, 8, 8, 8, 9, 9, 9, 10, 10, 10, 11, 11, 11)
w1<-c(0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, -1, -1, -1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1)
w2<-c(-1, -1, -1, 0, 0, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 0, -1, -1, -1, 0, 0, 0, 1, 1, 1)
s1<-c(-0.24132, -0.74571, 0.85059, 0.77328, 0.15347, -0.94949, -0.97068, 0.26282, 0.57142, 0.66691,
0.92762, -0.99334, -0.02421, -0.50774, -0.51018, -0.93534, 0.13127, 0.66763, -1.00000, 0.99705,
0.99621, 0.63952, 0.33773, -1.00000, 0.68169, -0.99524, 0.91475, -0.51733, -0.50540, 1.00000,
-0.86947, -0.86400, 0.49390)
s2<-c(0.91061, -0.87649, -0.15561, -0.20217, -0.71439, 0.98894, -0.43842, 0.98889, -0.67196, -0.80857,
0.99032, 0.45553, -1.00000, 1.00000, 0.99999, 0.09592, -1.00000, 0.78259, -0.98569, -0.95685,
-0.96119, -0.01310, -0.83059, 0.91607, -0.80435, 0.44163, 1.00000, 0.96132, -0.12044, -0.76850,
1.00000, 1.00000, -1.00000)
EE<-data.frame(WP,w1,w2,s1,s2)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:11),each=3))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:33)
#this randomizes subplots
s<-3
w<-11
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE35R7WP') {v <- c(379, 505, 127, 1, 253, 520, 142, 457, 205, 16, 45, 423, 297, 549, 171, 557, 494, 242, 179, 53, 516, 12,
390, 327, 75, 440, 566, 62, 125, 377, 42, 168, 357, 420, 483)
Full <-expand.grid(WP=c(1:7),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:35)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:7),each=5))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:35)
#this randomizes subplots
s<-5
w<-7
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE36R9WP') {v <- c(1, 487, 649, 163, 164, 650, 488, 2, 219, 543, 705, 57, 274, 679, 355, 112, 239, 725, 77, 563, 420, 15,
339, 582, 70, 475, 637, 394, 674, 26, 512, 188, 135, 297, 702, 378)
Full <-expand.grid(WP=c(1:9),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:36)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:9),each=4))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:36)
#this randomizes subplots
s<-4
w<-9
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE40R8WP') {v <- c(593, 161, 449, 305, 17, 50, 626, 482, 194, 338, 579, 147, 219, 435, 75, 620, 404, 44, 116, 476, 637, 277,
61, 565, 205, 646, 142, 214, 286, 502, 231, 519, 159, 591, 15, 608, 464, 104, 392, 32)
Full <-expand.grid(WP=c(1:8),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:40)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:8),each=5))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:40)
#this randomizes subplots
s<-5
w<-8
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE40R10WP') {WP<-c(1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 6, 7, 7, 7, 7, 8,
8, 8, 8, 9, 9, 9, 9, 10, 10, 10, 10)
w1<-c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1)
w2<-c(1, 1, 1, 1, 0, 0, 0, 0, -1, -1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, -1, -1, -1, -1, 1, 1, 1, 1, -1,
-1, -1, -1, 0, 0, 0, 0, -1, -1, -1, -1)
s1<-c(1.00000, -0.77966, -0.76995, 1.00000, -0.99675, -0.99944, -0.30072, 0.97375, -0.31686,
-0.65113, -0.12419, 1.00000, -0.31122, -0.06396, 0.97780, -0.69480, -0.99438, 0.97565,
-0.99830, -0.30613, 1.00000, -0.99993, -0.99898, 0.99929, 0.98534, 0.97935, -0.10962,
0.99028, -0.99895, 0.99941, 0.99988, -0.99996, 0.60424, 0.57595, 0.55075, 0.57399,
-0.99975, 1.00000, 0.99929, -0.99916)
s2<-c(-0.93578, 1.00000, 1.00000, -0.94037, -0.98676, 0.41478, 0.85200, -0.30739, -0.94247,
-0.93255, 0.65211, 0.44627, -0.73337, -0.95381, 0.96526, -0.05473, -0.99411, 0.00400,
0.98853, -0.02579, -0.99975, 1.00000, -0.99918, 0.99931, -1.00000, -1.00000, 1.00000,
-0.99382, -0.99897, 0.99943, -0.99996, 0.99988, 0.09254, 0.12663, 0.99971, 1.00000,
1.00000, -0.99992, 0.99931, -0.99901)
EE<-data.frame(WP,w1,w2,s1,s2)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:10),each=4))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:40)
#this randomizes subplots
s<-4
w<-10
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE42R7WP') {v <- c(43, 421, 547, 232, 169, 106, 2, 380, 2, 443, 317, 128, 290, 38, 290, 416, 164, 542, 375, 438, 60, 501,
186, 60, 432, 180, 558, 117, 54, 243, 391, 13, 202, 517, 139, 76, 525, 147, 210, 399, 21, 84)
Full <-expand.grid(WP=c(1:7),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:42)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:7),each=6))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:42)
#this randomizes subplots
s<-6
w<-7
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE44R11WP') {WP<-c(1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 6, 7, 7, 7, 7, 8, 8, 8, 8, 9, 9, 9,
9, 10, 10, 10, 10, 11, 11, 11, 11)
w1<-c(0, 0, 0, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, -1, -1, 1, 1, 1, 1, 0, 0, 0,
0, -1, -1, -1, -1, 1, 1, 1, 1)
w2<-c(1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 1, 1, -1, -1, -1, -1)
s1<-c(0.18532, -0.84429, 0.96867, -0.94704, 0.99218, 0.98611, -0.99838, -0.98064, 0.98252, 0.99579, -0.99042,
-0.98862, 0.97823, -0.98030, 1.00000, -0.99866, -0.76197, 1.00000, 0.11848, -0.99385, 0.13043, -0.76756,
-0.99414, 0.99393, 0.98086, 0.99739, -0.97939, -0.99959, 0.49499, 0.70847, -0.94145, 0.45942, 0.29916,
-1.00000, 0.18823, 0.18850, -0.51531, 0.57717, 0.24671, 0.24820, -0.98461, 0.24222, 1.00000, -0.99714)
s2<-c(0.77708, 0.90810, 0.01244, -0.99856, 0.90220, -0.98479, -0.99650, 0.99940, 0.91658, -0.98187, -1.00000,
0.98560, 0.91392, 0.98818, -0.98345, -0.99834, 0.76129, -0.12240, 0.98649, -0.92632, -0.96962, 0.71739,
-0.03351, 0.98480, 0.90914, -0.98216, 0.99276, -0.99943, -0.97279, 0.80294, -0.72137, -0.98347, -0.96515,
-0.51817, -1.00000, -0.99406, -0.99344, -0.83605, -0.99950, -0.99656, -0.30913, 0.90542, -0.99123, -0.25718)
EE<-data.frame(WP,w1,w2,s1,s2)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:11),each=4))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:44)
#this randomizes subplots
s<-4
w<-11
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE45R9WP') {v <- c(271, 595, 28, 676, 190, 218, 56, 542, 704, 380, 183, 669, 345, 507, 21, 616, 211, 697, 49, 292, 482, 563,
644, 77, 239, 150, 69, 717, 555, 474, 97, 664, 16, 421, 502, 224, 710, 62, 386, 548, 9, 495, 414, 576,
171)
Full <-expand.grid(WP=c(1:9),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:45)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:9),each=5))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:45)
#this randomizes subplots
s<-5
w<-9
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE48R8WP') {v <- c(625, 49, 337, 481, 193, 193, 450, 18, 450, 162, 594, 306, 147, 3, 579, 147, 219, 507, 476, 44, 620, 188,
260, 548, 445, 373, 517, 85, 373, 13, 30, 606, 462, 102, 174, 390, 207, 495, 351, 63, 639, 351, 648, 288,
504, 72, 216, 576)
Full <-expand.grid(WP=c(1:8),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:48)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:8),each=6))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:48)
#this randomizes subplots
s<-6
w<-8
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE50R10WP') {v <- c(61, 421, 241, 601, 781, 482, 662, 212, 32, 572, 643, 283, 733, 193, 13, 134, 584, 764, 44, 494, 565, 475,
205, 25, 655, 726, 186, 546, 276, 96, 787, 247, 67, 607, 427, 268, 628, 448, 88, 808, 119, 569, 749, 299,
209, 90, 630, 270, 810, 450)
Full <-expand.grid(WP=c(1:10),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:50)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:10),each=5))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:50)
#this randomizes subplots
s<-5
w<-10
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE54R9WP') {v <- c(505, 667, 181, 19, 505, 343, 704, 137, 56, 299, 542, 218, 3, 570, 246, 165, 651, 165, 13, 661, 94, 499,
418, 499, 545, 626, 59, 221, 59, 464, 213, 618, 51, 375, 699, 294, 565, 241, 79, 727, 565, 241, 8, 413,
494, 170, 89, 656, 540, 702, 135, 54, 459, 378)
Full <-expand.grid(WP=c(1:9),w1=c(-1,0,1),w2=c(-1,0,1),s1=c(-1,0,1),s2=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:54)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:9),each=6))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:54)
#this randomizes subplots
s<-6
w<-9
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE55R11WP') {WP<-c(1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 6, 6, 6, 6,
6, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11)
w1<-c(1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, -1,
-1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, -1, -1, -1, -1, -1)
w2<-c(-1, -1, -1, -1, -1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1,
1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0)
s1<-c(0.98429, -0.82507, -0.98451, 0.99750, 0.04071, -0.87034, -0.99475, -0.01204, 0.98776,
0.99778, 0.99995, -0.91951, 0.09887, -1.00000, 0.92910, 0.43305, -0.94788, 0.72799,
-1.00000, 0.99976, -0.87261, 0.09876, 1.00000, 0.95365, -0.96688, -0.87350, -0.02119,
-0.98713, 0.99174, 0.99849, -0.59677, 0.81325, 0.77006, -1.00000, -1.00000, 0.66729,
1.00000, 1.00000, 0.93112, 0.14840, 0.00000, -1.00000, 1.00000, 1.00000, -1.00000,
0.26815, -0.91270, -0.96677, 0.21395, -0.97235, 0.00000, 1.00000, 1.00000, -1.00000,
-1.00000)
s2<-c(0.91021, 0.97557, -0.98402, -0.99577, -0.02061, 0.88528, -0.94831, -0.02324, -0.99669,
0.99857, -0.33466, -0.74973, -1.00000, 1.00000, 1.00000, -0.98809, -0.95946, 0.91483,
1.00000, -0.08190, -1.00000, -0.12157, 0.99462, -0.92560, 0.93793, 0.87980, 0.02329,
-0.97363, -0.99471, 0.98086, -0.91112, -0.94322, -0.95956, 0.88470, 0.87411, 0.74084,
-0.95482, -0.95997, -1.00000, 0.78796, -1.00000, 0.00000, 1.00000, -1.00000, 1.00000,
-0.89773, 0.88177, 0.89081, -0.88393, 0.88409, -1.00000, -1.00000, 1.00000, 0.00000,
1.00000)
EE<-data.frame(WP,w1,w2,s1,s2)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:11),each=5))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:55)
#this randomizes subplots
s<-5
w<-11
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE60R10WP') {WP<-c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8,
8, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10)
w1<-c(-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 0)
w2<-c(0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1,
1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 0)
s1<-c(-0.25288, -0.72584, 0.83486, 0.36873, -0.98195, 0.90698, -0.31684, -1.00000, 1.00000, -1.00000, 1.00000, -0.04849, -0.80187, 0.93467, -0.93572,
0.99053, 0.15517, -0.99421, -0.92930, 1.00000, 0.99966, -0.98391, -0.99357, 1.00000, -0.73483, 0.10336, 0.99975, -0.97198, 0.94820, -0.99593,
0.97861, -0.11000, -0.97978, -0.01820, 0.99953, -1.00000, 0.99962, 1.00000, -0.92972, 1.00000, -0.97702, -1.00000, -1.00000, -0.54625, 0.08232,
0.48161, 0.98423, -0.95217, -1.00000, -0.20007, 0.85249, 1.00000, -1.00000, 1.00000, 0.96057, -1.00000, 0.99870, -1.00000, -0.09485, 0.00575)
s2<-c(-0.71040, -0.79763, 0.54092, 0.86236, 0.90678, -0.99326, -0.38804, -0.79542, 0.29917, 0.95800, -0.98796, 1.00000, 0.29070, 0.75794, -0.99059,
-0.20564, -1.00000, 0.99095, 0.99724, -0.99715, 0.97250, 0.07376, -0.99933, -1.00000, 0.90590, -0.09903, -0.99837, 0.12662, 0.90729, -0.99905,
0.98608, -0.15295, 0.98263, 0.02622, -0.99723, -1.00000, -0.96872, -1.00000, -1.00000, 1.00000, 1.00000, 0.01574, -0.71426, 0.90832, 0.99987,
-0.90413, -0.94444, -0.30111, -0.71224, 0.94995, -0.95141, -0.42612, -0.88513, 1.00000, -0.11086, 0.97592, 0.99732, -0.01926, -0.99837, -1.00000)
EE<-data.frame(WP,w1,w2,s1,s2)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:10),each=6))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:60)
#this randomizes subplots
s<-6
w<-10
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else if(des=='EE66R11WP') {WP<-c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8,
8, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11)
w1<-c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0,
0, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1)
w2<-c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1,
1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1, 1)
s1<-c(-0.08536, 0.04859, -0.97256, 0.96469, -0.99999, 1.00000, -0.97123, -0.06605, 0.98627, -0.98369, -0.00895, 0.99901, -0.03921, 0.97856, -0.99437,
-0.97307, 0.99560, -0.01215, -0.03334, 0.96042, -0.99434, 0.99813, -0.98866, 0.01316, 0.45629, -0.48804, 0.84232, -0.85490, -1.00000, 0.99969,
0.04630, 0.07748, -0.31119, -0.19387, 0.46002, -0.31104, -1.00000, 0.70299, 1.00000, 0.38519, 1.00000, -0.96757, -0.79257, -0.76546, -0.88832,
0.52015, -0.70245, -0.80991, 1.00000, -0.84982, 0.97308, -0.95245, -0.05020, 1.00000, 0.89930, 0.88143, -0.95362, -1.00000, 0.86108, 0.84376,
-0.82466, -1.00000, -0.82337, -0.98441, 0.67493, -0.98599)
s2<-c(-0.41563, 0.26767, -0.91084, -0.95726, 0.97809, 1.00000, 0.98496, -0.01963, -0.99834, -0.99104, -0.01267, 0.99875, -0.03929, 1.00000, -1.00000,
0.98629, -0.98627, 0.00130, -0.99050, -0.03640, -0.99788, 0.00250, 0.98432, 0.99998, -0.22692, 0.99952, 0.99991, -0.99631, 0.12799, -0.94216,
0.77437, 0.79950, -0.95416, -0.96031, -0.13838, -0.99194, -0.99358, -0.94996, -0.99574, 0.01866, 0.99006, 0.81201, -0.33748, 0.23941, -0.45839,
0.06763, -0.02703, -0.29384, -1.00000, 0.08272, 0.76364, -0.96491, 1.00000, -1.00000, 0.81487, 0.79824, -0.76645, -0.57502, 0.76975, 0.79419,
0.82274, -0.94014, 0.86408, -0.80296, 0.50191, -0.83181)
EE<-data.frame(WP,w1,w2,s1,s2)
if (randomize==TRUE) {#This orders the whole plots randomly
o<-c(rep(sample(1:11),each=6))
EE<-EE[order(order(o)), ]
rownames(EE)<-c(1:66)
#this randomizes subplots
s<-6
w<-11
rowp<-c(sample(1:s))
for (i in 1:(w-1)) {rowp<-c(rowp,i*s+sample(1:s))}
EE<-EE[order(rowp), ]
}
return(EE)
} else cat(" Design name misspelled-Enter EEw2s2( ) to display list of names","\n")
}
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.