All_Cutting.R

### calculate reliability per failure
library(readxl)
Auto <- read_xlsx('C:\\Users\\itais\\Documents\\Lab\\Results_nAuto_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='EASTMAN'
  if (i==2) b='BULMER'
  a=Auto[Auto$machine==b,]
  fa.sum=read.csv(file = paste("Cutting",b,".csv"))
  fa.sum=fa.sum[,2:3]
  fa.len=length(fa.sum[[1]])
  q<-data.frame(NULL)
  q=data.frame(ftype=a$fail,itype=a$item, event=a$`40179`)
  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[11])), origin = "1900-01-01") )
    print(M)
  }
  
}

## [1] "EASTMAN"
## [1] "מחשב"
## [1] NA
## [1] "2010-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.3824798   1.3008858 168.2717239
## 
## $haz4
## [1] 0.4624319 1.0716917 3.9302147

##       t0 
## 2850.467 
## [1] "EASTMAN"
## [1] "מיסבים"
## [1] NA
## [1] "2010-01-03"
## 3P weibull: Bad weibull 3p solution 
## 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.04878256 0.21667784 0.22069479 0.45957090
## 
## $haz4
## [1] 0.05980655 0.19404520 0.19568815 0.25476630

##       t0 
## 2577.606 
## [1] "EASTMAN"
## [1] "ראש"
## [1] "כתיבה"
## [1] "2010-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.095450e+00 7.369619e+00 5.762138e+14
## 
## $haz4
## [1]   2.533474   6.071217 121.227643

##      t0 
## 2833.07 
## [1] "EASTMAN"
## [1] "ראש"
## [1] "סכין"
## [1] "2010-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.02964791 0.06241300 0.07257982 0.08658294 0.11948201 0.13483476
## [7] 0.16292169
## 
## $haz4
## [1] 0.03228205 0.07752946 0.08353583 0.08836132 0.09132800 0.09074909
## [7] 0.08837259
## 
## $haz3
## [1] 0.02266458 0.05476819 0.06839608 0.09031178 0.15743399 0.19690945
## [7] 0.28379612

##       t0 
## 1790.764 
## [1] "BULMER"
## [1] "כבל"
## [1] NA
## [1] "2010-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.2855217    0.3014080    0.4572876    0.6635753    2.7266857
## [6]    4.2424580 4762.0719061
## 
## $haz4
## [1] 0.3784283 0.3982425 0.5691684 0.7468187 1.5666085 1.8496175 6.0555767
## 
## $haz3
## [1] 5.1923655 2.5375129 1.0039555 0.7643071 0.4766879 0.4367375 0.2266405

##       t0 
## 2742.963 
## [1] "BULMER"
## [1] "מחשב"
## [1] NA
## [1] "2010-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.1954571   0.7521400   1.1978318 108.3032159
## 
## $haz4
## [1] 0.2456676 0.7392616 0.9519901 3.0355640

##       t0 
## 2729.363 
## [1] "BULMER"
## [1] "מנוע"
## [1] NA
## [1] "2010-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] 4.879178e+00 1.723168e+01 3.616637e+36
## 
## $haz4
## [1]  5.899102e+00  1.419575e+01 1.797693e+307

##       t0 
## 2784.444 
## [1] "BULMER"
## [1] "ראש"
## [1] "סכין"
## [1] "2010-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.1874733 0.6056145 9.0742991
## 
## $haz4
## [1] 0.2266620 0.4989155 1.1465830

##       t0 
## 2632.915