### calculate reliability per failure
library(readxl)
Auto <- read_xlsx('C:\\Users\\itais\\Documents\\Lab\\Results_Autoclave_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:4) {
  a=Auto[Auto$X1==i,]
  fa.sum=read.csv(file = paste("auto",i,".csv"))
  fa.sum=fa.sum[,4:8]
  fa.len=length(fa.sum[[1]])
  while (is.na(fa.sum$common[fa.len])) fa.len=fa.len-1
  fa.sum=fa.sum[1:fa.len,]
  l=length(a[[1]])
  q<-data.frame(NULL)
  q=data.frame(ftype=a$fail,event=rep(0,l))
  if (i==1) {q$event=a$`30682`}
  if (i==2) {q$event=a$`32143`}
  if (i==3) {q$event=a$`35065`}
  if (i==4) {q$event=a$`41640`}
  for (j in 1:fa.len) {
    if (fa.sum$scount[j]<2) next
    q1 = q[q$ftype %in% c(as.character(fa.sum$x1[j]),as.character(fa.sum$x2[j]),as.character(fa.sum$x3[j]) ),]
    #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)
    
    M = f(tm = q1$event, pch = 5, title = paste('Autoclave #',i,' Reliability: ',as.character(fa.sum$x1[j])), 
          date = as.Date.numeric(as.integer(names(a[10+i])), origin = "1900-01-01") )
    print(M)
  }
  
}
## [1] "1984-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.03894323  0.04443799  0.06768001  0.08128300  0.08335210
##  [6]  0.09388456  0.12905382  0.37162989  2.53967437 16.70060169
## 
## $haz4
##  [1] 0.05328541 0.06015280 0.08498397 0.09706674 0.09877945 0.10705407
##  [7] 0.13040544 0.21669075 0.38102347 0.52749135
## 
## $haz3
##  [1] 0.1771396 0.1663520 0.1548564 0.1522721 0.1519656 0.1506339 0.1477778
##  [8] 0.1418472 0.1359852 0.1323721

##       t0 
## 10313.81 
## [1] "1984-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.01347220 0.01370338 0.01493880 0.04015387 0.07181628 0.09346019
##  [7] 0.10388070 0.11946659 0.17093742 0.18126626 0.19417764 0.23462234
## [13] 0.27057416 0.31239917 1.02805839
## 
## $haz4
##  [1] 0.01516034 0.01554440 0.01759897 0.05432762 0.08603721 0.10214161
##  [7] 0.10882843 0.11783252 0.14153842 0.14548272 0.15012384 0.16294841
## [13] 0.17264643 0.18242928 0.26153348
## 
## $haz3
##  [1] 0.01267887 0.01288891 0.01401306 0.03772470 0.06950470 0.09234451
##  [7] 0.10363120 0.12084197 0.18023986 0.19259676 0.20823534 0.25853787
## [13] 0.30482550 0.36040818 1.52993360

##       t0 
## 10955.63 
## [1] "1984-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] 2.363367e-01 2.566548e-01 2.652568e-01 2.652568e-01 7.797324e+00
## [6] 1.329105e+09
## 
## $haz4
## [1] 0.2946652 0.3141162 0.3220696 0.3220696 1.4736402 7.8996393
## 
## $haz3
## [1] 37.4182584  2.2933144  1.8635952  1.8635952  0.2210637  0.0644433

##       t0 
## 9941.451 
## [1] "1984-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.01909579 0.02108085 0.02631679 0.03391444 0.12533408 0.13931001
##  [7] 0.14053698 0.27678734 0.28054099 0.83538956
## 
## $haz4
##  [1] 0.02513737 0.02811706 0.03551152 0.04510211 0.10786810 0.11347724
##  [7] 0.11394373 0.15011779 0.15083198 0.20641406
## 
## $haz3
##  [1] 0.13144836 0.10443882 0.08561958 0.07627404 0.05890881 0.05817015
##  [7] 0.05811088 0.05424641 0.05418115 0.04983712

##       t0 
## 10905.33 
## [1] "1984-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.01559174 0.05151307 0.05973702 0.05984203 0.09363526 0.17416742
##  [7] 0.23307216 0.40348973 0.40541391 2.21509057
## 
## $haz4
##  [1] 0.01703697 0.06820970 0.07738178 0.07749379 0.10826784 0.15655484
##  [7] 0.18065607 0.22726331 0.22767077 0.37016922
## 
## $haz3
##  [1] 0.01732836 0.05809315 0.06649432 0.06659978 0.09856786 0.16383718
##  [7] 0.20541468 0.30868162 0.30974262 0.96233511

