Bayesian approach using nonhomogeneous Poisson process is considered for modelling software reliability problem. The usefulness of the iterative sampling-based method increases greatly as the dimension of a problem increases. Maximum likelihood estimator and Gibbs estimator are derived. Model selection based on a predictive likelihood is studied. A numerical example is given.