English

Maximum entropy Edgeworth estimates of the number of integer points in polytopes

Methodology 2010-08-10 v2 Computation

Abstract

Abstract: The number of points x=(x1,x2,...xn)x=(x_1 ,x_2 ,...x_n) that lie in an integer cube CC in RnR^n and satisfy the constraints jhij(xj)=si,1id\sum_j h_{ij}(x_j )=s_i ,1\le i\le d is approximated by an Edgeworth-corrected Gaussian formula based on the maximum entropy density pp on xCx \in C, that satisfies Ejhij(xj)=si,1idE\sum_j h_{ij}(x_j )=s_i ,1\le i\le d. Under pp, the variables X1,X2,...XnX_1 ,X_2 ,...X_n are independent with densities of exponential form. Letting SiS_i denote the random variable jhij(Xj)\sum_j h_{ij}(X_j ), conditional on S=s,XS=s, X is uniformly distributed over the integers in CC that satisfy S=sS=s. The number of points in CC satisfying S=sS=s is p{S=s}exp(I(p))p \{S=s\}\exp (I(p)) where I(p)I(p) is the entropy of the density pp. We estimate p{S=s}p \{S=s\} by pZ(s)p_Z(s), the density at ss of the multivariate Gaussian ZZ with the same first two moments as SS; and when dd is large we use in addition an Edgeworth factor that requires the first four moments of SS under pp. The asymptotic validity of the Edgeworth-corrected estimate is proved and demonstrated for counting contingency tables with given row and column sums as the number of rows and columns approaches infinity, and demonstrated for counting the number of graphs with a given degree sequence, as the number of vertices approaches infinity.

Keywords

Cite

@article{arxiv.0910.2497,
  title  = {Maximum entropy Edgeworth estimates of the number of integer points in polytopes},
  author = {Alexander Barvinok and J. A. Hartigan},
  journal= {arXiv preprint arXiv:0910.2497},
  year   = {2010}
}

Comments

29 pages 3 tables Revision updates references,and sharpens statement and proof of theorem 2