##       t0 
## 11003.11 
## [1] "1984-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.01210520 0.01669964 0.01669964 0.01768564 0.02330111 0.02979265
##  [7] 0.03742212 0.04324174 0.04357748 0.05777737 0.06102025 0.08976638
## [13] 0.08989917 0.08989917 0.09438596 0.11687551 0.13664495 0.23730578
## [19] 0.28211337 0.32222495 0.34213009 0.35780927 0.35975490 0.36913286
## [25] 0.42866446 0.43444799 1.53433804
## 
## $haz4
##  [1] 0.01397309 0.02166321 0.02166321 0.02328492 0.03217469 0.04161929
##  [7] 0.05161214 0.05852178 0.05890341 0.07358653 0.07659000 0.09907091
## [13] 0.09916072 0.09916072 0.10213349 0.11546706 0.12547037 0.16191484
## [19] 0.17353080 0.18247876 0.18651428 0.18953032 0.18989526 0.19162676
## [25] 0.20167419 0.20257334 0.28463361
## 
## $haz3
##  [1] 0.01982761 0.03183195 0.03183195 0.03393128 0.04391813 0.05272258
##  [7] 0.06084650 0.06598527 0.06626010 0.07628048 0.07822118 0.09196413
## [13] 0.09201688 0.09201688 0.09375535 0.10139977 0.10700823 0.12697261
## [19] 0.13328614 0.13815972 0.14036331 0.14201317 0.14221299 0.14316165
## [25] 0.14868836 0.14918499 0.19696646

##       t0 
## 10878.85 
## [1] "1984-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.02154986 0.02376129 0.07004259 0.07277054 0.07447969 0.14994812
##  [7] 0.32632308 0.32632308 0.39719428 3.34414199
## 
## $haz4
##  [1] 0.02719211 0.03067820 0.08654102 0.08902259 0.09054406 0.14035953
##  [7] 0.20103807 0.20103807 0.21669041 0.37908994
## 
## $haz3
##  [1] 0.02948016 0.03284979 0.08438232 0.08671323 0.08814753 0.13803426
##  [7] 0.20979996 0.20979996 0.23090040 0.54941568

##       t0 
## 10573.74 
## [1] "1984-01-03"

## 3P weibull: 3p optimization did not converge 
## 3P log-normal: 
## Weibull 2-paprameters and 3-parameters likelyhood

## Log-normal 2-paprameters and 3-parameters likelyhood

## Failure prediction 
## $haz2
## [1] 0.02496415 0.04616419 0.07297718 0.13804729 0.30775696 0.31237184
## [7] 2.30511958
## 
## $haz4
## [1] 0.03154167 0.06056089 0.08856825 0.13412685 0.19730728 0.19850781
## [7] 0.35558991

##       t0 
## 10817.91 
## [1] "1984-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.08431424  0.28331988 46.43827683
## 
## $haz4
## [1] 0.1019390 0.2334037 0.8451141

##      t0 
## 10029.5 
## [1] "1984-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.06426777 0.21327627 3.37148156
## 
## $haz4
## [1] 0.07770206 0.17570060 0.42894447

##       t0 
## 10705.43 
## [1] "1984-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.04074896 0.04131234 0.04473449 0.52642746 1.28269605
## 
## $haz4
## [1] 0.05021253 0.05075751 0.05396498 0.17086334 0.21079660
## 
## $haz3
## [1] 37.68656010  1.16130519  0.27402822  0.02078392  0.01611414

##      t0 
## 13925.7 
## [1] "1984-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.02239141 0.02707985 0.05192740 0.06493526 0.24638318 0.25264415
## [7] 0.65136043
## 
## $haz4
## [1] 0.02953426 0.03582059 0.06170999 0.07180929 0.13709987 0.13833679
## [7] 0.18331571
## 
## $haz3
## [1] 0.09524121 0.08416852 0.07235620 0.07031350 0.06320911 0.06311770
## [7] 0.06021322

##       t0 
## 10901.79 
## [1] "1984-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.01647236 0.04915209 0.06577820 0.15190313 0.18282738 0.44293363
## 
## $haz4
## [1] 0.01995622 0.06000274 0.07377409 0.11651273 0.12618593 0.17123993
## 
## $haz3
## [1] 0.01554157 0.04730320 0.06455679 0.16328176 0.20194855 0.57790269

