With the advancement of multi-constellation and multi-frequency global navigation satellite systems (GNSSs), more observations are available for high precision positioning applications. Although there is a lot of progress in the GNSS world, achieving realistic precision of the solution (neither too optimistic nor too pessimistic) is still an open problem. Weighting among different GNSS systems requires a realistic stochastic model for all observations to achieve the best linear unbiased estimation (BLUE) of unknown parameters in multi-GNSS data processing mode. In addition, the correct integer ambiguity resolution (IAR) becomes crucial in shortening the Time-To-Fix (TTF) in RTK, especially in challenging environmental conditions. In general, it is required to estimate various variances for observation types, consider the correlation between different observables, and compensate for the satellite elevation dependence of the observable precision. Quality control of GNSS signals, such as GPS, GLONASS, Galileo, and BeiDou can be performed by processing a zero or short baseline double difference pseudorange and carrier phase observations using the least-squares variance component estimation (LS-VCE). The efficacy of this method is investigated using real multi-GNSS data sets collected by the Trimble NETR9, SEPT POLARX5, and LEICA GR30 receivers. The results show that the standard deviation of observations depends on the system and the observable type in which a particular receiver could have the best performance. We also note that the estimated variances and correlations among different observations are also dependent on the receiver type. It is because the approaches utilized for the recovery techniques differ from one type of receiver to another kind. The reliability of IAR will improve if a realistic stochastic model is applied in single or multi-GNSS data processing. According to the results, for the data sets considered, a realistic stochastic model can increase the computed empirical success rate to 100% in multi-GNSS as well as a single system. As mentioned previously, the realistic precision of the solution can be achieved with a realistic stochastic model. However, using the estimated stochastic model, in fact, leads to better precision and accuracy for the estimated baseline components, up to 39% in multi-GNSS.