Positive harmonic function on $\mathbb{R}^n$ is a constant?

Since it is not clear whether the Wikipedia proof uses boundedness or not, please allow me to give a detailed proof that only uses nonnegativity.

Let $u$ be a nonnegative harmonic function in $\mathbb{R}^n$, and let $x,y\in\mathbb{R}^n$. Denote by $B_R(y)$ the ball of radius $R>0$ centred at $y$, and similarly by $B_r(x)$ the ball of radius $r>0$ centred at $x$. We consider such balls with $R>0$ large and $r=R-|x-y|$, that is, $B_r(x)\subset B_R(y)$ and the boundary of $B_r(x)$ is tangent to the boundary of $B_R(y)$. Then we have $$ u(x) = \frac1{|B_r|}\int_{B_r(x)}u \leq \frac1{|B_r|}\int_{B_R(y)}u = \frac{|B_R|}{|B_r|}u(y), $$ where we have used the mean value property in the first and last equalities, and the nonnegativity of $u$ in the second inequality. Now we send $R\to\infty$ and get $$ u(x) \leq u(y). $$ As $x$ and $y$ were arbitrary, we conclude that $u$ is constant.

Note that the Harnack inequality is lurking behind this proof. Also, it was crucial to consider balls with different radii in order to make use of nonnegativity. The Wikipedia proof uses balls of equal radii.

Taking $n=2$ just to simplify:

Note that, given $y\in\mathbb{R}^2$, let $r=||y||$. Then, $\forall n\in\mathbb{N}$: $$B(y,n)\subset B(0,2n)$$ Thus $$\int_{B(y,n)}u(z)dz\leq\int_{B(0,n+r)}u(z)dz$$ $$u(y) = \frac{1}{\pi n^2}\int_{B(y,n)}u(z)dz\leq\frac{(n+r)^2}{n^2}\frac{1}{\pi(n+r)^2}\int_{B(0,n+r)}u(z)dz=\left(\frac{n+r}{n}\right)^2u(0)$$ That is $$u(y)\leq\left(\frac{n+r}{n}\right)^2u(0),\qquad\forall n \in\mathbb{N}$$ Takin $lim_n\rightarrow \infty$, we have that $u(y)\leq u(0)$. But since $y\in\mathbb{R}^2$ was arbitrary, we conclude that $u(z)\leq u(0),\forall z\in\mathbb{R}^2$.
But then, given $R>0$, we have that $\max_{z\in\overline{B(0,R)}}u(z) \leq u(0)$, thus, $\max_{z\in\overline{B(0,R)}}u(z) = u(0)$. Since $0$ is an interior point, by the strong maximum principle, $u$ is constant equal to $u(0)$ in the ball. Since $R$ was arbitrary, it follows that $u(z) = u(0),\forall z\in\mathbb{R}^2$.