We look for one PhD student who will forecast the capability of such advanced higher-order statistics for HI intensity mapping, in the view of the upcoming SKA Observatory. They will apply these new statistical methods to simulated HI intensity mapping data for the first time. To this purpose, it will also be necessary to modify the existing techniques to account for the peculiarities of the HI intensity mapping signal, like its superb redshift resolution, opposite to its rather poor angular resolution (compared to optical surveys). Moreover, the telescope beam plays an important role, and it will have to be modelled accurately to extract the most of the information. All of this will be done with the most advanced techniques of machine learning and data analysis of cosmological maps. The goal is to forecast constraints on cosmological parameters as per the planned SKA Observatory HI intensity mapping survey, which will cover 5000 to 15,000 square degrees.