An approach to the selection of multistratum fractional factorial designs

Size: px
Start display at page:

Download "An approach to the selection of multistratum fractional factorial designs"

Transcription

1 An approach to the selection of multistratum fractional factorial designs Ching-Shui Cheng & Pi-Wen Tsai 11 August, 2008, DEMA2008 at Isaac Newton Institute Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 1 / 32

2 Multistratum fractional factorial designs Split-plot designs Huang, Chen & Voelkel (1998) Bingham & Sitter (1999a,b; 2001) Bingham, Schoen & Sitter (2004) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 2 / 32

3 Multistratum fractional factorial designs Split-plot designs Huang, Chen & Voelkel (1998) Bingham & Sitter (1999a,b; 2001) Bingham, Schoen & Sitter (2004) Blocked split-plot designs McLeod & Brewster (2004) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 2 / 32

4 Multistratum fractional factorial designs Split-plot designs Huang, Chen & Voelkel (1998) Bingham & Sitter (1999a,b; 2001) Bingham, Schoen & Sitter (2004) Blocked split-plot designs McLeod & Brewster (2004) Strip-plot designs (two processing stages) Miller (1997) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 2 / 32

5 Multistratum fractional factorial designs Experiments with multiple processing stages Split-lot designs: Mee & Bates (1998) Bingham, Sitter, Kelly, Moore & Olivas (2008) Ranjan, Bingham & Dean (2008) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 3 / 32

6 Multistratum fractional factorial designs Treatment structure factorial structure Unit (block, plot) structure simple orthogonal block structure (Nelder, 1965) Tjur block structure (Tjur, 1984) orthogonal block structure (Speed and Bailey, 1982; Bailey, 1985) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 4 / 32

7 Simple orthogonal block structures (Nelder 1965) Nesting: block/unit Crossing: row column Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 5 / 32

8 Blocked split-plot designs: McLeod & Brewster (2004) Block structure: 32 runs, 4 weeks/4 days/2 runs. block/whole-plot/subplot Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 6 / 32

9 Blocked split-plot designs: McLeod & Brewster (2004) Block structure: 32 runs, 4 weeks/4 days/2 runs. block/whole-plot/subplot Treatment structure: A B p q r s. Each of the six factors has two levels. A, B: hard to vary (whole-plot factors) p, q, r, s: easy to vary (subplot factors) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 6 / 32

10 Strip-plot designs: Miller (1997) Block structure: 32 runs, 2 block/(4 rows 4 columns) Treatment structure: 2 10 A, B, C, D, E, F: rst-stage factors a,b,c,d: second-stage factors Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 7 / 32

11 Strip-plot designs: Miller (1997) Block structure: 32 runs, 2 block/(4 rows 4 columns) Treatment structure: 2 10 A, B, C, D, E, F: rst-stage factors a,b,c,d: second-stage factors Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 7 / 32

12 Multiple processing stages Split-lot designs: Mee & Bates (1998) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 8 / 32

13 Multiple processing stages Split-lot designs: Mee & Bates (1998) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 8 / 32

14 Multiple processing stages Split-lot designs: Mee & Bates (1998) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 8 / 32

15 Multiple processing stages Split-lot designs: Mee & Bates (1998) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 8 / 32

16 Multiple processing stages Bingham et al (2008): 32 runs; two processing stages First-stage randomization: 8/4 Second-stage randomization: 8/4 Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 9 / 32

17 Multiple processing stages Bingham et al (2008): 32 runs; two processing stages First-stage randomization: 8/4 Second-stage randomization: 8/ v v v v 2 v v v v 3 v v v v 4 v v v v 5 v v v v 6 v v v v 7 v v v v 8 v v v v Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 9 / 32

18 Multiple processing stages Bingham et al (2008): 32 runs; two processing stages First-stage randomization: 8/4 Second-stage randomization: 8/ v v v v 2 v v v v 3 v v v v 4 v v v v 5 v v v v 6 v v v v 7 v v v v 8 v v v v Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 9 / 32

19 Multistratum fractional factorial designs E(y) = Xα, cov(y) = s ξ i S i, i=0 Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 10 / 32

20 Multistratum fractional factorial designs s E(y) = Xα, cov(y) = ξ i S i, where ξ i is an eigenvalue of cov(y), and S i is the orthogonal projection matrix onto the eigenspace S i of cov(y) with eigenvalue ξ i. i=0 Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 10 / 32

21 Multistratum fractional factorial designs s E(y) = Xα, cov(y) = ξ i S i, where ξ i is an eigenvalue of cov(y), and S i is the orthogonal projection matrix onto the eigenspace S i of cov(y) with eigenvalue ξ i. i=0 S 0 = 1 N J N=run size and J is the matrix of all 1's. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 10 / 32

22 Multistratum fractional factorial designs s E(y) = Xα, cov(y) = ξ i S i, where ξ i is an eigenvalue of cov(y), and S i is the orthogonal projection matrix onto the eigenspace S i of cov(y) with eigenvalue ξ i. i=0 S 0 = 1 N J N=run size and J is the matrix of all 1's. Each of the eigenspaces with i 1 is called a stratum. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 10 / 32

23 Multistratum fractional factorial designs s E(y) = Xα, cov(y) = ξ i S i, where ξ i is an eigenvalue of cov(y), and S i is the orthogonal projection matrix onto the eigenspace S i of cov(y) with eigenvalue ξ i. i=0 S 0 = 1 N J N=run size and J is the matrix of all 1's. Each of the eigenspaces with i 1 is called a stratum. ξ i : the ith stratum variance. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 10 / 32

24 Multistratum fractional factorial designs c S i var(c y) = ξ i c 2. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 11 / 32

25 Multistratum fractional factorial designs c S i var(c y) = ξ i c 2. c S i, d S i, i i cov(c y, d y) = 0. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 11 / 32

26 Multistratum fractional factorial designs c S i var(c y) = ξ i c 2. c S i, d S i, i i cov(c y, d y) = 0. Orthogonal designs: Im(X) = j T j, where each T j S i for some i. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 11 / 32

27 Multistratum fractional factorial designs c S i var(c y) = ξ i c 2. c S i, d S i, i i cov(c y, d y) = 0. Orthogonal designs: Im(X) = j T j, where each T j S i for some i. Chapter 10 of Bailey (2008) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 11 / 32

28 Multistratum fractional factorial designs Each stratum corresponds to a partition of the units into subsets of equal size. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 12 / 32

29 Multistratum fractional factorial designs Each stratum corresponds to a partition of the units into subsets of equal size. If the partition corresponding to stratum i is coarser than that of stratum i, then we say that the latter is nested in the former. We assume that in this case, ξ i ξ i. This is true for the usual mixed-eect models. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 12 / 32

30 Blocked split-plot designs: McLeod & Brewster (2004) Block structure: 32 runs, 4 weeks/4 days/2 runs. block/whole-plot/subplot Treatment structure: A B p q r s. Each of the six factors has two levels. A, B: hard to vary (whole-plot factors) p, q, r, s: easy to vary (subplot factors) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 13 / 32

31 Blocked split-plot designs: McLeod & Brewster (2004) Block structure: 32 runs, 4 weeks/4 days/2 runs. block/whole-plot/subplot Treatment structure: A B p q r s. Each of the six factors has two levels. A, B: hard to vary (whole-plot factors) p, q, r, s: easy to vary (subplot factors) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 13 / 32

32 Blocked split-plot designs: McLeod & Brewster (2004) Block structure: 32 runs, 4 weeks/4 days/2 runs. block/whole-plot/subplot Treatment structure: A B p q r s. Each of the six factors has two levels. A, B: hard to vary (whole-plot factors) p, q, r, s: easy to vary (subplot factors) Null ANOVA: Three strata block whole-plot subplots (units) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 13 / 32

