I have seen that the argument of why neutrino oscillations exist is that their mass matrix is nondiagonal and complex and so one needs to transform it into diagonal by unitary rotations:$$\nu_\alpha=\sum_i U_{\alpha i}\nu_i$$ $\alpha $ for flavor eigenstates and $i$ for mass eigenstates. At a later time the state of the neutrino $|\nu(t)>$ is given by a linear combination of the $\nu_i$ by their respective time evolution operator and so there is a probability to find then in different eigenstates at a later time by $|<\nu_\alpha|\nu(t)>|^2$.
Why not using the same argument to say that leptons or quarks change of flavor as they propagate? I don't know much of this, any wisdom on the subject would be appreciated.