##       t0 
## 10907.34 
## [1] "1984-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.04458389  0.04543455  0.07909287  0.26611336  0.48883347 10.47259442
## 
## $haz4
## [1] 0.05828583 0.05932568 0.09384005 0.18808962 0.23950967 0.48075158
## 
## $haz3
## [1] 7.26653434 0.98647857 0.09912799 0.04471027 0.03620407 0.01946223

##       t0 
## 12069.61 
## [1] "1984-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.02876821  0.03162030  0.03690260  0.06766441  0.07130365
##  [6]  0.07260615  0.07573379  0.08123278  0.08556108  0.09082748
## [11]  0.43777937  0.44019943  0.48777587  2.56176676 13.62782395
## 
## $haz4
##  [1] 0.04005285 0.04412414 0.05128533 0.08508702 0.08839432 0.08954996
##  [7] 0.09226742 0.09685932 0.10031882 0.10435894 0.22559541 0.22604721
## [13] 0.23446683 0.36859629 0.49104363
## 
## $haz3
##  [1] 0.1034737 0.1094255 0.1150036 0.1254380 0.1260236 0.1262198 0.1266650
##  [8] 0.1273719 0.1278710 0.1284221 0.1381282 0.1381529 0.1386046 0.1445365
## [15] 0.1489718

##       t0 
## 10256.87 
## [1] "1984-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.01029512 0.01932797 0.01981478 0.02314982 0.09434636 0.09449979
##  [7] 0.10231333 0.10247896 0.10533367 0.10533367 0.10722063 0.24979108
## [13] 0.25798149 0.27136578 0.27554922 0.37035616 0.51857174 2.10100607
## 
## $haz4
##  [1] 0.009423569 0.024262654 0.025063469 0.030469670 0.108212963
##  [6] 0.108327307 0.113965437 0.114081192 0.116053032 0.116053032
## [11] 0.117332861 0.181984169 0.184544872 0.188566715 0.189784796
## [16] 0.213429551 0.240408922 0.348967887
## 
## $haz3
##  [1] 0.01054018 0.02007446 0.02058536 0.02407627 0.09541584 0.09556440
##  [7] 0.10310612 0.10326550 0.10600946 0.10600946 0.10782004 0.23901273
## [13] 0.24628443 0.25811778 0.26180416 0.34393980 0.46793585 1.63385433

##       t0 
## 10884.61 
## [1] "1984-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.02808453 0.08488004 0.19086412
## 
## $haz4
## [1] 0.03395521 0.06992562 0.09576361

##       t0 
## 10561.44 
## [1] "1984-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.07081828  0.11380088  0.16848049  0.84832775 25.70103718
## 
## $haz4
## [1] 0.0912683 0.1367774 0.1808766 0.4001631 0.8922405
## 
## $haz3
## [1] 0.2560939 0.2182848 0.2060276 0.1834609 0.1650031

##       t0 
## 11160.67 
## [1] "1988-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] 1.292294e-01 1.292294e-01 1.344642e-01 1.399071e-01 1.417691e-01
##  [6] 1.617795e-01 1.845569e-01 1.945249e-01 3.503551e-01 3.690563e-01
## [11] 9.156418e-01 1.196864e+00 5.105194e+00 8.948605e+00 3.763150e+06
## 
## $haz4
##  [1] 0.1777959 0.1777959 0.1839294 0.1901753 0.1922821 0.2140205 0.2369445
##  [8] 0.2464346 0.3645468 0.3759472 0.5946165 0.6645579 1.0663478 1.2267846
## [15] 4.3856245
## 
## $haz3
##  [1] 4.0472945 4.0472945 1.6484584 1.2410946 1.1623260 0.7885069 0.6411319
##  [8] 0.6024933 0.4014144 0.3921604 0.2940726 0.2770621 0.2190210 0.2047808
## [15] 0.1022120

##       t0 
## 8642.746 
## [1] "1988-01-03"

## 3P weibull: Bad weibull 3p solution 
## 3P log-normal: 
## Weibull 2-paprameters and 3-parameters likelyhood

## Log-normal 2-paprameters and 3-parameters likelyhood

