Description Details Author(s) References Examples
Permite realizar los graficos de control multivariante T2 de Hotelling, de doble dimension y de dimension variable. Ademas calcula los ARL bajo y fuera de control.
Package: | Gcmulti |
Type: | Package |
Version: | 1.0 |
Date: | 2016-08-29 |
License: | GPL (>= 2) |
Jhonny Barahona
Omar Ruiz
Maintainer: Jhonny Barahona <jfbaraho@espol.edu.ec>
[1] Aparisi, F.; Epprecht, E.; Ruiz, O.; Veiga A., "Reducing Sampling Costs of Multivariate SPC with a Double-Dimension T^2 Control Chart",” International Journal of Production Research, vol. 1, no. 1, pag. 1-15, 2013.
[2] Aparisi, F.; Epprecht, E.; Ruiz, O., "T2 Control Charts with Variable Dimension", Journal of Quality Technology, vol. 44, no. 4, pp. 375-393, 2012.
[3] Ruiz, O., “Gráficos de control de calidad multivariantes con dimensión variable” Tesis doctoral, Dept. de Estadística e Investigación Operativa, Univ. Pol. de Valencia, España, 2013.
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 | ###Simulaci?n de datos normales multivariados###
mu=c(4.5,7,8.45)
sigma=matrix(c(2,1.5,2.4,1.5,3,3.1,2.4,3.1,4),ncol=3)
#la matriz de covarianzas debe ser definida positiva
#comprobaci?n por medio de los valores propios de la matriz (>0)
eigen(sigma)
library(MASS)
#datos multivariantes
dmulti=mvrnorm(250,mu,sigma)
colnames(dmulti)=c("va1","va2","va3")
View(dmulti)
#Datos a utilizar
#n1----->50 dmulti[1:50,]
#n2----->100 dmulti[1:100,]
#n5----->250 dmulti
dmulti_n1=dmulti[1:50,]
dmulti_n2=dmulti[1:100,]
#vector de medias bajo control
mu0=c(5.4,6.8,8.5)
#uso de la funci?n T2 de Hotelling
#para p1
grafica_T2(mu0[1:2],sigma[1:2,1:2],dmulti_n1[,1:2],k=1,LCS=11.983,d=1) #n=1
#para p
grafica_T2(mu0,sigma,dmulti_n1,k=1,LCS=14.321,d=1.2)#n=1
#uso de la funci?n T2 con dimensi?n variable
grafica_T2DV(p1=2,p=3,w=3.83,CLp1=17.09,CLp=10.62,disp1=1,disp=1.2,mu=mu0,sigma=sigma,X=dmulti_n1,k=1)
#Uso de la funci?n T2 con doble dimensi?n
grafica_T2DD(p1=2,p=3,w=2.89,CLP1=14.07,CLP=14.11,disp1=1,disp=1.2,mu=mu0,sigma=sigma,X=dmulti_n1,k=1)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.