33 Blocked split-plot designs: McLeod & Brewster (2004) Block structure: 32 runs, 4 weeks/4 days/2 runs. block/whole-plot/subplot Treatment structure: A B p q r s. Each of the six factors has two levels. A, B: hard to vary (whole-plot factors) p, q, r, s: easy to vary (subplot factors) Null ANOVA: Three strata block whole-plot subplots (units) U B W S Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 13 / 32

34 Blocked split-plot designs: McLeod & Brewster (2004) Block structure: 32 runs, 4 weeks/4 days/2 runs. block/whole-plot/subplot Treatment structure: A B p q r s. Each of the six factors has two levels. A, B: hard to vary (whole-plot factors) p, q, r, s: easy to vary (subplot factors) Null ANOVA: Three strata block whole-plot subplots (units) U 1 B 3 W 12 S 16 Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 13 / 32

35 Blocked split-plot designs: McLeod & Brewster (2004) Block structure: 32 runs, 4 weeks/4 days/2 runs. block/whole-plot/subplot Treatment structure: A B p q r s. Each of the six factors has two levels. A, B: hard to vary (whole-plot factors) p, q, r, s: easy to vary (subplot factors) Null ANOVA: Three strata block whole-plot subplots (units) U 1 B 3 W 12 S 16 ξ B ξ W ξ S Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 13 / 32

36 Strip-plot designs: Miller (1997) Block structure: 32 runs, 2 block/(4 rows 4 columns) Treatment structure: 2 10 A, B, C, D, E, F: rst-stage factors a,b,c,d: second-stage factors Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 14 / 32

37 Strip-plot designs: Miller (1997) Block structure: 32 runs, 2 block/(4 rows 4 columns) Treatment structure: 2 10 A, B, C, D, E, F: rst-stage factors a,b,c,d: second-stage factors U 1 B 1 6 R C E 18 Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 14 / 32

38 Multiple processing stages Split-lot designs: Mee & Bates (1998) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 15 / 32

39 Multiple processing stages Split-lot designs: Mee & Bates (1998) U 1 3 S 1 S S 3 S 4 3 E 3 Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 15 / 32

40 Multiple processing stages Bingham et al (2008): 32 runs; two processing stages First-stage randomization: 8/4 Second-stage randomization: 8/ v v v v 2 v v v v 3 v v v v 4 v v v v 5 v v v v 6 v v v v 7 v v v v 8 v v v v Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 16 / 32

41 Multiple processing stages Bingham et al (2008): 32 runs; two processing stages First-stage randomization: 8/4 Second-stage randomization: 8/ v v v v 2 v v v v 3 v v v v 4 v v v v 5 v v v v 6 v v v v 7 v v v v 8 v v v v U 1 S 1 S S 1 S 2 6 E 18 Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 16 / 32

42 Optimality criteria Existing works use modications of the minimum aberration criterion. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 17 / 32

43 Optimality criteria Existing works use modications of the minimum aberration criterion. Block eects are treated as xed. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 17 / 32

44 Optimality criteria Existing works use modications of the minimum aberration criterion. Block eects are treated as xed. For split-plot designs, the usual minimum aberration criterion for unstructured units is used. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 17 / 32

45 Optimality criteria Existing works use modications of the minimum aberration criterion. Block eects are treated as xed. For split-plot designs, the usual minimum aberration criterion for unstructured units is used. Stratum variances are not taken into account. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 17 / 32

46 New criterion Cheng & Tsai (2008) propose a unied approach to the selection of blocked and split-plot fractional factorial designs. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 18 / 32

47 New criterion Cheng & Tsai (2008) propose a unied approach to the selection of blocked and split-plot fractional factorial designs. This approach can be extended to experiments with more complicated block structures. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 18 / 32

48 Assumptions The three-factor and higher-order interactions are negligible. The main eects must be estimated. The two-factor interactions are equally important. Only two-level regular and orthogonal designs will be considered. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 19 / 32

49 Information capacity In justifying the minimum aberration criterion for designs with unstructured units, Cheng, Steinberg and Sun (1999) showed that it is a good surrogate for maximizing the number of estimable models containing all the main eects and a given number of two-factor interactions. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 20 / 32

50 Information capacity In justifying the minimum aberration criterion for designs with unstructured units, Cheng, Steinberg and Sun (1999) showed that it is a good surrogate for maximizing the number of estimable models containing all the main eects and a given number of two-factor interactions. In the multistratum situation, dierent estimable models may have dierent eciencies. Thus we take the average eciency over the candidate models. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 20 / 32

51 Information capacity Suppose the information for the main eects is specied to lie in certain strata. Then it is enough to consider the eciencies of the estimators of two-factor interactions, measured by D 1/k, where D is the determinant of the information matrix for the two-factor interactions, and k is the number of two-factor interactions in the model. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 21 / 32

52 Information capacity Suppose the information for the main eects is specied to lie in certain strata. Then it is enough to consider the eciencies of the estimators of two-factor interactions, measured by D 1/k, where D is the determinant of the information matrix for the two-factor interactions, and k is the number of two-factor interactions in the model. For a given design d, dene the information capacity I k (d) as the average of D 1/k over the set of models containing all the main eects and k two-factor interactions. The objective is to maximize I k (d). Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 21 / 32

53 Information capacity Suppose the information for the main eects is specied to lie in certain strata. Then it is enough to consider the eciencies of the estimators of two-factor interactions, measured by D 1/k, where D is the determinant of the information matrix for the two-factor interactions, and k is the number of two-factor interactions in the model. For a given design d, dene the information capacity I k (d) as the average of D 1/k over the set of models containing all the main eects and k two-factor interactions. The objective is to maximize I k (d). Sun (1993), Tsai, Gilmour and Mead (2000) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 21 / 32

54 Information capacity I k (d) is time-consuming to calculate and will be replaced by a good surrogate. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 22 / 32

55 Information capacity I k (d) is time-consuming to calculate and will be replaced by a good surrogate. I k (d), or its surrogate, depends on the stratum variances ξ 1,..., ξ s. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 22 / 32

56 Information capacity I k (d) is time-consuming to calculate and will be replaced by a good surrogate. I k (d), or its surrogate, depends on the stratum variances ξ 1,..., ξ s. The partial order among the ξ's can be used to simplify conditions for optimality. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 22 / 32

57 A surrogate Suppose the run size is 2 n p. Then there are a total of 2 n p 1 alias sets. Let f = 2 n p 1 n. Then f is the number of alias sets that do not contain main eects. Suppose that h i of these alias sets fall in the i th stratum. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 23 / 32

58 A surrogate Suppose the run size is 2 n p. Then there are a total of 2 n p 1 alias sets. Let f = 2 n p 1 n. Then f is the number of alias sets that do not contain main eects. Suppose that h i of these alias sets fall in the i th stratum. For each 1 j h i, let m (i) j contained in the jth of these alias sets. be the number of two-factor interactions Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 23 / 32

59 A surrogate Suppose the run size is 2 n p. Then there are a total of 2 n p 1 alias sets. Let f = 2 n p 1 n. Then f is the number of alias sets that do not contain main eects. Suppose that h i of these alias sets fall in the i th stratum. For each 1 j h i, let m (i) j contained in the jth of these alias sets. Let m be the vector (ξ 1/k m (1),..., ξ 1/k m (1) h 1,..., ξs 1/k be the number of two-factor interactions m (s),..., ξ 1/k 1 s m (s) h s ). Then it can be shown that I k (d) is a Schur concave function of m, and is nondecreasing in each component of m. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 23 / 32

60 A surrogate It follws that a good surrogate for maximizing I k (d) is to Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 24 / 32

