Abstract. Using graphs to find optimal block designs. Problem 1: Factorial design. Conference on Design of Experiments, Tianjin, June 2006
|
|
- Claire Waters
- 5 years ago
- Views:
Transcription
1 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, and either can be used. Using graphs to find optimal block designs R. A. Bailey University of St Andrews / QMUL (emerita) A block design is D-optimal if it maximizes the number of spanning trees; it is A-optimal if it minimizes the total of the pairwise resistances when the graph is thought of as an electrical network. Special session for Ching-Shui Cheng, 20 June 20 I shall report on some surprising results about optimal designs when replication is very low. 2/ / Conference on Design of Experiments, Tianjin, June 2006 Problem : Factorial design Problem There are several treatment factors, with various numbers of levels. There may not be room to include all combinations of these levels. There are several inherent factors (also called block factors) on the experimental units. The inherent factors may pose some constraints on how treatment factors can be applied. So how do we set about designing the experiment? Solution Use Desmond Patterson s design key. / Problem 2: Optimal incomplete-block design / Two papers about the second problem Problem The design must be in blocks, but each block is too small to contain every treatment. We want variance to be small, as measured by the D-criterion (minimizing the volume of the ellpsoid of confidence). How do we set about finding a D-optimal block design? Solution Represent the incomplete-block design as a graph (with vertices and edges), and maximize the number of spanning trees in the graph. 5/ I Ch ing Shui Che ng: Maximizing the total number of spanning trees in a graph: two related problems in graph theory and optimum design theory. Journal of Combinatorial Theory, Series B (98), I Ching-Shui Cheng: Graph and optimum design theories some connections and examples. Bulletin of the International Statistical Institute 9 (98), part, /
2 Some early papers Year RAB (design key) (design key) ; Two memorial conferences for R. C. Bose, India, 988 CSC optimal optimal optimal ; optimal and graphs Calcutta meeting on Combinatorial Mathematics and Applications: CSC spoke on On the optimality of (M.S)-optimal designs in large systems. and graphs; with trend with trend Delhi meeting on Probability, Statistics and Design of Experiments: RAB spoke on Cyclic designs and designs. / Two joint papers I C.-S. Cheng and R. A. Bailey: Optimality of some two-associate-class partially balanced incomplete-block designs. Annals of Statistics 9 (99), I R. A. Bailey, Ching-Shui Cheng and Patricia Kipnis: Construction of trend-resistant designs. Statistica Sinica 2 (992), 9. 8/ Some later papers and books Year RAB CSC ; optimal optimal optimal and graphs optimal and graphs optimal and graphs ; B (design key) 9/ With Jo Kunert, Dortmund, April / What makes a block design good for experiments? I have v treatments that I want to compare. I have b blocks. Each block has space for k treatments (not necessarily distinct). How should I choose a block design? / 2/
3 Two designs with v = 5, b =, k = : which is better? Two designs with v = 5, b =, k = : which is better? Conventions: columns are blocks; order of treatments within each block is irrelevant; order of blocks is irrelevant replications differ by queen-bee design binary non-binary A design is binary if no treatment occurs more than once in any block. The replication of a treatment is its number of occurrences. A design is a queen-bee design if there is a treatment that occurs in every block. / / Experimental units and incidence matrix There are bk experimental units. If ω is an experimental unit, put f (ω) = treatment on ω g(ω) = block containing ω. For i =,..., v and j =,..., b, let n ij = {ω : f (ω) = i and g(ω) = j} = number of experimental units in block j which have treatment i. The v b incidence matrix N has entries n ij. The G of a block design has one vertex for each treatment, one vertex for each block, one edge for each experimental unit, with edge ω joining vertex f (ω) to vertex g(ω). It is a bipartite graph, with n ij edges between treatment-vertex i and block-vertex j. 5/ 6/ Example : v =, b = k = Example 2: v = 8, b =, k = / 8/
4 Concurrence graph Example : v =, b = k = The concurrence graph G of a block design has one vertex for each treatment, one edge for each unordered pair α, ω, with α = ω, g(α) = g(ω) and f (α) = f (ω): this edge joins vertices f (α) and f (ω). There are no loops. If i = j then the number of edges between vertices i and j is λ ij = b s= n is n js ; this is called the concurrence of i and j, and is the (i, j)-entry of Λ = NN. 9/ 2 2 concurrence graph can recover design may have more symmetry more vertices more edges if k = 2 more edges if k 20/ Example 2: v = 8, b =, k = Laplacian matrices concurrence graph 2/ The Laplacian matrix L of the concurrence graph G is a v v matrix with (i, j)-entry as follows: if i = j then L ij = (number of edges between i and j) = λ ij ; L ii = valency of i = λ ij. j =i The Laplacian matrix L of the G is a (v + b) (v + b) matrix with (i, j)-entry as follows: L ii = valency { of i k if i is a block = replication r i of i if i is a treatment if i = j then L ij = (number of edges between i and j) 0 if i and j are both treatments = 0 if i and j are both blocks if i is a treatment and j is a block, or vice versa. n ij 22/ Connectivity Generalized inverse All row-sums of L and of L are zero, so both matrices have 0 as eigenvalue on the appropriate all- vector. Theorem The following are equivalent.. 0 is a simple eigenvalue of L; 2. G is a connected graph;. G is a connected graph;. 0 is a simple eigenvalue of L; 5. the design is connected in the sense that all differences between treatments can be estimated. From now on, assume connectivity. Call the remaining eigenvalues non-trivial. They are all non-negative. 2/ Under the assumption of connectivity, the Moore Penrose generalized inverse L of L is defined by ( L = L + ) v J v v J v, where J v is the v v all- matrix. (The matrix v J v is the orthogonal projector onto the null space of L.) The Moore Penrose generalized inverse L of L is defined similarly. 2/
5 Estimation and variance How do we calculate variance? We measure the response Y ω on each experimenal unit ω. If experimental unit ω has treatment i and is in block m (f (ω) = i and g(ω) = m), then we assume that Y ω = τ i + β m + random noise. We want to estimate contrasts i x i τ i with i x i = 0. In particular, we want to estimate all the simple differences τ i τ j. Put V ij = variance of the best linear unbiased estimator for τ i τ j. Theorem (Standard linear model theory) Assume that all the noise is independent, with variance σ 2. If i x i = 0, then the variance of the best linear unbiased estimator of i x i τ i is equal to (x L x)kσ 2. In particular, the variance of the best linear unbiased estimator of the simple difference τ i τ j is ( ) V ij = L ii + L jj 2L ij kσ 2. We want all the V ij to be small. 25/ 26/... Or we can use the Electrical networks Theorem The variance of the best linear unbiased estimator of the simple difference τ i τ j is ) V ij = ( L ii + L jj 2 L ij σ 2. 2/ We can consider the concurrence graph G as an electrical network with a -ohm resistance in each edge. Connect a -volt battery between vertices i and j. Current flows in the network, according to these rules.. Ohm s Law: In every edge, voltage drop = current resistance = current. 2. Kirchhoff s Voltage Law: The total voltage drop from one vertex to any other vertex is the same no matter which path we take from one to the other.. Kirchhoff s Current Law: At every vertex which is not connected to the battery, the total current coming in is equal to the total current going out. Find the total current I from i to j, then use Ohm s Law to define the effective resistance R ij between i and j as /I. 28/ Electrical networks: variance Example calculation: v = 2, b = 6, k = Theorem The effective resistance R ij between vertices i and j in G is So ( ) R ij = L ii + L jj 2L ij. V ij = R ij kσ 2. Effective resistances are easy to calculate without matrix inversion if the graph is sparse. 29/ 5 5 [20] [0] [] [0] [28] [0] [] V = I = 6 R = 6 0/
6 ... Or we can use the Example 2 yet again: v = 8, b =, k = V = 2 I = 8 R = If i and j are treatment vertices in the G and R ij is the effective resistance between them in G then V ij = R ij σ [0] [5] [2] concurrence graph / 2/ Average pairwise variance Optimality The variance of the best linear unbiased estimator of the simple difference τ i τ j is ( ) V ij = L ii + L jj 2L ij kσ 2. Put V = average value of the V ij. Then V = 2kσ2 Tr(L ) v = 2kσ 2, harmonic mean of θ,..., θ v where θ,..., θ v are the nontrivial eigenvalues of L. The design is called A-optimal if it minimizes the average of the variances V ij ; equivalently, it maximizes the harmonic mean of the non-trivial eigenvalues of the Laplacian matrix L; D-optimal if it minimizes the volume of the confidence ellipsoid for (τ,..., τ v ); equivalently, it maximizes the geometric mean of the non-trivial eigenvalues of the Laplacian matrix L; over all block designs with block size k and the given v and b. / / D-optimality: spanning trees What about the? A spanning tree for the graph is a collection of edges of the graph which form a tree (connected graph with no cycles) and which include every vertex. Cheng (98), after Gaffke (98), after Kirchhoff (8): product of non-trivial eigenvalues of L = v number of spanning trees. So a design is D-optimal if and only if its concurrence graph G has the maximal number of spanning trees. This is easy to calculate by hand when the graph is sparse. Theorem (Gaffke) Let G and G be the concurrence graph and for a connected incomplete-block design for v treatments in b blocks of size k. Then the number of spanning trees for G is equal to k b v+ times the number of spanning trees for G. So a block design is D-optimal if and only if its maximizes the number of spanning trees. If v > b it is easier to count spanning trees in the than in the concurrence graph. 5/ 6/
7 Example 2 one last time: v = 8, b =, k = concurrence graph 8 spanning trees 26 spanning trees Large blocks; many unreplicated treatments b blocks k plots v varieties. n plots. bn varieties all single replication whole design Whole design has v + bn varieties in b blocks of size k + n; the subdesign Γ has v core varieties in b blocks of size k; call the remaining varieties orphans. / 8/ : 0 + 5n treatments in 5 blocks of + n plots Pairwise resistance A n. A 2 A A n 5 6 B B n C C n D D n 9 0 E E n B B n 2 5 C. 6 C n 8 9 E n D. 0. E D n 9/ A n. A Resistance(A, A 2 ) = 2 B B n 2 5 C. 6 C n 8 9 E n D. 0. E D n Resistance(A, B ) = 2 + Resistance(block A, block B) in Γ Resistance(A, 8) = + Resistance(block A, 8) in Γ Resistance(, 8) = Resistance(, 8) in Γ 0/ Sum of the pairwise variances A less obvious consequence Theorem (cf Herzberg and Jarrett, 200) The sum of the variances of treatment differences in where = constant + V + nv + n 2 V 2, V = the sum of the variances of treatment differences in Γ V 2 = the sum of the variances of block differences in Γ V = the sum of the variances of sums of one treatment and one block in Γ. Consequence If n or b is large, and we want an A-optimal design, it may be best to make Γ a complete block design for k controls, even if there is no interest in comparisons between new treatments and controls, or between controls. (If Γ is equi-replicate then V 2 and V are both increasing functions of V.) Consequence For a given choice of k, make Γ as efficient as possible. / 2/
8 Spanning trees D-optimality under very low replication A spanning tree for the is a collection edges which provides a unique path between every pair of vertices. A n. A B B n 2 5 C. 6 C n 8 9 E n D. 0. E D n Consequence The whole design is D-optimal for v + bn treatments in b blocks of size k + n if and only if the core design Γ is D-optimal for v treatments in b blocks of size k. Consequence Even when n or b is large, D-optimal designs do not include uninteresting controls. The orphans make no difference to the number of spanning trees for the. / /
Abstract. Designs for variety trials with very low replication. How do we allow for variation between the plots? Variety Testing
Abstract Designs for variety trials with very low replication R A Bailey University of St Andrews QMUL (emerita) TU Dortmund, June 05 In the early stages of testing new varieties, it is common that there
More informationCombinatorics of optimal designs
Combinatorics of optimal designs R A Bailey and Peter J Cameron pjcameron@qmulacuk British Combinatorial Conference, St Andrews, July 2009 Mathematicians and statisticians There is a very famous joke about
More informationLaplacian eigenvalues and optimality: II. The Laplacian of a graph. R. A. Bailey and Peter Cameron
Laplacian eigenvalues and optimality: II. The Laplacian of a graph R. A. Bailey and Peter Cameron London Taught Course Centre, June 2012 The Laplacian of a graph This lecture will be about the Laplacian
More informationCombinatorics of optimal designs
Combinatorics of optimal designs R. A. Bailey and Peter J. Cameron Abstract To a combinatorialist, a design is usually a 2-design or balanced incompleteblock design. However, 2-designs do not necessarily
More informationAbstract. A substitute for the non-existent affine plane of order 6. Chapter 1. Outline. Square lattice designs for n 2 treatments in rn blocks of n
A substitute for the non-existent affine plane of order R. A. Bailey University of St Andrews QMUL (emerita) Algebra and Combinatorics seminar, University of St Andrews, 1 February 018 Joint work with
More informationLinear estimation in models based on a graph
Linear Algebra and its Applications 302±303 (1999) 223±230 www.elsevier.com/locate/laa Linear estimation in models based on a graph R.B. Bapat * Indian Statistical Institute, New Delhi 110 016, India Received
More informationCONSTRUCTION OF REGULAR GRAPH DESIGNS AND ITS GRAPHICAL REPRESENTATION
CHAPTER 3 CONSTRUCTION OF REGULAR GRAPH DESIGNS AND ITS GRAPHICAL REPRESENTATION 3.1 Introduction 3.2 Historical Review 3.3 Preliminary Results 3.4 Construction of RG Designs using a BIB design with b
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 informationThe Matrix-Tree Theorem
The Matrix-Tree Theorem Christopher Eur March 22, 2015 Abstract: We give a brief introduction to graph theory in light of linear algebra. Our results culminates in the proof of Matrix-Tree Theorem. 1 Preliminaries
More informationOn some matrices related to a tree with attached graphs
On some matrices related to a tree with attached graphs R. B. Bapat Indian Statistical Institute New Delhi, 110016, India fax: 91-11-26856779, e-mail: rbb@isid.ac.in Abstract A tree with attached graphs
More informationLaplacians of Graphs, Spectra and Laplacian polynomials
Laplacians of Graphs, Spectra and Laplacian polynomials Lector: Alexander Mednykh Sobolev Institute of Mathematics Novosibirsk State University Winter School in Harmonic Functions on Graphs and Combinatorial
More informationEfficiency of NNBD and NNBIBD using autoregressive model
2018; 3(3): 133-138 ISSN: 2456-1452 Maths 2018; 3(3): 133-138 2018 Stats & Maths www.mathsjournal.com Received: 17-03-2018 Accepted: 18-04-2018 S Saalini Department of Statistics, Loyola College, Chennai,
More information11 Block Designs. Linear Spaces. Designs. By convention, we shall
11 Block Designs Linear Spaces In this section we consider incidence structures I = (V, B, ). always let v = V and b = B. By convention, we shall Linear Space: We say that an incidence structure (V, B,
More information3.5 Efficiency factors
3.5. EFFICIENCY FACTORS 63 3.5 Efficiency factors For comparison we consider a complete-block design where the variance of each response is σ CBD. In such a design, Λ = rj Θ and k = t, so Equation (3.3)
More informationOn the adjacency matrix of a block graph
On the adjacency matrix of a block graph R. B. Bapat Stat-Math Unit Indian Statistical Institute, Delhi 7-SJSS Marg, New Delhi 110 016, India. email: rbb@isid.ac.in Souvik Roy Economics and Planning Unit
More informationLatin squares: Equivalents and equivalence
Latin squares: Equivalents and equivalence 1 Introduction This essay describes some mathematical structures equivalent to Latin squares and some notions of equivalence of such structures. According to
More information18.312: Algebraic Combinatorics Lionel Levine. Lecture 19
832: Algebraic Combinatorics Lionel Levine Lecture date: April 2, 20 Lecture 9 Notes by: David Witmer Matrix-Tree Theorem Undirected Graphs Let G = (V, E) be a connected, undirected graph with n vertices,
More informationLaplacian Energy of Graphs
Laplacian Energy of Graphs Sigma Xi Research Exhibition Claire T. Carney Library Uniersity of Massachusetts Dartmouth April 28 th & 29 th, 2009 Leanne R. Silia & Gary E. Dais Department of Mathematics
More informationConstruction of Pair-wise Balanced Design
Journal of Modern Applied Statistical Methods Volume 15 Issue 1 Article 11 5-2016 onstruction of Pair-wise Balanced Design Rajarathinam Arunachalam Manonmaniam Sundaranar University, arrathinam@yahoo.com
More informationLecture: Modeling graphs with electrical networks
Stat260/CS294: Spectral Graph Methods Lecture 16-03/17/2015 Lecture: Modeling graphs with electrical networks Lecturer: Michael Mahoney Scribe: Michael Mahoney Warning: these notes are still very rough.
More informationResistance distance in wheels and fans
Resistance distance in wheels and fans R B Bapat Somit Gupta February 4, 009 Abstract: The wheel graph is the join of a single vertex and a cycle, while the fan graph is the join of a single vertex and
More informationFinding normalized and modularity cuts by spectral clustering. Ljubjana 2010, October
Finding normalized and modularity cuts by spectral clustering Marianna Bolla Institute of Mathematics Budapest University of Technology and Economics marib@math.bme.hu Ljubjana 2010, October Outline Find
More informationHadamard and conference matrices
Hadamard and conference matrices Peter J. Cameron University of St Andrews & Queen Mary University of London Mathematics Study Group with input from Rosemary Bailey, Katarzyna Filipiak, Joachim Kunert,
More informationData Analysis and Manifold Learning Lecture 9: Diffusion on Manifolds and on Graphs
Data Analysis and Manifold Learning Lecture 9: Diffusion on Manifolds and on Graphs Radu Horaud INRIA Grenoble Rhone-Alpes, France Radu.Horaud@inrialpes.fr http://perception.inrialpes.fr/ Outline of Lecture
More informationHadamard and conference matrices
Hadamard and conference matrices Peter J. Cameron University of St Andrews & Queen Mary University of London Mathematics Study Group with input from Rosemary Bailey, Katarzyna Filipiak, Joachim Kunert,
More informationCombining Association Schemes
Chapter 4 Combining Association Schemes 4.1 Tensor products In this chapter we take an association scheme on Ω 1 and an association scheme on Ω, and combine them to obtain an association scheme on Ω 1
More informationE(s 2 )-OPTIMAL SUPERSATURATED DESIGNS
Statistica Sinica 7(1997), 929-939 E(s 2 )-OPTIMAL SUPERSATURATED DESIGNS Ching-Shui Cheng University of California, Berkeley Abstract: Tang and Wu (1997) derived a lower bound on the E(s 2 )-value of
More informationNotes 6 : First and second moment methods
Notes 6 : First and second moment methods Math 733-734: Theory of Probability Lecturer: Sebastien Roch References: [Roc, Sections 2.1-2.3]. Recall: THM 6.1 (Markov s inequality) Let X be a non-negative
More informationLecture 13: Spectral Graph Theory
CSE 521: Design and Analysis of Algorithms I Winter 2017 Lecture 13: Spectral Graph Theory Lecturer: Shayan Oveis Gharan 11/14/18 Disclaimer: These notes have not been subjected to the usual scrutiny reserved
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 informationSpanning Trees of Shifted Simplicial Complexes
Art Duval (University of Texas at El Paso) Caroline Klivans (Brown University) Jeremy Martin (University of Kansas) Special Session on Extremal and Probabilistic Combinatorics University of Nebraska, Lincoln
More informationThe QR Factorization
The QR Factorization How to Make Matrices Nicer Radu Trîmbiţaş Babeş-Bolyai University March 11, 2009 Radu Trîmbiţaş ( Babeş-Bolyai University) The QR Factorization March 11, 2009 1 / 25 Projectors A projector
More informationThe maximum size of a partial spread in H(4n + 1,q 2 ) is q 2n+1 + 1
The maximum size of a partial spread in H(4n +, 2 ) is 2n+ + Frédéric Vanhove Dept. of Pure Mathematics and Computer Algebra, Ghent University Krijgslaan 28 S22 B-9000 Ghent, Belgium fvanhove@cage.ugent.be
More informationLINEAR SPACES. Define a linear space to be a near linear space in which any two points are on a line.
LINEAR SPACES Define a linear space to be a near linear space in which any two points are on a line. A linear space is an incidence structure I = (P, L) such that Axiom LS1: any line is incident with at
More informationOn sums of eigenvalues, and what they reveal about graphs and manifolds
On sums of eigenvalues, and what they reveal about graphs and manifolds Evans Harrell Georgia Tech Université F. Rabelais www.math.gatech.edu/~harrell Tours le 30 janvier, 2014 Copyright 2014 by Evans
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 informationRoot systems and optimal block designs
Root systems and optimal block designs Peter J. Cameron School of Mathematical Sciences Queen Mary, University of London Mile End Road London E1 4NS, UK p.j.cameron@qmul.ac.uk Abstract Motivated by a question
More informationON MULTI-AVOIDANCE OF RIGHT ANGLED NUMBERED POLYOMINO PATTERNS
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 4 (2004), #A21 ON MULTI-AVOIDANCE OF RIGHT ANGLED NUMBERED POLYOMINO PATTERNS Sergey Kitaev Department of Mathematics, University of Kentucky,
More informationRESISTANCE DISTANCE IN WHEELS AND FANS
Indian J Pure Appl Math, 41(1): 1-13, February 010 c Indian National Science Academy RESISTANCE DISTANCE IN WHEELS AND FANS R B Bapat 1 and Somit Gupta Indian Statistical Institute, New Delhi 110 016,
More informationOn the inverse matrix of the Laplacian and all ones matrix
On the inverse matrix of the Laplacian and all ones matrix Sho Suda (Joint work with Michio Seto and Tetsuji Taniguchi) International Christian University JSPS Research Fellow PD November 21, 2012 Sho
More informationStrongly Regular Decompositions of the Complete Graph
Journal of Algebraic Combinatorics, 17, 181 201, 2003 c 2003 Kluwer Academic Publishers. Manufactured in The Netherlands. Strongly Regular Decompositions of the Complete Graph EDWIN R. VAN DAM Edwin.vanDam@uvt.nl
More informationLecture 24: Element-wise Sampling of Graphs and Linear Equation Solving. 22 Element-wise Sampling of Graphs and Linear Equation Solving
Stat260/CS294: Randomized Algorithms for Matrices and Data Lecture 24-12/02/2013 Lecture 24: Element-wise Sampling of Graphs and Linear Equation Solving Lecturer: Michael Mahoney Scribe: Michael Mahoney
More informationOn the sum of two largest eigenvalues of a symmetric matrix
On the sum of two largest eigenvalues of a symmetric matrix Javad Ebrahimi B. Department of Mathematics Simon Fraser University Burnaby, B.C. V5A 1S6 Vladimir Nikiforov Department of Mathematical Sci.
More informationHanna Furmańczyk EQUITABLE COLORING OF GRAPH PRODUCTS
Opuscula Mathematica Vol. 6 No. 006 Hanna Furmańczyk EQUITABLE COLORING OF GRAPH PRODUCTS Abstract. A graph is equitably k-colorable if its vertices can be partitioned into k independent sets in such a
More informationSpectral Clustering. Zitao Liu
Spectral Clustering Zitao Liu Agenda Brief Clustering Review Similarity Graph Graph Laplacian Spectral Clustering Algorithm Graph Cut Point of View Random Walk Point of View Perturbation Theory Point of
More informationLinear graph theory. Basic definitions of linear graphs
Linear graph theory Linear graph theory, a branch of combinatorial mathematics has proved to be a useful tool for the study of large or complex systems. Leonhard Euler wrote perhaps the first paper on
More informationCombinatorial and physical content of Kirchhoff polynomials
Combinatorial and physical content of Kirchhoff polynomials Karen Yeats May 19, 2009 Spanning trees Let G be a connected graph, potentially with multiple edges and loops in the sense of a graph theorist.
More informationKevin James. MTHSC 3110 Section 2.1 Matrix Operations
MTHSC 3110 Section 2.1 Matrix Operations Notation Let A be an m n matrix, that is, m rows and n columns. We ll refer to the entries of A by their row and column indices. The entry in the i th row and j
More informationPreliminaries and Complexity Theory
Preliminaries and Complexity Theory Oleksandr Romanko CAS 746 - Advanced Topics in Combinatorial Optimization McMaster University, January 16, 2006 Introduction Book structure: 2 Part I Linear Algebra
More informationConstruction of Partially Balanced Incomplete Block Designs
International Journal of Statistics and Systems ISS 0973-675 Volume, umber (06), pp. 67-76 Research India Publications http://www.ripublication.com Construction of Partially Balanced Incomplete Block Designs
More informationDeterminant of the distance matrix of a tree with matrix weights
Determinant of the distance matrix of a tree with matrix weights R. B. Bapat Indian Statistical Institute New Delhi, 110016, India fax: 91-11-26856779, e-mail: rbb@isid.ac.in Abstract Let T be a tree with
More informationApplying block intersection polynomials to study graphs and designs
Applying block intersection polynomials to study graphs and designs Leonard Soicher Queen Mary University of London CoCoA15, Colorado State University, Fort Collins, July 2015 Leonard Soicher (QMUL) Block
More informationBLOCK DESIGNS WITH FACTORIAL STRUCTURE
BLOCK DESIGNS WITH ACTORIAL STRUCTURE V.K. Gupta and Rajender Parsad I.A.S.R.I., Library Avenue, New Delhi - 0 0 The purpose of this talk is to expose the participants to the interesting work on block
More informationWalks, Springs, and Resistor Networks
Spectral Graph Theory Lecture 12 Walks, Springs, and Resistor Networks Daniel A. Spielman October 8, 2018 12.1 Overview In this lecture we will see how the analysis of random walks, spring networks, and
More informationThis appendix provides a very basic introduction to linear algebra concepts.
APPENDIX Basic Linear Algebra Concepts This appendix provides a very basic introduction to linear algebra concepts. Some of these concepts are intentionally presented here in a somewhat simplified (not
More informationData Analysis and Manifold Learning Lecture 3: Graphs, Graph Matrices, and Graph Embeddings
Data Analysis and Manifold Learning Lecture 3: Graphs, Graph Matrices, and Graph Embeddings Radu Horaud INRIA Grenoble Rhone-Alpes, France Radu.Horaud@inrialpes.fr http://perception.inrialpes.fr/ Outline
More informationGraphs and Their Applications (7)
11111111 I I I Graphs and Their Applications (7) by K.M. Koh* Department of Mathematics National University of Singapore, Singapore 117543 F.M. Dong and E.G. Tay Mathematics and Mathematics Education National
More informationLecture 10. Sublinear Time Algorithms (contd) CSC2420 Allan Borodin & Nisarg Shah 1
Lecture 10 Sublinear Time Algorithms (contd) CSC2420 Allan Borodin & Nisarg Shah 1 Recap Sublinear time algorithms Deterministic + exact: binary search Deterministic + inexact: estimating diameter in a
More informationFactorial designs. Experiments
Chapter 5: Factorial designs Petter Mostad mostad@chalmers.se Experiments Actively making changes and observing the result, to find causal relationships. Many types of experimental plans Measuring response
More informationPage 52. Lecture 3: Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 2008/10/03 Date Given: 2008/10/03
Page 5 Lecture : Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 008/10/0 Date Given: 008/10/0 Inner Product Spaces: Definitions Section. Mathematical Preliminaries: Inner
More information1 Counting spanning trees: A determinantal formula
Math 374 Matrix Tree Theorem Counting spanning trees: A determinantal formula Recall that a spanning tree of a graph G is a subgraph T so that T is a tree and V (G) = V (T ) Question How many distinct
More informationOn minors of the compound matrix of a Laplacian
On minors of the compound matrix of a Laplacian R. B. Bapat 1 Indian Statistical Institute New Delhi, 110016, India e-mail: rbb@isid.ac.in August 28, 2013 Abstract: Let L be an n n matrix with zero row
More informationDissertation Defense
Clustering Algorithms for Random and Pseudo-random Structures Dissertation Defense Pradipta Mitra 1 1 Department of Computer Science Yale University April 23, 2008 Mitra (Yale University) Dissertation
More informationOn Euclidean distance matrices
On Euclidean distance matrices R. Balaji and R. B. Bapat Indian Statistical Institute, New Delhi, 110016 November 19, 2006 Abstract If A is a real symmetric matrix and P is an orthogonal projection onto
More informationNotes on the Matrix-Tree theorem and Cayley s tree enumerator
Notes on the Matrix-Tree theorem and Cayley s tree enumerator 1 Cayley s tree enumerator Recall that the degree of a vertex in a tree (or in any graph) is the number of edges emanating from it We will
More informationUSING GRAPHS AND GAMES TO GENERATE CAP SET BOUNDS
USING GRAPHS AND GAMES TO GENERATE CAP SET BOUNDS JOSH ABBOTT AND TREVOR MCGUIRE Abstract. Let F 3 be the field with 3 elements and consider the k- dimensional affine space, F k 3, over F 3. A line of
More informationBlock designs and statistics
Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent
More informationCombinatorial optimization problems
Combinatorial optimization problems Heuristic Algorithms Giovanni Righini University of Milan Department of Computer Science (Crema) Optimization In general an optimization problem can be formulated as:
More informationCSC Linear Programming and Combinatorial Optimization Lecture 10: Semidefinite Programming
CSC2411 - Linear Programming and Combinatorial Optimization Lecture 10: Semidefinite Programming Notes taken by Mike Jamieson March 28, 2005 Summary: In this lecture, we introduce semidefinite programming
More informationTHE idea of network coding over error-free networks,
Path Gain Algebraic Formulation for the Scalar Linear Network Coding Problem Abhay T. Subramanian and Andrew Thangaraj, Member, IEEE arxiv:080.58v [cs.it] 9 Jun 00 Abstract In the algebraic view, the solution
More informationβ α : α β We say that is symmetric if. One special symmetric subset is the diagonal
Chapter Association chemes. Partitions Association schemes are about relations between pairs of elements of a set Ω. In this book Ω will always be finite. Recall that Ω Ω is the set of ordered pairs of
More informationOptimal and Efficient Designs for Functional Magnetic Resonance Imaging (fmri) Experiments
Optimal and Efficient Designs for Functional Magnetic Resonance Imaging (fmri) Experiments Ming-Hung (Jason) Kao Arizona State University with Ching-Shui Cheng and Federick Kin Hing Phoa Introduction fmri
More informationDesigning Information Devices and Systems I Discussion 13B
EECS 6A Fall 7 Designing Information Devices and Systems I Discussion 3B. Orthogonal Matching Pursuit Lecture Orthogonal Matching Pursuit (OMP) algorithm: Inputs: A set of m songs, each of length n: S
More informationX -1 -balance of some partially balanced experimental designs with particular emphasis on block and row-column designs
DOI:.55/bile-5- Biometrical Letters Vol. 5 (5), No., - X - -balance of some partially balanced experimental designs with particular emphasis on block and row-column designs Ryszard Walkowiak Department
More information1. Ignoring case, extract all unique words from the entire set of documents.
CS 378 Introduction to Data Mining Spring 29 Lecture 2 Lecturer: Inderjit Dhillon Date: Jan. 27th, 29 Keywords: Vector space model, Latent Semantic Indexing(LSI), SVD 1 Vector Space Model The basic idea
More informationCOUNTING AND ENUMERATING SPANNING TREES IN (di-) GRAPHS
COUNTING AND ENUMERATING SPANNING TREES IN (di-) GRAPHS Xuerong Yong Version of Feb 5, 9 pm, 06: New York Time 1 1 Introduction A directed graph (digraph) D is a pair (V, E): V is the vertex set of D,
More informationDISCRETIZED CONFIGURATIONS AND PARTIAL PARTITIONS
DISCRETIZED CONFIGURATIONS AND PARTIAL PARTITIONS AARON ABRAMS, DAVID GAY, AND VALERIE HOWER Abstract. We show that the discretized configuration space of k points in the n-simplex is homotopy equivalent
More informationAlgebraic Methods in Combinatorics
Algebraic Methods in Combinatorics Po-Shen Loh 27 June 2008 1 Warm-up 1. (A result of Bourbaki on finite geometries, from Răzvan) Let X be a finite set, and let F be a family of distinct proper subsets
More informationOn some incidence structures constructed from groups and related codes
On some incidence structures constructed from groups and related codes Dean Crnković Department of Mathematics University of Rijeka Croatia Algebraic Combinatorics and Applications The first annual Kliakhandler
More information1 Basic Combinatorics
1 Basic Combinatorics 1.1 Sets and sequences Sets. A set is an unordered collection of distinct objects. The objects are called elements of the set. We use braces to denote a set, for example, the set
More informationTree sets. Reinhard Diestel
1 Tree sets Reinhard Diestel Abstract We study an abstract notion of tree structure which generalizes treedecompositions of graphs and matroids. Unlike tree-decompositions, which are too closely linked
More informationChapter 10 Combinatorial Designs
Chapter 10 Combinatorial Designs BIBD Example (a,b,c) (a,b,d) (a,c,e) (a,d,f) (a,e,f) (b,c,f) (b,d,e) (b,e,f) (c,d,e) (c,d,f) Here are 10 subsets of the 6 element set {a, b, c, d, e, f }. BIBD Definition
More informationTopics in Graph Theory
Topics in Graph Theory September 4, 2018 1 Preliminaries A graph is a system G = (V, E) consisting of a set V of vertices and a set E (disjoint from V ) of edges, together with an incidence function End
More informationBBM402-Lecture 20: LP Duality
BBM402-Lecture 20: LP Duality Lecturer: Lale Özkahya Resources for the presentation: https://courses.engr.illinois.edu/cs473/fa2016/lectures.html An easy LP? which is compact form for max cx subject to
More informationKernels of Directed Graph Laplacians. J. S. Caughman and J.J.P. Veerman
Kernels of Directed Graph Laplacians J. S. Caughman and J.J.P. Veerman Department of Mathematics and Statistics Portland State University PO Box 751, Portland, OR 97207. caughman@pdx.edu, veerman@pdx.edu
More informationRoberto s Notes on Linear Algebra Chapter 10: Eigenvalues and diagonalization Section 3. Diagonal matrices
Roberto s Notes on Linear Algebra Chapter 10: Eigenvalues and diagonalization Section 3 Diagonal matrices What you need to know already: Basic definition, properties and operations of matrix. What you
More informationCombinatorial Optimisation, Problems I, Solutions Term /2015
/0/205 Combinatorial Optimisation, Problems I, Solutions Term 2 204/205 Tomasz Tkocz, t (dot) tkocz (at) warwick (dot) ac (dot) uk 3. By Problem 2, any tree which is not a single vertex has at least 2
More informationApproximate Counting and Markov Chain Monte Carlo
Approximate Counting and Markov Chain Monte Carlo A Randomized Approach Arindam Pal Department of Computer Science and Engineering Indian Institute of Technology Delhi March 18, 2011 April 8, 2011 Arindam
More informationThe crossdes Package
Date 2007-04-16 The crossdes Package April 16, 2007 Title Design and Randomization in Crossover Studies Version 1.0-7 Author Martin Oliver Sailer Depends AlgDesign, gtools,
More informationAn Algorithmist s Toolkit September 10, Lecture 1
18.409 An Algorithmist s Toolkit September 10, 2009 Lecture 1 Lecturer: Jonathan Kelner Scribe: Jesse Geneson (2009) 1 Overview The class s goals, requirements, and policies were introduced, and topics
More informationCurrent; Forest Tree Theorem; Potential Functions and their Bounds
April 13, 2008 Franklin Kenter Current; Forest Tree Theorem; Potential Functions and their Bounds 1 Introduction In this section, we will continue our discussion on current and induced current. Review
More informationLanguage Acquisition and Parameters: Part II
Language Acquisition and Parameters: Part II Matilde Marcolli CS0: Mathematical and Computational Linguistics Winter 205 Transition Matrices in the Markov Chain Model absorbing states correspond to local
More informationAn Introduction to Spectral Graph Theory
An Introduction to Spectral Graph Theory Mackenzie Wheeler Supervisor: Dr. Gary MacGillivray University of Victoria WheelerM@uvic.ca Outline Outline 1. How many walks are there from vertices v i to v j
More informationAnswers in blue. If you have questions or spot an error, let me know. 1. Find all matrices that commute with A =. 4 3
Answers in blue. If you have questions or spot an error, let me know. 3 4. Find all matrices that commute with A =. 4 3 a b If we set B = and set AB = BA, we see that 3a + 4b = 3a 4c, 4a + 3b = 3b 4d,
More informationTheorems of Erdős-Ko-Rado type in polar spaces
Theorems of Erdős-Ko-Rado type in polar spaces Valentina Pepe, Leo Storme, Frédéric Vanhove Department of Mathematics, Ghent University, Krijgslaan 28-S22, 9000 Ghent, Belgium Abstract We consider Erdős-Ko-Rado
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 informationUse 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[y i α βx i ] 2 (2) Q = i=1
Least squares fits This section has no probability in it. There are no random variables. We are given n points (x i, y i ) and want to find the equation of the line that best fits them. We take the equation
More informationEigenvectors Via Graph Theory
Eigenvectors Via Graph Theory Jennifer Harris Advisor: Dr. David Garth October 3, 2009 Introduction There is no problem in all mathematics that cannot be solved by direct counting. -Ernst Mach The goal
More informationEquivalent relaxations of optimal power flow
Equivalent relaxations of optimal power flow 1 Subhonmesh Bose 1, Steven H. Low 2,1, Thanchanok Teeraratkul 1, Babak Hassibi 1 1 Electrical Engineering, 2 Computational and Mathematical Sciences California
More information