Some interpretations, like the many-worlds interpretation, treat the wavefunction (modulo an overall phase factor) as objective and fundamental.
But consider the following example for a qubit: a classical probability distribution over wavefunctions with a 1/2 probability of $|0\rangle$ and a 1/2 probability for $|1\rangle$. Then, consider another classical probability distribution with a 1/2 probability for $\frac{1}{\sqrt 2}\left(|0\rangle+|1\rangle\right)$ and a 1/2 probability for $\frac{1}{\sqrt 2}\left(|0\rangle-|1\rangle\right)$.
Both examples are described by the same density matrix $\left(\begin{array}{cc} \frac{1}{2} &0\\0&\frac{1}{2}\end{array}\right)$ and can't be distinguished empirically by any experiment. If wavefunctions are objective and fundamental, why can't we distinguish between both examples?