The power function distribution is a flexible waiting time model that may provide better fit for some failure data. This paper presents the comparison of the maximum likelihood estimates and the Bayes estimates of two-parameter power function distribution. The Bayes estimates are obtained, using conjugate priors, under five loss functions consist of square error, precautionary, weighted, LINEX and DeGroot loss function. The Gibbs sampling algorithm is proposed to generate samples from posterior distributions and in result the Bayes estimates are computed. The comparison of the maximum likelihood estimates and the Bayes estimates are done through the root mean squared errors. One real-life data set is analyzed to illustrate the evaluation of proposed methods of estimation. Finally, results from the simulation are discussed to assess the performance behavior of the maximum likelihood estimates and the Bayes estimates.