He explains this by citing the fact that the square of the wave function which gives the probability density is maximum at the origin.
Not exactly.
The answer by @dsva explains why this is wrong and I'll just expand that in a minute, but first note that it's easy to see why this is wrong.
If the electron was most likely to be at the origin it would be likely to interact with the nucleus, something we expect not to happen. In fact there is a small chance of this because the nucleus is a non-zero size.
Essentially he's missing the need to sum over the shell, and the shell's volume is zero at the origin because of the volume element $dV = 4\pi r^2dr$ with $r=0$.
The probability of finding the electron in a region of zero radial size is zero, but we can evaluate the relative probability for two radii. For 1-s orbitals that is :
$$\frac {r_1^2 e^{-{r_1}/a}}{r_1^2 e^{-{r_1}/a}+r_2^2 e^{-{r_2}/a}}$$
and
$$\frac {r_2^2 e^{-{r_2}/a}}{r_1^2 e^{-{r_1}/a}+r_2^2 e^{-{r_2}/a}}$$
For $r_1=0$ we clearly get a relative probability of zero compared to any $r_2\neq 0$.
At the same time, we all agree that the Bohr radius is the distance at which probability of finding the electron is maximum for 1s orbital.
This is because of what an expectation value is and how it is calculated in QM.
An expectation value is an average over the entire space of the effect of an operator on a wavefunction. It is not simply the operator multiplied by the probability density.
In short an expectation value can be considered an average over all space, whereas the probability density is not an average over all space.