Data di ambil dari buku Montgomery, D.C., 2003. Applied Statistics and Probability for Engineers. John Wiley & Sonshalaman 316 table 7.10.
utss <- read_excel("C:/Users/ASUS/Downloads/utss.xlsx") View(utss)
selanjutnya memanggil packages qcc
library(qcc)
kemudian membuat peta kendali u Phase I
qcc(data = utss$Errors[1:10], sizes = 50, type = "u", plot = TRUE)
This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.
## List of 11
## $ call : language qcc(data = utss$Errors[1:10], type = "u", sizes = 50, plot = TRUE)
## $ type : chr "u"
## $ data.name : chr "utss$Errors[1:10]"
## $ data : num [1:10, 1] 2 3 8 1 1 4 1 4 5 1
## ..- attr(*, "dimnames")=List of 2
## $ statistics: Named num [1:10] 0.04 0.06 0.16 0.02 0.02 0.08 0.02 0.08 0.1 0.02
## ..- attr(*, "names")= chr [1:10] "1" "2" "3" "4" ...
## $ sizes : num [1:10] 50 50 50 50 50 50 50 50 50 50
## $ center : num 0.06
## $ std.dev : num 0.245
## $ nsigmas : num 3
## $ limits : num [1, 1:2] 0 0.164
## ..- attr(*, "dimnames")=List of 2
## $ violations:List of 2
## - attr(*, "class")= chr "qcc"
terlihat bahwa pada phase I tidak terdapat titik out control
selanjutnya lanjut ke phase II
qcc(data = utss\(`Errors`[1:10], sizes = 50 , type = "u", newdata = utss\)Errors[11:20], newsizes = 50 , plot = TRUE)
## List of 15
## $ call : language qcc(data = utss$Errors[1:10], type = "u", sizes = 50, newdata = utss$Errors[11:20], newsizes = 50, plot = TRUE)
## $ type : chr "u"
## $ data.name : chr "utss$Errors[1:10]"
## $ data : num [1:10, 1] 2 3 8 1 1 4 1 4 5 1
## ..- attr(*, "dimnames")=List of 2
## $ statistics : Named num [1:10] 0.04 0.06 0.16 0.02 0.02 0.08 0.02 0.08 0.1 0.02
## ..- attr(*, "names")= chr [1:10] "1" "2" "3" "4" ...
## $ sizes : num [1:10] 50 50 50 50 50 50 50 50 50 50
## $ center : num 0.06
## $ std.dev : num 0.245
## $ newstats : Named num [1:10] 0.16 0.04 0.08 0.06 0.08 0.02 0.16 0.06 0.14 0.08
## ..- attr(*, "names")= chr [1:10] "11" "12" "13" "14" ...
## $ newdata : num [1:10, 1] 8 2 4 3 4 1 8 3 7 4
## $ newsizes : num [1:10] 50 50 50 50 50 50 50 50 50 50
## $ newdata.name: chr "utss$Errors[11:20]"
## $ nsigmas : num 3
## $ limits : num [1, 1:2] 0 0.164
## ..- attr(*, "dimnames")=List of 2
## $ violations :List of 2
## - attr(*, "class")= chr "qcc"
Dari phase II terlihat tidak terdapat titik out control
You can also embed plots, for example:
Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.