# Dane
dane <- read.table(header = FALSE, text = "
  0.731
  0.734
  0.733
  0.734
  0.733
  0.734
  0.734
  0.733
  0.734
  0.715
  0.734
  0.732
  0.735
  0.734
  0.734
  0.731
  0.731
  0.735
  0.713
  0.731
  0.715
  0.736
  0.733
  0.733
  0.732
  0.732
  0.733
  0.715
  0.734
  0.734
  ")
library(qcc)
# The dataset
qcc(data = dane, 
# The chart type (in this case it lets qcc know that n = 1)    
type = "xbar.one", 
# Display the plot
    plot = TRUE) 
#Dataset n>1
chartdata2 <- read.table(header = FALSE, text = "
  0.07705     0.076030001 0.076959997 0.075489998 0.076920003
  0.078160003 0.077119999 0.076250002 0.076690003 0.076470003
  0.077550001 0.077079996 0.077129997 0.076839998 0.077440001
  0.07581     0.07728     0.076549999 0.077760004 0.077540003
  0.077459998 0.079109997 0.077830002 0.07711     0.078670003
  0.0766      0.078039996 0.076240003 0.077660002 0.075900003
  0.077919997 0.074440002 0.078419998 0.077370003 0.074500002
  0.077299997 0.076370001 0.075539999 0.076800004 0.078280002
  0.075460002 0.076229997 0.077959999 0.07491     0.078620002
  0.07931     0.077420004 0.07519     0.07784     0.07649
  0.076609999 0.076860003 0.077720001 0.075790003 0.078380004
  0.078050002 0.075450003 0.078400001 0.076810002 0.077770002
  0.077289999 0.078170002 0.077260002 0.076849997 0.075120002
  0.077639997 0.076820001 0.077200003 0.079180002 0.078299999
  0.077610001 0.075630002 0.079350002 0.076920003 0.075800002
  0.076509997 0.076970004 0.077919997 0.078189999 0.077660002
  0.07863     0.077200003 0.077310003 0.076530002 0.075680003
  0.075730003 0.07807     0.075580001 0.074589998 0.076499999
  0.076849997 0.076140001 0.077370003 0.075790003 0.07705
  0.076370001 0.076700002 0.075620003 0.077399999 0.076559998
  0.076180004 0.078189999 0.07626     0.076580003 0.077459998
  0.076059997 0.075350001 0.077310003 0.07553     0.077990003
  0.076889999 0.075240001 0.076609999 0.077090003 0.079120003
  0.07762     0.078919999 0.078730002 0.076849997 0.077660002
  0.078299999 0.075479999 0.076619998 0.078780003 0.07728")

# karta x-bar
xbar_chart2 <- qcc(data = chartdata2,
                   type = "xbar",
                   sizes = 5,
                   title = "Sample X-bar Chart Title", # Replacement title
                   digits = 2, # Limit the signifciant figures
                   plot = TRUE)

# Charakterystyki xbar_chart2

xbar_chart2

# Karta R-chart

xbar_chart2 <- qcc(data = chartdata2,
                   type = "R",
                   sizes = 5,
                   title = "Sample R Chart Title",
                   digits = 2,
                   plot = TRUE)
# Karta S-chart

xbar_chart2 <- qcc(data = chartdata2,
                   type = "stv",
                   sizes = 5,
                   title = "Sample R Chart Title",
                   digits = 2,
                   plot = TRUE)

#Process Capability Analysis

qq = qcc(chartdata2, type = "xbar", nsigmas = 3)  

process.capability(qq, spec.limits = c(0.071,0.075))

process.capability(qq, spec.limits = c(0.074,0.080))

qq_1 = qcc(chartdata2, type = "xbar.one", nsigmas = 3)  

process.capability(qq_1, spec.limits = c(0.071,0.075))

process.capability(qq_1, spec.limits = c(0.074,0.080))