The reference that suggests that $|z|$ is differentiable at $z=0$ is incorrect. Note that if $f(z)$ is differentiable, then $$f'(z)=\lim_{\Delta z\to 0}\frac{f(z+\Delta z)-f(z)}{\Delta z}$$
Taking $f(z)=|z|$ and $z=0$, we see that $\displaystyle \lim_{\Delta z\to 0}\frac{|\Delta z|-0}{\Delta z}$ fails to exist.
Similarly, if $f(z)=u(x,y)+iv(x,y)$, and $g(z)=|f(z)|^2=u^2(x,y)+v^2(x,y)$, then $g(z)$ is purely real. The only purely real function that is complex differentiable in an open neighborhood of a point is a function that is constant. So, $g$ is differentiable in a neighborhood of $z$ only if $f$ is constant there.
To show this, we appeal to the Cauchy-Riemann equations. Note that if $h(z)=\phi(x,y)$, where $\phi$ is a purely real-valued function, then the Cauchy-Riemann equations reveal
$$\frac{\partial \phi(x,y)}{\partial x}=\frac{\partial \phi(x,y)}{\partial y}=0\implies \phi(x,y)\,\,\text{is a constant}$$
Now, it is possible that $|f|^2$ might be differentiable at an isolated point $z$, but not in an open neighborhood of that point. Note that if $f$ is analytic in the open domain $O$, and $z\in O$ and $z+\Delta z\in O$, then
$$\begin{align}
\frac{|f(z+\Delta z)|^2-|f(z)|^2}{\Delta z}&=\frac{|f(z)+f'(z)\Delta z+O\left((\Delta z)^2\right)|^2-|f(z)|^2}{\Delta z}\\\\
&=\frac{2\text{Re}\left(\overline{f(z)}f'(z)\Delta z\right)}{\Delta z}+O\left(\Delta z\right)
\end{align}$$
Hence, if the limit $\lim_{\Delta z\to 0}\frac{2\text{Re}\left(\overline{f(z)}f'(z)\Delta z\right)}{\Delta z}$ exists, then $|f|^2$ is differentiable at $z$. Again, $|f|^2$ cannot be analytic at any point unless $f$ is a constant.
Example 3.3.1. Consider the function f(z) = |z|^2 In our usual notation, we clearly have: u = x^2 + y^2 and v = 0. The Cauchy-Riemann equations 2x = 0 and 2y = 0 can only be satisfied at z = 0. It follows that the function is differentiable only at the point z = 0, and is therefore analytic nowhere. Is this incorrect?
– Hass Aug 25 '17 at 13:39