diff --git a/scripts/01_xbar_chart.R b/scripts/01_xbar_chart.R
index 6934fe5474e72f15a667649126990b2ba041047a..098d312408d17c0deb1f0e1ae8603226183f00bd 100644
--- a/scripts/01_xbar_chart.R
+++ b/scripts/01_xbar_chart.R
@@ -19,3 +19,44 @@ pntw <- t(unstack(x = pnt, form = x ~ i))
 pntw
 
 rownames(pntw) <- NULL
+
+#-----------------------------------------------------------------------
+
+library(qcc)
+
+qm <- qcc(data = pntw, type = "xbar")
+qr <- qcc(data = pntw, type = "R")
+
+plot(qm)
+plot(qr)
+
+#=======================================================================
+# As amostras do mesmo processo com mais 10 coletas de 5 elementos.
+
+url <- paste0("http://www.leg.ufpr.br/~walmes/",
+              "data/business_economics_dataset",
+              "/EXAMPLES/PAINT122.DAT")
+
+pnt <- read.table(file = url)
+names(pnt) <- c("x", "i")
+str(pnt)
+
+xtabs(~i, data = pnt)
+
+# Passar para o formato wide.
+pntw <- t(unstack(x = pnt, form = x ~ i))
+pntw
+
+rownames(pntw) <- NULL
+
+#-----------------------------------------------------------------------
+
+library(qcc)
+
+qm <- qcc(data = pntw[1:25, ], type = "xbar", newdata = pntw[26:35, ])
+qr <- qcc(data = pntw[1:25, ], type = "R", newdata = pntw[26:35, ])
+
+plot(qm)
+plot(qr)
+
+#-----------------------------------------------------------------------