A new method for assessing the reliability of small computing systems that allows the generation of only small samples of operational data is proposed in the article. The result of applying the technique is the posterior failure distribution density, on the basis of which various reliability indicators can be calculated. The methodology consists of two stages: the first is the preparation of operational data, including detection of failures using machine learning methods, and the second is the construction of the failure distribution density using the adapted Rosenblatt–Parzen method. Increasing the efficiency of estimates using the proposed method is achieved by taking into account censored data, compensating for the shift of failure distribution densities and finding the optimal smoothing parameter.