u-chart

Description

This function builds a u-chart for the Poisson-based count data statistic.

Usage

cchart.u(x1 = NULL, n1 = NULL, type = "norm", u1 = NULL, x2 = NULL, 
n2 = NULL, lambda = NULL, u2 = NULL)

Arguments

Details

For a phase I u-chart, n1 must be specified and either x1 or u1. For a phase II u-chart, n2 must be specified, plus x2 or u2 and either phat, x1 and n1, or u1 and n1. It is important to note that the normal approximation used in the Shewhart u-chart is valid only for n*u large. For small n*p , it should be used an "improved u chart" (with non-normal correction) given by using the argument "CF".

Value

Returns a u-chart.

Author(s)

Daniela R. Recchia, Emanuel P. Barbosa

Examples

data(moonroof)
attach(moonroof)
cchart.u(x1 = yi[1:17], n1 = ni[1:17])
cchart.u(x1 = yi[1:17], n1 = ni[1:17], type = "CF", x2 = yi[18:34], n2 = ni[18:34])
cchart.u(type = "std", u2 = ui[18:34], n2 = ni[18:34], lambda = 1.4)

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(IQCC)
Loading required package: qcc
Package 'qcc', version 2.6
Type 'citation("qcc")' for citing this R package in publications.
Loading required package: MASS
Loading required package: micEcon

If you have questions, suggestions, or comments regarding one of the 'micEcon' packages, please use a forum or 'tracker' at micEcon's R-Forge site:
https://r-forge.r-project.org/projects/micecon/
Loading required package: miscTools
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/IQCC/cchart.u.Rd_%03d_medium.png", width=480, height=480)
> ### Name: cchart.u
> ### Title: u-chart
> ### Aliases: cchart.u
> 
> ### ** Examples
> 
> data(moonroof)
> attach(moonroof)
> cchart.u(x1 = yi[1:17], n1 = ni[1:17])
Warning message:
In qcc(x1, type = "u", n1, limits = c(l, u), center = lambda, title = "Shewhart u-chart (phase I)") :
  'std.dev' is not used when limits is given
> cchart.u(x1 = yi[1:17], n1 = ni[1:17], type = "CF", x2 = yi[18:34], n2 = ni[18:34])
Warning message:
In qcc(x2, type = "u", n2, limits = c(l, u), center = lambda, title = "Cornish-Fisher u-exact (phase II)") :
  'std.dev' is not used when limits is given
> cchart.u(type = "std", u2 = ui[18:34], n2 = ni[18:34], lambda = 1.4)
List of 11
 $ call      : language qcc(data = std, type = "u", sizes = n2, center = 0, limits = c(-3, 3),      title = "Stardardized u-chart (phase II)")
 $ type      : chr "u"
 $ data.name : chr "std"
 $ data      : num [1:17, 1] -27.04 7.53 -18.39 8.72 -4.87 ...
  ..- attr(*, "dimnames")=List of 2
 $ statistics: Named num [1:17] -1.229 0.538 -1.149 0.396 -0.304 ...
  ..- attr(*, "names")= chr [1:17] "1" "2" "3" "4" ...
 $ sizes     : int [1:17] 22 14 16 22 16 14 5 13 19 10 ...
 $ center    : num 0
 $ std.dev   : num NaN
 $ nsigmas   : num 3
 $ limits    : num [1, 1:2] -3 3
  ..- attr(*, "dimnames")=List of 2
 $ violations:List of 2
 - attr(*, "class")= chr "qcc"
Warning messages:
1: In sqrt(sum(sizes * data)/sum(sizes)) : NaNs produced
2: In qcc(std, type = "u", n2, center = 0, limits = c(-3, 3), title = "Stardardized u-chart (phase II)") :
  'std.dev' is not used when limits is given
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>