A number theory paper I wrote was recently rejected from a journal due to "working under the untenable hypothesis that the natural density behaves like a probability measure (as it is not $\sigma$-additive and not stable under finite intersections)".
My paper was about providing a lattice-theoretic proof that the probability two Gaussian integers is relatively prime is $\frac{1}{\zeta_{\mathbb{Q}(i)}(2)}$, where $\zeta_{\mathbb{Q}(i)}(2)$ is the Dedekind zeta function over the Gaussian integers (as the known proof by Collins and Johnson from 1989 used difficult analytic number theory techniques). The main objection from the referee was to an argument of the following form:
Theorem: Let $x,y \in \mathbb{Z}$ be randomly chosen. Then, $$P(\gcd(x,y)=1) = \frac{1}{\zeta(2)}$$ Proof: We have that $P(\gcd(x,y)=1) = P(p \nmid x \cap p \nmid y, \forall p)$, where $p$ is prime. For a fixed prime $p$, the probability that $p \mid x$ and $p \mid y$ is $\frac{1}{p^2}$ and so we have that the probability this does not occur is $1 - \frac{1}{p^2}$. Therefore, taking the product over all primes we have that $$P(\gcd(x,y)=1) = \prod_{p \; \text{prime}} \left(1 - \frac{1}{p^2}\right) = \left(\prod_{p \; \text{prime}} \left(\frac{1}{1-p^{-2}}\right) \right)^{-1} = \frac{1}{\zeta(2)}$$
Why is this invalid reasoning?