Cosmological correlators encode statistical properties of the initial conditions of our universe. Mathematically, they can often be written as Mellin integrals of a certain rational function associated to graphs, namely the flat space wavefunction. The singularities of these cosmological integrals are parameterized by binary hyperplane arrangements. Using different algebraic tools, we shed light on the differential and difference equations satisfied by these integrals. Moreover, we study a multivariate version of partial fractioning of the flat space wavefunction, and propose a graph-based algorithm to compute this decomposition.