## Failure prediction 
## $haz2
##  [1] 0.01977878 0.02019072 0.02025020 0.02025020 0.04855576 0.11570734
##  [7] 0.15281370 0.20961105 0.21846582 0.22033103 0.25898360 0.30885720
## [13] 0.89494704
## 
## $haz4
##  [1] 0.02591295 0.02656100 0.02665430 0.02665430 0.06282420 0.11100127
##  [7] 0.12772768 0.14699651 0.14952609 0.15004584 0.15991961 0.17063783
## [13] 0.23289904

##       t0 
## 9484.875 
## [1] "1988-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.01911383 0.02029534 0.02110841 0.02110841 0.03751772 0.16590056
##  [7] 0.21506137 0.21916031 0.21916031 0.26359091 0.26359091 0.26359091
## [13] 0.30184060 1.62439456
## 
## $haz4
##  [1] 0.02362307 0.02555832 0.02688056 0.02688056 0.05120151 0.14722517
##  [7] 0.16657200 0.16798968 0.16798968 0.18189501 0.18189501 0.18189501
## [13] 0.19213085 0.31469645
## 
## $haz3
##  [1] 0.01801599 0.01910782 0.01986048 0.01986048 0.03529197 0.16932618
##  [7] 0.22544649 0.23022435 0.23022435 0.28291798 0.28291798 0.28291798
## [13] 0.32953207 2.39348883

##       t0 
## 9444.868 
## [1] "1988-01-03"

## 3P weibull: 3p optimization did not converge 
## 3P log-normal: 
## Weibull 2-paprameters and 3-parameters likelyhood

## Log-normal 2-paprameters and 3-parameters likelyhood

## Failure prediction 
## $haz2
##  [1] 0.02883141 0.04654283 0.05731564 0.05758928 0.07384326 0.13290056
##  [7] 0.17745845 0.24300829 0.43246932 0.61436857 0.88081045 7.60493188
## 
## $haz4
##  [1] 0.03780322 0.06351677 0.07718085 0.07751026 0.09569608 0.14547390
##  [7] 0.17274579 0.20390036 0.26371832 0.30108736 0.33966545 0.56170566

##       t0 
## 9587.083 
## [1] "1988-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.02575429 0.03393348 0.03654384 0.17510469 0.23977839 0.25829630
## [7] 1.07910978
## 
## $haz4
## [1] 0.03377875 0.04479891 0.04803733 0.13422528 0.15325216 0.15774513
## [7] 0.23969038
## 
## $haz3
## [1] 0.17190599 0.10814534 0.10242966 0.06595709 0.06316274 0.06256862
## [7] 0.05418383

##       t0 
## 9511.698 
## [1] "1988-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.009382861 0.014268913 0.015516164 0.018226548 0.018638747
##  [6] 0.018638747 0.018929183 0.019190879 0.019792063 0.027901777
## [11] 0.039151163 0.056549536 0.063290805 0.063491112 0.063591486
## [16] 0.065733789 0.067729780 0.103405407 0.117318712 0.130977014
## [21] 0.160689049 0.162123700 0.162845490 0.185942960 0.196008440
## [26] 0.198313046 0.207483829 0.242610366 0.250053218 0.268638589
## [31] 0.268638589 0.293053688 0.299791295 0.308853393 0.350656982
## [36] 0.367262580 0.420410250 0.435224294 0.481972109 1.591005724
## 
## $haz4
##  [1] 0.009087203 0.017256883 0.019370967 0.023918473 0.024601643
##  [6] 0.024601643 0.025081409 0.025512532 0.026498571 0.039114086
## [11] 0.054414501 0.073891550 0.080374379 0.080559399 0.080651956
## [16] 0.082602797 0.084379266 0.110856183 0.119156972 0.126513769
## [21] 0.140401775 0.141011225 0.141315964 0.150457884 0.154111365
## [26] 0.154922753 0.158064739 0.168974438 0.171087335 0.176104515
## [31] 0.176104515 0.182196513 0.183788624 0.185874382 0.194759693
## [36] 0.197994276 0.207422750 0.209832593 0.216914534 0.296447168
## 
## $haz3
##  [1] 0.01254046 0.02460756 0.02720375 0.03234571 0.03307548 0.03307548
##  [7] 0.03358220 0.03403365 0.03505289 0.04685563 0.05927881 0.07359953
## [13] 0.07815577 0.07828474 0.07834924 0.07970564 0.08093604 0.09893005
## [19] 0.10451057 0.10945831 0.11883971 0.11925334 0.11946025 0.12569363
## [25] 0.12820095 0.12875920 0.13092605 0.13851932 0.14000362 0.14354774
## [31] 0.14354774 0.14789046 0.14903293 0.15053453 0.15699679 0.15937702
## [37] 0.16640606 0.16822543 0.17362905 0.24181400

