Nothing
EEw2s3 <- function(des='',randomize=FALSE) {
if (des=='') {
cat(" ", "\n")
cat("Catalog of D-efficient Estimation Equivalent RS","\n")
cat(" Designs for (2 wp factors and 3 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("EE24R8WP", width=11),format(" 8",width=11),format(" 3",width=21),"\n")
cat(format("EE28R7WP", width=11),format(" 7",width=11),format(" 4",width=21),"\n")
cat(format("EE32R8WP", width=11),format(" 8",width=11),format(" 4",width=21),"\n")
cat(format("EE35R7WP", width=11),format(" 7",width=11),format(" 5",width=21),"\n")
cat(format("EE40R8WP", width=11),format(" 8",width=11),format(" 5",width=21),"\n")
cat(format("EE42R7WP", width=11),format(" 7",width=11),format(" 6",width=21),"\n")
cat(format("EE48R8WP", width=11),format(" 8",width=11),format(" 6",width=21),"\n")
cat(" ","\n")
cat("==> to retrieve a design type EEw2s3('EE24R8WP') etc.","\n")
} else if (des=='EE21R7WP') {v <- c(1191, 183, 435, 492, 1311, 807, 969, 150, 1158, 1278, 1530, 522, 47, 1118, 1559, 1427, 1679, 41, 637, 196, 1645)
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),s3=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=='EE24R8WP') { v <- c(161, 1889, 449, 1298, 362, 1226, 11, 1451, 1739, 284, 860, 1868, 981, 1341, 621, 1902, 462, 174, 1351, 1351, 631, 144, 1728, 1152)
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),s3=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=='EE28R7WP') { v <- c(1569, 1317, 561, 57, 170, 1682, 422, 1178, 66, 1452, 507, 948, 1579, 193, 697, 508, 1153, 1657, 145, 397, 1287, 1539, 846, 27, 1694, 497, 623, 371)
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),s3=c(-1,0,1))
EE <- Full[v, ]
rownames(EE)<- c(1:28)
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=='EE32R8WP') {v <- c(1129, 1417,1921, 409, 1442, 2,578, 1730, 1315, 163, 1891, 451, 844, 1564, 1924, 556, 1701, 1269, 117, 1773, 438, 1878, 150, 1302, 1511, 647, 1871, 71, 1360, 208, 1000, 496)
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),s3=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=='EE35R7WP') {v <- c(1639, 1135, 127, 379, 820, 1129, 1381, 184, 1255, 436, 990, 1557, 45, 549, 1305, 508, 1516,
1264, 4, 823, 1489, 1615, 355, 607, 418, 1532, 1280, 83, 209, 524, 245, 1694, 749, 1190,
497 )
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),s3=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=='EE40R8WP') {v <- c(1433, 1721, 1793, 65, 641, 1490, 122, 1562, 626, 1130, 1875, 147, 651, 1587, 435, 1748, 1460,
884, 596, 20, 1357, 1213, 1933, 277, 205, 1894, 166, 1318, 454, 886, 295, 1231, 1447, 7,
1735, 1832, 1328, 392, 464, 824)
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),s3=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=='EE42R7WP') {v <- c(1310, 1562, 50, 932, 1373, 491, 1682, 611, 1241, 548, 422, 170, 1641, 1011, 192, 129, 1137,
1263, 1152, 1656, 396, 333, 711, 81, 1552, 40, 1300, 1426, 796, 544, 1518, 510, 6, 1266,
636, 1455, 1197, 1323, 1701, 1008, 315, 189)
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),s3=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=='EE48R8WP') {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)
w1<-c(-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, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 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, 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)
s1<-c(0.56084, -0.06357, 1.00000, 0.51562, -0.06737, 0.70429, -0.11015, -0.95967, -0.20421,
-0.35938, 0.99579, -0.99950, 1.00000, -0.86437, -0.85703, -0.09420, -0.92222, 0.88506,
0.56986, -0.21668, -0.07369, 0.91509, 0.91019, 0.86681, 1.00000, 1.00000, 1.00000,
0.09854, 0.15793, 0.17613, -0.12690, -0.90367, -0.66383, 0.10124, 0.92932, -0.97328,
-0.78345, -1.00000, 0.15928, 0.75899, -0.97078, 0.98320, -0.48317, 0.68254, -0.29762,
-0.43140, 1.00000, 0.57906)
s2<-c(-0.92854, 0.10102, -1.00000, -1.00000, 0.18066, -0.90389, -0.92912, -0.15766, -0.95595,
0.90964, 0.63698, -0.23680, 0.91050, -0.59467, -0.97021, -0.55536, 0.74053, -0.89746,
-1.00000, -1.00000, -0.99894, 0.82650, 0.86108, 1.00000, -1.00000, -1.00000, -1.00000,
0.99876, -1.00000, -1.00000, 0.96027, -0.84218, 0.62294, -0.45932, -0.09132, -0.92330,
-0.78227, 0.84761, -0.99991, 0.73327, -0.92855, -0.23682, -0.42720, 0.98957, 0.54908,
-0.42988, 0.87858, 0.99915)
s3<-c(-0.92486, 0.97572, 0.07996, -0.56142, 0.97940, -0.93836, 0.96591, -0.24995, -0.91156,
-0.89274, 0.99626, 0.99349, -0.99128, -0.72718, -0.97457, -0.32284, 0.90005, 0.75505,
-0.92204, 0.58001, 0.52779, -0.88350, -0.84440, 0.95054, -0.00586, 0.06026, 0.02739,
1.00000, 0.32016, 0.10662, -1.00000, 0.90604, 0.87733, 0.11895, 0.99909, -1.00000,
0.51880, -1.00000, -0.89626, 1.00000, -0.05022, -0.93309, 1.00000, -0.75242, 0.97294,
0.99664, -1.00000, -0.75579)
EE<-data.frame(WP,w1,w2,s1,s2,s3)
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 cat(" Design name misspelled-Enter EEw2s3( ) 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.