README.md

In Kuan-cheng-da/EATME: EWMA-p control charts with correction of measurement error

EATME

Overview

This package aims to address measurement error effects when constructing EWMA-p control charts. The method primarily focuses on binary random variables, but it can be applied to any continuous random variables by using sign statistic to transform them to discrete ones. With the correction of measurement error effects, we can obtain the corrected control limits of EWMA-p control chart and reasonably adjusted EWMA-p control charts. Two folders “R” and “man” contain the R code and the manual for the package EATME. The folder named “Real_data_analysis” contains IC and OC real datasets, additional R code for our estimation method (Real_data_analysis.R) and implementation of the existing R package (implement_qcc_rda.R), and two figures created by Real_data_analysis.R for the comparison.

Installation

You can install the released version of EATME from github

install.packages('EATME')
library(EATME)

Examples

This is an example which shows how find the control limits for two sided EWMA-p chart:

library(EATME)
#0.2 is the in-control probability
#0.05 is the smoothing parameter for EWMA chart
#5 is sample size
#1 is the correctness of true and false
#200 is ARL0
#10000 is the number of simulation for Monte Carlo method
#1 is the maximum error what we can allow between hat_ARL0 and ARL0
EWMA_p_two(0.2,0.05,5,1,1,200,10000,1)
#> $L1
#> [1] 2.320362
#> 
#> $L2
#> [1] 2.219673
#> 
#> $hat_ARL0
#> [1] 199.6776
#> 
#> $hat_MRL
#> [1] 131
#> 
#> $hat_SDRL
#> [1] 219.8172
#> 
#> $UCL
#> [1] 0.2664658
#> 
#> $LCL
#> [1] 0.1364184

This is an example which shows how plot the two sided adjusted EWMA-p chart :

set.seed(6)
#We generate the discrete IC data with measurement error.
IC = ME_data_generate(0.4,50,50,0.9,0.9)
#We generate the discrete IC data with measurement error.
OC = ME_data_generate(0.5,50,50,0.9,0.9)
#The real EWMA-p chart
EWMA_p_chart_two(IC$real_data,OC$real_data,0.05,50,1,1)

#The unadjusted EWMA-p chart
EWMA_p_chart_two(IC$obs_data,OC$obs_data,0.05,50,1,1)

#The adjusted EWMA-p chart
EWMA_p_chart_two(IC$obs_data,OC$obs_data,0.05,50,0.9,0.9)

Kuan-cheng-da/EATME documentation built on May 18, 2024, 3:06 a.m.