( ) (, ) arxiv: 1510.02269 hgm OpenXM search. d n A = (a ij ). A i a i Z d, Z d. i a ij > 0. β N 0 A = N 0 a 1 + + N 0 a n Z A (β; p) = Au=β,u N n 0 A-. u! = n i=1 u i!, p u = n i=1 pu i i. Z = Z A Au = β ( ). p u u! (1)
Z A (β; p) = 22(2 2 ): A = : û = u = ( 11 0 25 12 1 p25 3 p12 ( u1 u 2 u 3 u 4 Au=β,u N d 0 0 0 1 1 1 0 1 0 0 1 0 1 ) ( 4 7,..., u = 32 5 p u u!, β = (37, 36, 12) T. u ). Au = β u ) ( 0 11,..., u = 36 1 Z A (β; p) = p11 4 11!25!12! 2 F 1 ( 11, 12, 26; y), y = p 2p 3 p 1 p 4 (odds ), ( ) 1 1 Ker A = Zā, ā =, u = û + vā y v /(û+vā)! 1 1 2F 1 (y)/û!. ).
p(ξ) = (exp ξ 1,..., exp ξ n ), ψ(ξ) = log Z(β; p(ξ)). i = p i. E[U i ] = p i i Z Z p=p(ξ) = ψ(ξ) ξ i, ψ(ξ) ξ i = Au=β u i exp(u ξ)/u! Z ( u ) E[U i ] u i. MLE (Maximal likelihood estimation) : ξ E[U] (Au = β u MLE ). 1 2 MLE? E[U]? 3. 4?
( :,, 2014, ) η i = E[U i ], η = (η i ). ξ- η- ( ) moment map E[U]. : ψ(ξ) (strictly ). ξ ψ(ξ) =: ξ η ψ (ξ ) = maxarg ξ (ξ ξ ψ(ξ)) ψ (ξ ) = ξ. ξ- ξ (η ), dual flat line. D[ξ : ξ ] = ψ(ξ) ψ(ξ ) ψ(ξ )(ξ ξ ) (divergence), D[ξ : ξ ] = D [ξ, ξ ]. : ψ(ξ) = log Z(β; p(ξ)).
U = (U 1,..., U n ), Fisher (Fisher s maximal likelihood estimation, MLE), U ( ) pu U! /Z(β, p) p p. Proposition ( ) p Fisher s MLE p E i (p) = p i Z p i /Z = U i.. pu U! /Z(β, p) Ui log p i log U! log Z F. p, F p i = U i p i Z/ p i Z = 0., U i = p i Z/ p i Z.
: M.Michalek, B.Sturmfels, C.Uhler, P.Zwiernik, Exponential Varieties, arxiv:1412.6185 : ξ positive definite m m symmetric matrix. x m-, Z(ξ) = exp R m ( 1 2 x T ξx log Z = m 2 log(2π) 1 2 log det(ξ), log Z = 1 2 ξ 1 ) dx = (2π)m/2 det(ξ) 1/2. Positive defifinite symmmetric matrix PD m convex cone. PM m convex cone. ξ ij log Z = 1 2 ξ (maximal likelihood estimation). log R n x i x j exp X ( 1 ) 2 xt ξx dx = 1 2 E[X i X j ] ( exp 1 ) 2 X T ξx /Z
( MLE) 4 7 32 5 (4, 7, 32, 5) β fix ψ., p 2(ξ)p 3 (ξ) p 1 (ξ)p 4 (ξ), 10.4167. MLE odds. ( ) ξ, MLE odds.. 1,,, :, 39. 71-100 (2010). 2, Algebraic Statistical Methods for Conditional Inference of Discrete Statistical Models, 2015/3, ( ).
R (Z n unimodular ), S (Z d unimodular ) α 1 O SAR =... O, α i 0, α i α i+1. O α d Ker(A : Z n Z d ) Z- {Re d+1,..., Re n }, (Re d+i ) T ā i, Ā = λ = Āξ. exp(λ i) = pāi. 22: Ā = ( 1, 1, 1, 1), p 2p 3 p 1 p 4 ā 1. ā n d (n d) n = n j=1 pāij j, odds ( odds = odds )
MLE? E[U]? ( ). ξ ξ Im A T. E[U](ξ) = E[U](ξ ). Theorem Newton New(Z) ( Z ξ p ) n d., ( odds log ) moment map. E[U] : R n /Im A T relint(new(z)) R n relint., f (ξ) = η ξ log Z(β, p(ξ)). 22: 25 < E[U 21 ] < 36. (1) MLE..
moment map E. 2 3. 4 7 2 32 5 6, generalized odds 2. Z Newton polytope 2.,. p R 2.
moment map E[U](p)..
log Z A.. 1 M.Saito, B.Sturmfels, N.Takayama, Hypergeometric polynomials and Integer Programming, Compositio Mathematica, 115 (1999), 185 204... 2 Holonomic Gradient Method (HGM) (2011 ). Z. HGM= ( ) Z. 3 M.Ogawa, A.Takemura, N.Takayama, An Application of A-hypergeometric Equations to Conditional Maximal Likelihood Estimation of 2 m Contingency Tables, in preparation. D (2015/3 ). 4 Y.Goto, Contiguity Relations of Lauricella s F D Revisited, arxiv:1412.3256. F D 2 m. 5 K.Ohara, N.Takayama, Pfaffian Systems of A-Hypergeometric Systems II Holonomic Gradient Method, arxiv:1505.02947. A-. Macaulay. E[U] before HGM.
E[U](ξ)? (k β ). P k (u, ξ) = exp(u ξ) u!z k (ξ), u S k, k = 1, 2,..., S k = {u N n Au = kβ}, Z k (ξ) = exp(u ξ), u! u S k : k? 1 J.Cornfield, A Statistical Problem Arising from Retrospective Studies, Proceedings of 3rd Berkeley Symposium on Mathematical Statistics and Probability 4 (1956), 135 148. 2,,, 1982. 3 R.L.Plackett, Analysis of Categorical Data, 2nd ed, Griffin, 1981. (pp. 41 (2 2 table), pp. 65 66 (r s table))
: k pu /u! Z(kβ;p)? λ = Āξ, exp(ξ) = p. m = m(λ) ( ) ( IPS (iterative proportional scaling) ). { β = Am, λ = Ā log m. (2) : km. 22, p 2p 3 p 1 p 4 = 1/3 3. 0.0 0.1 0.2 0.3 0.00 0.05 0.10 0.15 0.20 0.25 0 2 4 6 8 10 h2 0 2 4 6 8 10 h3
Theorem β N 0 A int(r 0 A). M = diag(m i ). ( ) P k (u, ξ) = det(ām 1 Ā T ) 1/2 n (u i km i ) 2 (2πk) (n d)/2 exp 2km i, sup i u i km i <φ(k) i=1 P k (u, ξ) P k (u, ξ) 1 0 (k ), φ(k) : φ(k) = o(k 2/3 ), k/φ(k) 2 = o(1).
MLE 1 p = u/ u ( u = u 1 + + u n ). (m odds λ. 2 E[U](p) m.) 2 ξ dual flat line p. (Ė(U)(ξ) ): 1 Y.Goto, K.Matsumoto, Pfaffian equations and contiguity relations of the hypergeometric function of type (k + 1, k + n + 2) and their applications, arxiv:1602.01637. - r 1 r 2 contingency tables HGM (Aomoto-Gel fand E(k, n) twisted cohomology group ). 2 modular method, - -.
. A = 0 0 0 1 1 1 1 1 0 0 1 0 1 0 0 1 1 0 1 0 1 1 1 0 1 1 0 0 2 2 2. p 1 p 2 p 4 p 5, 0 p 3 p 6 p 7. p 4 + p 5 + p 6 + p 7, p 1 + p 4 + p 6, p 2 + p 3 + p 5 + p 7, p 1 + p 2 + p 4 + p 5. cell p = exp(ξ)? (MLE, ). 19 132 0 9 11 52 6 97 p 1, p 2, p 3, p 4. (p 5, p 6, p 7 ). η = (E[U 5 ], E[U 6 ], E[U 7 ]) p η P 0 = (19, 132, 9, 11, 52, 6, 97)/326 (51.9194, 5.99193, 97.0891) P 0 + (0, 0, 0, 0, h 1, h 2, h 3 ) h = (0.000256154, 0.000152585, 0.00310983) (52.0006, 6.00006, 96.9993)