qcc: Quality Control Charts
In qcc: Quality Control Charts
Description
Usage
Arguments
Value
Note
Author(s)
References
See Also
Examples
Description
Create an object of class 'qcc' to perform statistical quality control. This object may then be used to plot Shewhart charts, drawing OC curves, computes capability indices, and more.
Usage
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
qcc(data, type, sizes, center, std.dev, limits,
data.name, labels, newdata, newsizes, newdata.name,
newlabels, nsigmas = 3, confidence.level,
rules = shewhart.rules, plot = TRUE, ...)
## S3 method for class 'qcc'
print(x, ...)
## S3 method for class 'qcc'
summary(object, digits = getOption("digits"), ...)
## S3 method for class 'qcc'
plot(x, add.stats = TRUE, chart.all = TRUE,
label.limits = c("LCL ", "UCL"), title, xlab, ylab, ylim,
axes.las = 0, digits = getOption("digits"),
restore.par = TRUE, ...)
Arguments
data
a data frame, a matrix or a vector containing observed data for
the variable to chart. Each row of a data frame or a matrix, and each
value of a vector, refers to a sample or ”rationale group”.
type
a character string specifying the group statistics to compute.
Available methods are:
Statistic charted Chart description
"xbar" mean means of a continuous process variable
"R" range ranges of a continuous process variable
"S" standard deviation standard deviations of a continuous variable
"xbar.one" mean one-at-time data of a continuous process variable
"p" proportion proportion of nonconforming units
"np" count number of nonconforming units
"c" count nonconformities per unit
"u" count average nonconformities per unit
"g" count number of non-events between events
Furthermore, a user specified type of chart, say "newchart", can be provided. This requires the definition of "stats.newchart", "sd.newchart", and "limits.newchart". As an example, see stats.xbar.
sizes
a value or a vector of values specifying the sample sizes associated with each group. For continuous data provided as data frame or matrix the sample sizes are obtained counting the non-NA elements of each row. For "p", "np" and "u" charts the argument sizes is required.
center
a value specifying the center of group statistics or the ”target” value of the process.
std.dev
a value or an available method specifying the within-group standard deviation(s) of the process.
Several methods are available for estimating the standard deviation in case of a continuous process variable; see sd.xbar, sd.xbar.one, sd.R, sd.S.
limits
a two-values vector specifying control limits.
data.name
a string specifying the name of the variable which appears on the plots. If not provided is taken from the object given as data.
labels
a character vector of labels for each group.
newdata
a data frame, matrix or vector, as for the data argument, providing further data to plot but not included in the computations.
newsizes
a vector as for the sizes argument providing further data sizes to plot but not included in the computations.
newdata.name
a string specifying the name of the variable which appears on the plots. If not provided is taken from the object given as newdata.
newlabels
a character vector of labels for each new group defined in the argument newdata.
nsigmas
a numeric value specifying the number of sigmas to use for computing control limits. It is ignored when the confidence.level argument is provided.
confidence.level
a numeric value between 0 and 1 specifying the confidence level of the computed probability limits.
rules
a function of rules to apply to the chart. By default, the shewhart.rules function is used.
plot
logical. If TRUE a Shewhart chart is plotted.
add.stats
a logical value indicating whether statistics and other information should be printed at the bottom of the chart.
chart.all
a logical value indicating whether both statistics for data and for newdata (if given) should be plotted.
label.limits
a character vector specifying the labels for control limits.
title
a string giving the label for the main title.
xlab
a string giving the label for the x-axis.
ylab
a string giving the label for the y-axis.
ylim
a numeric vector specifying the limits for the y-axis.
axes.las
numeric in {0,1,2,3} specifying the style of axis labels. See help(par).
digits
the number of significant digits to use.
restore.par
a logical value indicating whether the previous par settings must be restored. If you need to add points, lines, etc. to a control chart set this to FALSE.
object
an object of class 'qcc'.
x
an object of class 'qcc'.
...
additional arguments to be passed to the generic function.
Value
Returns an object of class 'qcc'.
Note
For a nice blog post discussing the qcc package, in particular how to implement the Western Eletric Rules (WER), see http://blog.yhathq.com/posts/quality-control-in-r.html.
Author(s)
Luca Scrucca
References
Mason, R.L. and Young, J.C. (2002) Multivariate Statistical Process Control with Industrial Applications, SIAM.
Montgomery, D.C. (2005) Introduction to Statistical Quality Control, 5th ed. New York: John Wiley & Sons.
Ryan, T. P. (2000), Statistical Methods for Quality Improvement, 2nd ed. New York: John Wiley & Sons, Inc.
Scrucca, L. (2004). qcc: an R package for quality control charting and statistical process control. R News 4/1, 11-17.
Wetherill, G.B. and Brown, D.W. (1991) Statistical Process Control. New York: Chapman & Hall.
See Also
shewhart.rules, cusum, ewma, oc.curves, process.capability, qcc.groups
Examples
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
##
## Continuous data
##
data(pistonrings)
attach(pistonrings)
diameter <- qcc.groups(diameter, sample)
qcc(diameter[1:25,], type="xbar")
qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,])
q <- qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,], plot=FALSE)
plot(q, chart.all=FALSE)
qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,], nsigmas=2)
qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,], confidence.level=0.99)
qcc(diameter[1:25,], type="R")
qcc(diameter[1:25,], type="R", newdata=diameter[26:40,])
qcc(diameter[1:25,], type="S")
qcc(diameter[1:25,], type="S", newdata=diameter[26:40,])
# add warning limits at 2 std. deviations
q <- qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,], plot=FALSE)
(warn.limits <- limits.xbar(q$center, q$std.dev, q$sizes, 2))
plot(q, restore.par = FALSE)
abline(h = warn.limits, lty = 3, col = "chocolate")
# variable control limits
out <- c(9, 10, 30, 35, 45, 64, 65, 74, 75, 85, 99, 100)
diameter <- qcc.groups(pistonrings$diameter[-out], sample[-out])
qcc(diameter[1:25,], type="xbar")
qcc(diameter[1:25,], type="R")
qcc(diameter[1:25,], type="S")
qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,])
qcc(diameter[1:25,], type="R", newdata=diameter[26:40,])
qcc(diameter[1:25,], type="S", newdata=diameter[26:40,])
detach(pistonrings)
##
## Attribute data
##
data(orangejuice)
attach(orangejuice)
qcc(D[trial], sizes=size[trial], type="p")
# remove out-of-control points (see help(orangejuice) for the reasons)
inc <- setdiff(which(trial), c(15,23))
q1 <- qcc(D[inc], sizes=size[inc], type="p")
qcc(D[inc], sizes=size[inc], type="p", newdata=D[!trial], newsizes=size[!trial])
detach(orangejuice)
data(orangejuice2)
attach(orangejuice2)
names(D) <- sample
qcc(D[trial], sizes=size[trial], type="p")
q2 <- qcc(D[trial], sizes=size[trial], type="p", newdata=D[!trial], newsizes=size[!trial])
detach(orangejuice2)
# put on the same graph the two orange juice samples
oldpar <- par(no.readonly = TRUE)
par(mfrow=c(1,2), mar=c(5,5,3,0))
plot(q1, title="First samples", ylim=c(0,0.5), add.stats=FALSE, restore.par=FALSE)
par("mar"=c(5,0,3,3), yaxt="n")
plot(q2, title="Second samples", add.stats=FALSE, ylim=c(0,0.5))
par(oldpar)
data(circuit)
attach(circuit)
qcc(x[trial], sizes=size[trial], type="c")
# remove out-of-control points (see help(circuit) for the reasons)
inc <- setdiff(which(trial), c(6,20))
qcc(x[inc], sizes=size[inc], type="c", labels=inc)
qcc(x[inc], sizes=size[inc], type="c", labels=inc,
newdata=x[!trial], newsizes=size[!trial], newlabels=which(!trial))
qcc(x[inc], sizes=size[inc], type="u", labels=inc,
newdata=x[!trial], newsizes=size[!trial], newlabels=which(!trial))
detach(circuit)
data(pcmanufact)
attach(pcmanufact)
qcc(x, sizes=size, type="u")
detach(pcmanufact)
data(dyedcloth)
attach(dyedcloth)
qcc(x, sizes=size, type="u")
# standardized control chart
q <- qcc(x, sizes=size, type="u", plot=FALSE)
z <- (q$statistics - q$center)/sqrt(q$center/q$size)
plot(z, type="o", ylim=range(z,3,-3), pch=16)
abline(h=0, lty=2)
abline(h=c(-3,3), lty=2)
detach(dyedcloth)
##
## Continuous one-at-time data
##
# viscosity data (Montgomery, pag. 242)
x <- c(33.75, 33.05, 34, 33.81, 33.46, 34.02, 33.68, 33.27, 33.49, 33.20,
33.62, 33.00, 33.54, 33.12, 33.84)
qcc(x, type="xbar.one")
qcc(x, type="xbar.one", std.dev = "SD")
Example output
Package 'qcc' version 2.7
Type 'citation("qcc")' for citing this R package in publications.
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "xbar")
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00979
$ nsigmas : num 3
$ limits : num [1, 1:2] 74 74
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "xbar", newdata = diameter[26:40, ])
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00979
$ newstats : Named num [1:15] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1, 1:2] 74 74
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "xbar", newdata = diameter[26:40, ], nsigmas = 2)
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00979
$ newstats : Named num [1:15] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 2
$ limits : num [1, 1:2] 74 74
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "xbar", newdata = diameter[26:40, ], confidence.level = 0.99)
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00979
$ newstats : Named num [1:15] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name : chr "diameter[26:40, ]"
$ confidence.level: num 0.99
$ limits : num [1, 1:2] 74 74
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "R")
$ type : chr "R"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 0.038 0.019 0.036 0.022 0.026 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.0228
$ std.dev : num 0.00979
$ nsigmas : num 3
$ limits : num [1, 1:2] 0 0.0481
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "R", newdata = diameter[26:40, ])
$ type : chr "R"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 0.038 0.019 0.036 0.022 0.026 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.0228
$ std.dev : num 0.00979
$ newstats : Named num [1:15] 0.044 0.025 0.015 0.019 0.017 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1, 1:2] 0 0.0481
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "S")
$ type : chr "S"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 0.01477 0.0075 0.01475 0.00908 0.01222 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.00924
$ std.dev : num 0.00983
$ nsigmas : num 3
$ limits : num [1, 1:2] 0 0.0193
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "S", newdata = diameter[26:40, ])
$ type : chr "S"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 0.01477 0.0075 0.01475 0.00908 0.01222 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.00924
$ std.dev : num 0.00983
$ newstats : Named num [1:15] 0.01655 0.01033 0.00691 0.0075 0.00673 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1, 1:2] 0 0.0193
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
LCL UCL
73.99242 74.00993
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "xbar")
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00986
$ nsigmas : num 3
$ limits : num [1:25, 1:2] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "R")
$ type : chr "R"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 0.038 0.009 0.036 0.022 0.026 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.0223
$ std.dev : num 0.00986
$ nsigmas : num 3
$ limits : num [1:25, 1:2] 0 0 0 0 0 0 0 0 0 0 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "S")
$ type : chr "S"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 0.01477 0.00458 0.01475 0.00908 0.01222 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.00939
$ std.dev : num 0.00993
$ nsigmas : num 3
$ limits : num [1:25, 1:2] 0 0 0 0 0 0 0 0 0 0 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "xbar", newdata = diameter[26:40, ])
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00986
$ newstats : Named num [1:15] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1:40, 1:2] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "R", newdata = diameter[26:40, ])
$ type : chr "R"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 0.038 0.009 0.036 0.022 0.026 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.0223
$ std.dev : num 0.00986
$ newstats : Named num [1:15] 0.044 0.025 0.015 0.019 0.017 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1:40, 1:2] 0 0 0 0 0 0 0 0 0 0 ...
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "S", newdata = diameter[26:40, ])
$ type : chr "S"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 0.01477 0.00458 0.01475 0.00908 0.01222 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.00939
$ std.dev : num 0.00993
$ newstats : Named num [1:15] 0.01655 0.01033 0.00691 0.0075 0.00673 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1:40, 1:2] 0 0 0 0 0 0 0 0 0 0 ...
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = D[trial], type = "p", sizes = size[trial])
$ type : chr "p"
$ data.name : chr "D[trial]"
$ data : int [1:30, 1] 12 15 8 10 4 7 16 9 14 10 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:30] 0.24 0.3 0.16 0.2 0.08 0.14 0.32 0.18 0.28 0.2 ...
..- attr(*, "names")= chr [1:30] "1" "2" "3" "4" ...
$ sizes : int [1:30] 50 50 50 50 50 50 50 50 50 50 ...
$ center : num 0.231
$ std.dev : num 0.422
$ nsigmas : num 3
$ limits : num [1:30, 1:2] 0.0524 0.0524 0.0524 0.0524 0.0524 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = D[inc], type = "p", sizes = size[inc], newdata = D[!trial], newsizes = size[!trial])
$ type : chr "p"
$ data.name : chr "D[inc]"
$ data : int [1:28, 1] 12 15 8 10 4 7 16 9 14 10 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:28] 0.24 0.3 0.16 0.2 0.08 0.14 0.32 0.18 0.28 0.2 ...
..- attr(*, "names")= chr [1:28] "1" "2" "3" "4" ...
$ sizes : int [1:28] 50 50 50 50 50 50 50 50 50 50 ...
$ center : num 0.215
$ std.dev : num 0.411
$ newstats : Named num [1:24] 0.18 0.12 0.24 0.1 0.12 0.08 0.12 0.06 0.14 0.12 ...
..- attr(*, "names")= chr [1:24] "29" "30" "31" "32" ...
$ newdata : int [1:24, 1] 9 6 12 5 6 4 6 3 7 6 ...
$ newsizes : int [1:24] 50 50 50 50 50 50 50 50 50 50 ...
$ newdata.name: chr "D[!trial]"
$ nsigmas : num 3
$ limits : num [1:52, 1:2] 0.0407 0.0407 0.0407 0.0407 0.0407 ...
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = D[trial], type = "p", sizes = size[trial])
$ type : chr "p"
$ data.name : chr "D[trial]"
$ data : int [1:24, 1] 9 6 12 5 6 4 6 3 7 6 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:24] 0.18 0.12 0.24 0.1 0.12 0.08 0.12 0.06 0.14 0.12 ...
..- attr(*, "names")= chr [1:24] "31" "32" "33" "34" ...
$ sizes : int [1:24] 50 50 50 50 50 50 50 50 50 50 ...
$ center : num 0.111
$ std.dev : num 0.314
$ nsigmas : num 3
$ limits : num [1:24, 1:2] 0 0 0 0 0 0 0 0 0 0 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x[trial], type = "c", sizes = size[trial])
$ type : chr "c"
$ data.name : chr "x[trial]"
$ data : int [1:26, 1] 21 24 16 12 15 5 28 20 31 25 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named int [1:26] 21 24 16 12 15 5 28 20 31 25 ...
..- attr(*, "names")= chr [1:26] "1" "2" "3" "4" ...
$ sizes : int [1:26] 100 100 100 100 100 100 100 100 100 100 ...
$ center : num 19.8
$ std.dev : num 4.45
$ nsigmas : num 3
$ limits : num [1, 1:2] 6.48 33.21
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x[inc], type = "c", sizes = size[inc], labels = inc)
$ type : chr "c"
$ data.name : chr "x[inc]"
$ data : int [1:24, 1] 21 24 16 12 15 28 20 31 25 20 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named int [1:24] 21 24 16 12 15 28 20 31 25 20 ...
..- attr(*, "names")= chr [1:24] "1" "2" "3" "4" ...
$ sizes : int [1:24] 100 100 100 100 100 100 100 100 100 100 ...
$ center : num 19.7
$ std.dev : num 4.43
$ nsigmas : num 3
$ limits : num [1, 1:2] 6.36 32.97
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = x[inc], type = "c", sizes = size[inc], labels = inc, newdata = x[!trial], newsizes = size[!trial]| __truncated__
$ type : chr "c"
$ data.name : chr "x[inc]"
$ data : int [1:24, 1] 21 24 16 12 15 28 20 31 25 20 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named int [1:24] 21 24 16 12 15 28 20 31 25 20 ...
..- attr(*, "names")= chr [1:24] "1" "2" "3" "4" ...
$ sizes : int [1:24] 100 100 100 100 100 100 100 100 100 100 ...
$ center : num 19.7
$ std.dev : num 4.43
$ newstats : Named int [1:20] 16 18 12 15 24 21 28 20 25 19 ...
..- attr(*, "names")= chr [1:20] "27" "28" "29" "30" ...
$ newdata : int [1:20, 1] 16 18 12 15 24 21 28 20 25 19 ...
$ newsizes : int [1:20] 100 100 100 100 100 100 100 100 100 100 ...
$ newdata.name: chr "x[!trial]"
$ nsigmas : num 3
$ limits : num [1, 1:2] 6.36 32.97
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = x[inc], type = "u", sizes = size[inc], labels = inc, newdata = x[!trial], newsizes = size[!trial]| __truncated__
$ type : chr "u"
$ data.name : chr "x[inc]"
$ data : int [1:24, 1] 21 24 16 12 15 28 20 31 25 20 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:24] 0.21 0.24 0.16 0.12 0.15 0.28 0.2 0.31 0.25 0.2 ...
..- attr(*, "names")= chr [1:24] "1" "2" "3" "4" ...
$ sizes : int [1:24] 100 100 100 100 100 100 100 100 100 100 ...
$ center : num 0.197
$ std.dev : num 0.443
$ newstats : Named num [1:20] 0.16 0.18 0.12 0.15 0.24 0.21 0.28 0.2 0.25 0.19 ...
..- attr(*, "names")= chr [1:20] "27" "28" "29" "30" ...
$ newdata : int [1:20, 1] 16 18 12 15 24 21 28 20 25 19 ...
$ newsizes : int [1:20] 100 100 100 100 100 100 100 100 100 100 ...
$ newdata.name: chr "x[!trial]"
$ nsigmas : num 3
$ limits : num [1, 1:2] 0.0636 0.3297
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x, type = "u", sizes = size)
$ type : chr "u"
$ data.name : chr "x"
$ data : int [1:20, 1] 10 12 8 14 10 16 11 7 10 15 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:20] 2 2.4 1.6 2.8 2 3.2 2.2 1.4 2 3 ...
..- attr(*, "names")= chr [1:20] "1" "2" "3" "4" ...
$ sizes : int [1:20] 5 5 5 5 5 5 5 5 5 5 ...
$ center : num 1.93
$ std.dev : num 1.39
$ nsigmas : num 3
$ limits : num [1, 1:2] 0.0661 3.7939
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x, type = "u", sizes = size)
$ type : chr "u"
$ data.name : chr "x"
$ data : int [1:10, 1] 14 12 20 11 7 10 21 16 19 23
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:10] 1.4 1.5 1.538 1.1 0.737 ...
..- attr(*, "names")= chr [1:10] "1" "2" "3" "4" ...
$ sizes : num [1:10] 10 8 13 10 9.5 10 12 10.5 12 12.5
$ center : num 1.42
$ std.dev : num 1.19
$ nsigmas : num 3
$ limits : num [1:10, 1:2] 0.291 0.158 0.431 0.291 0.262 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x, type = "xbar.one")
$ type : chr "xbar.one"
$ data.name : chr "x"
$ data : num [1:15, 1] 33.8 33 34 33.8 33.5 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:15] 33.8 33 34 33.8 33.5 ...
..- attr(*, "names")= chr [1:15] "1" "2" "3" "4" ...
$ sizes : int [1:15] 1 1 1 1 1 1 1 1 1 1 ...
$ center : num 33.5
$ std.dev : num 0.426
$ nsigmas : num 3
$ limits : num [1, 1:2] 32.2 34.8
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x, type = "xbar.one", std.dev = "SD")
$ type : chr "xbar.one"
$ data.name : chr "x"
$ data : num [1:15, 1] 33.8 33 34 33.8 33.5 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:15] 33.8 33 34 33.8 33.5 ...
..- attr(*, "names")= chr [1:15] "1" "2" "3" "4" ...
$ sizes : int [1:15] 1 1 1 1 1 1 1 1 1 1 ...
$ center : num 33.5
$ std.dev : num 0.342
$ nsigmas : num 3
$ limits : num [1, 1:2] 32.5 34.5
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
qcc documentation built on May 2, 2019, 9:15 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Description Usage Arguments Value Note Author(s) References See Also Examples
Description
Create an object of class 'qcc' to perform statistical quality control. This object may then be used to plot Shewhart charts, drawing OC curves, computes capability indices, and more.
Usage
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | qcc(data, type, sizes, center, std.dev, limits,
data.name, labels, newdata, newsizes, newdata.name,
newlabels, nsigmas = 3, confidence.level,
rules = shewhart.rules, plot = TRUE, ...)
## S3 method for class 'qcc'
print(x, ...)
## S3 method for class 'qcc'
summary(object, digits = getOption("digits"), ...)
## S3 method for class 'qcc'
plot(x, add.stats = TRUE, chart.all = TRUE,
label.limits = c("LCL ", "UCL"), title, xlab, ylab, ylim,
axes.las = 0, digits = getOption("digits"),
restore.par = TRUE, ...)
|
Arguments
data |
a data frame, a matrix or a vector containing observed data for the variable to chart. Each row of a data frame or a matrix, and each value of a vector, refers to a sample or ”rationale group”. | ||||||||||||||||||||||||||||||
type |
a character string specifying the group statistics to compute.
Furthermore, a user specified type of chart, say | ||||||||||||||||||||||||||||||
sizes |
a value or a vector of values specifying the sample sizes associated with each group. For continuous data provided as data frame or matrix the sample sizes are obtained counting the non- | ||||||||||||||||||||||||||||||
center |
a value specifying the center of group statistics or the ”target” value of the process. | ||||||||||||||||||||||||||||||
std.dev |
a value or an available method specifying the within-group standard deviation(s) of the process. | ||||||||||||||||||||||||||||||
limits |
a two-values vector specifying control limits. | ||||||||||||||||||||||||||||||
data.name |
a string specifying the name of the variable which appears on the plots. If not provided is taken from the object given as data. | ||||||||||||||||||||||||||||||
labels |
a character vector of labels for each group. | ||||||||||||||||||||||||||||||
newdata |
a data frame, matrix or vector, as for the | ||||||||||||||||||||||||||||||
newsizes |
a vector as for the | ||||||||||||||||||||||||||||||
newdata.name |
a string specifying the name of the variable which appears on the plots. If not provided is taken from the object given as newdata. | ||||||||||||||||||||||||||||||
newlabels |
a character vector of labels for each new group defined in the argument | ||||||||||||||||||||||||||||||
nsigmas |
a numeric value specifying the number of sigmas to use for computing control limits. It is ignored when the | ||||||||||||||||||||||||||||||
confidence.level |
a numeric value between 0 and 1 specifying the confidence level of the computed probability limits. | ||||||||||||||||||||||||||||||
rules |
a function of rules to apply to the chart. By default, the | ||||||||||||||||||||||||||||||
plot |
logical. If | ||||||||||||||||||||||||||||||
add.stats |
a logical value indicating whether statistics and other information should be printed at the bottom of the chart. | ||||||||||||||||||||||||||||||
chart.all |
a logical value indicating whether both statistics for | ||||||||||||||||||||||||||||||
label.limits |
a character vector specifying the labels for control limits. | ||||||||||||||||||||||||||||||
title |
a string giving the label for the main title. | ||||||||||||||||||||||||||||||
xlab |
a string giving the label for the x-axis. | ||||||||||||||||||||||||||||||
ylab |
a string giving the label for the y-axis. | ||||||||||||||||||||||||||||||
ylim |
a numeric vector specifying the limits for the y-axis. | ||||||||||||||||||||||||||||||
axes.las |
numeric in {0,1,2,3} specifying the style of axis labels. See | ||||||||||||||||||||||||||||||
digits |
the number of significant digits to use. | ||||||||||||||||||||||||||||||
restore.par |
a logical value indicating whether the previous | ||||||||||||||||||||||||||||||
object |
an object of class 'qcc'. | ||||||||||||||||||||||||||||||
x |
an object of class 'qcc'. | ||||||||||||||||||||||||||||||
... |
additional arguments to be passed to the generic function. |
Value
Returns an object of class 'qcc'.
Note
For a nice blog post discussing the qcc package, in particular how to implement the Western Eletric Rules (WER), see http://blog.yhathq.com/posts/quality-control-in-r.html.
Author(s)
Luca Scrucca
References
Mason, R.L. and Young, J.C. (2002) Multivariate Statistical Process Control with Industrial Applications, SIAM.
Montgomery, D.C. (2005) Introduction to Statistical Quality Control, 5th ed. New York: John Wiley & Sons.
Ryan, T. P. (2000), Statistical Methods for Quality Improvement, 2nd ed. New York: John Wiley & Sons, Inc.
Scrucca, L. (2004). qcc: an R package for quality control charting and statistical process control. R News 4/1, 11-17.
Wetherill, G.B. and Brown, D.W. (1991) Statistical Process Control. New York: Chapman & Hall.
See Also
shewhart.rules, cusum, ewma, oc.curves, process.capability, qcc.groups
Examples
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 | ##
## Continuous data
##
data(pistonrings)
attach(pistonrings)
diameter <- qcc.groups(diameter, sample)
qcc(diameter[1:25,], type="xbar")
qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,])
q <- qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,], plot=FALSE)
plot(q, chart.all=FALSE)
qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,], nsigmas=2)
qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,], confidence.level=0.99)
qcc(diameter[1:25,], type="R")
qcc(diameter[1:25,], type="R", newdata=diameter[26:40,])
qcc(diameter[1:25,], type="S")
qcc(diameter[1:25,], type="S", newdata=diameter[26:40,])
# add warning limits at 2 std. deviations
q <- qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,], plot=FALSE)
(warn.limits <- limits.xbar(q$center, q$std.dev, q$sizes, 2))
plot(q, restore.par = FALSE)
abline(h = warn.limits, lty = 3, col = "chocolate")
# variable control limits
out <- c(9, 10, 30, 35, 45, 64, 65, 74, 75, 85, 99, 100)
diameter <- qcc.groups(pistonrings$diameter[-out], sample[-out])
qcc(diameter[1:25,], type="xbar")
qcc(diameter[1:25,], type="R")
qcc(diameter[1:25,], type="S")
qcc(diameter[1:25,], type="xbar", newdata=diameter[26:40,])
qcc(diameter[1:25,], type="R", newdata=diameter[26:40,])
qcc(diameter[1:25,], type="S", newdata=diameter[26:40,])
detach(pistonrings)
##
## Attribute data
##
data(orangejuice)
attach(orangejuice)
qcc(D[trial], sizes=size[trial], type="p")
# remove out-of-control points (see help(orangejuice) for the reasons)
inc <- setdiff(which(trial), c(15,23))
q1 <- qcc(D[inc], sizes=size[inc], type="p")
qcc(D[inc], sizes=size[inc], type="p", newdata=D[!trial], newsizes=size[!trial])
detach(orangejuice)
data(orangejuice2)
attach(orangejuice2)
names(D) <- sample
qcc(D[trial], sizes=size[trial], type="p")
q2 <- qcc(D[trial], sizes=size[trial], type="p", newdata=D[!trial], newsizes=size[!trial])
detach(orangejuice2)
# put on the same graph the two orange juice samples
oldpar <- par(no.readonly = TRUE)
par(mfrow=c(1,2), mar=c(5,5,3,0))
plot(q1, title="First samples", ylim=c(0,0.5), add.stats=FALSE, restore.par=FALSE)
par("mar"=c(5,0,3,3), yaxt="n")
plot(q2, title="Second samples", add.stats=FALSE, ylim=c(0,0.5))
par(oldpar)
data(circuit)
attach(circuit)
qcc(x[trial], sizes=size[trial], type="c")
# remove out-of-control points (see help(circuit) for the reasons)
inc <- setdiff(which(trial), c(6,20))
qcc(x[inc], sizes=size[inc], type="c", labels=inc)
qcc(x[inc], sizes=size[inc], type="c", labels=inc,
newdata=x[!trial], newsizes=size[!trial], newlabels=which(!trial))
qcc(x[inc], sizes=size[inc], type="u", labels=inc,
newdata=x[!trial], newsizes=size[!trial], newlabels=which(!trial))
detach(circuit)
data(pcmanufact)
attach(pcmanufact)
qcc(x, sizes=size, type="u")
detach(pcmanufact)
data(dyedcloth)
attach(dyedcloth)
qcc(x, sizes=size, type="u")
# standardized control chart
q <- qcc(x, sizes=size, type="u", plot=FALSE)
z <- (q$statistics - q$center)/sqrt(q$center/q$size)
plot(z, type="o", ylim=range(z,3,-3), pch=16)
abline(h=0, lty=2)
abline(h=c(-3,3), lty=2)
detach(dyedcloth)
##
## Continuous one-at-time data
##
# viscosity data (Montgomery, pag. 242)
x <- c(33.75, 33.05, 34, 33.81, 33.46, 34.02, 33.68, 33.27, 33.49, 33.20,
33.62, 33.00, 33.54, 33.12, 33.84)
qcc(x, type="xbar.one")
qcc(x, type="xbar.one", std.dev = "SD")
|
Example output
Package 'qcc' version 2.7
Type 'citation("qcc")' for citing this R package in publications.
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "xbar")
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00979
$ nsigmas : num 3
$ limits : num [1, 1:2] 74 74
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "xbar", newdata = diameter[26:40, ])
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00979
$ newstats : Named num [1:15] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1, 1:2] 74 74
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "xbar", newdata = diameter[26:40, ], nsigmas = 2)
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00979
$ newstats : Named num [1:15] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 2
$ limits : num [1, 1:2] 74 74
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "xbar", newdata = diameter[26:40, ], confidence.level = 0.99)
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00979
$ newstats : Named num [1:15] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name : chr "diameter[26:40, ]"
$ confidence.level: num 0.99
$ limits : num [1, 1:2] 74 74
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "R")
$ type : chr "R"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 0.038 0.019 0.036 0.022 0.026 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.0228
$ std.dev : num 0.00979
$ nsigmas : num 3
$ limits : num [1, 1:2] 0 0.0481
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "R", newdata = diameter[26:40, ])
$ type : chr "R"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 0.038 0.019 0.036 0.022 0.026 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.0228
$ std.dev : num 0.00979
$ newstats : Named num [1:15] 0.044 0.025 0.015 0.019 0.017 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1, 1:2] 0 0.0481
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "S")
$ type : chr "S"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 0.01477 0.0075 0.01475 0.00908 0.01222 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.00924
$ std.dev : num 0.00983
$ nsigmas : num 3
$ limits : num [1, 1:2] 0 0.0193
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "S", newdata = diameter[26:40, ])
$ type : chr "S"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 0.01477 0.0075 0.01475 0.00908 0.01222 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.00924
$ std.dev : num 0.00983
$ newstats : Named num [1:15] 0.01655 0.01033 0.00691 0.0075 0.00673 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1, 1:2] 0 0.0193
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
LCL UCL
73.99242 74.00993
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "xbar")
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00986
$ nsigmas : num 3
$ limits : num [1:25, 1:2] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "R")
$ type : chr "R"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 0.038 0.009 0.036 0.022 0.026 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.0223
$ std.dev : num 0.00986
$ nsigmas : num 3
$ limits : num [1:25, 1:2] 0 0 0 0 0 0 0 0 0 0 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = diameter[1:25, ], type = "S")
$ type : chr "S"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:25] 0.01477 0.00458 0.01475 0.00908 0.01222 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.00939
$ std.dev : num 0.00993
$ nsigmas : num 3
$ limits : num [1:25, 1:2] 0 0 0 0 0 0 0 0 0 0 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "xbar", newdata = diameter[26:40, ])
$ type : chr "xbar"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 74
$ std.dev : num 0.00986
$ newstats : Named num [1:15] 74 74 74 74 74 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1:40, 1:2] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "R", newdata = diameter[26:40, ])
$ type : chr "R"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 0.038 0.009 0.036 0.022 0.026 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.0223
$ std.dev : num 0.00986
$ newstats : Named num [1:15] 0.044 0.025 0.015 0.019 0.017 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1:40, 1:2] 0 0 0 0 0 0 0 0 0 0 ...
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = diameter[1:25, ], type = "S", newdata = diameter[26:40, ])
$ type : chr "S"
$ data.name : chr "diameter[1:25, ]"
$ data : num [1:25, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:25] 0.01477 0.00458 0.01475 0.00908 0.01222 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ sizes : Named int [1:25] 5 3 5 5 5 4 4 5 4 5 ...
..- attr(*, "names")= chr [1:25] "1" "2" "3" "4" ...
$ center : num 0.00939
$ std.dev : num 0.00993
$ newstats : Named num [1:15] 0.01655 0.01033 0.00691 0.0075 0.00673 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata : num [1:15, 1:5] 74 74 74 74 74 ...
..- attr(*, "dimnames")=List of 2
$ newsizes : Named int [1:15] 5 5 5 5 5 5 5 5 5 5 ...
..- attr(*, "names")= chr [1:15] "26" "27" "28" "29" ...
$ newdata.name: chr "diameter[26:40, ]"
$ nsigmas : num 3
$ limits : num [1:40, 1:2] 0 0 0 0 0 0 0 0 0 0 ...
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = D[trial], type = "p", sizes = size[trial])
$ type : chr "p"
$ data.name : chr "D[trial]"
$ data : int [1:30, 1] 12 15 8 10 4 7 16 9 14 10 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:30] 0.24 0.3 0.16 0.2 0.08 0.14 0.32 0.18 0.28 0.2 ...
..- attr(*, "names")= chr [1:30] "1" "2" "3" "4" ...
$ sizes : int [1:30] 50 50 50 50 50 50 50 50 50 50 ...
$ center : num 0.231
$ std.dev : num 0.422
$ nsigmas : num 3
$ limits : num [1:30, 1:2] 0.0524 0.0524 0.0524 0.0524 0.0524 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = D[inc], type = "p", sizes = size[inc], newdata = D[!trial], newsizes = size[!trial])
$ type : chr "p"
$ data.name : chr "D[inc]"
$ data : int [1:28, 1] 12 15 8 10 4 7 16 9 14 10 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:28] 0.24 0.3 0.16 0.2 0.08 0.14 0.32 0.18 0.28 0.2 ...
..- attr(*, "names")= chr [1:28] "1" "2" "3" "4" ...
$ sizes : int [1:28] 50 50 50 50 50 50 50 50 50 50 ...
$ center : num 0.215
$ std.dev : num 0.411
$ newstats : Named num [1:24] 0.18 0.12 0.24 0.1 0.12 0.08 0.12 0.06 0.14 0.12 ...
..- attr(*, "names")= chr [1:24] "29" "30" "31" "32" ...
$ newdata : int [1:24, 1] 9 6 12 5 6 4 6 3 7 6 ...
$ newsizes : int [1:24] 50 50 50 50 50 50 50 50 50 50 ...
$ newdata.name: chr "D[!trial]"
$ nsigmas : num 3
$ limits : num [1:52, 1:2] 0.0407 0.0407 0.0407 0.0407 0.0407 ...
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = D[trial], type = "p", sizes = size[trial])
$ type : chr "p"
$ data.name : chr "D[trial]"
$ data : int [1:24, 1] 9 6 12 5 6 4 6 3 7 6 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:24] 0.18 0.12 0.24 0.1 0.12 0.08 0.12 0.06 0.14 0.12 ...
..- attr(*, "names")= chr [1:24] "31" "32" "33" "34" ...
$ sizes : int [1:24] 50 50 50 50 50 50 50 50 50 50 ...
$ center : num 0.111
$ std.dev : num 0.314
$ nsigmas : num 3
$ limits : num [1:24, 1:2] 0 0 0 0 0 0 0 0 0 0 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x[trial], type = "c", sizes = size[trial])
$ type : chr "c"
$ data.name : chr "x[trial]"
$ data : int [1:26, 1] 21 24 16 12 15 5 28 20 31 25 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named int [1:26] 21 24 16 12 15 5 28 20 31 25 ...
..- attr(*, "names")= chr [1:26] "1" "2" "3" "4" ...
$ sizes : int [1:26] 100 100 100 100 100 100 100 100 100 100 ...
$ center : num 19.8
$ std.dev : num 4.45
$ nsigmas : num 3
$ limits : num [1, 1:2] 6.48 33.21
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x[inc], type = "c", sizes = size[inc], labels = inc)
$ type : chr "c"
$ data.name : chr "x[inc]"
$ data : int [1:24, 1] 21 24 16 12 15 28 20 31 25 20 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named int [1:24] 21 24 16 12 15 28 20 31 25 20 ...
..- attr(*, "names")= chr [1:24] "1" "2" "3" "4" ...
$ sizes : int [1:24] 100 100 100 100 100 100 100 100 100 100 ...
$ center : num 19.7
$ std.dev : num 4.43
$ nsigmas : num 3
$ limits : num [1, 1:2] 6.36 32.97
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = x[inc], type = "c", sizes = size[inc], labels = inc, newdata = x[!trial], newsizes = size[!trial]| __truncated__
$ type : chr "c"
$ data.name : chr "x[inc]"
$ data : int [1:24, 1] 21 24 16 12 15 28 20 31 25 20 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named int [1:24] 21 24 16 12 15 28 20 31 25 20 ...
..- attr(*, "names")= chr [1:24] "1" "2" "3" "4" ...
$ sizes : int [1:24] 100 100 100 100 100 100 100 100 100 100 ...
$ center : num 19.7
$ std.dev : num 4.43
$ newstats : Named int [1:20] 16 18 12 15 24 21 28 20 25 19 ...
..- attr(*, "names")= chr [1:20] "27" "28" "29" "30" ...
$ newdata : int [1:20, 1] 16 18 12 15 24 21 28 20 25 19 ...
$ newsizes : int [1:20] 100 100 100 100 100 100 100 100 100 100 ...
$ newdata.name: chr "x[!trial]"
$ nsigmas : num 3
$ limits : num [1, 1:2] 6.36 32.97
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 15
$ call : language qcc(data = x[inc], type = "u", sizes = size[inc], labels = inc, newdata = x[!trial], newsizes = size[!trial]| __truncated__
$ type : chr "u"
$ data.name : chr "x[inc]"
$ data : int [1:24, 1] 21 24 16 12 15 28 20 31 25 20 ...
..- attr(*, "dimnames")=List of 2
$ statistics : Named num [1:24] 0.21 0.24 0.16 0.12 0.15 0.28 0.2 0.31 0.25 0.2 ...
..- attr(*, "names")= chr [1:24] "1" "2" "3" "4" ...
$ sizes : int [1:24] 100 100 100 100 100 100 100 100 100 100 ...
$ center : num 0.197
$ std.dev : num 0.443
$ newstats : Named num [1:20] 0.16 0.18 0.12 0.15 0.24 0.21 0.28 0.2 0.25 0.19 ...
..- attr(*, "names")= chr [1:20] "27" "28" "29" "30" ...
$ newdata : int [1:20, 1] 16 18 12 15 24 21 28 20 25 19 ...
$ newsizes : int [1:20] 100 100 100 100 100 100 100 100 100 100 ...
$ newdata.name: chr "x[!trial]"
$ nsigmas : num 3
$ limits : num [1, 1:2] 0.0636 0.3297
..- attr(*, "dimnames")=List of 2
$ violations :List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x, type = "u", sizes = size)
$ type : chr "u"
$ data.name : chr "x"
$ data : int [1:20, 1] 10 12 8 14 10 16 11 7 10 15 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:20] 2 2.4 1.6 2.8 2 3.2 2.2 1.4 2 3 ...
..- attr(*, "names")= chr [1:20] "1" "2" "3" "4" ...
$ sizes : int [1:20] 5 5 5 5 5 5 5 5 5 5 ...
$ center : num 1.93
$ std.dev : num 1.39
$ nsigmas : num 3
$ limits : num [1, 1:2] 0.0661 3.7939
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x, type = "u", sizes = size)
$ type : chr "u"
$ data.name : chr "x"
$ data : int [1:10, 1] 14 12 20 11 7 10 21 16 19 23
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:10] 1.4 1.5 1.538 1.1 0.737 ...
..- attr(*, "names")= chr [1:10] "1" "2" "3" "4" ...
$ sizes : num [1:10] 10 8 13 10 9.5 10 12 10.5 12 12.5
$ center : num 1.42
$ std.dev : num 1.19
$ nsigmas : num 3
$ limits : num [1:10, 1:2] 0.291 0.158 0.431 0.291 0.262 ...
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x, type = "xbar.one")
$ type : chr "xbar.one"
$ data.name : chr "x"
$ data : num [1:15, 1] 33.8 33 34 33.8 33.5 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:15] 33.8 33 34 33.8 33.5 ...
..- attr(*, "names")= chr [1:15] "1" "2" "3" "4" ...
$ sizes : int [1:15] 1 1 1 1 1 1 1 1 1 1 ...
$ center : num 33.5
$ std.dev : num 0.426
$ nsigmas : num 3
$ limits : num [1, 1:2] 32.2 34.8
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
List of 11
$ call : language qcc(data = x, type = "xbar.one", std.dev = "SD")
$ type : chr "xbar.one"
$ data.name : chr "x"
$ data : num [1:15, 1] 33.8 33 34 33.8 33.5 ...
..- attr(*, "dimnames")=List of 2
$ statistics: Named num [1:15] 33.8 33 34 33.8 33.5 ...
..- attr(*, "names")= chr [1:15] "1" "2" "3" "4" ...
$ sizes : int [1:15] 1 1 1 1 1 1 1 1 1 1 ...
$ center : num 33.5
$ std.dev : num 0.342
$ nsigmas : num 3
$ limits : num [1, 1:2] 32.5 34.5
..- attr(*, "dimnames")=List of 2
$ violations:List of 2
- attr(*, "class")= chr "qcc"
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.