Dipartimento di Matematica
|
|
- Chrystal Griffith
- 6 years ago
- Views:
Transcription
1 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, Torino-Italia
2 On the applicability of algebraic statistics to regression models with errors in variables Suela Ruffa Politecnico di Torino Abstract In this paper we study the applicability of algebraic statistics methods to regression linear models with errors in variables (EIV models). In fact the advantage in the algebraic statistics methods is that we can find saturated linear models identifiable with designs which are not regular fractions of full factorial designs. First we introduce the main concepts of algebraic statistics applied to the design of experiments (see [1]) and next we formulate the model of interest: the Berkson s model (see []) in the case of a polynomial linear regression problem. By mean of simple designs we illustrate the main issues met in the generalization of algebraic statistic theory to EIV models. Mainly these are due to the fact that the design points are unknown realizations of random variables. For the application we have in mind, the assumption of normality is natural. The generalization from these simple designs to higher dimensional and non regular fractional designs is straightforward. The paper concludes with a discussion on the applicability of currently available software for algebraic statistics analysis to designs with errors in variables. Keywords: EIV models, design of experiments, random design. 1 Algebraic statistics in design of experiments Algebraic statistics is a recent trend in mathematical and applied statistics for which concepts of algebraic commutative algebra and algebraic geometry are applied to the study of statistical models. The result of algebraic statistics of most practical relevance is in design of experiments. For any design D, that is a finite set of distinct points in IR m, it returns an IR-vector space basis of the set of real functions defined over D in a time polynomial in m and in n (see [1]), where n is the number of design points. That is it returns a saturated polynomial model which is identifiable by D. Computationally harder but still possible is to obtain the family of all saturated linear polynomial identifiable models with a hierarchical structure (see below). The bases considered have two main properties: 1. they consist of n monomials, where n is the number of points in D;. they are hierarchical, that is if the monomial x α = x α 1 1 x αm m is included in the basis and x β divides x α, then also x β = x β 1 1 x βm m is included in the basis. 1
3 The necessary computations can be performed with softwares for computational commutative algebra, like CoCoA with the function IdealOfPoints (available on the site dima. unige.it). This function works only for no replication designs even if the theory can been extended also for design with replicated points (see [4]). The theory on which this computation is based uses concepts of algebraic commutative algebra and in particular the concept of Gröbner basis (see any computational commutative algebra book as [5]). The following is technical but is the key point of all algebraic statistics theory. Over a field k (like R or C) we consider the polynomial ring k[x 1,..., x m ], the ring of m-dimensional polynomials with coefficients in k. Given a design D we define the design ideal, denoted by I =Ideal(D) as the set of all polynomials in k[x 1,..., x m ] whose zeros are design points. We say that a subset G = {g 1,..., g s } of the ideal I k[x 1,..., x m ] is a Gröbner basis if 1. g i I for all i = 1,..., s and. LT (g 1 ),..., LT (g s ) = LT (I) where LT (g i ) is the leading term (with respect to the term ordering that we fixed) of the polynomial g i and LT (I) is the ideal generated by the leading terms of the ideal I. For the notion of term ordering we refer to e.g. [5]. Here we simply note that a term ordering is a way to order monomials x α = x α 1 1 x αm m with α i Z 0, i = 1,..., n by ordering their exponents. A term ordering then corresponds to an order on the m-dimensional vectors with non-negative entries Z m 0. Given a design D we can prove that a monomial basis for the polynomials whose zeros include the design points is Est τ (D) = {x α : x α is not divisible by any LT of the elements of Gröbner basis of Ideal(D)} = {x α : x α LT (g) : g Ideal(D) } where τ is the chosen term ordering and Est τ (D) has n elements. It has been shown [1] that, given the design D the model β α x α (1) x α Est τ (D) is unambiguously identifiable with the design D. As consequence with a n point design we can identify n distinct monomials. Moreover Est τ (D) satisfies the two properties of monomial basis on the previous page. We can define also the fan of the design D as the set of all possible Est τ (D) when the term ordering τ varies over all term orderings. This is a finite set (see [10]). With an example we will clarify the steps of the process leading from a design to a saturated identifiable hierarchical model. Consider the full factorial design. Its design ideal is Ideal(D) = x 1 1, x 1 = {f 1 (x 1 1) + f (x 1) : f 1, f R[x 1, x ]}. The monomial basis with respect to any term ordering is Est(D) = {1, x 1, x, x 1 x } and the linear hierarchical saturated model identified with the design is then f(x 1, x ) = β 0 + β 1 x 1 + β x + β 3 x 1 x. Algebraic statistics methods work well if the points are fixed, but can we say something if the points are realizations of random variables? In this paper we answer this question.
4 EIV Berkson s model A regression model describes the relationship between a response variable and a set of variables called predictors or factors. Usually the response is modeled as a random variable while the predictors are considered as fixed constants. Very often the errors associated to factor settings are ignored. However there are situations where we cannot ignore that the levels of the factors are set with some error. In this case to describe the relationship between the response variables and the predictors we can use an errors in variables model (EIV model). In the case of simple linear regression (only one predictor, m = 1) the EIV model is: Y j = β 1 + β X j + ε j X j = t j + e j () where j runs over the n observations, e j and ε j are the errors associated to the variables X and Y respectively in the j-th observation and t j = E(X j ). If we assume that we can observe the value x j from the random variable X j whose mean is unknown, we can formulate a structural or a functional relationship. The models that arise from these two relationships are the most used in the literature (see e.g. [3, 8]). But in the application to design of experiments it seems to us more reasonable to assume that the t j with j = 1,..., n are known and that the x j are unknown. In fact they are the values in which we would like to observe our response. This situation is known as Berkson s model (see [, 7]). We cannot observe the true values that generate the data but in fact we suppose that the means of the random variables are known. If the planned value of the ith predictor in the jth observation is t j a value X j = t j + e j rather than the value t j generates the data point. The general form of a polynomial linear regression EIV Berkson s model is given in Equation (3). Y = β 1 f 1 (X) + + β p f p (X) + ε (3) X = t + e where X = (X 1,..., X m ) T is the column vector whose entries are predictors, f 1,..., f p are known polynomial functions of the predictors X 1,..., X m, ε is the error associated to the response variable Y, e = (e 1,..., e m ) T is the error associated to the factors (X 1,..., X m ) T, t = E(X) and β i are the regression coefficients. In the Equation (3) we omitted the subindex indicating the observation j to simplify the notations. In many cases (e.g. [7]) the assumptions on the errors are the following: e i N(0, σ i ); e i are stochastically independent, e.g. E(e i e j ) = 0 i j,..., E( m i=1 e i) = 0. Thus X i is a random variable with distribution Norm(t i, σ i ) for each i = 1,..., m and the vector X has a normal multivariate distribution Norm((t 1,..., t m ) T, Σ) with diagonal variance covariance matrix Σ =diag(σ 1,..., σ m). The previous independence assumptions on the e j mean that, from a practical viewpoint, the settings on one factor are independent of the rest of the factors. In the following sections we also assume that the random vector e i is independent of the random vector e j for any observation i j, i, j (1,..., n). We next use EIV Berkson s model to define a random design. 3
5 3 Random designs Under Berkson s model the factor settings are unknown realizations of random variables. Thus we define a random design as the set of these random variables. In the application of algebraic statistics to random designs a problem is that the environment where we have to work changes. Now we have random polynomials, that is polynomials whose coefficients are random variables thus we cannot work in a field as R. To apply algebraic statistics methods it is necessary that the set of polynomial coefficients is a field. So we have to prove that the normal independent random variables are in a field. Theorem 1 shows that this holds and thus we can perform the usual computations with polynomials but now their coefficients include also random variables. In fact we prove that the following theorem holds true. Theorem 1. Given the probability space (Ω, F, P), the set of random variables A = {z : Ω R such that P(z(ω) = 0) = 0 or P(z(ω) = 0) = 1} (4) with usual sum and product operations between functions in R is a field. Proof. The proof follows from arguments in undergraduate probability and algebra (see [13]). We want only to note that without the condition P(z(ω) = 0) = 0 or P(z(ω) = 0) = 1 the set A = {z : Ω R measurable} is not a field. In fact consider the two random variables { 1 ω B z 1 (ω) = 0 ω B { 0 ω B z (ω) = 1 ω B where B is a measurable set and B is its complementary set. These random variables are in A but z 1 (ω) z (ω) = 0 A with z 1 (ω) 0 A and z (ω) 0 A, this means that we obtain zero divisors and that A is not a field. We want to work with no replication designs (even if the theory can been extended also for design with replicated points (see [4])) because efficient software is available only in this case. If the design points are realizations of normal random variables we cannot exclude the possibility of replications but we can study the conditions under which the probability of such an event is small enough. In the next two subsections we find these conditions for random designs in one dimension and for a full factorial random design. Next we shall see how that can be extended to general designs. 3.1 One dimensional case µ 1 µ µ3 Figure 1: Means of one dimensional random design Consider the application of Berkson s model to a design with one factor and three levels. We have three independent random variables X i Norm(µ i, σ i ) for i = 1,, 3 4
6 and, without lack of generality, we can consider that the means are ordered as in Figure 1. Given ε > 0, 0 < α < 1, three random variables X i with i = 1,, 3 we say that the random variables composing a random design are sufficiently distant if P( X j X i < ε) < α i j. (5) If α is small then with high probability the realizations of these random variables are different from each others. In this case we can consider that with high probability the design has no replication. Once the distributions of X i are known one between α and ε determines the other one. With some standard calculation we can obtain conditions involving means and variances of random variables that assure the Equation (5) is satisfied. Under the normality assumption for X i Norm(µ i, σ i ) for i = 1,, 3, the conditions are the following: σ j + σ i < ε + µ j µ i for all j = i + 1, i = 1, (6) where is the αth percentile of a standard normal distribution. Thus ε reflects the variance of the factor settings and α is used to make sure that the design has no replication with high probability. We can extend the arguments to one dimensional designs with n points. In fact for such a design it is sufficient to verify that Equation (6) is satisfied for each set of tree adjacent points. Thus for a design with n points we obtain (n ) conditions of the type (6). 3. The full factorial random design We consider the design in Figure where the dots are the means of the random variables (µ 31, µ 3 ) (µ 41, µ 4 ) x (µ 1, µ ) x 1 (µ 11, µ 1 ) Figure : Means of a full factorial random design X i Norm((µ i1, µ i ), Σ i ) where Σ i is diag(σ i1, σ i), for i = 1,, 3, 4. As in the case of one dimensional designs we have to verify that the points of the random design are sufficiently distant. In this case we have to find conditions for which: P( X i X j < ε) < α i j where for X = (X 1, X ) T then X = X 1 + X. The following proposition gives sufficient conditions for that to happen. 5
7 Proposition. Let X = (X 1, X ) T Norm((µ X1,µ X ), Σ X ) and Y = (Y 1, Y ) T Norm((µ Y 1,µ Y ), Σ Y ) be two independent random vectors whose variance covariance matrices Σ X and Σ Y are diagonal. For ε > 0 if then P( X 1 Y 1 < ε ) < α or P( X Y < P( X Y < ε) < α. ε ) < α (7) To prove this proposition we use the independence assumption within and between random vectors. The last hypothesis in the case of independent random variables follows from the diagonal variance covariance matrix (see [13]). Now we go back to our example and the application of Proposition 9 implies the following conditions involving means and variances: σ 11 + σ 1 < ε + µ 11 µ 1 σ 1 + σ 31 < ε + µ 1 µ 31 σ 3 + σ 4 < ε + µ 3 µ 4 (8) σ 11 + σ 41 < ε + µ 11 µ 1 σ 1 + σ 41 < ε + µ 1 µ 41 σ 31 + σ 11 < ε + µ 31 µ 11 where σ i1 and σ i are the first and the second element of the diagonal matrix Σ i, is the α-th percentile of a standard normal distribution as in (6) and ε is a small positive number (see Proposition ). 3.3 Generalization to higher dimensional designs In general if we work in an m-dimensional space we can generalize the Proposition as follow. Proposition 3. Let X= (X 1,..., X m ) T Norm((µ X1,..., µ Xm ), Σ X ) and Y = (Y 1,..., Y m ) T Norm((µ Y 1,..., µ Y m ), Σ Y ) be two random independent random variables whose variance covariance matrices Σ X and Σ Y are diagonal. For ε > 0 if then P( X i Y i < ε m ) < α for at least one i (1,..., m) (9) P( X Y < ε) < α. Thus in the case of higher dimension designs we can find conditions similar to those in Equation (8). For any l m full factorial design with m factors and l levels each, the total number of conditions is ( l m ) if we impose that each point has to be sufficiently 6
8 distant from any other point (i.e. one condition for each pair of points). This number can be further reduced by considering only conditions relative to adjacent points, as in the one dimensional case. For example for m = we need only (l 1)(l 1) conditions, reducing the order of magnitude of the number of conditions in Proposition 3 from l m to l. In general the number of necessary conditions it easily seem to be smaller than m 1 ( m 1)(l 1) m. The conditions that guarantee the distance between the points for an l m full factorial design also guarantee this for any of its fractions. In fact any design is a fraction of a full factorial design with level sets of minimal size. 4 Applicability of algebraic statistics methods Now we return to the simple case of a random design. If this design satisfies the conditions given in Equation (8) then Theorem 4 allows us to apply algebraic statistics methods to the random design. Theoretically we can construct the design ideal as the set of polynomials whose zeros are random variables X 11,..., X 4 and with the method illustrated in [1] we can obtain the Gröbner bases of the design. In fact when we work with a random design we do not know precisely which are the factor settings. All we can say is that the fan of the random design is maximal with probability 1 (see [1]). For the full factorial random design, the fan is composed of five bases with four monomials. These are obtained varying the monomial ordering over all possible monomial orderings and are the following ones: {1, x 1, x, x 1 x } {1, x 1, x 1, x 3 1} {1, x 1, x, x 1} (10) {1, x 1, x, x } {1, x, x, x 3 }. Now we can define when with high probability the random design can be associated to a non random design whose points are the means of random variables and thus reduce our random design to a determinist design and make direct use of e.g. IdealOfPoints. For 0 < α < 1 we define ε α of a random design as a value for which P( X i µ i < ε α ) = 1 α. (11) Recall t i = µ i for i = 1,..., n. The value ε α can be seen as the precision of the instrument used in the measurement process. So Equation (11) can be interpreted in the following two ways: if we fix the probability level 1 α with α small, we can compute the value ε α such that Equation (11) is satisfied. The computed value can be compared with the precision of the measurement instrument η: if η < ε α then we can say with high (1 α) probability that the monomial basis for the random design is {1, x 1, x, x 1 x }. This is the monomial basis of the full factorial non-random design whose points are the means of the random variables. If η > ε α then we can use some classical statistical methods to evaluate which of the identifiable models in (10) is best associated with the design; 7
9 if we fix the value of ε α = η, then we compute the probability P( X i µ i < η). If this probability is sufficiently high we can affirm that with probability 1 α the monomial basis for the random design is {1, x 1, x, x 1 x }, and the other four models in (10) are not identifiable, that is there is collinearity among their elements. If the probability is low, we have again to use some other classical statistical methods to find which of the linear model in (10) better suits our design. Similar arguments to those in Section 3.3 generalize this to high dimensional design. Acknowledgements The author thanks Prof. Eva Riccomagno for her help during this work and Hugo Maruri- Aguilar for his useful suggestions. References [1] J. Abbott, A. Bigatti, M. Kreuzer, and L. Robbiano. Computing ideals of points. J. Symbolic Comput., 30: , 000. [] J. Berkson. Are there two regression? J. Am. Statist. Ass., 45: , [3] G. Casella and R. L. Berger. Statistical Inference. Duxbury, 00. [4] A. M. Cohen, A. Di Bucchianico, and E. Riccomagno. Replications with Gröbner bases. In Proceedings of the 6th International Workshop on Model-Oriented Data Analysis (Puchberg/Schneeberg, 001), pages Physica, Heidelberg, 001. [5] D. Cox, J. Little, and D. O Shea. Using Algebraic Geometry. Springer Verlag, [6] A. N. Donev. Dealing with errors in variables in response surface exploration. Commun. Statist.-Theory methods, 9: , 000. [7] A. N. Donev. Design of experiment in presence of errors in factor levels. Jour. Statist. Plan. Inf., 16: , 004. [8] W. A. Fuller. Measurement Error Models. Wiley, New York, [9] M. Kreuzer and L. Robbiano. Computational Commutative Algebra 1. Springer Verlag, 000. [10] T. Mora and L. Robbiano. The Gröbner fan of an ideal. J. Symbolic Comput., 6(-3):183 08, Computational aspects of commutative algebra. [11] T. Nummi and J. Möttönen. Estimation and prediction for low degree polynomial models under measurement errors with an application to forest harversters. Appl. Statist., 53:Part , 004. [1] G. Pistone, E. Riccomagno, and H. P. Wynn. Algebraic Statistics. Chapmann & hall/crc, 001. [13] S. Ruffa. La Statistica algebrica dei modelli di regressione con errori nelle variabili. Politecnico di Torino,
Dipartimento 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 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 informationMCS 563 Spring 2014 Analytic Symbolic Computation Friday 31 January. Quotient Rings
Quotient Rings In this note we consider again ideals, but here we do not start from polynomials, but from a finite set of points. The application in statistics and the pseudo code of the Buchberger-Möller
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 informationAlexander duality in experimental designs
Ann Inst Stat Math (2013) 65:667 686 DOI 10.1007/s10463-012-0390-9 Alexander duality in experimental designs Hugo Maruri-Aguilar Eduardo Sáenz-de-Cabezón Henry P. Wynn Received: 11 April 2012 / Revised:
More informationPOLYNOMIAL DIVISION AND GRÖBNER BASES. Samira Zeada
THE TEACHING OF MATHEMATICS 2013, Vol. XVI, 1, pp. 22 28 POLYNOMIAL DIVISION AND GRÖBNER BASES Samira Zeada Abstract. Division in the ring of multivariate polynomials is usually not a part of the standard
More informationGroebner Bases, Toric Ideals and Integer Programming: An Application to Economics. Tan Tran Junior Major-Economics& Mathematics
Groebner Bases, Toric Ideals and Integer Programming: An Application to Economics Tan Tran Junior Major-Economics& Mathematics History Groebner bases were developed by Buchberger in 1965, who later named
More informationS-Polynomials and Buchberger s Algorithm
S-Polynomials and Buchberger s Algorithm J.M. Selig Faculty of Business London South Bank University, London SE1 0AA, UK 1 S-Polynomials As we have seen in previous talks one of the problems we encounter
More informationConsistent Bivariate Distribution
A Characterization of the Normal Conditional Distributions MATSUNO 79 Therefore, the function ( ) = G( : a/(1 b2)) = N(0, a/(1 b2)) is a solu- tion for the integral equation (10). The constant times of
More informationGroebner Bases and Applications
Groebner Bases and Applications Robert Hines December 16, 2014 1 Groebner Bases In this section we define Groebner Bases and discuss some of their basic properties, following the exposition in chapter
More informationAlgebraic Representations of Gaussian Markov Combinations
Submitted to the Bernoulli Algebraic Representations of Gaussian Markov Combinations M. SOFIA MASSA 1 and EVA RICCOMAGNO 2 1 Department of Statistics, University of Oxford, 1 South Parks Road, Oxford,
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 information1 xa 2. 2 xan n. + c 2 x α 2
Operations Research Seminar: Gröbner Bases and Integer Programming Speaker: Adam Van Tuyl Introduction In this talk I will discuss how to use some of the tools of commutative algebra and algebraic geometry
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 informationPREMUR Seminar Week 2 Discussions - Polynomial Division, Gröbner Bases, First Applications
PREMUR 2007 - Seminar Week 2 Discussions - Polynomial Division, Gröbner Bases, First Applications Day 1: Monomial Orders In class today, we introduced the definition of a monomial order in the polyomial
More informationOil Fields and Hilbert Schemes
Oil Fields and Hilbert Schemes Lorenzo Robbiano Università di Genova Dipartimento di Matematica Lorenzo Robbiano (Università di Genova) Oil Fields and Hilbert Schemes March, 2008 1 / 35 Facts In the realm
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 informationNEW CLASSES OF SET-THEORETIC COMPLETE INTERSECTION MONOMIAL IDEALS
NEW CLASSES OF SET-THEORETIC COMPLETE INTERSECTION MONOMIAL IDEALS M. R. POURNAKI, S. A. SEYED FAKHARI, AND S. YASSEMI Abstract. Let be a simplicial complex and χ be an s-coloring of. Biermann and Van
More informationGRÖBNER BASES AND POLYNOMIAL EQUATIONS. 1. Introduction and preliminaries on Gróbner bases
GRÖBNER BASES AND POLYNOMIAL EQUATIONS J. K. VERMA 1. Introduction and preliminaries on Gróbner bases Let S = k[x 1, x 2,..., x n ] denote a polynomial ring over a field k where x 1, x 2,..., x n are indeterminates.
More informationPolynomials, Ideals, and Gröbner Bases
Polynomials, Ideals, and Gröbner Bases Notes by Bernd Sturmfels for the lecture on April 10, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra We fix a field K. Some examples of fields
More informationDESIGN OF EXPERIMENTS AND BIOCHEMICAL NETWORK INFERENCE
DESIGN OF EXPERIMENTS AND BIOCHEMICAL NETWORK INFERENCE REINHARD LAUBENBACHER AND BRANDILYN STIGLER Abstract. Design of experiments is a branch of statistics that aims to identify efficient procedures
More informationAlgebraic Statistics
Algebraic Statistics Part 1: Eva Riccomagno Department of Mathematics, Università di Genova riccomagno@dima.unige.it Part 2: Henry P. Wynn Department of Statistics, London School of Economics H.Wynn@lse.ac.uk...
More informationLecture 2: Gröbner Basis and SAGBI Basis
Lecture 2: Gröbner Basis and SAGBI Basis Mohammed Tessema Suppose we have a graph. Suppose we color the graph s vertices with 3 colors so that if the vertices are adjacent they are not the same colors.
More informationNotes on Linear Algebra and Matrix Theory
Massimo Franceschet featuring Enrico Bozzo Scalar product The scalar product (a.k.a. dot product or inner product) of two real vectors x = (x 1,..., x n ) and y = (y 1,..., y n ) is not a vector but a
More informationComputing Minimal Polynomial of Matrices over Algebraic Extension Fields
Bull. Math. Soc. Sci. Math. Roumanie Tome 56(104) No. 2, 2013, 217 228 Computing Minimal Polynomial of Matrices over Algebraic Extension Fields by Amir Hashemi and Benyamin M.-Alizadeh Abstract In this
More informationExtremal Behaviour in Sectional Matrices
Extremal Behaviour in Sectional Matrices Elisa Palezzato 1 joint work with Anna Maria Bigatti 1 and Michele Torielli 2 1 University of Genova, Italy 2 Hokkaido University, Japan arxiv:1702.03292 Ph.D.
More informationBinomial Ideals from Graphs
Binomial Ideals from Graphs Danielle Farrar University of Washington Yolanda Manzano St. Mary s University Juan Manuel Torres-Acevedo University of Puerto Rico Humacao August 10, 2000 Abstract This paper
More informationSHUHONG GAO, VIRGÍNIA M. RODRIGUES, AND JEFFREY STROOMER
GRÖBNER BASIS STRUCTURE OF FINITE SETS OF POINTS SHUHONG GAO, VIRGÍNIA M. RODRIGUES, AND JEFFREY STROOMER Abstract. We study the relationship between certain Gröbner bases for zerodimensional radical ideals,
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 informationLIMITING CASES OF BOARDMAN S FIVE HALVES THEOREM
Proceedings of the Edinburgh Mathematical Society Submitted Paper Paper 14 June 2011 LIMITING CASES OF BOARDMAN S FIVE HALVES THEOREM MICHAEL C. CRABB AND PEDRO L. Q. PERGHER Institute of Mathematics,
More informationCHAPTER 0 PRELIMINARY MATERIAL. Paul Vojta. University of California, Berkeley. 18 February 1998
CHAPTER 0 PRELIMINARY MATERIAL Paul Vojta University of California, Berkeley 18 February 1998 This chapter gives some preliminary material on number theory and algebraic geometry. Section 1 gives basic
More informationAlgebra 2. Rings and fields. Rings. A.M. Cohen, H. Cuypers, H. Sterk
and fields 2 and fields A.M. Cohen, H. Cuypers, H. Sterk A.M. Cohen, H. Cuypers, H. Sterk 2 September 25, 2006 1 / 31 and fields Multiplication turns each of the sets Z, Q, R, C, Z[X ], Q[X ], R[X ], C[X
More informationJournal of Algebra 226, (2000) doi: /jabr , available online at on. Artin Level Modules.
Journal of Algebra 226, 361 374 (2000) doi:10.1006/jabr.1999.8185, available online at http://www.idealibrary.com on Artin Level Modules Mats Boij Department of Mathematics, KTH, S 100 44 Stockholm, Sweden
More informationOn the minimal free resolution of a monomial ideal.
On the minimal free resolution of a monomial ideal. Caitlin M c Auley August 2012 Abstract Given a monomial ideal I in the polynomial ring S = k[x 1,..., x n ] over a field k, we construct a minimal free
More informationJournal of Symbolic Computation
Journal of Symbolic Computation 50 (013) 55 Contents lists available at SciVerse ScienceDirect Journal of Symbolic Computation www.elsevier.com/locate/jsc On the evaluation of multivariate polynomials
More informationarxiv: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 informationTHE CONE OF BETTI TABLES OVER THREE NON-COLLINEAR POINTS IN THE PLANE
JOURNAL OF COMMUTATIVE ALGEBRA Volume 8, Number 4, Winter 2016 THE CONE OF BETTI TABLES OVER THREE NON-COLLINEAR POINTS IN THE PLANE IULIA GHEORGHITA AND STEVEN V SAM ABSTRACT. We describe the cone of
More informationA decoding algorithm for binary linear codes using Groebner bases
A decoding algorithm for binary linear codes using Groebner bases arxiv:1810.04536v1 [cs.it] 9 Oct 2018 Harinaivo ANDRIATAHINY (1) e-mail : hariandriatahiny@gmail.com Jean Jacques Ferdinand RANDRIAMIARAMPANAHY
More informationSYMMETRIC POLYNOMIALS
SYMMETRIC POLYNOMIALS KEITH CONRAD Let F be a field. A polynomial f(x 1,..., X n ) F [X 1,..., X n ] is called symmetric if it is unchanged by any permutation of its variables: for every permutation σ
More informationMINIMAL POLYNOMIALS AND CHARACTERISTIC POLYNOMIALS OVER RINGS
JP Journal of Algebra, Number Theory and Applications Volume 0, Number 1, 011, Pages 49-60 Published Online: March, 011 This paper is available online at http://pphmj.com/journals/jpanta.htm 011 Pushpa
More information1 Mixed effect models and longitudinal data analysis
1 Mixed effect models and longitudinal data analysis Mixed effects models provide a flexible approach to any situation where data have a grouping structure which introduces some kind of correlation between
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 informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,500 08,000.7 M Open access books available International authors and editors Downloads Our authors
More informationOn The Belonging Of A Perturbed Vector To A Subspace From A Numerical View Point
Applied Mathematics E-Notes, 7(007), 65-70 c ISSN 1607-510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ On The Belonging Of A Perturbed Vector To A Subspace From A Numerical View
More informationLecture 1. (i,j) N 2 kx i y j, and this makes k[x, y]
Lecture 1 1. Polynomial Rings, Gröbner Bases Definition 1.1. Let R be a ring, G an abelian semigroup, and R = i G R i a direct sum decomposition of abelian groups. R is graded (G-graded) if R i R j R i+j
More informationAlgebraic Generation of Orthogonal Fractional Factorial Designs
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
More informationINITIAL COMPLEX ASSOCIATED TO A JET SCHEME OF A DETERMINANTAL VARIETY. the affine space of dimension k over F. By a variety in A k F
INITIAL COMPLEX ASSOCIATED TO A JET SCHEME OF A DETERMINANTAL VARIETY BOYAN JONOV Abstract. We show in this paper that the principal component of the first order jet scheme over the classical determinantal
More informationComputing syzygies with Gröbner bases
Computing syzygies with Gröbner bases Steven V Sam July 2, 2008 1 Motivation. The aim of this article is to motivate the inclusion of Gröbner bases in algebraic geometry via the computation of syzygies.
More informationLecture 20: Linear model, the LSE, and UMVUE
Lecture 20: Linear model, the LSE, and UMVUE Linear Models One of the most useful statistical models is X i = β τ Z i + ε i, i = 1,...,n, where X i is the ith observation and is often called the ith response;
More information4 Hilbert s Basis Theorem and Gröbner basis
4 Hilbert s Basis Theorem and Gröbner basis We define Gröbner bases of ideals in multivariate polynomial rings and see how they work in tandem with the division algorithm. We look again at the standard
More informationLecture 7: Schwartz-Zippel Lemma, Perfect Matching. 1.1 Polynomial Identity Testing and Schwartz-Zippel Lemma
CSE 521: Design and Analysis of Algorithms I Winter 2017 Lecture 7: Schwartz-Zippel Lemma, Perfect Matching Lecturer: Shayan Oveis Gharan 01/30/2017 Scribe: Philip Cho Disclaimer: These notes have not
More informationThe Geometry of Statistical Models for Two-Way Contingency Tables with Fixed Odds Ratios
Rend. Istit. Mat. Univ. Trieste Vol. XXXVII, 7 84 (2005) The Geometry of Statistical Models for Two-Way Contingency Tables with Fixed Odds Ratios Enrico Carlini and Fabio Rapallo ( ) Contribution to School
More informationOn the Computation of the Adjoint Ideal of Curves with Ordinary Singularities
Applied Mathematical Sciences Vol. 8, 2014, no. 136, 6805-6812 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.49697 On the Computation of the Adjoint Ideal of Curves with Ordinary Singularities
More information1. Introduction to commutative rings and fields
1. Introduction to commutative rings and fields Very informally speaking, a commutative ring is a set in which we can add, subtract and multiply elements so that the usual laws hold. A field is a commutative
More informationGENERATING IDEALS IN SUBRINGS OF K[[X]] VIA NUMERICAL SEMIGROUPS
GENERATING IDEALS IN SUBRINGS OF K[[X]] VIA NUMERICAL SEMIGROUPS SCOTT T. CHAPMAN Abstract. Let K be a field and S be the numerical semigroup generated by the positive integers n 1,..., n k. We discuss
More informationProblem Set 1 Solutions
Math 918 The Power of Monomial Ideals Problem Set 1 Solutions Due: Tuesday, February 16 (1) Let S = k[x 1,..., x n ] where k is a field. Fix a monomial order > σ on Z n 0. (a) Show that multideg(fg) =
More information1. Introduction to commutative rings and fields
1. Introduction to commutative rings and fields Very informally speaking, a commutative ring is a set in which we can add, subtract and multiply elements so that the usual laws hold. A field is a commutative
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 information1 Commutative Rings with Identity
1 Commutative Rings with Identity The first-year courses in (Abstract) Algebra concentrated on Groups: algebraic structures where there is basically one algebraic operation multiplication with the associated
More informationComputation of the Minimal Associated Primes
Computation of the Minimal Associated Primes Santiago Laplagne Departamento de Matemática, Universidad de Buenos Aires Buenos Aires, Argentina slaplagn@dm.uba.ar Abstract. We propose a new algorithm for
More informationLebesgue-Stieltjes measures and the play operator
Lebesgue-Stieltjes measures and the play operator Vincenzo Recupero Politecnico di Torino, Dipartimento di Matematica Corso Duca degli Abruzzi, 24, 10129 Torino - Italy E-mail: vincenzo.recupero@polito.it
More informationKey words. n-d systems, free directions, restriction to 1-D subspace, intersection ideal.
ALGEBRAIC CHARACTERIZATION OF FREE DIRECTIONS OF SCALAR n-d AUTONOMOUS SYSTEMS DEBASATTAM PAL AND HARISH K PILLAI Abstract In this paper, restriction of scalar n-d systems to 1-D subspaces has been considered
More informationAlgebraic Varieties. Notes by Mateusz Micha lek for the lecture on April 17, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra
Algebraic Varieties Notes by Mateusz Micha lek for the lecture on April 17, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra Algebraic varieties represent solutions of a system of polynomial
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationInferring the Causal Decomposition under the Presence of Deterministic Relations.
Inferring the Causal Decomposition under the Presence of Deterministic Relations. Jan Lemeire 1,2, Stijn Meganck 1,2, Francesco Cartella 1, Tingting Liu 1 and Alexander Statnikov 3 1-ETRO Department, Vrije
More informationGeneral error locator polynomials for binary cyclic codes with t 2 and n < 63
General error locator polynomials for binary cyclic codes with t 2 and n < 63 April 22, 2005 Teo Mora (theomora@disi.unige.it) Department of Mathematics, University of Genoa, Italy. Emmanuela Orsini (orsini@posso.dm.unipi.it)
More informationInteresting Examples on Maximal Irreducible Goppa Codes
Interesting Examples on Maximal Irreducible Goppa Codes Marta Giorgetti Dipartimento di Fisica e Matematica, Universita dell Insubria Abstract. In this paper a full categorization of irreducible classical
More informationIntegral Extensions. Chapter Integral Elements Definitions and Comments Lemma
Chapter 2 Integral Extensions 2.1 Integral Elements 2.1.1 Definitions and Comments Let R be a subring of the ring S, and let α S. We say that α is integral over R if α isarootofamonic polynomial with coefficients
More informationPrimary Decomposition of Ideals Arising from Hankel Matrices
Primary Decomposition of Ideals Arising from Hankel Matrices Paul Brodhead University of Wisconsin-Madison Malarie Cummings Hampton University Cora Seidler University of Texas-El Paso August 10 2000 Abstract
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 informationCombining Linear Equation Models via Dempster s Rule
Combining Linear Equation Models via Dempster s Rule Liping Liu Abstract. This paper proposes a concept of imaginary extreme numbers, which are like traditional complex number a + bi but with i = 1 being
More informationWeighted Least Squares
Weighted Least Squares The standard linear model assumes that Var(ε i ) = σ 2 for i = 1,..., n. As we have seen, however, there are instances where Var(Y X = x i ) = Var(ε i ) = σ2 w i. Here w 1,..., w
More informationRational Normal Curves as Set-Theoretic Complete Intersections of Quadrics
Rational Normal Curves as Set-Theoretic Complete Intersections of Quadrics Maria-Laura Torrente Dipartimento di Matematica, Università di Genova, Via Dodecaneso 35, I-16146 Genova, Italy torrente@dimaunigeit
More informationMIT Algebraic techniques and semidefinite optimization February 16, Lecture 4
MIT 6.972 Algebraic techniques and semidefinite optimization February 16, 2006 Lecture 4 Lecturer: Pablo A. Parrilo Scribe: Pablo A. Parrilo In this lecture we will review some basic elements of abstract
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 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 informationGeometry of Phylogenetic Inference
Geometry of Phylogenetic Inference Matilde Marcolli CS101: Mathematical and Computational Linguistics Winter 2015 References N. Eriksson, K. Ranestad, B. Sturmfels, S. Sullivant, Phylogenetic algebraic
More informationA Generalization of Wilson s Theorem
A Generalization of Wilson s Theorem R. Andrew Ohana June 3, 2009 Contents 1 Introduction 2 2 Background Algebra 2 2.1 Groups................................. 2 2.2 Rings.................................
More informationChange of Ordering for Regular Chains in Positive Dimension
Change of Ordering for Regular Chains in Positive Dimension X. Dahan, X. Jin, M. Moreno Maza, É. Schost University of Western Ontario, London, Ontario, Canada. École polytechnique, 91128 Palaiseau, France.
More informationControl of Hybrid Petri Nets using Max-Plus Algebra
Control of Hybrid Petri Nets using Max-Plus Algebra FABIO BALDUZZI*, ANGELA DI FEBBRARO*, ALESSANDRO GIUA, SIMONA SACONE^ *Dipartimento di Automatica e Informatica Politecnico di Torino Corso Duca degli
More informationD-MATH Algebra I HS 2013 Prof. Brent Doran. Exercise 11. Rings: definitions, units, zero divisors, polynomial rings
D-MATH Algebra I HS 2013 Prof. Brent Doran Exercise 11 Rings: definitions, units, zero divisors, polynomial rings 1. Show that the matrices M(n n, C) form a noncommutative ring. What are the units of M(n
More informationarxiv: v1 [math.lo] 12 Dec 2017
MODAL OPERATORS AND TORIC IDEALS RICCARDO CAMERLO, GIOVANNI PISTONE, AND FABIO RAPALLO arxiv:1712.04191v1 [math.lo] 12 Dec 2017 Abstract. In the present paper we consider modal propositional logic and
More informationPrinciples of Real Analysis I Fall I. The Real Number System
21-355 Principles of Real Analysis I Fall 2004 I. The Real Number System The main goal of this course is to develop the theory of real-valued functions of one real variable in a systematic and rigorous
More informationStandard Bases for Linear Codes over Prime Fields
Standard Bases for Linear Codes over Prime Fields arxiv:1708.05490v1 cs.it] 18 Aug 2017 Jean Jacques Ferdinand RANDRIAMIARAMPANAHY 1 e-mail : randriamiferdinand@gmail.com Harinaivo ANDRIATAHINY 2 e-mail
More informationSmooth supersaturated models
Smooth supersaturated models Ron A. Bates a, Hugo Maruri-Aguilar a,, Henry P. Wynn a a Department of Statistics, London School of Economics, London WC2A 2AE, UK Abstract In areas such as kernel smoothing
More informationON CERTAIN CLASSES OF CURVE SINGULARITIES WITH REDUCED TANGENT CONE
ON CERTAIN CLASSES OF CURVE SINGULARITIES WITH REDUCED TANGENT CONE Alessandro De Paris Università degli studi di Napoli Federico II Dipartimento di Matematica e Applicazioni R. Caccioppoli Complesso Monte
More informationCOS 424: Interacting with Data
COS 424: Interacting with Data Lecturer: Rob Schapire Lecture #14 Scribe: Zia Khan April 3, 2007 Recall from previous lecture that in regression we are trying to predict a real value given our data. Specically,
More informationAbstract Algebra for Polynomial Operations. Maya Mohsin Ahmed
Abstract Algebra for Polynomial Operations Maya Mohsin Ahmed c Maya Mohsin Ahmed 2009 ALL RIGHTS RESERVED To my students As we express our gratitude, we must never forget that the highest appreciation
More informationThe Multivariate Gaussian Distribution [DRAFT]
The Multivariate Gaussian Distribution DRAFT David S. Rosenberg Abstract This is a collection of a few key and standard results about multivariate Gaussian distributions. I have not included many proofs,
More informationSplitting sets and weakly Matlis domains
Commutative Algebra and Applications, 1 8 de Gruyter 2009 Splitting sets and weakly Matlis domains D. D. Anderson and Muhammad Zafrullah Abstract. An integral domain D is weakly Matlis if the intersection
More informationComputing Free Resolutions in Macaulay2
Computing Free Resolutions in Macaulay2 Madeline Brandt October 6, 2015 Introduction We will let R = k[x 1,..., x r ]. Definition 1. A free resolution of an R-module M is a complex F : F n φ n φ 1 F1 F0
More informationTHE FUNDAMENTAL THEOREM ON SYMMETRIC POLYNOMIALS. Hamza Elhadi S. Daoub
THE TEACHING OF MATHEMATICS 2012, Vol. XV, 1, pp. 55 59 THE FUNDAMENTAL THEOREM ON SYMMETRIC POLYNOMIALS Hamza Elhadi S. Daoub Abstract. In this work we are going to extend the proof of Newton s theorem
More informationTropical Islands. Jan Verschelde
Tropical Islands Jan Verschelde University of Illinois at Chicago Department of Mathematics, Statistics, and Computer Science http://www.math.uic.edu/ jan jan@math.uic.edu Graduate Computational Algebraic
More informationThe generic Gröbner walk
The generic Gröbner walk K. Fukuda, A. N. Jensen, N. Lauritzen, R. Thomas arxiv:math/0501345v1 [math.ac] 21 Jan 2005 February 28, 2008 Abstract The Gröbner walk is an algorithm for conversion between Gröbner
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 informationConstructive algebra. Thierry Coquand. May 2018
Constructive algebra Thierry Coquand May 2018 Constructive algebra is algebra done in the context of intuitionistic logic 1 H. Lombardi, C. Quitté Commutative Algebra: Constructive Methods, 2016 I. Yengui
More informationPrimary Decompositions of Powers of Ideals. Irena Swanson
Primary Decompositions of Powers of Ideals Irena Swanson 2 Department of Mathematics, University of Michigan Ann Arbor, Michigan 48109-1003, USA iswanson@math.lsa.umich.edu 1 Abstract Let R be a Noetherian
More informationMCS 563 Spring 2014 Analytic Symbolic Computation Monday 27 January. Gröbner bases
Gröbner bases In this lecture we introduce Buchberger s algorithm to compute a Gröbner basis for an ideal, following [2]. We sketch an application in filter design. Showing the termination of Buchberger
More informationOn the stable set of associated prime ideals of monomial ideals and square-free monomial ideals
On the stable set of associated prime ideals of monomial ideals and square-free monomial ideals Kazem Khashyarmanesh and Mehrdad Nasernejad The 10th Seminar on Commutative Algebra and Related Topics, 18-19
More information