R/mcs_delay.R
mcs_delays.Rd
mcs_delays()
is no longer under active development, switching to mcs_delay
is recommended.
mcs_delays(
date_prod,
date_register,
date_repair,
date_report,
time,
status,
distribution = "lognormal",
details = FALSE
)
A vector of class character
or Date
, in the
format "yyyy-mm-dd", indicating the date of production of a unit.
Use NA
for missing elements.
A vector of class character
or Date
, in
the format "yyyy-mm-dd", indicating the date of registration of a unit.
Use NA
for missing elements.
a vector of class character
or Date
, in the
format "yyyy-mm-dd", indicating the date of repair of a failed unit.
Use NA
for missing elements.
a vector of class character
or Date
, in the
format "yyyy-mm-dd", indicating the date of report of a failed unit.
Use NA
for missing elements.
A numeric vector of operating times.
A vector of binary data (0 or 1) indicating whether unit i is a right censored observation (= 0) or a failure (= 1).
Supposed distribution of the random variable. Only
"lognormal"
is implemented.
A logical. If FALSE
the output consists of a vector with
corrected operating times for the censored units and the input operating
times for the failed units. If TRUE
the output consists of a detailed
list, i.e the same vector as described before, simulated random numbers and
estimated distribution parameters.
A numerical vector of corrected operating times for the censored units
and the input operating times for the failed units if
details = FALSE
. If details = TRUE
the output is a list which
consists of the following elements:
time
: A numeric vector of corrected operating times for the censored
observations and input operating times for failed units.
x_sim_regist
: Simulated random numbers of specified distribution with
estimated parameters for delay in registration. The length of x_sim_regist
is equal to the number of censored observations.
x_sim_report
: Simulated random numbers of specified distribution with
estimated parameters for delay in report. The length of x_sim_report
is
equal to the number of censored observations.
coefficients_regist
: Estimated coefficients of supposed distribution for
delay in registration.
coefficients_report
: Estimated coefficients of supposed distribution for
delay in report
This function is a wrapper that combines both, mcs_delay_register and mcs_delay_report functions for the adjustment of operating times of censored units.
date_of_production <- c("2014-07-28", "2014-02-17", "2014-07-14",
"2014-06-26", "2014-03-10", "2014-05-14",
"2014-05-06", "2014-03-07", "2014-03-09",
"2014-04-13", "2014-05-20", "2014-07-07",
"2014-01-27", "2014-01-30", "2014-03-17",
"2014-02-09", "2014-04-14", "2014-04-20",
"2014-03-13", "2014-02-23", "2014-04-03",
"2014-01-08", "2014-01-08")
date_of_registration <- c("2014-08-17", "2014-03-29", "2014-12-06",
"2014-09-09", "2014-05-14", "2014-07-01",
"2014-06-16", "2014-04-03", "2014-05-23",
"2014-05-09", "2014-05-31", "2014-08-12",
"2014-04-13", "2014-02-15", "2014-07-07",
"2014-03-12", "2014-05-27", "2014-06-02",
"2014-05-20", "2014-03-21", "2014-06-19",
"2014-02-12", "2014-03-27")
date_of_repair <- c(NA, "2014-09-15", "2015-07-04", "2015-04-10", NA,
NA, "2015-04-24", NA, "2015-04-25", "2015-04-24",
"2015-06-12", NA, "2015-05-04", NA, NA,
"2015-05-22", NA, "2015-09-17", NA, "2015-08-15",
"2015-11-26", NA, NA)
date_of_report <- c(NA, "2014-10-09", "2015-08-28", "2015-04-15", NA,
NA, "2015-05-16", NA, "2015-05-28", "2015-05-15",
"2015-07-11", NA, "2015-08-14", NA, NA,
"2015-06-05", NA, "2015-10-17", NA, "2015-08-21",
"2015-12-02", NA, NA)
op_time <- rep(1000, length(date_of_repair))
status <- c(0, 1, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0)
# Example 1 - Simplified vector output:
x_corrected <- mcs_delays(
date_prod = date_of_production,
date_register = date_of_registration,
date_repair = date_of_repair,
date_report = date_of_report,
time = op_time,
status = status,
distribution = "lognormal",
details = FALSE
)
#> Warning: `mcs_delays()` was deprecated in weibulltools 2.0.0.
#> Please use `mcs_delay()` instead.
# Example 2 - Detailed list output:
list_detail <- mcs_delays(
date_prod = date_of_production,
date_register = date_of_registration,
date_repair = date_of_repair,
date_report = date_of_report,
time = op_time,
status = status,
distribution = "lognormal",
details = TRUE
)
#> Warning: `mcs_delays()` was deprecated in weibulltools 2.0.0.
#> Please use `mcs_delay()` instead.