##      t0 
## 9474.73 
## [1] "1988-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.07259944  0.23947600 12.06186651
## 
## $haz4
## [1] 0.08777534 0.19728438 0.57715320

##       t0 
## 8700.709 
## [1] "1988-01-03"

## 3P weibull: 3p optimization did not converge 
## 3P log-normal: 
## Weibull 2-paprameters and 3-parameters likelyhood

## Log-normal 2-paprameters and 3-parameters likelyhood

## Failure prediction 
## $haz2
## [1] 0.01813227 0.03259828 0.06330007 0.18737977 0.18863710 0.18884741
## [7] 0.77576481
## 
## $haz4
## [1] 0.02171439 0.04309274 0.07625929 0.14130933 0.14172372 0.14179276
## [7] 0.22595359

##       t0 
## 9435.029 
## [1] "1988-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.05414749 0.17542984 1.22259794
## 
## $haz4
## [1] 0.06546627 0.14452207 0.28478377

##       t0 
## 9400.455 
## [1] "1988-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.03423788 0.05255366 0.05859374 0.05921149 0.06482646 0.58954335
## [7] 0.70035261 3.27472844
## 
## $haz4
## [1] 0.04604263 0.06742330 0.07348351 0.07408004 0.07931949 0.23222027
## [7] 0.24458635 0.34889661
## 
## $haz3
## [1] 0.1711151 0.1370500 0.1339505 0.1336843 0.1315716 0.1110085 0.1102097
## [8] 0.1045953

##       t0 
## 9265.902 
## [1] "1988-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.01104337 0.01912679 0.02345921 0.10016119 0.11072665 0.12395830
##  [7] 0.14207666 0.18270748 0.19163386 0.19625030 0.29808468 0.31506016
## [13] 0.33831818 0.38162728 1.83516961
## 
## $haz4
##  [1] 0.008537457 0.021363593 0.028563133 0.121269327 0.129926261
##  [6] 0.139912270 0.152287634 0.175846497 0.180406240 0.182691279
## [11] 0.223637161 0.229151012 0.236261989 0.248330704 0.404267478
## 
## $haz3
##  [1] 0.01060185 0.01803380 0.02205565 0.09822280 0.10935438 0.12348197
##  [7] 0.14314731 0.18849096 0.19866831 0.20396073 0.32532369 0.34634583
## [13] 0.37548565 0.43074436 2.79253262

##       t0 
## 9887.389 
## [1] "1988-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.1155912  0.3934329 43.3773757
## 
## $haz4
## [1] 0.1397539 0.3241167 1.1653948

##       t0 
## 9599.594 
## [1] "1988-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.02974737  0.04134267  0.08220585  0.09791217  0.09813091
##  [6]  0.10352196  0.47323877  0.50129632  0.71580869 17.21676144
## 
## $haz4
##  [1] 0.03863589 0.05578676 0.10345203 0.11780856 0.11799680 0.12254427
##  [7] 0.27077625 0.27677434 0.31396970 0.62079774
## 
## $haz3
##  [1] 0.06180941 0.08519847 0.12273313 0.13099272 0.13109579 0.13354934
##  [7] 0.19430327 0.19637150 0.20890835 0.30890158

##       t0 
## 9009.874 
## [1] "1988-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] 1.110547e-01 2.497379e-01 7.378099e-01 2.439984e+04
## 
## $haz4
## [1] 0.1404867 0.2730479 0.4997155 2.6113188
## 
## $haz3
## [1] 0.1938769 0.2877622 0.3868171 1.0232992

##       t0 
## 8293.343 
## [1] "1988-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] 9.193145e-02 1.558021e-01 1.721902e-01 1.222793e+00 2.325559e+03
## 
## $haz4
## [1] 0.1182321 0.1829725 0.1967476 0.5249591 1.7484891
## 
## $haz3
## [1] 0.4405707 0.2914739 0.2817024 0.2073046 0.1481185

