Suppose that we have a Hilbert space $\mathcal{H}$, an orthogonal projection $P$ (i.e. $P=|\Psi\rangle\langle\Psi|$ for some $\Psi\in\mathcal{H}$ with $\|\Psi\|=1$) and another non-negative bounded linear operator (not necessarily a projection!) $R$ satisfying $R\leq P$.
In my lecture notes, a lemma states that $R=RP=PR$, which is used to prove a theorem about joint measurability of POVMs in quantum mechanics. The proof is as follows:
Let $I$ denote the identity operator on $\mathcal{H}$, and $R\leq P$ implies $(I-P)R(I-P)=0$;
Because non-negative bounded operators have square roots due to the spectral theorem, we have $\sqrt{R}(I-P)=0$. Applying $\sqrt{R}$ to the left gives $R(I-P)=0$, hence $R=RP$ and similarly $R=PR$.
What I do not understand: why does $R\leq P$ imply $(I-P)R(I-P)=0$ ?
I consulted the following questions, but I am still stuck: