The following code may be used for fitting an accelerated life test model in R.
Currently the R code implements the lognormal and Weibull distribution for failure times. However, the code can be easily adapted to implement other distributions as well (such as the Gumbel distribution).
As an example, temperature (on the Arrhenius scale) is used as acceleration factor in the code. But it also possible to use other acceleration factors (for instance, voltage stress). Moreover, the code can be extended to include multiple accelerators in the accelerated life test model, and even their interaction.