##       t0 
## 8622.778 
## [1] "1988-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.02223812 0.02382723 0.03605663 0.03705940 0.09673277 0.09782434
##  [7] 0.09782434 0.10230637 0.11563772 0.42530041 0.56661136 0.57064864
## [13] 2.90559437
## 
## $haz4
##  [1] 0.02950104 0.03195815 0.04925391 0.05054856 0.10539979 0.10612930
##  [7] 0.10612930 0.10905653 0.11717349 0.20870110 0.22899859 0.22949785
## [13] 0.33699404
## 
## $haz3
##  [1] 0.05665787 0.06086789 0.07869192 0.07961784 0.10401098 0.10423909
##  [7] 0.10423909 0.10514126 0.10754211 0.12894762 0.13297896 0.13307650
## [13] 0.15331858

##       t0 
## 9164.702 
## [1] "1988-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.0285726 0.0844777 0.1747630
## 
## $haz4
## [1] 0.03454530 0.06959416 0.09188640

##       t0 
## 9229.169 
## [1] "1996-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.02743119 0.02849217 0.03388747 0.03437098 0.05009551 0.05130297
##  [7] 0.05862013 0.17239432 0.20795513 0.32923184 0.36465700 0.42979628
## [13] 0.53382088 2.34298507
## 
## $haz4
##  [1] 0.03746776 0.03903142 0.04663467 0.04728818 0.06631224 0.06761333
##  [7] 0.07508678 0.14329878 0.15577212 0.18609939 0.19273909 0.20329864
## [13] 0.21696250 0.29947722
## 
## $haz3
##  [1] 0.23383724 0.19239829 0.14129461 0.13932496 0.11122752 0.11019679
##  [7] 0.10516189 0.08407802 0.08194142 0.07754626 0.07669490 0.07540286
## [13] 0.07382661 0.06565317

##       t0 
## 6415.805 
## [1] "1996-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.04107888 0.05058801 0.17765993 0.20607841 0.45128013 3.82438224
## 
## $haz4
## [1] 0.05345768 0.06618503 0.16725414 0.18089003 0.25432942 0.43771596
## 
## $haz3
## [1] 0.22063529 0.16376806 0.10696868 0.10439103 0.09416641 0.07880947

##       t0 
## 6742.588 
## [1] "1996-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.02027025 0.03950192 0.05613172 0.13855175 0.18054324 0.23658957
## [7] 0.27300183 1.30477942
## 
## $haz4
## [1] 0.02398095 0.05260041 0.07239009 0.13230518 0.15091844 0.16989973
## [7] 0.17986144 0.27840486
## 
## $haz3
## [1] 0.02312449 0.04494629 0.06197035 0.13198295 0.16213260 0.19880327
## [7] 0.22094604 0.64174331

##       t0 
## 6465.421 
## [1] "1996-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.02157030 0.02328423 0.04304815 0.15396340 0.16035707 0.17282110
## [7] 0.38814062
## 
## $haz4
## [1] 0.02799668 0.03050858 0.05414617 0.11053419 0.11226352 0.11541418
## [7] 0.14636623
## 
## $haz3
## [1] 0.33817754 0.15884537 0.06166503 0.03661959 0.03622193 0.03551846
## [7] 0.02961192

##       t0 
## 7286.803 
## [1] "1996-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.04120947 0.05778355 0.05909559 0.05945810 0.06167607 0.09247353
##  [7] 0.13313974 0.63622566 1.66832794 7.92407026
## 
## $haz4
##  [1] 0.05625365 0.07493888 0.07626770 0.07663140 0.07882500 0.10466636
##  [7] 0.12985036 0.24363976 0.30851797 0.39545958
## 
## $haz3
##  [1] 0.2278489 0.1759894 0.1747332 0.1744042 0.1725398 0.1589306 0.1517341
##  [8] 0.1356565 0.1299772 0.1233083

##       t0 
## 5956.705 
## [1] "1996-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.03264152 0.04831113 0.04837908 0.05283173 0.23771485 0.51408243
## [7] 1.42709263
## 
## $haz4
## [1] 0.04338716 0.06082604 0.06089279 0.06512515 0.14399917 0.18184335
## [7] 0.22368899
## 
## $haz3
## [1] 0.18997299 0.11968402 0.11959754 0.11474382 0.08305731 0.07634965
## [7] 0.07006490

##       t0 
## 6148.078 
## [1] "1996-01-03"

## 3P weibull: Bad weibull 3p solution 
## 3P log-normal: 
## Weibull 2-paprameters and 3-parameters likelyhood