61 A surrogate It follws that a good surrogate for maximizing I k (d) is to (a) maximize s i=1 ξ 1/k i h i j=1 m (i) j, and then Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 24 / 32

62 A surrogate It follws that a good surrogate for maximizing I k (d) is to (a) maximize (b) minimize s i=1 s i=1 ξ 1/k i ξ 2/k i h i j=1 h i j=1 m (i) j, and then (m (i) j ) 2 among those which maximize (a) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 24 / 32

63 A useful theorem The following is a sucient condition for a design d to be optimal under (a) and (b) with respect to all k and all ξ 1,..., ξ s such that ξ i ξ i if the partition corresponding to stratum i is nested in that of stratum i: Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 25 / 32

64 A useful theorem The following is a sucient condition for a design d to be optimal under (a) and (b) with respect to all k and all ξ 1,..., ξ s such that ξ i ξ i if the partition corresponding to stratum i is nested in that of stratum i: (i) d maximizes i T hi m(i) j=1 j, and Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 25 / 32

65 A useful theorem The following is a sucient condition for a design d to be optimal under (a) and (b) with respect to all k and all ξ 1,..., ξ s such that ξ i ξ i if the partition corresponding to stratum i is nested in that of stratum i: (i) d maximizes i T hi (ii) minimizes hi i T hi i T m(i) j=1 j. m(i) j=1 j (m(i) j=1 j, and ) 2 among those which maximize for all subsets T of {1,..., s} such that if i T, then i T for all i such that the partition corresponding to stratum i is nested in that of stratum i. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 25 / 32

66 A useful theorem The following is a sucient condition for a design d to be optimal under (a) and (b) with respect to all k and all ξ 1,..., ξ s such that ξ i ξ i if the partition corresponding to stratum i is nested in that of stratum i: (i) d maximizes i T hi (ii) minimizes hi i T hi i T m(i) j=1 j. m(i) j=1 j (m(i) j=1 j, and ) 2 among those which maximize for all subsets T of {1,..., s} such that if i T, then i T for all i such that the partition corresponding to stratum i is nested in that of stratum i. This result can be used to rule out inadmissible designs. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 25 / 32

67 Blocked split-plot designs (McLeod & Brewster, 2004) I = pqrs, δ 1 = ABpq, δ 2 = Bpr Block Stratum Source of variation df Bpr 1 Aqr 1 ABpq 1 WP Stratum A 1 B 1 AB 1 qr = ps 1 pr = qs 1 pq = rs 1 Residuals 6 Unit Stratum p 1 q 1 r 1 s 1 Ap 1 Aq 1 Ar 1 As 1 Bp 1 Bq 1 Br 1 Bs 1 Residuals 4 B 3 W 12 S 16 Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 26 / 32

68 Blocked split-plot designs (McLeod & Brewster, 2004) I = pqrs, δ 1 = ABpq, δ 2 = Bpr Block Stratum Source of variation df Bpr 1 Aqr 1 ABpq 1 WP Stratum A 1 B 1 AB 1 qr = ps 1 pr = qs 1 pq = rs 1 Residuals 6 Unit Stratum p 1 q 1 r 1 s 1 Ap 1 Aq 1 Ar 1 As 1 Bp 1 Bq 1 Br 1 Bs 1 Residuals 4 B 3 m (1) i : 000 W 12 m (2) i : S 16 m (3) 1 : Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 26 / 32

69 Blocked split-plot designs (McLeod & Brewster, 2004) I = pqrs, δ 1 = ABpq, δ 2 = Bpr Block Stratum Source of variation df Bpr 1 Aqr 1 ABpq 1 WP Stratum A 1 B 1 AB 1 qr = ps 1 pr = qs 1 pq = rs 1 Residuals 6 Unit Stratum p 1 q 1 r 1 s 1 Ap 1 Aq 1 Ar 1 As 1 Bp 1 Bq 1 Br 1 Bs 1 Residuals 4 B 3 m (1) i : 000 W 12 m (2) i : S 16 m (3) 1 : mi m 2 i (1)+(2)+(3) (2)+(3) (3) 8 8 Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 26 / 32

70 Blocked split-plot designs: 4/4/2 Two optimal 4/4/2 designs by McLeod & Brewster (2004) The m i values for d 1 and d 2 : d 1 m i (1) : 000 d 2 m i (1) : 000 m i (2) : m i (2) : m i (3) : m i (3) : Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 27 / 32

71 Blocked split-plot designs: 4/4/2 Two optimal 4/4/2 designs by McLeod & Brewster (2004) The m i values for d 1 and d 2 : d 1 m i (1) : 000 d 2 m i (1) : 000 m i (2) : m i (2) : m i (3) : m i (3) : d 1 mi m 2 i (1)+(2)+(3) (2)+(3) (3) 8 8 d 2 mi m 2 i Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 27 / 32

72 Blocked split-plot designs: 4/4/2 Two optimal 4/4/2 designs by McLeod & Brewster (2004) The m i values for d 1 and d 2 : d 1 m i (1) : 000 d 2 m i (1) : 000 m i (2) : m i (2) : m i (3) : m i (3) : d 1 mi m 2 i (1)+(2)+(3) (2)+(3) (3) 8 8 d 2 mi m 2 i d 1 dominates d 2 Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 27 / 32

73 Blocked split-plot designs: 4/4/2 d 1 mi m 2 i (1)+(2)+(3) (2)+(3) (3) 8 8 d 2 mi m 2 i Another admissible design d 3 m i (1): 100 m i (2): m i (3): d 3 mi m 2 i Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 28 / 32

74 Blocked split-plot designs: 4/4/2 d 1 mi m 2 i (1)+(2)+(3) (2)+(3) (3) 8 8 d 2 mi m 2 i Another admissible design d 3 m i (1): 100 m i (2): m i (3): d 3 mi m 2 i The optimality between d 1 and d 3 depends on the variance ratios. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 28 / 32

75 Blocked split-plot designs A search of all 32-run designs with various numbers of factors and block structures shows that there is a single admissible design in most of the cases. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 29 / 32

76 Strip-plot designs (Miller 97) Two processing-stage experiment block/(row column) B R C E Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 30 / 32

77 Strip-plot designs (Miller 97) Two processing-stage experiment block/(row column) all strata B R C E Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 30 / 32

78 Strip-plot designs (Miller 97) Two processing-stage experiment block/(row column) all strata without the Block stratum B R C E Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 30 / 32

79 Strip-plot designs (Miller 97) Two processing-stage experiment block/(row column) all strata without the Block stratum without the Block & Column strata B R C E Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 30 / 32

80 Strip-plot designs (Miller 97) Two processing-stage experiment block/(row column) all strata without the Block stratum without the Block & Column strata without the Block & Row strata B R C E Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 30 / 32

81 Strip-plot designs (Miller 97) Two processing-stage experiment block/(row column) all strata without the Block stratum without the Block & Column strata without the Block & Row strata Unit stratum (without Block, Row B R C E & Column strata) Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 30 / 32

82 Strip-plot designs (Miller 97) The design constructed by Miller is one of two admissible designs. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 31 / 32

83 Strip-plot designs (Miller 97) The design constructed by Miller is one of two admissible designs. There is a single admissible design in most of the 32-run cases. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 31 / 32

84 Summary Introduce a unied approach to the selection of multistratum factorial designs. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 32 / 32

85 Summary Introduce a unied approach to the selection of multistratum factorial designs. Provide a useful theorem for nding optimal designs. Cheng & Tsai (DEMA2008) Multi-stratum FF designs Isaac Newton Institute 32 / 32

Multistratum Fractional Factorial Designs. Ching-Shui Cheng and Pi-Wen Tsai. University of California, Berkeley and National Taiwan Normal University

