R code for a statistical algorithm - repost

Hi, I need to convert an existing statistical algorithm into a functional R code. The algorithm deals with, generally speaking, a version of survival analysis. The input for the analysis is simple distribution of "time since the event" dates (in our case here - times since forest fire) for forest patches of particular size. In summary, the algorithm does three things: (1) it identifies various parameter of this distribution, (2) it determines the epochs (periods) with different survival probabilities (in our case here - hazards of burning), using the Bayes Information Criterion to find the optimal number of breakpoints in the distribution, and finally (3) assess parameters of sub-distributions. The project will require good skills in statistics and R. I would like to stress the importance of the statistical experience in particular. In the context of this project, I am a "user", not a statistician, who developed the method. I attached a text with full details on the algorithm. The attached paper contains two example datasets which should be helpful in testing the code. These will be available to the coder as data files.I also have other texts written by this author with description of the same protocol. Additional texts will be available to the coder selected for the job, however it is my understanding that the attached text should be sufficient. I would like to have the code well documented and, preferably, written as a function so I can easily plug it in into a larger R program. The requested code can take advantage of existing R functions (e.g. BIC estimation) as soon as it can generate the required range of output variables. Finally, I would like to mention that I will need help with similar projects in the future and I am interested in finding a skillful and responsible professional for long-term cooperation.