Diseño Box Behnken
#Box-Behnken design
library(rsm)
e1 <- bbd(3, n0 = 3,randomize=FALSE,coding =list(x1~(Time-10)/5,
x2~(Amplitude-40)/10,
x3~(Temperature - 30)/10))
e1 run.order std.order Time Amplitude Temperature
1 1 1 5 30 30
2 2 2 15 30 30
3 3 3 5 50 30
4 4 4 15 50 30
5 5 5 5 40 20
6 6 6 15 40 20
7 7 7 5 40 40
8 8 8 15 40 40
9 9 9 10 30 20
10 10 10 10 50 20
11 11 11 10 30 40
12 12 12 10 50 40
13 13 13 10 40 30
14 14 14 10 40 30
15 15 15 10 40 30
Data are stored in coded form using these coding formulas ...
x1 ~ (Time - 10)/5
x2 ~ (Amplitude - 40)/10
x3 ~ (Temperature - 30)/10write.csv(x = e1, file = "e1.csv")
e2 <- read.csv(file = "e2.csv", header = T)
e2
coded.data(e2, formulas =list(x1~(Time-10)/5, x2~(Amplitude-40)/10, x3~(Temperature - 30)/10)) run.order std.order Time Amplitude Temperature Y
1 1 1 5 30 30 NA
2 2 2 15 30 30 NA
3 3 3 5 50 30 NA
4 4 4 15 50 30 NA
5 5 5 5 40 20 NA
6 6 6 15 40 20 NA
7 7 7 5 40 40 NA
8 8 8 15 40 40 NA
9 9 9 10 30 20 NA
10 10 10 10 50 20 NA
11 11 11 10 30 40 NA
12 12 12 10 50 40 NA
13 13 13 10 40 30 NA
14 14 14 10 40 30 NA
15 15 15 10 40 30 NA
16 1 1 5 30 30 NA
17 2 2 15 30 30 NA
18 3 3 5 50 30 NA
19 4 4 15 50 30 NA
20 5 5 5 40 20 NA
21 6 6 15 40 20 NA
22 7 7 5 40 40 NA
23 8 8 15 40 40 NA
24 9 9 10 30 20 NA
25 10 10 10 50 20 NA
26 11 11 10 30 40 NA
27 12 12 10 50 40 NA
28 13 13 10 40 30 NA
29 14 14 10 40 30 NA
30 15 15 10 40 30 NA
31 1 1 5 30 30 NA
32 2 2 15 30 30 NA
33 3 3 5 50 30 NA
34 4 4 15 50 30 NA
35 5 5 5 40 20 NA
36 6 6 15 40 20 NA
37 7 7 5 40 40 NA
38 8 8 15 40 40 NA
39 9 9 10 30 20 NA
40 10 10 10 50 20 NA
41 11 11 10 30 40 NA
42 12 12 10 50 40 NA
43 13 13 10 40 30 NA
44 14 14 10 40 30 NA
45 15 15 10 40 30 NA
Data are stored in coded form using these coding formulas ...
x1 ~ (Time - 10)/5
x2 ~ (Amplitude - 40)/10
x3 ~ (Temperature - 30)/10library (rsm)
ccf <- ccd(basis = 3, randomize = F, n0 = c(3,3), alpha="faces",
coding =list(x1~(Time-10)/5,
x2~(Amplitude-40)/10,
x3~(Temperature - 30)/10),
)
ccf run.order std.order Time Amplitude Temperature Block
1 1 1 5 30 20 1
2 2 2 15 30 20 1
3 3 3 5 50 20 1
4 4 4 15 50 20 1
5 5 5 5 30 40 1
6 6 6 15 30 40 1
7 7 7 5 50 40 1
8 8 8 15 50 40 1
9 9 9 10 40 30 1
10 10 10 10 40 30 1
11 11 11 10 40 30 1
12 1 1 5 40 30 2
13 2 2 15 40 30 2
14 3 3 10 30 30 2
15 4 4 10 50 30 2
16 5 5 10 40 20 2
17 6 6 10 40 40 2
18 7 7 10 40 30 2
19 8 8 10 40 30 2
20 9 9 10 40 30 2
Data are stored in coded form using these coding formulas ...
x1 ~ (Time - 10)/5
x2 ~ (Amplitude - 40)/10
x3 ~ (Temperature - 30)/10##CCF
library (rsm)
rq.ccf <- ccd(basis = 3, randomize = F, n0 = c(3,3), alpha="faces",
coding =list(x1~(Time-10)/5,
x2~(Amplitude-40)/10,
x3~(Temperature - 30)/10),
wbreps=c(3,3))
rq.ccf run.order std.order Time Amplitude Temperature Block
1 1 1 5 30 20 1
1.1 2 2 5 30 20 1
1.2 3 3 5 30 20 1
2 4 4 15 30 20 1
2.1 5 5 15 30 20 1
2.2 6 6 15 30 20 1
3 7 7 5 50 20 1
3.1 8 8 5 50 20 1
3.2 9 9 5 50 20 1
4 10 10 15 50 20 1
4.1 11 11 15 50 20 1
4.2 12 12 15 50 20 1
5 13 13 5 30 40 1
5.1 14 14 5 30 40 1
5.2 15 15 5 30 40 1
6 16 16 15 30 40 1
6.1 17 17 15 30 40 1
6.2 18 18 15 30 40 1
7 19 19 5 50 40 1
7.1 20 20 5 50 40 1
7.2 21 21 5 50 40 1
8 22 22 15 50 40 1
8.1 23 23 15 50 40 1
8.2 24 24 15 50 40 1
11 25 25 10 40 30 1
21 26 26 10 40 30 1
31 27 27 10 40 30 1
110 1 1 5 40 30 2
22 2 2 15 40 30 2
32 3 3 10 30 30 2
41 4 4 10 50 30 2
51 5 5 10 40 20 2
61 6 6 10 40 40 2
71 7 7 5 40 30 2
81 8 8 15 40 30 2
9 9 9 10 30 30 2
10 10 10 10 50 30 2
111 11 11 10 40 20 2
12 12 12 10 40 40 2
13 13 13 5 40 30 2
14 14 14 15 40 30 2
15 15 15 10 30 30 2
16 16 16 10 50 30 2
17 17 17 10 40 20 2
18 18 18 10 40 40 2
19 19 19 10 40 30 2
20 20 20 10 40 30 2
211 21 21 10 40 30 2
Data are stored in coded form using these coding formulas ...
x1 ~ (Time - 10)/5
x2 ~ (Amplitude - 40)/10
x3 ~ (Temperature - 30)/10