Multistratum Fractional Factorial Designs. Ching-Shui Cheng and Pi-Wen Tsai. University of California, Berkeley and National Taiwan Normal University Multistratum Fractional Factorial Designs Ching-Shui Cheng and Pi-Wen Tsai University of California, Berkeley and National Taiwan Normal University Abstract: Recent work on multistratum fractional factorial

More information

Optimal Two-Level Regular Fractional Factorial Block and. Split-Plot Designs

Optimal Two-Level Regular Fractional Factorial Block and. Split-Plot Designs Optimal Two-Level Regular Fractional Factorial Block and Split-Plot Designs BY CHING-SHUI CHENG Department of Statistics, University of California, Berkeley, California 94720, U.S.A. cheng@stat.berkeley.edu

More information

By Ming-Chung Chang and Ching-Shui Cheng Academia Sinica and University of California, Berkeley

By Ming-Chung Chang and Ching-Shui Cheng Academia Sinica and University of California, Berkeley Submitted to the Annals of Statistics A BAYESIAN APPROACH TO THE SELECTION OF TWO-LEVEL MULTI-STRATUM FACTORIAL DESIGNS By Ming-Chung Chang and Ching-Shui Cheng Academia Sinica and University of California,

More information

Introduction. Chapter 1

Introduction. Chapter 1 Chapter Introduction Many of the fundamental ideas and principles of experimental design were developed by Sir R. A. Fisher at the Rothamsted Experimental Station (Fisher, 926). This agricultural background

More information

Maximal Rank - Minimum Aberration Regular Two-Level Split-Plot Fractional Factorial Designs

Maximal Rank - Minimum Aberration Regular Two-Level Split-Plot Fractional Factorial Designs Sankhyā : The Indian Journal of Statistics 2007, Volume 69, Part 2, pp. 344-357 c 2007, Indian Statistical Institute Maximal Rank - Minimum Aberration Regular Two-Level Split-Plot Fractional Factorial

More information

Automatic analysis of variance for orthogonal plot structures

Automatic analysis of variance for orthogonal plot structures Automatic analysis of variance for orthogonal plot structures Heiko Grossmann Queen Mary, University of London Isaac Newton Institute Cambridge, 13 October 2011 Outline 1. Introduction 2. Theory 3. AutomaticANOVA

More information

FRACTIONAL FACTORIAL SPLIT-PLOT DESIGNS WITH MINIMUM ABERRATION AND MAXIMUM ESTIMATION CAPACITY

FRACTIONAL FACTORIAL SPLIT-PLOT DESIGNS WITH MINIMUM ABERRATION AND MAXIMUM ESTIMATION CAPACITY Statistica Sinica 12(2002), 885-903 FRACTIONAL FACTORIAL SPLIT-PLOT DESIGNS WITH MINIMUM ABERRATION AND MAXIMUM ESTIMATION CAPACITY Rahul Mukerjee and Kai-Tai Fang Indian Institute of Management and Hong

More information

Minimum Aberration and Related Criteria for Fractional Factorial Designs

Minimum Aberration and Related Criteria for Fractional Factorial Designs Minimum Aberration and Related Criteria for Fractional Factorial Designs Hegang Chen Division of Biostatistics and Bioinformatics 660 West Redwood Street University of Maryland School of Medicine Baltimore,

More information

Abstract. The design key in single- and multi-phase experiments. Introduction of the design key. H. Desmond Patterson

Abstract. The design key in single- and multi-phase experiments. Introduction of the design key. H. Desmond Patterson The design key in single- and multi-phase experiments R. A. Bailey University of St Andrews QMUL (emerita) Biometrics by the Harbour: Meeting of the Australasian Region of the International Biometric Society,

More information

SCTT The pqr-method august 2016

SCTT The pqr-method august 2016 SCTT The pqr-method august 2016 A. Doledenok, M. Fadin, À. Menshchikov, A. Semchankau Almost all inequalities considered in our project are symmetric. Hence if plugging (a 0, b 0, c 0 ) into our inequality

More information

A General Criterion for Factorial Designs Under Model Uncertainty

A General Criterion for Factorial Designs Under Model Uncertainty A General Criterion for Factorial Designs Under Model Uncertainty Steven Gilmour Queen Mary University of London http://www.maths.qmul.ac.uk/ sgg and Pi-Wen Tsai National Taiwan Normal University Fall

More information

Some optimal criteria of model-robustness for two-level non-regular fractional factorial designs

Some optimal criteria of model-robustness for two-level non-regular fractional factorial designs Some optimal criteria of model-robustness for two-level non-regular fractional factorial designs arxiv:0907.052v stat.me 3 Jul 2009 Satoshi Aoki July, 2009 Abstract We present some optimal criteria to

More information

A UNIFIED APPROACH TO FACTORIAL DESIGNS WITH RANDOMIZATION RESTRICTIONS

A UNIFIED APPROACH TO FACTORIAL DESIGNS WITH RANDOMIZATION RESTRICTIONS Calcutta Statistical Association Bulletin Vol. 65 (Special 8th Triennial Symposium Proceedings Volume) 2013, Nos. 257-260 A UNIFIED APPROACH TO FACTORIAL DESIGNS WITH RANDOMIZATION RESTRICTIONS PRITAM

More information

A General Criterion for Factorial Designs Under Model Uncertainty

A General Criterion for Factorial Designs Under Model Uncertainty A General Criterion for Factorial Designs Under Model Uncertainty Pi-Wen Tsai Department of Mathematics National Taiwan Normal University Taipei 116, Taiwan, R.O.C. E-mail: pwtsai@math.ntnu.edu.tw Steven

More information

Design of Experiments: III. Hasse diagrams in designed experiments

Design of Experiments: III. Hasse diagrams in designed experiments Design of Experiments: III. Hasse diagrams in designed experiments R. A. Bailey r.a.bailey@qmul.ac.uk 57a Reunião Anual da RBras, May 202 Outline I: Factors II: Models I. Factors Part I Factors pen. Each

More information

Selecting an Orthogonal or Nonorthogonal Two-Level Design for Screening

Selecting an Orthogonal or Nonorthogonal Two-Level Design for Screening Selecting an Orthogonal or Nonorthogonal Two-Level Design for Screening David J. Edwards 1 (with Robert W. Mee 2 and Eric D. Schoen 3 ) 1 Virginia Commonwealth University, Richmond, VA 2 University of

More information

I-Optimal Versus D-Optimal Split-Plot Response Surface Designs

I-Optimal Versus D-Optimal Split-Plot Response Surface Designs I-Optimal Versus D-Optimal Split-Plot Response Surface Designs BRADLEY JONES SAS Institute and Universiteit Antwerpen PETER GOOS Universiteit Antwerpen and Erasmus Universiteit Rotterdam Response surface

More information

Optimal blocking of two-level fractional factorial designs

Optimal blocking of two-level fractional factorial designs Journal of Statistical Planning and Inference 91 (2000) 107 121 www.elsevier.com/locate/jspi Optimal blocking of two-level fractional factorial designs Runchu Zhang a;, DongKwon Park b a Department of

More information

Unit 6: Fractional Factorial Experiments at Three Levels

Unit 6: Fractional Factorial Experiments at Three Levels Unit 6: Fractional Factorial Experiments at Three Levels Larger-the-better and smaller-the-better problems. Basic concepts for 3 k full factorial designs. Analysis of 3 k designs using orthogonal components

More information

Stat 579: Generalized Linear Models and Extensions

Stat 579: Generalized Linear Models and Extensions Stat 579: Generalized Linear Models and Extensions Mixed models Yan Lu Feb, 2018, week 7 1 / 17 Some commonly used experimental designs related to mixed models Two way or three way random/mixed effects

More information

Obtaining Non-isomorphic Two-level Regular Fractional Factorial Designs

