Using the spectroscopic imaging capabilities of the Swedish Solar Telescope, we aim to provide the first investigation on the nature and dynamics of elongated absorption features (fibrils) observed in H$\alpha$ in the internetwork. We observe and identify a number of internetwork fibrils, which form away from the kilogauss, network magnetic flux, and we provide a synoptic view on their behaviour. The internetwork fibrils are found to support wave-like behaviour, which we interpret as Magnetohydrodynamic (MHD) kink waves. The properties of these waves, that is, amplitude, period, and propagation speed, are measured from time-distance diagrams and we attempt to exploit them via magneto-seismology in order to probe the variation of plasma properties along the wave-guides. We found that the Internetwork (IN) fibrils appear, disappear, and re-appear on timescales of tens of minutes, suggesting that they are subject to repeated heating. No clear photospheric footpoints for the fibrils are found in photospheric magnetograms or H$\alpha$ wing images. However, we suggest that they are magnetised features as the majority of them show evidence of supporting propagating MHD kink waves, with a modal period of $120$~s. Additionally, one IN fibril is seen to support a flow directed along its elongated axis, suggesting a guiding field. The wave motions are found to propagate at speeds significantly greater than estimates for typical chromospheric sound speeds. Through their interpretation as kink waves, the measured speeds provide an estimate for local average Alfv\'en speeds. Furthermore, the amplitudes of the waves are also found to vary as a function of distance along the fibrils, which can be interpreted as evidence of stratification of the plasma in the neighbourhood of the IN fibril.

{Observations of the solar atmosphere have shown that magnetohydrodynamic waves are ubiquitous throughout. Improvements in instrumentation and the techniques used for measurement of the waves now enables subtleties of competing theoretical models to be compared with the observed waves behaviour. Some studies have already begun to undertake this process. However, the techniques employed for model comparison have generally been unsuitable and can lead to erroneous conclusions about the best model. The aim here is to introduce some robust statistical techniques for model comparison to the solar waves community, drawing on the experiences from other areas of astrophysics. In the process, we also aim to investigate the physics of coronal loop oscillations. } {The methodology exploits least-squares fitting to compare models to observational data. We demonstrate that the residuals between the model and observations contain significant information about the ability for the model to describe the observations, and show how they can be assessed using various statistical tests. In particular we discuss the Kolmogorov-Smirnoff one and two sample tests, as well as the runs test. We also highlight the importance of including any observational trend line in the model-fitting process.} {To demonstrate the methodology, an observation of an oscillating coronal loop undergoing standing kink motion is used. The model comparison techniques provide evidence that a Gaussian damping profile provides a better description of the observed wave attenuation than the often used exponential profile. This supports previous analysis from Pascoe et al. (2016). Further, we use the model comparison to provide evidence of time-dependent wave properties of a kink oscillation, attributing the behaviour to the thermodynamic evolution of the local plasma.}