Non-parametric shrinkage mean estimation for quadratic loss functions with unknown covariance matrices