## Log-normal 2-paprameters and 3-parameters likelyhood

## Failure prediction 
## $haz2
##  [1] 0.01593654 0.01625147 0.13463265 0.19478943 0.19689643 0.20189485
##  [7] 0.24414494 0.25782904 0.28834678 1.72469807
## 
## $haz4
##  [1] 0.01400615 0.01449710 0.15130076 0.18772583 0.18880839 0.19133444
##  [7] 0.21061616 0.21618398 0.22763496 0.40407528

##      t0 
## 6954.36 
## [1] "1996-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.04497607 0.06522795 0.26819174 0.81693808
## 
## $haz4
## [1] 0.05619809 0.07588687 0.16042691 0.22226496
## 
## $haz3
## [1] 0.31598916 0.12079399 0.05877764 0.04568520

##       t0 
## 7180.992 
## [1] "1996-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.02592926 0.03911453 0.04795852 0.11748270 0.20835002 0.29115895
## [7] 1.07216706
## 
## $haz4
## [1] 0.03392814 0.05122070 0.06094833 0.10873192 0.14002255 0.15752098
## [7] 0.21584179
## 
## $haz3
## [1] 0.06629758 0.07684005 0.08024509 0.09052688 0.09524419 0.09765272
## [7] 0.10558989

##       t0 
## 6192.396 
## [1] "1996-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.04754154  0.07056195  0.07440942  0.17482859  0.19415821  1.01207203
## [7] 23.85345824
## 
## $haz4
## [1] 0.06323974 0.09035204 0.09434699 0.16698707 0.17666845 0.33182629
## [7] 0.56400745
## 
## $haz3
## [1] 0.1711888 0.1631678 0.1626023 0.1572325 0.1568098 0.1523628 0.1478844

##       t0 
## 5726.186 
## [1] "1996-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.05316565 0.16417505 0.97336898
## 
## $haz4
## [1] 0.06427919 0.13525019 0.24202301

##       t0 
## 6379.531 
## [1] "1996-01-03"

## 3P weibull: 3p optimization did not converge 
## 3P log-normal: 
## Weibull 2-paprameters and 3-parameters likelyhood

## Log-normal 2-paprameters and 3-parameters likelyhood

## Failure prediction 
## $haz2
##  [1]  0.03040932  0.06181496  0.06181496  0.12904216  0.12904216
##  [6]  0.13071888  0.13184781  0.21497710  2.45802071 16.43154124
## 
## $haz4
##  [1] 0.04045749 0.08147454 0.08147454 0.13796255 0.13796255 0.13904245
##  [7] 0.13976305 0.18208905 0.39186213 0.51879312

##       t0 
## 5853.445 
## [1] "1996-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.09672723   0.19447606   0.20002505   0.23670819   0.90248048
## [6]   2.08227499 358.08538258
## 
## $haz4
## [1] 0.1281296 0.2419117 0.2472922 0.2806590 0.6003400 0.8286683 2.2238290
## 
## $haz3
## [1] 0.2377079 0.3395638 0.3424478 0.3587006 0.4527923 0.4955211 0.6714742

##       t0 
## 7197.623 
## [1] "1996-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.1782335    0.6065699 4949.5235194
## 
## $haz4
## [1] 0.2154907 0.4997026 2.7452620

##      t0 
## 6206.63 
## [1] "1996-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.07589061   0.12168749   0.12168749   0.13806623   0.40215348
## [6]   1.98055309 476.57388754
## 
## $haz4
## [1] 0.1009357 0.1509711 0.1509711 0.1659965 0.3133119 0.5611956 1.2657928
## 
## $haz3
## [1] 0.3316379 0.2811249 0.2811249 0.2757797 0.2518848 0.2360707 0.2125025

##       t0 
## 5949.682 
## [1] "1996-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.01773046 0.03341763 0.05080854 0.08823798 0.12702600 0.19320645
## [7] 0.24448918 0.61023483
## 
## $haz4
## [1] 0.02158148 0.04444497 0.06377832 0.09214486 0.11147635 0.13336342
## [7] 0.14522588 0.18679438
## 
## $haz3
## [1] 0.02015112 0.03787745 0.05544008 0.08847688 0.11832981 0.16293276
## [7] 0.19378829 0.36703943