Obtaining Non-isomorphic Two-level Regular Fractional Factorial Designs Obtaining Non-isomorphic Two-level Regular Fractional Factorial Designs by Chunfang Lin B.Sc., 2002. University of Science and Technology of China. a project submitted in partial fulfillment of the requirements

More information

Bounds on the maximum numbers of clear two-factor interactions for 2 (n 1+n 2 ) (k 1 +k 2 ) fractional factorial split-plot designs

Bounds on the maximum numbers of clear two-factor interactions for 2 (n 1+n 2 ) (k 1 +k 2 ) fractional factorial split-plot designs 1816 Science in China: Series A Mathematics 2006 Vol. 49 No. 12 1816 1829 DOI: 10.1007/s11425-006-2032-2 Bounds on the maximum numbers of clear two-factor interactions for 2 (n 1+n 2 ) (k 1 +k 2 ) fractional

More information

OPTIMAL FRACTIONAL FACTORIAL SPLIT-PLOT DESIGNS FOR MODEL SELECTION

OPTIMAL FRACTIONAL FACTORIAL SPLIT-PLOT DESIGNS FOR MODEL SELECTION OPTIMAL FRACTIONAL FACTORIAL SPLIT-PLOT DESIGNS FOR MODEL SELECTION by Joslin Goh B.Math., University of Waterloo, 2007 a Project submitted in partial fulfillment of the requirements for the degree of

More information

SOME NEW THREE-LEVEL ORTHOGONAL MAIN EFFECTS PLANS ROBUST TO MODEL UNCERTAINTY

SOME NEW THREE-LEVEL ORTHOGONAL MAIN EFFECTS PLANS ROBUST TO MODEL UNCERTAINTY Statistica Sinica 14(2004), 1075-1084 SOME NEW THREE-LEVEL ORTHOGONAL MAIN EFFECTS PLANS ROBUST TO MODEL UNCERTAINTY Pi-Wen Tsai, Steven G. Gilmour and Roger Mead National Health Research Institutes, Queen

More information

A note on optimal foldover design

A note on optimal foldover design Statistics & Probability Letters 62 (2003) 245 250 A note on optimal foldover design Kai-Tai Fang a;, Dennis K.J. Lin b, HongQin c;a a Department of Mathematics, Hong Kong Baptist University, Kowloon Tong,

More information

ETH p. 1/21. Split Plot Designs. Large and small units Confounding main effects Repeated measures anova

ETH p. 1/21. Split Plot Designs. Large and small units Confounding main effects Repeated measures anova ETH p 1/21 Split Plot Designs Large and small units Confounding main effects Repeated measures anova ETH p 2/21 Study in Dental Medicine Can measurement of electric resistance help in detecting tooth decay?

More information

Some Nonregular Designs From the Nordstrom and Robinson Code and Their Statistical Properties

Some Nonregular Designs From the Nordstrom and Robinson Code and Their Statistical Properties Some Nonregular Designs From the Nordstrom and Robinson Code and Their Statistical Properties HONGQUAN XU Department of Statistics, University of California, Los Angeles, CA 90095-1554, U.S.A. (hqxu@stat.ucla.edu)

More information

Lecture 12: Diagonalization

Lecture 12: Diagonalization Lecture : Diagonalization A square matrix D is called diagonal if all but diagonal entries are zero: a a D a n 5 n n. () Diagonal matrices are the simplest matrices that are basically equivalent to vectors

More information

arxiv: v1 [stat.me] 16 Dec 2008

arxiv: v1 [stat.me] 16 Dec 2008 Recent Developments in Nonregular Fractional Factorial Designs Hongquan Xu, Frederick K. H. Phoa and Weng Kee Wong University of California, Los Angeles May 30, 2018 arxiv:0812.3000v1 [stat.me] 16 Dec

More information

Two-Way Factorial Designs

Two-Way Factorial Designs 81-86 Two-Way Factorial Designs Yibi Huang 81-86 Two-Way Factorial Designs Chapter 8A - 1 Problem 81 Sprouting Barley (p166 in Oehlert) Brewer s malt is produced from germinating barley, so brewers like

More information

Optimal Selection of Blocked Two-Level. Fractional Factorial Designs

Optimal Selection of Blocked Two-Level. Fractional Factorial Designs Applied Mathematical Sciences, Vol. 1, 2007, no. 22, 1069-1082 Optimal Selection of Blocked Two-Level Fractional Factorial Designs Weiming Ke Department of Mathematics and Statistics South Dakota State

More information

460 HOLGER DETTE AND WILLIAM J STUDDEN order to examine how a given design behaves in the model g` with respect to the D-optimality criterion one uses

460 HOLGER DETTE AND WILLIAM J STUDDEN order to examine how a given design behaves in the model g` with respect to the D-optimality criterion one uses Statistica Sinica 5(1995), 459-473 OPTIMAL DESIGNS FOR POLYNOMIAL REGRESSION WHEN THE DEGREE IS NOT KNOWN Holger Dette and William J Studden Technische Universitat Dresden and Purdue University Abstract:

More information

Stat/F&W Ecol/Hort 572 Review Points Ané, Spring 2010

Stat/F&W Ecol/Hort 572 Review Points Ané, Spring 2010 1 Linear models Y = Xβ + ɛ with ɛ N (0, σ 2 e) or Y N (Xβ, σ 2 e) where the model matrix X contains the information on predictors and β includes all coefficients (intercept, slope(s) etc.). 1. Number of

More information

Choosing initial values for the EM algorithm for nite mixtures

Choosing initial values for the EM algorithm for nite mixtures Computational Statistics & Data Analysis 41 (2003) 577 590 www.elsevier.com/locate/csda Choosing initial values for the EM algorithm for nite mixtures Dimitris Karlis, Evdokia Xekalaki Department of Statistics,

More information

Lec 10: Fractions of 2 k Factorial Design

Lec 10: Fractions of 2 k Factorial Design December 5, 2011 Fraction of 2 k experiments Screening: Some of the factors may influence the results. We want to figure out which. The number of combinations, 2 k, is too large for a complete investigation.

More information

USING REGULAR FRACTIONS OF TWO-LEVEL DESIGNS TO FIND BASELINE DESIGNS

USING REGULAR FRACTIONS OF TWO-LEVEL DESIGNS TO FIND BASELINE DESIGNS Statistica Sinica 26 (2016, 745-759 doi:http://dx.doi.org/10.5705/ss.202014.0099 USING REGULAR FRACTIONS OF TWO-LEVEL DESIGNS TO FIND BASELINE DESIGNS Arden Miller and Boxin Tang University of Auckland

More information

Orthogonal and nearly orthogonal designs for computer experiments

Orthogonal and nearly orthogonal designs for computer experiments Biometrika (2009), 96, 1,pp. 51 65 C 2009 Biometrika Trust Printed in Great Britain doi: 10.1093/biomet/asn057 Advance Access publication 9 January 2009 Orthogonal and nearly orthogonal designs for computer

More information

To reduce the cost of industrial experimentation one can forgo complete randomization. An example of a design with restricted

To reduce the cost of industrial experimentation one can forgo complete randomization. An example of a design with restricted Research Article (www.interscience.wiley.com) DOI: 10.1002/qre.1047 Published online 8 July 2009 in Wiley InterScience Design and Analysis of Industrial Strip-plot Experiments Heidi Arnouts, a Peter Goos

More information

CBSE Sample Paper-03 (Unsolved) SUMMATIVE ASSESSMENT II MATHEMATICS Class IX. Time allowed: 3 hours Maximum Marks: 90

CBSE Sample Paper-03 (Unsolved) SUMMATIVE ASSESSMENT II MATHEMATICS Class IX. Time allowed: 3 hours Maximum Marks: 90 CBSE Sample Paper-3 (Unsolved) SUMMATIVE ASSESSMENT II MATHEMATICS Class IX Time allowed: 3 hours Maximum Marks: 9 General Instructions: a) All questions are compulsory. b) The question paper consists

