Algebraic Generation of Orthogonal Fractional Factorial Designs
|
|
- Annabel Robbins
- 5 years ago
- Views:
Transcription
1 Algebraic Generation of Orthogonal Fractional Factorial Designs Roberto Fontana Dipartimento di Matematica, Politecnico di Torino Special Session on Advances in Algebraic Statistics AMS 2010 Spring Southeastern Sectional Meeting Lexington, March 27th-28th,2010 Joint paper with Giovanni Pistone, Politecnico di Torino [Fontana and Pistone(2010)]
2 Outline Aim of the work Full factorial design, fractions and counting functions Strata Generation of fractions Concluding remarks
3 Aim of the work To find a general method for fractional factorial design generation, where general means without any restriction for the number of factors/level
4 Full factorial design - complex coding of levels We refer to [Pistone and Rogantin(2008)] m factors, D 1,..., D m D j has n j levels coded with the n j -th roots of the unity: D j = {ω 0,..., ω nj 1} ( 1 2π ω k = exp n j D is the full factorial design with complex coding D = D 1 D j D m. ζ = (ζ 1,..., ζ m ) is a point of D #D is the cardinality of D, #D = m i=1 n j ) ( k exp i 2π n j ) k ;
5 Fractions and Counting functions A fraction F of D is a multiset (F, f ) whose underlying set of elements F is contained in D and f is the multiplicity function f : F N that for each element in F gives the number of times it belongs to the multiset F. Definition The counting function R of a fraction F is a response defined on D so that for each ζ D, R(ζ) equals the number of appearances of ζ in the fraction. We denote by c α the coefficients of the representation of R on D using the monomial basis {X α, α L = Z n1 Z nj Z nm }: R(ζ) = α L c α X α (ζ) ζ D c α C.
6 Full factorial design - complex coding of levels L is the full factorial design with integer coding α = (α 1,..., α m ) L = Z n1 Z nj Z nm, α j = 0,..., n j 1, j = 1,..., m is an element of L X j is the j-th component function, which maps a point to its i-th component: X j : D (ζ 1,..., ζ m ) ζ j D j ; the function X j is called simple term or, by abuse of terminology, factor. X α is the interaction term X α 1 1 Xm αm, i.e. the function X α : D (ζ 1,..., ζ m ) ζ α 1 1 ζαm m ;
7 Counting functions Proposition Let F be a fraction of a full factorial design D and R = α L c αx α be its counting function. 1 The coefficients c α are equal to 1 #D X α (ζ) ; ζ F 2 The term X α is centered on F, i.e. E F (X α ) = 1 if, and only if, c α = c [ α] = 0. #F ζ F X α (ζ) = 0, 3 The terms X α and X β are orthogonal on F, i.e. E F (X α X β ) = 0, if, and only if, c [α β] = 0.
8 Projectivity and Orthogonal Arrays Definition A fraction F factorially projects onto the I -factors, I {1,..., m}, if the projection is a multiple full factorial design, i.e. a full factorial design where each point appears equally often.
9 Projectivity and Orthogonal Arrays Definition A fraction F factorially projects onto the I -factors, I {1,..., m}, if the projection is a multiple full factorial design, i.e. a full factorial design where each point appears equally often. Definition A fraction F is a mixed orthogonal array of strength t if it factorially projects onto any I -factors with #I = t.
10 Projectivity and Orthogonal Arrays Definition A fraction F factorially projects onto the I -factors, I {1,..., m}, if the projection is a multiple full factorial design, i.e. a full factorial design where each point appears equally often. Definition A fraction F is a mixed orthogonal array of strength t if it factorially projects onto any I -factors with #I = t. Proposition (Projectivity)
11 Projectivity and Orthogonal Arrays Definition A fraction F factorially projects onto the I -factors, I {1,..., m}, if the projection is a multiple full factorial design, i.e. a full factorial design where each point appears equally often. Definition A fraction F is a mixed orthogonal array of strength t if it factorially projects onto any I -factors with #I = t. Proposition (Projectivity) 1 A fraction factorially projects onto the I -factors if, and only if, all the coefficients of the counting function involving only the I -factors are 0.
12 Projectivity and Orthogonal Arrays Definition A fraction F factorially projects onto the I -factors, I {1,..., m}, if the projection is a multiple full factorial design, i.e. a full factorial design where each point appears equally often. Definition A fraction F is a mixed orthogonal array of strength t if it factorially projects onto any I -factors with #I = t. Proposition (Projectivity) 1 A fraction factorially projects onto the I -factors if, and only if, all the coefficients of the counting function involving only the I -factors are 0. 2 A fraction is an orthogonal array of strength t if, and only if, all the coefficients of the counting function up to the order t are zero: c α = 0 for all α of order up to t, α (0, 0,..., 0).
13 Orthogonal array of strength t and indicator functions Let F = α L b αx α an indicator function of a fraction F of D Let C = {α L : 0 < α t}; F is an indicator function of an Orthogonal Array of strengt t if it is a solution of the following system { b α = β L b βb [α β] b α = 0, α C Remark F is an indicator function F (F 1) = 0 b α = β L b βb [α β]
14 Level set of X α Proposition (level set 1) Let X j the simple term with level set Ω nj = {ω 0,..., ω nj 1}. Let s consider the term Xj r and let s define { 1 r = 0 s j,r = n j / gcd(r, n j ) r > 0 Over D, the term X r j takes all the values of Ω sj,r equally often. Example n j = 4, D j = {ω 0, ω 1, ω 2, ω 3 } {1, i, 1, i} ζ j Xj 0 Xj 1 Xj 2 Xj ω 1 1 ω 1 ω 2 ω 3 ω 2 1 ω 2 1 ω 2 ω 3 1 ω 3 ω 2 ω 1
15 Level set of X α Proposition (level set 2) Let X α = X α 1 1 Xm αm an interaction. X α i i takes values in Ω si,αi where s i,αi is determined according to the previous Proposition level set 1. Let s define s α = lcm(s 1,α1,..., s m,αm ). Over D, the term X α takes all the values of Ω sα equally often. Example n j = n k = 4, X 2 j X 2 k. From Proposition level set 1 we have s j = s k = 2. We obtain s = 2. Indeed X 2 j (ζ j) {1, 1}, X 2 k (ζ k) {1, 1} X 2 j X 2 k (ζ j, ζ k ) {1, 1}
16 Strata Definition Given a term X α, α L = Z n1... Z nm the full design D is partitioned into the the following strata { } Dh α = ζ D : X α (ζ) = ω h where ω h Ω sα and s α is determined according to the previous Propositions level set 1 and level set 2. We use n α,h to denote the number of points of the fraction F that are in the stratum D α h, with h = 0,..., s α 1, n α,h = ζ D α h y ζ
17 Counting function as linear combination of indicator functions Let s consider the indicator functions 1 ζ of all the single points of D { 1 ζ = (ζ 1,..., ζ m ) 1 ζ : D (ζ 1,..., ζ m ) 0 ζ (ζ 1,..., ζ m ) The counting function R of a fraction F can be written as ζ D y ζ1 ζ with y ζ R(ζ) {0, 1,..., n,...}. The particular case in which R is an indicator function corresponds to y ζ {0, 1}. Proposition Let F be a fraction of D. Its counting fraction R can be expressed both as R = α c αx α and R = ζ D y ζ1 ζ. The relation between the coefficients c α and y ζ is c α = 1 y ζ X #D α (ζ) ζ D
18 c α and n α,h Proposition Let F be a fraction of D with counting fraction R = α L c αx α. Each c α, α L, depends on n α,h, h = 0,..., s α 1, as c α = 1 s α 1 n α,h ω h #D h=0 where s α is determined by X α (see Proposition level set 1 and Proposition level set 2 ). Viceversa, each n α,h, h = 0,..., s α 1, depends on c [ kα], k = 0,..., s 1 as n α,h = #D s α s α 1 k=0 c [ kα] ω [hk]
19 Conditions on n α,h to have X α centered Let F be a fraction of D with counting fraction R = 1 sα 1 #D α L h=0 n α,hω h X α. Proposition (*) Let X α be a term with level set Ω sα on full design D. s α is prime The term X α is centered on the fraction F if, and only if, its s α levels appear equally often: n α,0 = n α,1 =... = n α,s 1 = λ α s α is not prime Let P α (ζ) = s α 1 h=0 n α,hζ h and let s denote by Φ sα the cyclotomic polynomial of the s α -roots of the unity. The term X α is centered on the fraction F if, and only if, the remainder H α (ζ) = P α (ζ)modφ sα (ζ) whose coefficients are integer linear combinations of n α,h, h = 0,..., s α 1, is identically zero.
20 Conditions on R(ζ) to have X α centered n α,h are related to the values of the counting function R of a fraction F by the following relation n α,h = y ζ, ζ D α h Proposition * allows to express the condition X α is centered on F ( c α = 0) as integer linear combinations of the values y ζ of the counting function over the full design D.
21 Orthogonal Arrays For Orthogonal Arrays, using Proposition projectivity we have c α = 0 α C where C = {α L : 0 < α t} and α is the number of non null elements of α.
22 OA(n, s m, t), s prime, as a system of linear equations Let s consider OA(n, s m, t), s prime. We consider α C and, using Proposition *, we write the condition c α = 0 using strata as ζ D0 α y ζ = λ y ζ = λ Varying α C we get Considering λ as a variable, ζ D α 1... ζ Ds 1 α y ζ = λ where à = [A, 1] and Ỹ = (Y, λ). AY = λ1 ÃỸ = 0
23 Computation of generators for OA(n, s m, t), s prime The sum of two Orthogonal Arrays, Y 1 OA(n 1, s m, t) and Y 2 OA(n 2, s m, t), is an Orthogonal Array Y 1 + Y 2 OA(n 1 + n 2, s m, t). The Hilbert Basis [Schrijver(1986)] is a minimal set of generators such that any OA(n, s m, t) becomes a linear combination of the generators with positive or null integer coefficients. This approach extends that of [Carlini and Pistone(2007)] where homogeneous 2-level fractions were considered and the conditions ζ F X α (ζ) = 0 were used. c α = 1 #D The advantage of using strata is that we avoid computations with complex numbers (X α (ζ)). We give some examples. For the computation we use 4ti2 [4ti2 team(2008)].
24 Some examples of OA(n, s m, t), s prime Some classes of OA s: A(n, 2 5, 2) The matrix à has 30 rows and 33 columns. We find 26, 142 solutions A(n, 3 3, 2) The matrix à has has 54 rows and 28 columns. We find 66 solutions: 12 have 9 points, all different 54 have 18 points, 17 different. A(n, 3 4, 3) The matrix à has has 192 rows and 82 columns. We find 131, 892 solutions
25 One example of OA(18, 3 3, 2) Classification of OA(n, 3 4, 3) # support # of points max rep N , , , , , , , , , , , , , , ,184 Total 131,892
26 OA(n, s n sn k, t) as a system of linear equations k Let s consider mixed level orthogonal arrays. OA(n, s m, t), without the condition s prime, become a particular case. Using Proposition *, we express the condition c α = 0, α C using strata. Let s consider, as an example, OA(n, 2 3 3, 3). Let s take α (1, 1, 0, 0) C = {α Z 2 Z 3 Z 3 Z 3 : 1 α 3}. We get s α = 6. The condition c 1,1,0,0 = 0 is equivalent to X 1,1,0,0 is centered that is equivalent to H α (ζ) P α (ζ)modφ 6 (ζ) is identically zero. We have and therefore P α (ζ) = n (1,1,0,0),0 + n (1,1,0,0),1 ζ n (1,1,0,0),5 ζ 5 Φ 6 (ζ) = 1 ζ + ζ 2 H α (ζ) P α (ζ)modφ 6 (ζ) = n (1,1,0,0),0 n (1,1,0,0),2 n (1,1,0,0),3 + n (1,1,0,0),6 + (n (1,1,0,0),1 + n (1,1,0,0),2 n (1,1,0,0),5 n (1,1,0,0),6 )ζ
27 OA(n, s n sn k, t) as a system of linear equations k Varying α C and reminding that n α,h = ζ Dh α y ζ we get ÃỸ = 0. Some classes of OA s: The matrix à has 89 rows and 54 columns. We find 403 solutions OA(n, 2 3 3, 3) OA(n, 4 2, 2) OA(n, 6 2, 2) The matrix à has 10 rows and 16 columns. We find 24 solutions corresponding to all the LHD. The matrix à has 18 rows and 36 columns. We find 620 solutions corresponding to all the LHD. Classification of OA(n, 2 3 3, 3) # support # of points max rep N Total 403
28 Sampling - example find one single replicate OA(9, 3 3, 2);
29 Sampling - example find one single replicate OA(9, 3 3, 2); find one solution of AY = 1 with y ζ {0, 1}
30 Sampling - example find one single replicate OA(9, 3 3, 2); find one solution of AY = 1 with y ζ {0, 1} we used standard Simulated Annealing with objective function V (Y ) defined as the number of equations of AY = 1 that are satisfied; we implemented this algorithm using SAS/IML.
31 Sampling - example find one single replicate OA(9, 3 3, 2); find one solution of AY = 1 with y ζ {0, 1} we used standard Simulated Annealing with objective function V (Y ) defined as the number of equations of AY = 1 that are satisfied; we implemented this algorithm using SAS/IML. we found a solution in 2, 702 iterations (a couple of seconds on a standard laptop).
32 Sampling - example find one single replicate OA(9, 3 3, 2); find one solution of AY = 1 with y ζ {0, 1} we used standard Simulated Annealing with objective function V (Y ) defined as the number of equations of AY = 1 that are satisfied; we implemented this algorithm using SAS/IML. we found a solution in 2, 702 iterations (a couple of seconds on a standard laptop). we have also experimented the algorithm on
33 Sampling - example find one single replicate OA(9, 3 3, 2); find one solution of AY = 1 with y ζ {0, 1} we used standard Simulated Annealing with objective function V (Y ) defined as the number of equations of AY = 1 that are satisfied; we implemented this algorithm using SAS/IML. we found a solution in 2, 702 iterations (a couple of seconds on a standard laptop). we have also experimented the algorithm on OA(9, 3 4 ; 3); we found one solution in 2, 895 iterations;
34 Sampling - example find one single replicate OA(9, 3 3, 2); find one solution of AY = 1 with y ζ {0, 1} we used standard Simulated Annealing with objective function V (Y ) defined as the number of equations of AY = 1 that are satisfied; we implemented this algorithm using SAS/IML. we found a solution in 2, 702 iterations (a couple of seconds on a standard laptop). we have also experimented the algorithm on OA(9, 3 4 ; 3); we found one solution in 2, 895 iterations; 4 4 sudoku; we found one solution in 2, 852 iterations;
35 Sampling - example find one single replicate OA(9, 3 3, 2); find one solution of AY = 1 with y ζ {0, 1} we used standard Simulated Annealing with objective function V (Y ) defined as the number of equations of AY = 1 that are satisfied; we implemented this algorithm using SAS/IML. we found a solution in 2, 702 iterations (a couple of seconds on a standard laptop). we have also experimented the algorithm on OA(9, 3 4 ; 3); we found one solution in 2, 895 iterations; 4 4 sudoku; we found one solution in 2, 852 iterations; 9 9 sudoku; we did not find any solution in 100, 000 iterations;
36 Sampling example move from one single replicate OA(9, 3 3, 2) to another single replicate OA(9, 3 3, 2) ;
37 Sampling example move from one single replicate OA(9, 3 3, 2) to another single replicate OA(9, 3 3, 2) ; move from one solution of AY = 1 with y ζ {0, 1} to another solution of AY = 1 with y ζ {0, 1};
38 Sampling example move from one single replicate OA(9, 3 3, 2) to another single replicate OA(9, 3 3, 2) ; move from one solution of AY = 1 with y ζ {0, 1} to another solution of AY = 1 with y ζ {0, 1}; we used 4ti2 to generate M, the Markov Basis corresponding to the homogeneous system AX = 0.
39 Sampling example move from one single replicate OA(9, 3 3, 2) to another single replicate OA(9, 3 3, 2) ; move from one solution of AY = 1 with y ζ {0, 1} to another solution of AY = 1 with y ζ {0, 1}; we used 4ti2 to generate M, the Markov Basis corresponding to the homogeneous system AX = 0. M contains 81 different moves;
40 Sampling example move from one single replicate OA(9, 3 3, 2) to another single replicate OA(9, 3 3, 2) ; move from one solution of AY = 1 with y ζ {0, 1} to another solution of AY = 1 with y ζ {0, 1}; we used 4ti2 to generate M, the Markov Basis corresponding to the homogeneous system AX = 0. M contains 81 different moves; as an initial fraction we considered the nine-run regular fraction F 0 whose indicator function R 0 is R 0 = 1 3 (1 + X 1X 2 X 3 + X1 2X 2 2X 3 2 ). We run simulations generating R i as R i = R i 1 + ɛ MRi 1 M Ri 1 (and moving to F i if R i 0);
41 Sampling example move from one single replicate OA(9, 3 3, 2) to another single replicate OA(9, 3 3, 2) ; move from one solution of AY = 1 with y ζ {0, 1} to another solution of AY = 1 with y ζ {0, 1}; we used 4ti2 to generate M, the Markov Basis corresponding to the homogeneous system AX = 0. M contains 81 different moves; as an initial fraction we considered the nine-run regular fraction F 0 whose indicator function R 0 is R 0 = 1 3 (1 + X 1X 2 X 3 + X1 2X 2 2X 3 2 ). We run simulations generating R i as R i = R i 1 + ɛ MRi 1 M Ri 1 (and moving to F i if R i 0); we obtained all the 12 different 9-run fractions, each one with 9 different points as known in the literature and as shown previously.
42 Concluding remarks + no restrictions on the number of levels
43 Concluding remarks + no restrictions on the number of levels + no computation with complex numbers
44 Concluding remarks + no restrictions on the number of levels + no computation with complex numbers + Hilbert basis, Markov basis, MCMC available
45 Concluding remarks + no restrictions on the number of levels + no computation with complex numbers + Hilbert basis, Markov basis, MCMC available computational effort
46 Some references 4ti2 team, ti2 a software package for algebraic, geometric and combinatorial problems on linear spaces. Available at Carlini, E., Pistone, G., Hibert bases for orthogonal arrays. Journal of Statistical Theory and practice, Fontana, R., Pistone, G., Algebraic strata for non symmetrical orthogonal fractional factorial designs and application. La matematica e le sue applicazioni 2010/1, Politecnico di Torino DIMAT. Pistone, G., Rogantin, M., Indicator function and complex coding for mixed fractional factorial designs. J. Statist. Plann. Inference 138 (3), Schrijver, A., Theory of linear and integer programming. Wiley-Interscience Series in Discrete Mathematics. John Wiley & Sons Ltd., Chichester, a Wiley-Interscience Publication.
47 Thanks for your attention!
arxiv:math/ v1 [math.st] 13 Mar 2007
Indicator function and complex coding for mixed fractional factorial designs 1 Giovanni Pistone arxiv:math/0703365v1 [math.st] 13 Mar 2007 Abstract Department of Mathematics - Politecnico di Torino Maria-Piera
More informationIndicator function and complex coding for mixed fractional factorial designs 1
Indicator function and complex coding for mixed fractional factorial designs 1 Giovanni Pistone Department of Mathematics - Politecnico di Torino Maria-Piera Rogantin Department of Mathematics - Università
More informationDipartimento di Matematica
Dipartimento di Matematica G. PISTONE, M. P. ROGANTIN INDICATOR FUNCTION AND COMPLEX CODING FOR MIXED FRACTIONAL FACTORIAL DESIGNS (REVISED 2006) Rapporto interno N. 17, luglio 2006 Politecnico di Torino
More informationCharacterizations of indicator functions of fractional factorial designs
Characterizations of indicator functions of fractional factorial designs arxiv:1810.08417v2 [math.st] 26 Oct 2018 Satoshi Aoki Abstract A polynomial indicator function of designs is first introduced by
More informationDipartimento di Matematica
Dipartimento di Matematica G. PISTONE, E. RICCOMAGNO, M. P. ROGANTIN ALGEBRAIC STATISTICS METHODS FOR DOE Rapporto interno N. 18, settembre 2006 Politecnico di Torino Corso Duca degli Abruzzi, 24-10129
More informationCharacterization of multivariate Bernoulli distributions with given margins
arxiv:176.1357v1 [math.st] 5 Jun 217 Characterization of multivariate Bernoulli distributions with given margins Roberto Fontana 1 and Patrizia Semeraro 2 1 Department of Mathematical Sciences, Politecnico
More informationToric statistical models: parametric and binomial representations
AISM (2007) 59:727 740 DOI 10.1007/s10463-006-0079-z Toric statistical models: parametric and binomial representations Fabio Rapallo Received: 21 February 2005 / Revised: 1 June 2006 / Published online:
More informationOn the decomposition of orthogonal arrays
On the decomposition of orthogonal arrays Wiebke S. Diestelkamp Department of Mathematics University of Dayton Dayton, OH 45469-2316 wiebke@udayton.edu Jay H. Beder Department of Mathematical Sciences
More informationOptimal Fractional Factorial Plans for Asymmetric Factorials
Optimal Fractional Factorial Plans for Asymmetric Factorials Aloke Dey Chung-yi Suen and Ashish Das April 15, 2002 isid/ms/2002/04 Indian Statistical Institute, Delhi Centre 7, SJSS Marg, New Delhi 110
More informationABSTRACT. Department of Mathematics. interesting results. A graph on n vertices is represented by a polynomial in n
ABSTRACT Title of Thesis: GRÖBNER BASES WITH APPLICATIONS IN GRAPH THEORY Degree candidate: Angela M. Hennessy Degree and year: Master of Arts, 2006 Thesis directed by: Professor Lawrence C. Washington
More informationORTHOGONAL ARRAYS OF STRENGTH 3 AND SMALL RUN SIZES
ORTHOGONAL ARRAYS OF STRENGTH 3 AND SMALL RUN SIZES ANDRIES E. BROUWER, ARJEH M. COHEN, MAN V.M. NGUYEN Abstract. All mixed (or asymmetric) orthogonal arrays of strength 3 with run size at most 64 are
More informationClassification of three-word indicator functions of two-level factorial designs
AISM (2006) 58: 595 608 DOI 0.007/s063-006-0033-0 N. Balakrishnan Po Yang Classification of three-word indicator functions of two-level factorial designs Received: 2 July 200 / Revised: 2 April 2005 /
More informationOn the construction of asymmetric orthogonal arrays
isid/ms/2015/03 March 05, 2015 http://wwwisidacin/ statmath/indexphp?module=preprint On the construction of asymmetric orthogonal arrays Tianfang Zhang and Aloke Dey Indian Statistical Institute, Delhi
More informationBricklaying and the Hermite Normal Form
Bricklaying and the Hermite Normal Form William J. Gilbert Pure Mathematics Department University of Waterloo Waterloo, Ontario Canada N2L 3G1 AMS Classifications: Primary: 15A36 Appeared in the American
More informationParameter inequalities for orthogonal arrays with mixed levels
Parameter inequalities for orthogonal arrays with mixed levels Wiebke S. Diestelkamp Department of Mathematics University of Dayton Dayton, OH 45469-2316 wiebke@udayton.edu Designs, Codes and Cryptography,
More informationFactorial 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 informationIndicator Functions and the Algebra of the Linear-Quadratic Parametrization
Biometrika (2012), 99, 1, pp. 1 12 C 2012 Biometrika Trust Printed in Great Britain Advance Access publication on 31 July 2012 Indicator Functions and the Algebra of the Linear-Quadratic Parametrization
More informationThe decomposability of simple orthogonal arrays on 3 symbols having t + 1 rows and strength t
The decomposability of simple orthogonal arrays on 3 symbols having t + 1 rows and strength t Wiebke S. Diestelkamp Department of Mathematics University of Dayton Dayton, OH 45469-2316 USA wiebke@udayton.edu
More informationFILTERED RINGS AND MODULES. GRADINGS AND COMPLETIONS.
FILTERED RINGS AND MODULES. GRADINGS AND COMPLETIONS. Let A be a ring, for simplicity assumed commutative. A filtering, or filtration, of an A module M means a descending sequence of submodules M = M 0
More informationDipartimento di Matematica
Dipartimento di Matematica S. RUFFA ON THE APPLICABILITY OF ALGEBRAIC STATISTICS TO REGRESSION MODELS WITH ERRORS IN VARIABLES Rapporto interno N. 8, maggio 006 Politecnico di Torino Corso Duca degli Abruzzi,
More informationInformation Geometry and Algebraic Statistics
Matematikos ir informatikos institutas Vinius Information Geometry and Algebraic Statistics Giovanni Pistone and Eva Riccomagno POLITO, UNIGE Wednesday 30th July, 2008 G.Pistone, E.Riccomagno (POLITO,
More informationMCS 563 Spring 2014 Analytic Symbolic Computation Monday 14 April. Binomial Ideals
Binomial Ideals Binomial ideals offer an interesting class of examples. Because they occur so frequently in various applications, the development methods for binomial ideals is justified. 1 Binomial Ideals
More informationSome 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 informationN. Makhijani, R. K. Sharma and J. B. Srivastava
International Electronic Journal of Algebra Volume 18 2015 21-33 UNITS IN F q kc p r C q N. Makhijani, R. K. Sharma and J. B. Srivastava Received: 5 February 2014 Communicated by A. Çiğdem Özcan Abstract.
More informationFigure 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 informationFRACTIONAL FACTORIAL DESIGNS OF STRENGTH 3 AND SMALL RUN SIZES
FRACTIONAL FACTORIAL DESIGNS OF STRENGTH 3 AND SMALL RUN SIZES ANDRIES E. BROUWER, ARJEH M. COHEN, MAN V.M. NGUYEN Abstract. All mixed (or asymmetric) orthogonal arrays of strength 3 with run size at most
More informationConnections between the resolutions of general two-level factorial designs
AISM (2006) 58: 609 68 DOI 0.007/s0463-005-0020-x N. Balakrishnan Po Yang Connections between the resolutions of general two-level factorial designs Received: 28 February 2005 / Published online: 3 June
More informationForms of four-word indicator functions with implications to two-level factorial designs
Ann Inst Stat Math (20) 63:375 386 DOI 0.007/s0463-009-0222-8 Forms of four-word indicator functions with implications to two-level factorial designs N. Balakrishnan Po Yang Received: 9 March 2008 / Revised:
More informationOptimal Ramp Schemes and Related Combinatorial Objects
Optimal Ramp Schemes and Related Combinatorial Objects Douglas R. Stinson David R. Cheriton School of Computer Science University of Waterloo BCC 2017, Glasgow, July 3 7, 2017 1 / 18 (t, n)-threshold Schemes
More informationApplications of Algebraic Geometry IMA Annual Program Year Workshop Applications in Biology, Dynamics, and Statistics March 5-9, 2007
Applications of Algebraic Geometry IMA Annual Program Year Workshop Applications in Biology, Dynamics, and Statistics March 5-9, 2007 Information Geometry (IG) and Algebraic Statistics (AS) Giovanni Pistone
More informationAn Additive Characterization of Fibers of Characters on F p
An Additive Characterization of Fibers of Characters on F p Chris Monico Texas Tech University Lubbock, TX c.monico@ttu.edu Michele Elia Politecnico di Torino Torino, Italy elia@polito.it January 30, 2009
More informationALGEBRA I FORM I. Textbook: Algebra, Second Edition;Prentice Hall,2002
ALGEBRA I FORM I Textbook: Algebra, Second Edition;Prentice Hall,00 Prerequisites: Students are expected to have a knowledge of Pre Algebra and proficiency of basic math skills including: positive and
More informationA GENERAL CONSTRUCTION FOR SPACE-FILLING LATIN HYPERCUBES
Statistica Sinica 6 (016), 675-690 doi:http://dx.doi.org/10.5705/ss.0015.0019 A GENERAL CONSTRUCTION FOR SPACE-FILLING LATIN HYPERCUBES C. Devon Lin and L. Kang Queen s University and Illinois Institute
More informationSamurai Sudoku-Based Space-Filling Designs
Samurai Sudoku-Based Space-Filling Designs Xu Xu and Peter Z. G. Qian Department of Statistics University of Wisconsin Madison, Madison, WI 53706 Abstract Samurai Sudoku is a popular variation of Sudoku.
More informationKnots-quivers correspondence, lattice paths, and rational knots
Knots-quivers correspondence, lattice paths, and rational knots 1 CAMGSD, Departamento de Matemática, Instituto Superior Técnico, Portugal 2 Mathematical Institute SANU, Belgrade, Serbia Warszawa, September
More informationMarkov chain Monte Carlo methods for the Box-Behnken designs and centrally symmetric configurations
arxiv:1502.02323v1 [math.st] 9 Feb 2015 Markov chain Monte Carlo methods for the Box-Behnken designs and centrally symmetric configurations Satoshi Aoki, Takayuki Hibi and Hidefumi Ohsugi October 27, 2018
More informationA 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 informationA Sampling Kaczmarz-Motzkin Algorithm for Linear Feasibility
A Sampling Kaczmarz-Motzkin Algorithm for Linear Feasibility Jamie Haddock Graduate Group in Applied Mathematics, Department of Mathematics, University of California, Davis Copper Mountain Conference on
More informationOn the intersection of Hermitian surfaces
On the intersection of Hermitian surfaces Nicola Durante and Gary Ebert Abstract In [6] and [3] the authors determine the structure of the intersection of two Hermitian surfaces of PG(3, q 2 ) under the
More informationBulletin of the Transilvania University of Braşov Vol 7(56), No Series III: Mathematics, Informatics, Physics, 59-64
Bulletin of the Transilvania University of Braşov Vol 756), No. 2-2014 Series III: Mathematics, Informatics, Physics, 59-64 ABOUT QUATERNION ALGEBRAS AND SYMBOL ALGEBRAS Cristina FLAUT 1 and Diana SAVIN
More informationAn infinite family of Goethals-Seidel arrays
An infinite family of Goethals-Seidel arrays Mingyuan Xia Department of Mathematics, Central China Normal University, Wuhan, Hubei 430079, China Email: xiamy@ccnu.edu.cn Tianbing Xia, Jennifer Seberry
More informationAnale. Seria Informatică. Vol. XIII fasc Annals. Computer Science Series. 13 th Tome 1 st Fasc. 2015
24 CONSTRUCTION OF ORTHOGONAL ARRAY-BASED LATIN HYPERCUBE DESIGNS FOR DETERMINISTIC COMPUTER EXPERIMENTS Kazeem A. Osuolale, Waheed B. Yahya, Babatunde L. Adeleke Department of Statistics, University of
More informationA Likelihood Ratio Test
A Likelihood Ratio Test David Allen University of Kentucky February 23, 2012 1 Introduction Earlier presentations gave a procedure for finding an estimate and its standard error of a single linear combination
More informationHILBERT BASIS OF THE LIPMAN SEMIGROUP
Available at: http://publications.ictp.it IC/2010/061 United Nations Educational, Scientific and Cultural Organization and International Atomic Energy Agency THE ABDUS SALAM INTERNATIONAL CENTRE FOR THEORETICAL
More informationQuasi-reducible Polynomials
Quasi-reducible Polynomials Jacques Willekens 06-Dec-2008 Abstract In this article, we investigate polynomials that are irreducible over Q, but are reducible modulo any prime number. 1 Introduction Let
More informationMTH Linear Algebra. Study Guide. Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education
MTH 3 Linear Algebra Study Guide Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education June 3, ii Contents Table of Contents iii Matrix Algebra. Real Life
More informationMath 101 Study Session Spring 2016 Test 4 Chapter 10, Chapter 11 Chapter 12 Section 1, and Chapter 12 Section 2
Math 101 Study Session Spring 2016 Test 4 Chapter 10, Chapter 11 Chapter 12 Section 1, and Chapter 12 Section 2 April 11, 2016 Chapter 10 Section 1: Addition and Subtraction of Polynomials A monomial is
More informationLinear Models Review
Linear Models Review Vectors in IR n will be written as ordered n-tuples which are understood to be column vectors, or n 1 matrices. A vector variable will be indicted with bold face, and the prime sign
More informationThe Dual Lattice, Integer Linear Systems and Hermite Normal Form
New York University, Fall 2013 Lattices, Convexity & Algorithms Lecture 2 The Dual Lattice, Integer Linear Systems and Hermite Normal Form Lecturers: D. Dadush, O. Regev Scribe: D. Dadush 1 Dual Lattice
More informationMutually Orthogonal Latin Squares: Covering and Packing Analogues
Squares: Covering 1 1 School of Computing, Informatics, and Decision Systems Engineering Arizona State University Mile High Conference, 15 August 2013 Latin Squares Definition A latin square of side n
More information6 Markov Chain Monte Carlo (MCMC)
6 Markov Chain Monte Carlo (MCMC) The underlying idea in MCMC is to replace the iid samples of basic MC methods, with dependent samples from an ergodic Markov chain, whose limiting (stationary) distribution
More informationNormed & Inner Product Vector Spaces
Normed & Inner Product Vector Spaces ECE 174 Introduction to Linear & Nonlinear Optimization Ken Kreutz-Delgado ECE Department, UC San Diego Ken Kreutz-Delgado (UC San Diego) ECE 174 Fall 2016 1 / 27 Normed
More informationThe primitive root theorem
The primitive root theorem Mar Steinberger First recall that if R is a ring, then a R is a unit if there exists b R with ab = ba = 1. The collection of all units in R is denoted R and forms a group under
More informationD-bounded Distance-Regular Graphs
D-bounded Distance-Regular Graphs CHIH-WEN WENG 53706 Abstract Let Γ = (X, R) denote a distance-regular graph with diameter D 3 and distance function δ. A (vertex) subgraph X is said to be weak-geodetically
More informationMath 2331 Linear Algebra
6. Orthogonal Projections Math 2 Linear Algebra 6. Orthogonal Projections Jiwen He Department of Mathematics, University of Houston jiwenhe@math.uh.edu math.uh.edu/ jiwenhe/math2 Jiwen He, University of
More informationMarkov bases and subbases for bounded contingency tables
Markov bases and subbases for bounded contingency tables arxiv:0905.4841v2 [math.co] 23 Jun 2009 Fabio Rapallo Abstract Ruriko Yoshida In this paper we study the computation of Markov bases for contingency
More informationOptimal Designs for Stated Choice Experiments Generated from Fractional Factorial Designs
Optimal Designs for Stated Choice Experiments Generated from Fractional Factorial Designs Stephen Bush School of Mathematical Sciences, University of Technology Sydney, PO Box 2 Broadway NSW 2007, Australia
More informationPPs over Finite Fields defined by functional equations
Permutation Polynomials over Finite Fields defined by functional equations Neranga Fernando Joint work with Dr. Xiang-dong Hou Department of Mathematics and Statistics University of South Florida AMS Spring
More informationExplicit Methods in Algebraic Number Theory
Explicit Methods in Algebraic Number Theory Amalia Pizarro Madariaga Instituto de Matemáticas Universidad de Valparaíso, Chile amaliapizarro@uvcl 1 Lecture 1 11 Number fields and ring of integers Algebraic
More information8. Statistical Equilibrium and Classification of States: Discrete Time Markov Chains
8. Statistical Equilibrium and Classification of States: Discrete Time Markov Chains 8.1 Review 8.2 Statistical Equilibrium 8.3 Two-State Markov Chain 8.4 Existence of P ( ) 8.5 Classification of States
More information5 Non-unique factorizations
5 Non-unique factorizations In this chapter we briefly discuss some aspects of non-unique factorization, ending with an application of quadratic forms (and more generally ideals) to factorization problems
More informationLacunary Polynomials over Finite Fields Course notes
Lacunary Polynomials over Finite Fields Course notes Javier Herranz Abstract This is a summary of the course Lacunary Polynomials over Finite Fields, given by Simeon Ball, from the University of London,
More informationMATRICES. a m,1 a m,n A =
MATRICES Matrices are rectangular arrays of real or complex numbers With them, we define arithmetic operations that are generalizations of those for real and complex numbers The general form a matrix of
More informationAn infinite family of Goethals-Seidel arrays
University of Wollongong Research Online Faculty of Informatics - Papers (Archive) Faculty of Engineering and Information Sciences 2005 An infinite family of Goethals-Seidel arrays M. Xia Central China
More informationOrthogonal arrays from Hermitian varieties
Orthogonal arrays from Hermitian varieties A Aguglia L Giuzzi Abstract A simple orthogonal array OA(q 2n 1, q 2n 2, q, 2) is constructed by using the action of a large subgroup of P GL(n + 1, q 2 ) on
More informationMAT Linear Algebra Collection of sample exams
MAT 342 - Linear Algebra Collection of sample exams A-x. (0 pts Give the precise definition of the row echelon form. 2. ( 0 pts After performing row reductions on the augmented matrix for a certain system
More informationUsing Matrix Analysis to Approach the Simulated Annealing Algorithm
Using Matrix Analysis to Approach the Simulated Annealing Algorithm Dennis I. Merino, Edgar N. Reyes Southeastern Louisiana University Hammond, LA 7040 and Carl Steidley Texas A&M University- Corpus Christi
More informationComplete permutation polynomials of monomial type
Complete permutation polynomials of monomial type Giovanni Zini (joint works with D. Bartoli, M. Giulietti and L. Quoos) (based on the work of thesis of E. Franzè) Università di Perugia Workshop BunnyTN
More informationUniform semi-latin squares and their Schur-optimality
Uniform semi-latin squares and their Schur-optimality Leonard H. Soicher School of Mathematical Sciences Queen Mary University of London Mile End Road, London E1 4NS, UK L.H.Soicher@qmul.ac.uk September
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 7: More on Householder Reflectors; Least Squares Problems Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 15 Outline
More informationarxiv: v1 [math.ac] 16 Nov 2018
On Hilbert Functions of Points in Projective Space and Structure of Graded Modules Damas Karmel Mgani arxiv:1811.06790v1 [math.ac] 16 Nov 2018 Department of Mathematics, College of Natural and Mathematical
More informationSampling Contingency Tables
Sampling Contingency Tables Martin Dyer Ravi Kannan John Mount February 3, 995 Introduction Given positive integers and, let be the set of arrays with nonnegative integer entries and row sums respectively
More informationSHABNAM AKHTARI AND JEFFREY D. VAALER
ON THE HEIGHT OF SOLUTIONS TO NORM FORM EQUATIONS arxiv:1709.02485v2 [math.nt] 18 Feb 2018 SHABNAM AKHTARI AND JEFFREY D. VAALER Abstract. Let k be a number field. We consider norm form equations associated
More informationThe Gaussian coefficient revisited
The Gaussian coefficient revisited Richard EHRENBORG and Margaret A. READDY Abstract We give new -(1+)-analogue of the Gaussian coefficient, also now as the -binomial which, lie the original -binomial
More informationAn Infeasible Interior-Point Algorithm with full-newton Step for Linear Optimization
An Infeasible Interior-Point Algorithm with full-newton Step for Linear Optimization H. Mansouri M. Zangiabadi Y. Bai C. Roos Department of Mathematical Science, Shahrekord University, P.O. Box 115, Shahrekord,
More informationCONSTRUCTION OF SLICED SPACE-FILLING DESIGNS BASED ON BALANCED SLICED ORTHOGONAL ARRAYS
Statistica Sinica 24 (2014), 1685-1702 doi:http://dx.doi.org/10.5705/ss.2013.239 CONSTRUCTION OF SLICED SPACE-FILLING DESIGNS BASED ON BALANCED SLICED ORTHOGONAL ARRAYS Mingyao Ai 1, Bochuan Jiang 1,2
More informationON PERMUTATION POLYNOMIALS OF PRESCRIBED SHAPE
ON PERMUTATION POLYNOMIALS OF PRESCRIBED SHAPE AMIR AKBARY, DRAGOS GHIOCA, AND QIANG WANG Abstract. We count permutation polynomials of F q which are sums of m + 2 monomials of prescribed degrees. This
More informationarxiv: v1 [math.co] 27 Jul 2015
Perfect Graeco-Latin balanced incomplete block designs and related designs arxiv:1507.07336v1 [math.co] 27 Jul 2015 Sunanda Bagchi Theoretical Statistics and Mathematics Unit Indian Statistical Institute
More informationSolutions 2016 AB Exam
Solutions 206 AB Exam Texas A&M High School Math Contest October 22, 206. If (x, y) is a point on the circle x 2 + y 2 = and the distance from (x, y) to (0, ) is 6 5, what is the value of y? ANSWER: 7
More informationDimensional exploration techniques for photonics
Dimensional exploration techniques for photonics Philippe A.J.G. CHEVALIER Ghent University FACULTEIT INGENIEURSWETENSCHAPPEN EN Who is this physicist? 2/30 FACULTEIT INGENIEURSWETENSCHAPPEN EN Who is
More informationTHREE IDENTITIES BETWEEN STIRLING NUMBERS AND THE STABILIZING CHARACTER SEQUENCE
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 60, October 1976 THREE IDENTITIES BETWEEN STIRLING NUMBERS AND THE STABILIZING CHARACTER SEQUENCE MICHAEL GILPIN Abstract. Let x denote the stabilizing
More informationUSING 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 informationc Igor Zelenko, Fall
c Igor Zelenko, Fall 2017 1 18: Repeated Eigenvalues: algebraic and geometric multiplicities of eigenvalues, generalized eigenvectors, and solution for systems of differential equation with repeated eigenvalues
More informationAlgebra II- Comprehensive/ Part A
Algebra II- Comprehensive/ Part A COURSE DESCRIPTION: This course builds upon algebraic concepts covered in Algebra I and prepares students for advanced-level courses. Students extend their knowledge and
More information{ independent variable some property or restriction about independent variable } where the vertical line is read such that.
Page 1 of 5 Introduction to Review Materials One key to Algebra success is identifying the type of work necessary to answer a specific question. First you need to identify whether you are dealing with
More informationOctober 28, S4.4 Theorems about Zeros of Polynomial Functions
MAT 171 Precalculus Algebra Dr. Claude Moore Cape Fear Community College CHAPTER 4: Polynomial and Rational Functions 4.1 Polynomial Functions and Models 4.2 Graphing Polynomial Functions 4.3 Polynomial
More informationLecture 15: Algebraic Geometry II
6.859/15.083 Integer Programming and Combinatorial Optimization Fall 009 Today... Ideals in k[x] Properties of Gröbner bases Buchberger s algorithm Elimination theory The Weak Nullstellensatz 0/1-Integer
More informationQUASI-ORTHOGONAL ARRAYS AND OPTIMAL FRACTIONAL FACTORIAL PLANS
Statistica Sinica 12(2002), 905-916 QUASI-ORTHOGONAL ARRAYS AND OPTIMAL FRACTIONAL FACTORIAL PLANS Kashinath Chatterjee, Ashish Das and Aloke Dey Asutosh College, Calcutta and Indian Statistical Institute,
More informationCommuting birth-and-death processes
Commuting birth-and-death processes Caroline Uhler Department of Statistics UC Berkeley (joint work with Steven N. Evans and Bernd Sturmfels) MSRI Workshop on Algebraic Statistics December 18, 2008 Birth-and-death
More informationOn rational numbers associated with arithmetic functions evaluated at factorials
On rational numbers associated with arithmetic functions evaluated at factorials Dan Baczkowski (joint work with M. Filaseta, F. Luca, and O. Trifonov) (F. Luca) Fix r Q, there are a finite number of positive
More informationLINEAR ALGEBRA 1, 2012-I PARTIAL EXAM 3 SOLUTIONS TO PRACTICE PROBLEMS
LINEAR ALGEBRA, -I PARTIAL EXAM SOLUTIONS TO PRACTICE PROBLEMS Problem (a) For each of the two matrices below, (i) determine whether it is diagonalizable, (ii) determine whether it is orthogonally diagonalizable,
More informationUNDERSTANDING RULER AND COMPASS CONSTRUCTIONS WITH FIELD THEORY
UNDERSTANDING RULER AND COMPASS CONSTRUCTIONS WITH FIELD THEORY ISAAC M. DAVIS Abstract. By associating a subfield of R to a set of points P 0 R 2, geometric properties of ruler and compass constructions
More informationStochastic Models: Markov Chains and their Generalizations
Scuola di Dottorato in Scienza ed Alta Tecnologia Dottorato in Informatica Universita di Torino Stochastic Models: Markov Chains and their Generalizations Gianfranco Balbo e Andras Horvath Outline Introduction
More informationRings of Residues. S. F. Ellermeyer. September 18, ; [1] m
Rings of Residues S F Ellermeyer September 18, 2006 If m is a positive integer, then we obtain the partition C = f[0] m ; [1] m ; : : : ; [m 1] m g of Z into m congruence classes (This is discussed in
More informationLinear equations in linear algebra
Linear equations in linear algebra Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra Pearson Collections Samy T. Linear
More informationNumerical Optimization
Linear Programming - Interior Point Methods Computer Science and Automation Indian Institute of Science Bangalore 560 012, India. NPTEL Course on Example 1 Computational Complexity of Simplex Algorithm
More informationLinear Algebra I for Science (NYC)
Element No. 1: To express concrete problems as linear equations. To solve systems of linear equations using matrices. Topic: MATRICES 1.1 Give the definition of a matrix, identify the elements and the
More informationTranslation Symmetry, Space Groups, Bloch functions, Fermi energy
Translation Symmetry, Space Groups, Bloch functions, Fermi energy Roberto Orlando and Silvia Casassa Università degli Studi di Torino July 20, 2015 School Ab initio Modelling of Solids (UniTo) Symmetry
More informationA Robust APTAS for the Classical Bin Packing Problem
A Robust APTAS for the Classical Bin Packing Problem Leah Epstein 1 and Asaf Levin 2 1 Department of Mathematics, University of Haifa, 31905 Haifa, Israel. Email: lea@math.haifa.ac.il 2 Department of Statistics,
More informationStable Process. 2. Multivariate Stable Distributions. July, 2006
Stable Process 2. Multivariate Stable Distributions July, 2006 1. Stable random vectors. 2. Characteristic functions. 3. Strictly stable and symmetric stable random vectors. 4. Sub-Gaussian random vectors.
More information