Как получить значение надежности пакетом weibullR графика Reliabilty v / s Time в R? Я получаю значение eta и beta, которое отображается на графике. Я хочу получить значение надежности, которое отображается на графике. Ниже мой код.
library("WeibullR")
event <- dataset$censor
time <- dataset$RunDays
originaltest <- data.frame(event, time)
originaltest <- originaltest[!apply(is.na(originaltest) | originaltest == "", 1, all),]
weibl <- wblr(originaltest, col="darkgreen",label="censored dataset",
dist = "weibull2p", ylab = "Unreliability")
data <- wblr.fit(weibl, col = "Red", method.fit = "mle", dist = "weibull2p", is.plot.legend = TRUE)
#Calculate eta and Beta
eta <- data$fit[[1]]$eta
beta <- data$fit[[1]]$beta
# calculate the reliability on various days
reliability <- data$data$dpoints
plot(x =reliability$time, y = ((1-reliability$ppp)*100), xlab = "Rundays", ylab = "reliability",
xlim=c(0,1000), ylim=c(0,100), col="red", type = "l", grid = grid(nx =30, ny = 30, col = "yellow"))
legend("bottomright", legend=c(paste0("eta= ", eta),paste0("beta= ", beta)),
col=c("red", "blue"), box.col="green", pch = c(1,1))