More information

Eton College King's Scholarship Examination 2013

Eton College King's Scholarship Examination 2013 Eton College King's Scholarship Examination 2013 MATHEMATICS A Answer Question 1 and as many of the other five questions as you can. Question 1 is worth 50 marks. All other questions are worth 10 marks

More information

From Handout #1, the randomization model for a design with a simple block structure can be written as

From Handout #1, the randomization model for a design with a simple block structure can be written as Hout 4 Strata Null ANOVA From Hout 1 the romization model for a design with a simple block structure can be written as C œ.1 \ X α % (4.1) w where α œðα á α > Ñ E ÐÑœ % 0 Z œcov ÐÑ % with cov( % 3 % 4

More information

arxiv: v2 [math.ho] 10 May 2018

arxiv: v2 [math.ho] 10 May 2018 Constructive Harmonic Conjugates Mark Mandelkern arxiv:1805.00978v2 [math.ho] 10 May 2018 Abstract In the study of the real projective plane, harmonic conjugates have an essential role, with applications

More information

Florida State University Libraries

Florida State University Libraries Florida State University Libraries Electronic Theses, Treatises and Dissertations The Graduate School 2006 Efficient Mixed-Level Fractional Factorial Designs: Evaluation, Augmentation and Application Yong

More information

Abstract. Using graphs to find optimal block designs. Problem 1: Factorial design. Conference on Design of Experiments, Tianjin, June 2006

Abstract. Using graphs to find optimal block designs. Problem 1: Factorial design. Conference on Design of Experiments, Tianjin, June 2006 Abstract Ching-Shui Cheng was one of the pioneers of using graph theory to prove results about optimal incomplete-block designs. There are actually two graphs associated with an incomplete-block design,

More information

Question 1 ( 1.0 marks) places of decimals? Solution: Now, on dividing by 2, we obtain =

Question 1 ( 1.0 marks) places of decimals? Solution: Now, on dividing by 2, we obtain = Question 1 ( 1.0 marks) The decimal expansion of the rational number places of decimals? will terminate after how many The given expression i.e., can be rewritten as Now, on dividing 0.043 by 2, we obtain

More information

Moment Aberration Projection for Nonregular Fractional Factorial Designs

Moment Aberration Projection for Nonregular Fractional Factorial Designs Moment Aberration Projection for Nonregular Fractional Factorial Designs Hongquan Xu Department of Statistics University of California Los Angeles, CA 90095-1554 (hqxu@stat.ucla.edu) Lih-Yuan Deng Department

More information

EXISTENCE AND CONSTRUCTION OF RANDOMIZATION DEFINING CONTRAST SUBSPACES FOR REGULAR FACTORIAL DESIGNS

EXISTENCE AND CONSTRUCTION OF RANDOMIZATION DEFINING CONTRAST SUBSPACES FOR REGULAR FACTORIAL DESIGNS Submitted to the Annals of Statistics EXISTENCE AND CONSTRUCTION OF RANDOMIZATION DEFINING CONTRAST SUBSPACES FOR REGULAR FACTORIAL DESIGNS By Pritam Ranjan, Derek R. Bingham and Angela M. Dean, Acadia

More information

Factorial Designs and Harmonic Analysis on Finite Abelian Groups

Factorial Designs and Harmonic Analysis on Finite Abelian Groups Factorial Designs and Harmonic Analysis on Finite Abelian Groups P.M. van de Ven EURANDOM P.O. Box 513, 5600 MB, Eindhoven The Netherlands email: vandeven@eurandom.tue.nl A. Di Bucchianico EURANDOM and

More information

8. Quadrilaterals. If AC = 21 cm, BC = 29 cm and AB = 30 cm, find the perimeter of the quadrilateral ARPQ.

8. Quadrilaterals. If AC = 21 cm, BC = 29 cm and AB = 30 cm, find the perimeter of the quadrilateral ARPQ. 8. Quadrilaterals Q 1 Name a quadrilateral whose each pair of opposite sides is equal. Mark (1) Q 2 What is the sum of two consecutive angles in a parallelogram? Mark (1) Q 3 The angles of quadrilateral

More information

These outputs can be written in a more convenient form: with y(i) = Hc m (i) n(i) y(i) = (y(i); ; y K (i)) T ; c m (i) = (c m (i); ; c m K(i)) T and n

These outputs can be written in a more convenient form: with y(i) = Hc m (i) n(i) y(i) = (y(i); ; y K (i)) T ; c m (i) = (c m (i); ; c m K(i)) T and n Binary Codes for synchronous DS-CDMA Stefan Bruck, Ulrich Sorger Institute for Network- and Signal Theory Darmstadt University of Technology Merckstr. 25, 6428 Darmstadt, Germany Tel.: 49 65 629, Fax:

More information

STAT22200 Spring 2014 Chapter 8A

STAT22200 Spring 2014 Chapter 8A STAT22200 Spring 2014 Chapter 8A Yibi Huang May 13, 2014 81-86 Two-Way Factorial Designs Chapter 8A - 1 Problem 81 Sprouting Barley (p166 in Oehlert) Brewer s malt is produced from germinating barley,

More information

Ð"Ñ + Ð"Ñ, Ð"Ñ +, +, + +, +,,

ÐÑ + ÐÑ, ÐÑ +, +, + +, +,, Handout #11 Confounding: Complete factorial experiments in incomplete blocks Blocking is one of the important principles in experimental design. In this handout we address the issue of designing complete

More information

Mathematics Revision Guides Vectors Page 1 of 19 Author: Mark Kudlowski M.K. HOME TUITION. Mathematics Revision Guides Level: GCSE Higher Tier VECTORS

Mathematics Revision Guides Vectors Page 1 of 19 Author: Mark Kudlowski M.K. HOME TUITION. Mathematics Revision Guides Level: GCSE Higher Tier VECTORS Mathematics Revision Guides Vectors Page of 9 M.K. HOME TUITION Mathematics Revision Guides Level: GCSE Higher Tier VECTORS Version:.4 Date: 05-0-05 Mathematics Revision Guides Vectors Page of 9 VECTORS

More information

Construction and analysis of Es 2 efficient supersaturated designs

Construction and analysis of Es 2 efficient supersaturated designs Construction and analysis of Es 2 efficient supersaturated designs Yufeng Liu a Shiling Ruan b Angela M. Dean b, a Department of Statistics and Operations Research, Carolina Center for Genome Sciences,

More information

Automatic generation of generalised regular factorial designs. l Environnement, Jouy en Josas, France

Automatic generation of generalised regular factorial designs. l Environnement, Jouy en Josas, France Automatic generation of generalised regular factorial designs André Kobilinsky (1), Hervé Monod (1) and Rosemary A. Bailey (2) (1) INRA, UR MaIAGE Mathématiques et Informatique Appliqués du Génome à l

More information

Topic 2 [312 marks] The rectangle ABCD is inscribed in a circle. Sides [AD] and [AB] have lengths

Topic 2 [312 marks] The rectangle ABCD is inscribed in a circle. Sides [AD] and [AB] have lengths Topic 2 [312 marks] 1 The rectangle ABCD is inscribed in a circle Sides [AD] and [AB] have lengths [12 marks] 3 cm and (\9\) cm respectively E is a point on side [AB] such that AE is 3 cm Side [DE] is

More information

Use of α-resolvable designs in the construction of two-factor experiments of split-plot type

