We present results of the analysis of 70 RR Lyrae stars located in the bar of the Large Magellanic Cloud (LMC). Combining the spectroscopically determined metallicity of these stars from the literature with precise periods from the OGLE III catalog (Soszynski et al. 2009, J/AcA/59/1) and multi-epoch K_s_ photometry from the VISTA survey of the Magellanic Clouds system, we derive a new near-infrared period-luminosity-metallicity (PL_Ks_Z) relation for RR Lyrae variables. In order to fit the relation we use a fitting method developed specifically for this study. The zero-point of the relation is estimated two different ways: by assuming the value of the distance to the LMC and by using Hubble Space Telescope parallaxes of five RR Lyrae stars in the Milky Way (MW). The difference in distance moduli derived by applying these two approaches is ~0.2 mag. To investigate this point further we derive the PL_Ks_Z relation based on 23 MW RR Lyrae stars that had been analyzed in Baade-Wesselink studies. We compared the derived PL_Ks_Z relations for RR Lyrae stars in the MW and LMC. Slopes and zero-points are different, but still consistent within the errors. The shallow slope of the metallicity term is confirmed by both LMC and MW variables. The astrometric space mission Gaia is expected to provide a huge contribution to the determination of the RR Lyrae PL_Ks_Z relation; however, calculating an absolute magnitude from the trigonometric parallax of each star and fitting a PL_Ks_Z relation directly to period and absolute magnitude leads to biased results. We present a tool to achieve an unbiased solution by modeling the data and inferring the slope and zero-point of the relation via statistical methods.