In Auburngrads/SMRD: Statistical Methods For Reliability Data

SMRD:::vinny()
library(SMRD)

In this echapter...

Analyze the gaaslaser dataset

DT::datatable(gaaslaser, 
              options = list(pageLength=17),
              rownames = FALSE)

GaAsLaser.rmd <- 
  frame.to.rmd(gaaslaser, 
               response.column = 1, 
               unit.column = 2, 
               time.column = 3, 
               response.units = "Increase in Operating Current (%)")

summary(GaAsLaser.rmd)
plot(GaAsLaser.rmd)

Trellis Plots

trellis.plot(GaAsLaser.rmd, order.groups = T)

trellis.plot(GaAsLaser.rmd, order.groups = F)

Convert to `life.data` objects

GaAsLaser.ld <- SMRD:::rmd.to.ld(GaAsLaser.rmd, 
                          fail.level = 10, 
                          x.axis = "sqrt")

SMRD:::plot.rmd.residual(GaAsLaser.ld)

GaAsLaser.ld <- SMRD:::rmd.to.ld(GaAsLaser.rmd,
                          fail.level = 10)

SMRD:::plot.rmd.residual(GaAsLaser.ld)

summary(GaAsLaser.ld)

ML Probability plots

mleprobplot(GaAsLaser.ld, dist = "Weibull")
mleprobplot(GaAsLaser.ld, dist = "Lognormal")
mleprobplot(GaAsLaser.ld, dist = "normal")

Model with censoring times

GaAsLaser.censor.ld <- rmd.to.ld(GaAsLaser.rmd,
                                 fail.level = 8,
                                 censor.time = 3000)

summary(GaAsLaser.censor.ld)

mleprobplot(GaAsLaser.censor.ld, dist = "Weibull")
mleprobplot(GaAsLaser.censor.ld, dist = "Lognormal")
mleprobplot(GaAsLaser.censor.ld, dist = "normal")