Use of α-resolvable designs in the construction of two-factor experiments of split-plot type Biometrical Letters Vol. 53 206, No. 2, 05-8 DOI: 0.55/bile-206-0008 Use of α-resolvable designs in the construction of two-factor experiments of split-plot type Kazuhiro Ozawa, Shinji Kuriki 2, Stanisław

More information

DATE: MATH ANALYSIS 2 CHAPTER 12: VECTORS & DETERMINANTS

DATE: MATH ANALYSIS 2 CHAPTER 12: VECTORS & DETERMINANTS NAME: PERIOD: DATE: MATH ANALYSIS 2 MR. MELLINA CHAPTER 12: VECTORS & DETERMINANTS Sections: v 12.1 Geometric Representation of Vectors v 12.2 Algebraic Representation of Vectors v 12.3 Vector and Parametric

More information

Figure 9.1: A Latin square of order 4, used to construct four types of design

Figure 9.1: A Latin square of order 4, used to construct four types of design 152 Chapter 9 More about Latin Squares 9.1 Uses of Latin squares Let S be an n n Latin square. Figure 9.1 shows a possible square S when n = 4, using the symbols 1, 2, 3, 4 for the letters. Such a Latin

More information

Assignment 9 Answer Keys

Assignment 9 Answer Keys Assignment 9 Answer Keys Problem 1 (a) First, the respective means for the 8 level combinations are listed in the following table A B C Mean 26.00 + 34.67 + 39.67 + + 49.33 + 42.33 + + 37.67 + + 54.67

More information

Answer Keys to Homework#10

Answer Keys to Homework#10 Answer Keys to Homework#10 Problem 1 Use either restricted or unrestricted mixed models. Problem 2 (a) First, the respective means for the 8 level combinations are listed in the following table A B C Mean

More information

This model of the conditional expectation is linear in the parameters. A more practical and relaxed attitude towards linear regression is to say that

This model of the conditional expectation is linear in the parameters. A more practical and relaxed attitude towards linear regression is to say that Linear Regression For (X, Y ) a pair of random variables with values in R p R we assume that E(Y X) = β 0 + with β R p+1. p X j β j = (1, X T )β j=1 This model of the conditional expectation is linear

More information

Lecture 11: Eigenvalues and Eigenvectors

Lecture 11: Eigenvalues and Eigenvectors Lecture : Eigenvalues and Eigenvectors De nition.. Let A be a square matrix (or linear transformation). A number λ is called an eigenvalue of A if there exists a non-zero vector u such that A u λ u. ()

More information

Factorial Treatment Structure: Part I. Lukas Meier, Seminar für Statistik

Factorial Treatment Structure: Part I. Lukas Meier, Seminar für Statistik Factorial Treatment Structure: Part I Lukas Meier, Seminar für Statistik Factorial Treatment Structure So far (in CRD), the treatments had no structure. So called factorial treatment structure exists if

More information

Statistica Sinica Preprint No: SS

Statistica Sinica Preprint No: SS Preprint No: SS-2018-0120 Title Order-of-Addition Modeling Manuscript ID SS-2018-0120 URL http://www.stat.sinica.edu.tw/statistica/ DOI 10.5705/ss.202018.0120 Complete List of Authors Robert W Mee Corresponding

More information

Theorems. Least squares regression

Theorems. Least squares regression Theorems In this assignment we are trying to classify AML and ALL samples by use of penalized logistic regression. Before we indulge on the adventure of classification we should first explain the most

More information

Overcoming pseudoreplication in experimental designs

Overcoming pseudoreplication in experimental designs Overcoming pseudoreplication in experimental designs Ky Mathews, Nicole Cocks, David Hughes, Brian Cullis National Institute of Applied Statistics Research Australia University of Wollongong Australasian

More information

A Short Overview of Orthogonal Arrays

A Short Overview of Orthogonal Arrays A Short Overview of Orthogonal Arrays John Stufken Department of Statistics University of Georgia Isaac Newton Institute September 5, 2011 John Stufken (University of Georgia) Orthogonal Arrays September

More information

Generalized Shifted Inverse Iterations on Grassmann Manifolds 1

Generalized Shifted Inverse Iterations on Grassmann Manifolds 1 Proceedings of the Sixteenth International Symposium on Mathematical Networks and Systems (MTNS 2004), Leuven, Belgium Generalized Shifted Inverse Iterations on Grassmann Manifolds 1 J. Jordan α, P.-A.

More information

Split-plot Designs. Bruce A Craig. Department of Statistics Purdue University. STAT 514 Topic 21 1

Split-plot Designs. Bruce A Craig. Department of Statistics Purdue University. STAT 514 Topic 21 1 Split-plot Designs Bruce A Craig Department of Statistics Purdue University STAT 514 Topic 21 1 Randomization Defines the Design Want to study the effect of oven temp (3 levels) and amount of baking soda

More information

Constrained Optimization

Constrained Optimization Constrained Optimization Joshua Wilde, revised by Isabel Tecu, Takeshi Suzuki and María José Boccardi August 13, 2013 1 General Problem Consider the following general constrained optimization problem:

More information

Optimal Common-Centroid-Based Unit Capacitor Placements for Yield Enhancement of Switched-Capacitor Circuits

Optimal Common-Centroid-Based Unit Capacitor Placements for Yield Enhancement of Switched-Capacitor Circuits Optimal Common-Centroid-Based Unit Capacitor Placements for Yield Enhancement of Switched-Capacitor Circuits CHIEN-CHIH HUANG, National Central University CHIN-LONG WEY, National Chiao Tung University

More information

Analysis of variance using orthogonal projections

Analysis of variance using orthogonal projections Analysis of variance using orthogonal projections Rasmus Waagepetersen Abstract The purpose of this note is to show how statistical theory for inference in balanced ANOVA models can be conveniently developed

More information

Definitive Screening Designs

Definitive Screening Designs Definitive Screening Designs Bradley Jones September 2011 Copyright 2008, SAS Institute Inc. All rights reserved. Joint work with Chris Nachtsheim Outline 1. Motivation 2. Design Structure 3. Design Construction

More information

Re-sampling and exchangeable arrays University Ave. November Revised January Summary

Re-sampling and exchangeable arrays University Ave. November Revised January Summary Re-sampling and exchangeable arrays Peter McCullagh Department of Statistics University of Chicago 5734 University Ave Chicago Il 60637 November 1997 Revised January 1999 Summary The non-parametric, or

More information

One Variable Calculus. Izmir University of Economics Econ 533: Quantitative Methods and Econometrics

One Variable Calculus. Izmir University of Economics Econ 533: Quantitative Methods and Econometrics Izmir University of Economics Econ 533: Quantitative Methods and Econometrics One Variable Calculus Introduction Finding the best way to do a specic task involves what is called an optimization problem.

More information

Small Units inside Large Units

Small Units inside Large Units Chapter 8 Small Units inside Large Units 8.1 Experimental units bigger than observational units 8.1.1 The context Example 8.1 (Example 1.2 continued: Calf-feeding) Four feed treatments are compared on

More information

Algorithmic Construction of Efficient Fractional Factorial Designs With Large Run Sizes

Algorithmic Construction of Efficient Fractional Factorial Designs With Large Run Sizes Algorithmic Construction of Efficient Fractional Factorial Designs With Large Run Sizes Hongquan Xu Department of Statistics University of California Los Angeles, CA 90095-1554 (hqxu@stat.ucla.edu) February

More information

Solutions to Exercises

Solutions to Exercises 1 c Atkinson et al 2007, Optimum Experimental Designs, with SAS Solutions to Exercises 1. and 2. Certainly, the solutions to these questions will be different for every reader. Examples of the techniques

More information

Extended Bayesian Information Criteria for Model Selection with Large Model Spaces

