To be able to calculate a parameterized aeroacoustic liner impedance, a robust statistical metamodel is constructed as a function of the frequency and of the control parameters that are the percentage of open area and the sound pressure level. This construction is based on the use of simulated data generated with a computationally expensive aeroacoustic model, which translates to a very small training dataset. This means that the learning process has to be used and the probabilistic learning on manifolds algorithm is chosen. Although the aeroacoustic simulation is conducted on a large aeroacoustic computational model, some approximations are introduced, generating model errors that are taken into account by a probability model in the constructed training dataset. This probability model is calibrated using dimensionless experiments available from the open literature. Despite the fact that only a small amount of data is available, a novel statistical metamodel is successfully developed for which the predictions are consistent. This statistical framework allows for exhibiting a confidence region of the parameterized aeroacoustic liner impedance, which gives an information about the level of uncertainties as a function of the frequency and the control parameters.