Есть ли в пакете weibullR какая-либо функция для получения графика времени Reliabilty v / s в R? - PullRequest
0 голосов
/ 18 января 2019

Я строю график зависимости надежности от времени в R с использованием пакета WeibullR. Но по умолчанию это печать Ненадежность против времени. Есть ли способ сделать это надежность в зависимости от графика времени? Заранее спасибо

используемая библиотека:

библиотека (WeibullR)

библиотека (WeibullR) df <- data.frame (время = c (10000, 10000, 20000, 20000, 30000, 30000, 30000, 30000, 40000, 50000, 50000, 60000, 70000, 70000, 70000, 70000, 80000, 80000, 80000, 80000, 90000, 90000, 100000), событие = rep (1, 23) weibl <- 1- wblr (df, col = "darkgreen", label = "набор цензурированных данных", dist = "weibull2p", ylab = "check") weibl_fit <- wblr.fit (weR, col = "Red", method.fit = "rr") data <- wblr.conf (weibl_fit, col = "blue") участок (данные) </p>

По умолчанию это график ненадежности [%] против времени до отказа. Ожидаемая надежность - время.

1 Ответ

0 голосов
/ 18 марта 2019

Здесь df - это кадр данных, имеющий время и событие в качестве имен двух столбцов. Для библиотеки ниже необходимо, чтобы имена ваших столбцов были в этом формате (а именно: событие и время).

Чтобы построить график зависимости надежности от времени, в качестве входных данных я использую выше df.

библиотека (WeibullR)

weibl <- wblr (df, col = "darkgreen", label = "набор данных с цензурой", dist = "weibull2p", ylab = "Ненадежность") </p>

data <- wblr.fit (weibl, col = "Red", method.fit = "mle", dist = "weibull2p") надежность <- data $ data $ dpoints </p>

график (x = надежность $ время, y = ((1-надежность $ ppp) * 100), xlab = "время", ylab = "надежность", xlim = c (1000, 100000), ylim = c ( 0,100), col = "оранжевый")

Приведенный выше график даст вам надежность по сравнению с графиком времени.

...