Extended Bayesian Information Criteria for Model Selection with Large Model Spaces Extended Bayesian Information Criteria for Model Selection with Large Model Spaces Jiahua Chen, University of British Columbia Zehua Chen, National University of Singapore (Biometrika, 2008) 1 / 18 Variable

More information

5.) For each of the given sets of vectors, determine whether or not the set spans R 3. Give reasons for your answers.

5.) For each of the given sets of vectors, determine whether or not the set spans R 3. Give reasons for your answers. Linear Algebra - Test File - Spring Test # For problems - consider the following system of equations. x + y - z = x + y + 4z = x + y + 6z =.) Solve the system without using your calculator..) Find the

More information

The restarted QR-algorithm for eigenvalue computation of structured matrices

The restarted QR-algorithm for eigenvalue computation of structured matrices Journal of Computational and Applied Mathematics 149 (2002) 415 422 www.elsevier.com/locate/cam The restarted QR-algorithm for eigenvalue computation of structured matrices Daniela Calvetti a; 1, Sun-Mi

More information

Interaction balance in symmetrical factorial designs with generalized minimum aberration

Interaction balance in symmetrical factorial designs with generalized minimum aberration Interaction balance in symmetrical factorial designs with generalized minimum aberration Mingyao Ai and Shuyuan He LMAM, School of Mathematical Sciences, Peing University, Beijing 100871, P. R. China Abstract:

More information

Eigenvalue problems and optimization

Eigenvalue problems and optimization Notes for 2016-04-27 Seeking structure For the past three weeks, we have discussed rather general-purpose optimization methods for nonlinear equation solving and optimization. In practice, of course, we

More information

On the Application of Affine Resolvable Designs to Variety. Trials

On the Application of Affine Resolvable Designs to Variety. Trials On the Application of Affine Resolvable Designs to Variety. Trials Tadeusz Caliński, Stanisław Czajka, and Wiesław Pilarczyk Depatment of Mathematical and Statistical Methods Poznań University of Life

More information

A connection between number theory and linear algebra

A connection between number theory and linear algebra A connection between number theory and linear algebra Mark Steinberger Contents 1. Some basics 1 2. Rational canonical form 2 3. Prime factorization in F[x] 4 4. Units and order 5 5. Finite fields 7 6.

More information

Similarity of Triangle

Similarity of Triangle Similarity of Triangle 95 17 Similarity of Triangle 17.1 INTRODUCTION Looking around you will see many objects which are of the same shape but of same or different sizes. For examples, leaves of a tree

More information

BASIC ALGORITHMS IN LINEAR ALGEBRA. Matrices and Applications of Gaussian Elimination. A 2 x. A T m x. A 1 x A T 1. A m x

BASIC ALGORITHMS IN LINEAR ALGEBRA. Matrices and Applications of Gaussian Elimination. A 2 x. A T m x. A 1 x A T 1. A m x BASIC ALGORITHMS IN LINEAR ALGEBRA STEVEN DALE CUTKOSKY Matrices and Applications of Gaussian Elimination Systems of Equations Suppose that A is an n n matrix with coefficents in a field F, and x = (x,,

More information

Problem Description The problem we consider is stabilization of a single-input multiple-state system with simultaneous magnitude and rate saturations,

Problem Description The problem we consider is stabilization of a single-input multiple-state system with simultaneous magnitude and rate saturations, SEMI-GLOBAL RESULTS ON STABILIZATION OF LINEAR SYSTEMS WITH INPUT RATE AND MAGNITUDE SATURATIONS Trygve Lauvdal and Thor I. Fossen y Norwegian University of Science and Technology, N-7 Trondheim, NORWAY.

More information

Plan of Class 4. Radial Basis Functions with moving centers. Projection Pursuit Regression and ridge. Principal Component Analysis: basic ideas

Plan of Class 4. Radial Basis Functions with moving centers. Projection Pursuit Regression and ridge. Principal Component Analysis: basic ideas Plan of Class 4 Radial Basis Functions with moving centers Multilayer Perceptrons Projection Pursuit Regression and ridge functions approximation Principal Component Analysis: basic ideas Radial Basis

More information

Multigrid Techniques for Nonlinear Eigenvalue Problems; Solutions of a. Sorin Costiner and Shlomo Ta'asan

Multigrid Techniques for Nonlinear Eigenvalue Problems; Solutions of a. Sorin Costiner and Shlomo Ta'asan Multigrid Techniques for Nonlinear Eigenvalue Problems; Solutions of a Nonlinear Schrodinger Eigenvalue Problem in 2D and 3D Sorin Costiner and Shlomo Ta'asan Department of Applied Mathematics and Computer

More information

Analysis of variance. Gilles Guillot. September 30, Gilles Guillot September 30, / 29

Analysis of variance. Gilles Guillot. September 30, Gilles Guillot September 30, / 29 Analysis of variance Gilles Guillot gigu@dtu.dk September 30, 2013 Gilles Guillot (gigu@dtu.dk) September 30, 2013 1 / 29 1 Introductory example 2 One-way ANOVA 3 Two-way ANOVA 4 Two-way ANOVA with interactions

More information

Expanding brackets and factorising

Expanding brackets and factorising CHAPTER 8 Epanding brackets and factorising 8 CHAPTER Epanding brackets and factorising 8.1 Epanding brackets There are three rows. Each row has n students. The number of students is 3 n 3n. Two students

More information

Undercomplete Independent Component. Analysis for Signal Separation and. Dimension Reduction. Category: Algorithms and Architectures.

Undercomplete Independent Component. Analysis for Signal Separation and. Dimension Reduction. Category: Algorithms and Architectures. Undercomplete Independent Component Analysis for Signal Separation and Dimension Reduction John Porrill and James V Stone Psychology Department, Sheeld University, Sheeld, S10 2UR, England. Tel: 0114 222

More information

Step 2: Select Analyze, Mixed Models, and Linear.

Step 2: Select Analyze, Mixed Models, and Linear. Example 1a. 20 employees were given a mood questionnaire on Monday, Wednesday and again on Friday. The data will be first be analyzed using a Covariance Pattern model. Step 1: Copy Example1.sav data file

More information

IE 5531: Engineering Optimization I

IE 5531: Engineering Optimization I IE 5531: Engineering Optimization I Lecture 15: Nonlinear optimization Prof. John Gunnar Carlsson November 1, 2010 Prof. John Gunnar Carlsson IE 5531: Engineering Optimization I November 1, 2010 1 / 24

More information

The legacy of Sir Ronald A. Fisher. Fisher s three fundamental principles: local control, replication, and randomization.

The legacy of Sir Ronald A. Fisher. Fisher s three fundamental principles: local control, replication, and randomization. 1 Chapter 1: Research Design Principles The legacy of Sir Ronald A. Fisher. Fisher s three fundamental principles: local control, replication, and randomization. 2 Chapter 2: Completely Randomized Design

More information

PRACTICE PROBLEMS FOR THE FINAL

PRACTICE PROBLEMS FOR THE FINAL PRACTICE PROBLEMS FOR THE FINAL Here are a slew of practice problems for the final culled from old exams:. Let P be the vector space of polynomials of degree at most. Let B = {, (t ), t + t }. (a) Show

More information

Introduction and mathematical preliminaries

Introduction and mathematical preliminaries Chapter Introduction and mathematical preliminaries Contents. Motivation..................................2 Finite-digit arithmetic.......................... 2.3 Errors in numerical calculations.....................

More information

Stat 502 Design and Analysis of Experiments Large Factorial Designs

Stat 502 Design and Analysis of Experiments Large Factorial Designs 1 Stat 502 Design and Analysis of Experiments Large Factorial Designs Fritz Scholz Department of Statistics, University of Washington December 3, 2013 2 Exploration of Many Factors The typical setting

More information