##      t0 
## 6457.36 
## [1] "1996-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.01635138 0.03830760 0.05468265 0.10023572 0.18599447 0.18599447
## [7] 0.22098754 0.80439535
## 
## $haz4
## [1] 0.01834484 0.05091020 0.06934720 0.10439393 0.14147274 0.14147274
## [7] 0.15161542 0.21977483
## 
## $haz3
## [1] 0.01673563 0.03917985 0.05556794 0.09999565 0.18048772 0.18048772
## [7] 0.21248599 0.70912319

##   t0 
## 6396 
## [1] "1996-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.02635978 0.02693204 0.11473061 0.15119937 0.30427013
## 
## $haz4
## [1] 0.03400663 0.03471277 0.08676316 0.09582092 0.11544598
## 
## $haz3
## [1] 0.02246632 0.02295962 0.12503660 0.18276158 0.52494624

##       t0 
## 6071.364 
## [1] "2014-01-03"

## 3P weibull: 3p optimization did not converge 
## 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.05434004 0.21519176 0.27145083 0.33905719 0.36470291 0.42914777
## [7] 1.50962979
## 
## $haz4
## [1] 0.04380732 0.25963241 0.29404106 0.32166653 0.32948816 0.34468640
## [7] 0.37017908

##       t0 
## 784.5534 
## [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.09793478 0.27070270 0.27070270 0.33996837 0.76593591 2.74745052
## 
## $haz4
## [1] 0.1193022 0.3024544 0.3024544 0.3368151 0.4123173 0.4038458
## 
## $haz3
## [1] 0.1119087 0.2777906 0.2777906 0.3330172 0.6117235 1.4778564

##       t0 
## 697.0014 
## [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.1161518 0.1208329 0.1233232
## 
## $haz4
## [1] 0.14043168 0.09954421 0.07259142

##       t0 
## 848.8073 
## [1] "2014-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.1027106 0.1027106 0.4329741 0.4372613 0.4398406 0.4432877 1.8161158
## 
## $haz4
## [1] 0.1208585 0.1208585 0.3614325 0.3621363 0.3625489 0.3630876 0.3368287
## 
## $haz3
## [1]   0.06239392   0.06239392   0.64712447   0.66317676   0.67299645
## [6]   0.68631331 195.94910957

##       t0 
## 603.3656 
## [1] "2014-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.1780640 0.1820607 0.1836684 0.5997766 0.8383254 2.6790118
## 
## $haz4
## [1] 0.2308121 0.2350373 0.2367026 0.3751991 0.3789332 0.3181185
## 
## $haz3
## [1] 83.93990975  2.57508858  2.02080745  0.10592637  0.08124782  0.03902651

##       t0 
## 1446.968 
## [1] "2014-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] 5.079371e-01 2.961096e+00 2.961096e+00 8.880003e+07
## 
## $haz4
## [1]  0.6216648  2.6390554  2.6390554 16.5070006
## 
## $haz3
## [1]   0.6770723   2.2936435   2.2936435 142.2578742

##       t0 
## 850.8402 
## [1] "2014-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.04910198  0.09072687  0.09140483  0.09697132  0.13705871
##  [6]  0.28989824  0.40848745  0.43069223  0.43983200  0.49788065
## [11]  0.50297365  0.71167351  0.74924399  0.74924399  0.75273435
## [16]  1.03478226  1.03478226  1.03926689  1.20599074 18.68561986
## 
## $haz4
##  [1] 0.0398233 0.1112778 0.1124540 0.1220582 0.1872071 0.3659497 0.4543381
##  [8] 0.4678819 0.4732368 0.5045867 0.5071375 0.5913088 0.6032306 0.6032306
## [15] 0.6042996 0.6740717 0.6740717 0.6749707 0.7050452 0.9631222
## 
## $haz3
##  [1]   0.04308313   0.07489581   0.07543068   0.07984474   0.11285786
##  [6]   0.25916017   0.39476332   0.42226201   0.43377262   0.50948159
## [11]   0.51633811   0.82682290   0.88881858   0.88881858   0.89467244
## [16]   1.42080870   1.42080870   1.43002376   1.79159401 363.48630267

##       t0 
## 843.2448 
## [1] 0
## [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.3853534  1.0214055 88.0156520
## 
## $haz4
## [1] 0.4659062 0.8414512 1.0834827

##       t0 
## 455.3413 
## [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.2390509 0.5561514 4.9597762
## 
## $haz4
## [1] 0.2890212 0.4581669 0.4760149

##       t0 
## 553.8678