Data fitting to distribution and finding the probability & Reliability

I have 5 tables of data failure for 5 equipment. It show (running hour & number of failure) as we have maximum 4 types of failure. First, I want to these data to be fitted to a distribution like (normal, Beta, Gamma ext... ), we will select the best fit. Then we will check the probability of the equipment failure and its reliability. I have done some work in R studio for the distribution fitting for the first one. I need the work to be done by R studio or Python or MATLAB. Attached the tables and My coding attempt.