The following R code implements a fatigue-limit model. Details of the implemented fatigue-limit model can be found in this 1997 paper by Pascual and Meeker, and in Meeker and Escobar’s 1998 book Statistical Methods for Reliability Data (pp. 593-597).
The implemented model assumes that the fatigue-limit is fixed (i.e., identical fatigue-limit for all specimens). But in some practical cases it may be more realistic to assume that the fatigue-limits of individual specimens are considerably different. Under such circumstances a random fatigue-limit model should be used, which allows for a different fatigue-limit for each specimen. A random fatigue-limit model is described by Pascual and Meeker in this 1998 paper. However, the implemented fixed fatigue-limit model can still be applied when the differences in individual fatigue limits between specimens are small.
Continue reading R code for a Fatigue-Limit Model