Omschrijving van het project

Los van wat al in de Charter genoemd is toch even deze woorden.

Probleem: Plotselinge toename van klachten van externe klanten. Bijna een vertienvoudiging.

Oplossing: De band moet langzamer en dat kan alleen met behoud van capaciteit wanneer van 3 naar 4 banden/lijnen wordt overgeschakeld.

Uit de grafieken zal ook naar voren komen dat er een duidelijk verband is tussen uitval van loins en het voorkomen van product categorie 3 en de bandsnelheid.

De praktische onderliggendeoorzaak ligt in het vaardigheidsniveau van de nieuwe medewerkers aan de ‘lijn’.

Waarom dan toch proberen om via de weg van de statistiek te gaan, alleen maar om te ontdekken wat je toch al wist? Wenu, het nut zal mogelijk pas tot uiting komen wanneer meerdere invloeden worden gemeten, die ogenschijnlijk niets met elkaar van doen hebben.

Veel meten kan aan de andere kant een kostenpost zijn. Data moeten vergaard worden en bewaard worden. Als gegevens niet toch al routinematig deel uitmaken van het proces, kan voortdurend extra meten een belasting vormen voor het proces.

Kanttekening bij de gekozen oplossing: vier in plaats van drie banden gebruiken voor dezelfde productie is een vorm van verspilling.

In onderstaande mindmap wordt behalve Defects/Mankementen ook Underutilisation/Onderbenutting genoemd. De gekozen oplossing die de mankementen wil bestrijden, moet dan ook van tijdelijke aard zijn, althans voor zover men wil vasthouden aan de lean principes, ook in de vaak chaotische omstandigheden van alledag, die minder tot reflectie uitnodigen.

Nu is het de verwachting dat het inderdaad om een tijdelijke maatregel gaat. Nieuwe werkkrachten leren immers gaandeweg de nodige vaardigheden aan, waarna de snelheid naar verwachting, zo nodig stapsgewijs, naar een gemiddeld acceptabel niveau teruggebracht kan worden.

Voorlopig echter staat klanttevredenheid voorop.

Lees de metingen

Metingen inlezen en de kolomnamen veranderen. De gegevens staat, zoals gewoonlijk, in een excel bestand.

Kolom namen:

De werkuren worden doorgaans geïgnoreerd aangezien de waarden steeds hetzelfde zijn. Ze zijn wel zichtbaar in de lijst.

Hier de bovenste 6 rijen meetgegevens

Hier de 50 rijen meetreeks

Print categorie tegen de gemeten datums

Stats

Category 3

print(f41())

In deze boxplot is de bandsnelheid geen continu variabele. Let bijvoorbeeld op de sprong tussen 600 en 800.

print(fbox1(dt()))

Loins

Hier is de bandsnelheid discreet (gefactord, geen continu variabele). Van snelheid 800 bestaan slechts 3 meetpunten.

In deze boxplot is de bandsnelheid geen continu variabele. Let bijvoorbeeld op de sprong tussen 600 en 800.

print(fbox2(dt()))

Gecombineerd

print(f7())

T.Tests

stat_data_t1 <- function() {
  t  <-dt() %>% 
    select(c3_kg,v) %>% 
    filter(v==600) 
  
  return(t)
}  
stat_data_t2 <- function() {
  t  <-dt() %>% 
    select(c3_kg,v) %>% 
    filter(v==850) 
  
  return(t)
}  

t1 <- stat_data_t1()
t1
summary(t1)
##      c3_kg            v      
##  Min.   :1011   Min.   :600  
##  1st Qu.:1204   1st Qu.:600  
##  Median :1296   Median :600  
##  Mean   :1316   Mean   :600  
##  3rd Qu.:1370   3rd Qu.:600  
##  Max.   :1802   Max.   :600
var(t1)
##          c3_kg v
## c3_kg 32794.69 0
## v         0.00 0
sd(t1$c3_kg)
## [1] 181.0931
mean(t1$c3_kg)
## [1] 1316.2
max(t1$c3_kg)
## [1] 1802
t2 <- stat_data_t2()
t2
class(t2)
## [1] "data.frame"
t.test(t1,t2,var.equal=TRUE)
## 
##  Two Sample t-test
## 
## data:  t1 and t2
## t = -4.8184, df = 50, p-value = 1.387e-05
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -1641.4004  -675.5663
## sample estimates:
## mean of x mean of y 
##   958.100  2116.583

Gemiddelden van categorie 3 product in kg en loins uitval per bandsnelheid. Het aantal metingen is aangegeven met ‘n’.

## # A tibble: 5 × 4
##       v cat3_mean loins_mean     n
##   <dbl>     <dbl>      <dbl> <int>
## 1   600      1316         92    20
## 2   800      2846        418     3
## 3   850      3383        537     6
## 4   900      3326        510    10
## 5   950      3281        526    11

Data selectie op bandsnelheid 600

## # A tibble: 20 × 3
## # Groups:   v [1]
##    c3_kg ul_st     v
##    <dbl> <dbl> <dbl>
##  1  1802    80   600
##  2  1633    93   600
##  3  1300    60   600
##  4  1446    77   600
##  5  1507    60   600
##  6  1358    88   600
##  7  1309    70   600
##  8  1293    40   600
##  9  1408    87   600
## 10  1352   101   600
## 11  1288    89   600
## 12  1223   123   600
## 13  1011    78   600
## 14  1112    94   600
## 15  1205   105   600
## 16  1244   185   600
## 17  1300   139   600
## 18  1148   108   600
## 19  1185    79   600
## 20  1200    87   600

Facet diagrammen van categorie 3 product in kg bij de gemeten bandsnelheden. Het aantal meetmomenten bij snelheid 600 is het grootst.

Categorie 3 afgezet tegen de bandsnelheid (v). De bandsnelheid is hier continu weergegeven. Er zijn in de metingen maar 5 discrete bandsnelheden ingesteld geweest en in de grafiek verwerkt.

print(f5())

print(f6())
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

print(f7())

stat_data3 <- function() {
  t  <-dt() %>% 
    select(c3_kg,ul_st,v) %>% 
    group_by(v) 
  
  return(t)
}

print(fbox1(stat_data3()))

print(fbox2(stat_data3()))

Heatmap

library(rpart)
library(rpart.plot)

regres <- function(fl) {
  
  str(dt())
  
  dfCW <- dt()
  m1 <- rpart(c3_kg ~.,data=dfCW,method="anova")
  m1
  rpart.plot(m1,type=3,digits=4, fallen.leaves = fl)
  
  m2 <- rpart(ul_st ~.,data=dfCW,method="anova")
  m2
  
  rpart.plot(m2,type=2,digits=3, fallen.leaves = fl)
  
}
regres(fall <- TRUE)
## 'data.frame':    50 obs. of  5 variables:
##  $ dag  : POSIXct, format: "2022-08-22" "2022-08-23" ...
##  $ c3_kg: num  3290 2745 3745 3132 3732 ...
##  $ ul_st: num  600 400 550 620 570 450 510 480 490 520 ...
##  $ v    : num  950 950 950 950 950 950 950 950 950 950 ...
##  $ wu   : num  18 18 18 18 18 18 18 18 18 18 ...

?rpart

Uitgesplitste data (v: de bandsnelheid )

print(f4())