### calculate reliability per failure
library(readxl)
Auto <- read_xlsx('C:\\Users\\itais\\Documents\\Lab\\Results_cnc_reduced.xlsx')
Sys.setlocale(category = "LC_ALL", locale = "Hebrew")
## [1] "LC_COLLATE=Hebrew_Israel.1255;LC_CTYPE=Hebrew_Israel.1255;LC_MONETARY=Hebrew_Israel.1255;LC_NUMERIC=C;LC_TIME=Hebrew_Israel.1255"
library(WeibullR)
library(fitdistrplus)
## Loading required package: MASS
## Loading required package: survival
source('~/CNC Matrix/itais_functions.R')
for (i in 1:2) {
if (i==1) b='EAGLE'
if (i==2) b='MATRIX'
a=Auto[Auto$machine==b,]
fa.sum=read.csv(file = paste("CNC",b,".csv"))
fa.sum=fa.sum[,2:3]
fa.len=length(fa.sum[[1]])
l=length(a[[1]])
q<-data.frame(NULL)
q=data.frame(ftype=a$fail,itype=a$item, event=rep(0,l))
if (i==1) {q$event=a$`40909`}
if (i==2) {q$event=a$`41640`}
for (j in 1:fa.len) {
if (fa.sum$cnt1[j]<2) next
fa.ix=regexpr("+[.]",fa.sum$fa1[j])
ftype=NULL
itype=NULL
if (fa.ix>-1) {
ftype=substr(fa.sum$fa1[j],1,fa.ix-2)
itype=substr(fa.sum$fa1[j],fa.ix+2,nchar(as.character(fa.sum$fa1[j])) )
} else {
ftype=as.character(fa.sum$fa1[j])
itype=NA
}
if (is.na(itype)){
q1 = q[q$ftype==ftype,]
} else {
q1 = q[q$ftype==ftype & q$itype==itype,]
}
#qw=wblr.fit(wblr(q1$event, main=paste('Autoclave #',i,' Reliability: ',as.character(fa.sum$x1[j])),
# dist='weibull3p', ties.handler='lowest'), modify.by.t0 =T )
#plot.wblr(qw)
print(b)
print(ftype)
print(itype)
M = f(tm = q1$event, pch = 5, title = paste('CNC ',b,' Reliability: ',as.character(ftype),as.character(itype)),
date = as.Date.numeric(as.integer(names(a[9+i])), origin = "1900-01-01") )
print(M)
}
}
## [1] "EAGLE"
## [1] "אמולסיה"
## [1] "barrel"
## [1] "2012-01-03"
## 3P weibull:
## 3P log-normal: 3p optimization did not converge
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.1081437 0.1647723 0.2589658 0.3908713 1.2582837 13.1337013
##
## $haz4
## [1] 0.1418297 0.2081339 0.2876102 0.3624624 0.5575724 0.7819658
##
## $haz3
## [1] 0.2739725 0.3008366 0.3169908 0.3276935 0.3494998 0.3792313
## t0
## 1456.741
## [1] "EAGLE"
## [1] "אמולסיה"
## [1] "pump"
## [1] "2012-01-03"
## 3P weibull:
## 3P log-normal: 3p optimization did not converge
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.09976032 0.17130338 0.25791189 0.36261330 1.06821329
##
## $haz4
## [1] 0.1265267 0.2021663 0.2356164 0.2448094 0.2013571
##
## $haz3
## [1] 0.1034124 0.1738499 0.2545113 0.3477453 0.9181240
## t0
## 761.22
## [1] "EAGLE"
## [1] "אמולסיה"
## [1] "system"
## [1] "2012-01-03"
## 3P weibull:
## 3P log-normal:
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.1100926 0.1156352 0.1207246 0.1383047 0.1393169 0.1432330 0.1525400
##
## $haz4
## [1] 0.14001696 0.15268071 0.15215614 0.10642496 0.10327198 0.09142660
## [7] 0.06676112
##
## $haz3
## [1] 0.35027206 0.19662219 0.15202031 0.08681879 0.08458372 0.07672876
## [7] 0.06193215
## t0
## 995.6092
## [1] "EAGLE"
## [1] "גריז"
## [1] NA
## [1] "2012-01-03"
## 3P weibull:
## 3P log-normal: 3p optimization did not converge
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.1151015 0.1289521 0.3905384 0.5394576 3.5193718
##
## $haz4
## [1] 0.1486060 0.1652199 0.3100362 0.3378319 0.3637238
##
## $haz3
## [1] 2.46888631 0.54955126 0.12362530 0.10336907 0.05010646
## t0
## 1377.912
## [1] "EAGLE"
## [1] "החלפת כלים"
## [1] NA
## [1] "2012-01-03"
## 3P weibull:
## 3P log-normal:
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.1781497 0.2025239 0.3918245 0.6215413 1.9156930 208.1663426
##
## $haz4
## [1] 0.2334930 0.2627568 0.4300133 0.5534713 0.8299484 1.2256742
##
## $haz3
## [1] 1.4336244 0.8294425 0.4379202 0.3629831 0.2735659 0.1511774
## t0
## 1046.49
## [1] "EAGLE"
## [1] "ואקום"
## [1] "pump"
## [1] "2012-01-03"
## 3P weibull:
## 3P log-normal:
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.1095283 0.1220939 0.1771666 0.1946287 0.2411560 0.5250732 0.6544764
## [8] 1.8924828
##
## $haz4
## [1] 0.1462644 0.1631683 0.2164521 0.2279341 0.2501553 0.2816807 0.2785140
## [8] 0.2236044
##
## $haz3
## [1] 0.5242844 0.3811437 0.2635321 0.2494879 0.2248114 0.1719356 0.1617477
## [8] 0.1250577
## t0
## 766.5964
## [1] "EAGLE"
## [1] "ואקום"
## [1] "system"
## [1] "2012-01-03"
## 3P weibull:
## 3P log-normal: Bad log-normal 3p solution
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.03876548 0.05936286 0.07725009 0.15137217 0.18961982 0.22934910
## [7] 0.23495005 0.26271786 0.27298443 0.59565247
##
## $haz4
## [1] 0.03959906 0.07722584 0.10404317 0.16541301 0.17991729 0.18921956
## [7] 0.19020651 0.19422099 0.19539177 0.19943973
##
## $haz3
## [1] 0.02621079 0.03916444 0.05289987 0.14371442 0.21991936 0.32859916
## [7] 0.34679991 0.44922889 0.49270938 5.11816591
## t0
## 1222.36
## [1] "EAGLE"
## [1] "ספינדל"
## [1] NA
## [1] "2012-01-03"
## 3P weibull:
## 3P log-normal: 3p optimization did not converge
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.06604287 0.09002792 0.16512254 0.20107252 0.20978655 0.22298160
## [7] 0.27649239 0.46815985 0.90173422
##
## $haz4
## [1] 0.0768539 0.1193384 0.1930475 0.2081765 0.2107410 0.2139975 0.2215440
## [8] 0.2179783 0.1865434
##
## $haz3
## [1] 0.07877131 0.10534294 0.17101506 0.19734548 0.20339985 0.21235835
## [7] 0.24649970 0.34997142 0.53150116
## t0
## 824.9471
## [1] "EAGLE"
## [1] "פילטר"
## [1] NA
## [1] "2012-01-03"
## 3P weibull:
## 3P log-normal:
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.06274216 0.26148582 0.27128348 0.29282905 0.31557646 0.34608990
## [7] 1.67041789 15.03737344
##
## $haz4
## [1] 0.06872171 0.30236101 0.30922893 0.32340353 0.33712495 0.35379325
## [7] 0.55984053 0.59436385
##
## $haz3
## [1] 0.08046901 0.28177438 0.28893833 0.30419880 0.31964048 0.33939749
## [7] 0.82605486 2.28889635
## t0
## 924.9426
## [1] "EAGLE"
## [1] "ציר"
## [1] NA
## [1] "2012-01-03"
## 3P weibull: Bad weibull 3p solution
## 3P log-normal: Bad log-normal 3p solution
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.01561784 0.10931245 0.11993225 0.15185619 0.16388933 0.16388933
## [7] 0.16388933 0.16720074 0.31700544
##
## $haz4
## [1] 0.00297309 0.13496589 0.14182524 0.15632794 0.16006788 0.16006788
## [7] 0.16006788 0.16097187 0.17420790
## t0
## 1470.594
## [1] "MATRIX"
## [1] "אמולסיה"
## [1] "barrel"
## [1] "2014-01-03"
## 3P weibull: Bad weibull 3p solution
## 3P log-normal: Bad log-normal 3p solution
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.1524366 0.1747749 0.1831111 0.1858455 0.2019559
##
## $haz4
## [1] 0.17861272 0.19186229 0.16257457 0.15219370 0.09617105
## t0
## 476.6246
## [1] "MATRIX"
## [1] "אמולסיה"
## [1] "system"
## [1] "2014-01-03"
## 3P weibull: Bad weibull 3p solution
## 3P log-normal: Bad log-normal 3p solution
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.0782703 0.3162126 0.3208468 0.4505566 0.6938666 4.3516047
##
## $haz4
## [1] 0.08529558 0.34822725 0.35107089 0.41375060 0.48057475 0.58834047
## t0
## 826.3322
## [1] "MATRIX"
## [1] "החלפת כלים"
## [1] NA
## [1] "2014-01-03"
## 3P weibull:
## 3P log-normal:
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.09056023 0.18968912 0.22787275 0.24042295 0.27511313 0.44532964
## [7] 0.77315644 0.82755466 3.22011130
##
## $haz4
## [1] 0.1069055 0.2479372 0.2829501 0.2927368 0.3161864 0.3829345 0.4225593
## [8] 0.4248450 0.3876923
##
## $haz3
## [1] 0.1185905 0.2271447 0.2583143 0.2678091 0.2924953 0.3913679 0.5292096
## [8] 0.5483368 1.0626953
## t0
## 612.7695
## [1] "MATRIX"
## [1] "ואקום"
## [1] "pump"
## [1] "2014-01-03"
## 3P weibull: 3p optimization did not converge
## 3P log-normal: 3p optimization did not converge
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.78478936 0.05142443 0.02671331
##
## $haz4
## [1] 0.94883863 0.04236432 0.01898474
## t0
## 2912.119
## [1] "MATRIX"
## [1] "ואקום"
## [1] "system"
## [1] "2014-01-03"
## 3P weibull: 3p optimization did not converge
## 3P log-normal: 3p optimization did not converge
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 2.961408e+00 1.022122e+01 1.472628e+26
##
## $haz4
## [1] 3.580449e+00 8.420418e+00 1.797693e+307
## t0
## 570.8738
## [1] "MATRIX"
## [1] "פילטר"
## [1] NA
## [1] "2014-01-03"
## 3P weibull:
## 3P log-normal:
## Weibull 2-paprameters and 3-parameters likelyhood
## Log-normal 2-paprameters and 3-parameters likelyhood
## Failure prediction
## $haz2
## [1] 0.1647806 0.2191518 0.2426301 0.2588591 0.2930594 0.3415332 0.3701964
## [8] 0.4291866 0.4809005 0.5413987 0.5413987 0.5698238 0.5964482 0.6655701
## [15] 0.6673406 1.2915642 1.7118470 3.4542042
##
## $haz4
## [1] 0.1991918 0.3008362 0.3366183 0.3584750 0.3974587 0.4386948 0.4568282
## [8] 0.4831231 0.4970963 0.5060073 0.5060073 0.5081553 0.5092473 0.5088731
## [15] 0.5088142 0.4410698 0.3938434 0.2725688
##
## $haz3
## [1] 0.3055992 0.3671241 0.3822754 0.3910960 0.4067486 0.4244691 0.4333040
## [8] 0.4488881 0.4604721 0.4722715 0.4722715 0.4773032 0.4817674 0.4924077
## [15] 0.4926643 0.5560858 0.5835976 0.6553023
## t0
## 305.4429