Multipactor breakdown is a detrimental electromagnetic phenomenon caused by resonant secondary electron emissions synchronizing with field oscillation, which frequently takes place in powerful microwave devices and accelerating structures. Regarded as the principal failure mode of space microwave systems, multipactor may cause the performance to degenerate or even hardware operation to deteriorate catastrophically, thus multipactor becomes a major limitation in promoting the further development of space communication technology. Meanwhile, higher power capacity and volume integration accordingly lead to continuously growing multipactor hazard. In order to prevent multipactor from occurring, the accurate predictive technique to determine multipactor susceptibility has become a key issue for the mechanical design and performance optimization of microwave devices in the ground stage. Compared with the existing approaches to investigating the multipactor, statistical theories are able to conduct multipactor threshold calculation and mechanism analysis, with the stochastic nature of secondary emission fully considered from the probabilistic perspective. Currently, stationary statistical theory of multipactor has been developed for efficient multipactor threshold analysis of the parallel-plate geometry. However, it has not been further extended to the coaxial geometry which is commonly involved in radio frequency (RF) systems. For this reason, the stationary statistical modeling of the coaxial multipactor with all influencing factors considered is detailed in this paper. Due to the field nonuniformity and the secondary emission randomness, analytic equation of electron trajectories in the coaxial geometry is approximately derived by using the perturbation approach. Based on the implicit correlation between electron emission velocity and transit time, the joint probability density function is constructed for the calculation of the probability density distribution of electron transit time. Afterwards, a system of integral equations for depicting electron multiplication process in the coaxial geometry is formulated and solved with a novel and general iteration method. Finally, this stationary statistical theory is applied to the full multipactor susceptibility chart of coaxial transmission lines with typical coating materials in space engineering, such as silver, copper, alumina and alodine. A comparison shows that the calculation results are in reasonable agreement with the experimental measurements provided by the Europe Space Agent. What is more, there exists significant difference between multipactor susceptibility curves of the parallel-plate geometry and the coaxial geometry. This research is of great significance for optimizing the mechanism design and material selection of multipactor-free microwave devices.