Title: | Maxwell Control Charts |
---|---|
Description: | Computes Control limits, coefficients of control limits, various performance metrics and depicts control charts for monitoring Maxwell-distributed quality characteristics. |
Authors: | Zahid Khan [aut], Zsolt T. Kosztyan [aut, cre] |
Maintainer: | Zsolt T. Kosztyan <[email protected]> |
License: | GPL (>= 2) |
Version: | 0.0.1 |
Built: | 2024-10-31 05:50:59 UTC |
Source: | https://github.com/kzst/mxcc |
This dataset contains the failure times (in hours) of a vertical boring machine, used to illustrate the control chart for monitoring the Maxwell distribution parameter.The data was originally reported by Krishna and Malik (2012).
data("failure_time")
data("failure_time")
A data frame consisted of 8 sample batches each wiht 4 observations.
The failure times in this dataset are organized into 8 subgroups, each containing 4 observations. The failure times are measured in hours. These data are used to construct control charts for monitoring the scale parameter of the Maxwell distribution.
Krishna, H.and Malik M. Reliability estimation in Maxwell distribution with progressively Type-II censored data.Journal of Statistical Computation and Simulation. 2012;82(4):pp. 623–641. doi:10.1080/00949655.2010.550291
Hossain, M.P., Omar, M.H. and Riaz, M., 2017. New V control chart for the Maxwell distribution. Journal of Statistical Computation and Simulation, 87(3), pp.594-606. doi:10.1080/00949655.2016.1222391
data("failure_time") failure_time
data("failure_time") failure_time
It calculates the Average Run Length (ARL) for either the V or VSQ control charts, based on the specified sample size, shift constant, and false alarm probability. The user can choose between the two types of control charts.
mxarl(n_val = 1, delta = seq(1, 3, length.out = 100), alpha = 0.0027, type = "V")
mxarl(n_val = 1, delta = seq(1, 3, length.out = 100), alpha = 0.0027, type = "V")
n_val |
A numeric vector specifying the sample sizes. Default is |
delta |
A numeric vector specifying the shift constants for the control chart. Default is |
alpha |
A numeric value specifying the significance level (false alarm probability). Default is |
type |
A character string specifying the type of chart to be used. Can be either |
This function computes the Average Run Length (ARL) for both V and VSQ control charts by calculating the power and ARL values for the given sample sizes and shift constants
A plot showing the ARL curves for the specified sample sizes and shift constants
Zahid Khan
Hossain, M.P., Omar, M.H. and Riaz, M., 2017. New V control chart for the Maxwell distribution. Journal of Statistical Computation and Simulation, 87(3), pp.594-606. doi:10.1080/00949655.2016.1222391
Shah, F., Khan, Z., Aslam, M. and Kadry, S., 2021. Statistical Development of the VSQ‐Control Chart for Extreme Data with an Application to the Carbon Fiber Industry. Mathematical Problems in Engineering, 2021(1), p.9766986. doi:10.1155/2021/9766986
mxarl(n_val = c(5,10), delta = seq(1, 3, length.out = 100), alpha = 0.0027, type = "V")
mxarl(n_val = c(5,10), delta = seq(1, 3, length.out = 100), alpha = 0.0027, type = "V")
The function mxk
calculates the coefficients for V and VSQ control charts used to monitor the scale parameter of Maxwell-distributed quality characteristics. It computes two coefficients based on the chosen chart type. For the V chart, the function returns L1
and L2
, while for the VSQ chart, it returns P1
and P2
.
mxk(n = 1, alpha = 0.0027, type = "V")
mxk(n = 1, alpha = 0.0027, type = "V")
n |
Sample size. The number of observations in each subgroup (numeric). Defaults to |
alpha |
Probability of false alarm (Type I error). Defaults to |
type |
The type of chart. Accepts either |
The mxk
function calculates the coefficients used in V and VSQ control charts for monitoring the scale parameter of Maxwell-distributed data. The user must specify the sample size n
, the probability of a false alarm alpha
, and the type of chart ("V"
or "VSQ"
). By default, n
is set to 1, and alpha
is set to 0.0027.
Depending on the type
of chart:
L1 |
Coefficient L1 for the V chart. |
L2 |
Coefficient L2 for the V chart. |
P1 |
Coefficient P1 for the VSQ chart. |
P2 |
Coefficient P2 for the VSQ chart. |
Zahid Khan
Hossain, M.P., Omar, M.H. and Riaz, M., 2017. New V control chart for the Maxwell distribution. Journal of Statistical Computation and Simulation, 87(3), pp.594-606. doi:10.1080/00949655.2016.1222391
Shah, F., Khan, Z., Aslam, M. and Kadry, S., 2021. Statistical Development of the VSQ‐Control Chart for Extreme Data with an Application to the Carbon Fiber Industry. Mathematical Problems in Engineering, 2021(1), p.9766986. doi:10.1155/2021/9766986
mxk(n = 4, alpha = 0.0027, type = "VSQ")
mxk(n = 4, alpha = 0.0027, type = "VSQ")
This function calculates the L-sigma multiplier (L) for control chart based on the specified type: V chart or VSQ chart. The calculation is based on the sample size and the false alarm probability.
mxm(n = 1, alpha = 0.0027, type = "V")
mxm(n = 1, alpha = 0.0027, type = "V")
n |
Sample size used in the chart. Default is 1. |
alpha |
False alarm probability. Default is 0.0027. |
type |
The type of control chart. Can be "V" for V chart or "VSQ" for VSQ chart. Default is "V". |
The function computes the L-sigma multiplier for either V chart or VSQ chart based on the specified type. If the type is "V", it uses the gamma distribution with the specified parameters. If the type is "VSQ", it uses the ch distribution. The output is the L value that represent multiplier in the control limits for the respective chart.
Returns the L-sigma multiplier (L) as a numeric value
Zahid Khan
Hossain, M.P., Omar, M.H. and Riaz, M., 2017. New V control chart for the Maxwell distribution. Journal of Statistical Computation and Simulation, 87(3), pp.594-606. doi:10.1080/00949655.2016.1222391
Shah, F., Khan, Z., Aslam, M. and Kadry, S., 2021. Statistical Development of the VSQ‐Control Chart for Extreme Data with an Application to the Carbon Fiber Industry. Mathematical Problems in Engineering, 2021(1), p.9766986. doi:10.1155/2021/9766986
mxm(n = 5, alpha = 0.005, type = "V")
mxm(n = 5, alpha = 0.005, type = "V")
The mxp
function calculates the power of V chart and VSQ control chart for monitoring the Maxwell scale parameter. It computes the probability of detecting a shift in the process, depending on the specified sample size, significance level, and the shift magnitude
mxp(n = 1, alpha = 0.0027, delta = 1, type)
mxp(n = 1, alpha = 0.0027, delta = 1, type)
n |
The sample size for each subgroup (Integer). Default is |
alpha |
Probability of false alarm (type I error) for the control chart(numeric). Default is |
delta |
The shift constant representing the magnitude of the shift to detect(numeric). Default is |
type |
Specifies the type of control chart to be used. Options are |
The function calculates the power of a control chart based on the provided sample size (n
), false alarm probability (alpha
), and shift constant (delta
). The chart type, either "V"
or "VSQ"
, determines which chart is used for the calculations. Power is a critical metric that evaluates the sensitivity of the control chart to detecting process shifts, allowing users to monitor for deviations from the expected process behavior.
A numeric value representing the power of the control chart to detect the process shift.
Zahid Khan
Hossain, M.P., Omar, M.H. and Riaz, M., 2017. New V control chart for the Maxwell distribution. Journal of Statistical Computation and Simulation, 87(3), pp.594-606. doi:10.1080/00949655.2016.1222391
Shah, F., Khan, Z., Aslam, M. and Kadry, S., 2021. Statistical Development of the VSQ‐Control Chart for Extreme Data with an Application to the Carbon Fiber Industry. Mathematical Problems in Engineering, 2021(1), p.9766986. doi:10.1155/2021/9766986
result <- mxp(n = 5, alpha = 0.0027, delta = 2, type = "V") print(result)
result <- mxp(n = 5, alpha = 0.0027, delta = 2, type = "V") print(result)
This function computes the power curves for V or VSQ control charts based on the Maxwell distribution. It allows the user to analyze the performance of these charts under different sample sizes and shifts in the process parameter.
mxpw(n_val = 1, delta = seq(1, 3, length.out = 100), alpha = 0.0027, type = "V")
mxpw(n_val = 1, delta = seq(1, 3, length.out = 100), alpha = 0.0027, type = "V")
n_val |
A numeric vector specifying the sample sizes. Default is |
delta |
A numeric vector specifying the shift constants for the control chart. Default is |
alpha |
A numeric value specifying the significance level (false alarm probability). Default is |
type |
A character string specifying the type of chart to be used. Can be either |
This function calculates the power curves for either the V or VSQ control charts, depending on the specified type
parameter. It computes the power values for different sample sizes and shift constants. The function uses the Gamma and Chi distribution functions for the V and VSQ charts, respectively.
The function returns a plot of power curves for the specified control chart type
Zahid Khan
Hossain, M.P., Omar, M.H. and Riaz, M., 2017. New V control chart for the Maxwell distribution. Journal of Statistical Computation and Simulation, 87(3), pp.594-606. doi:10.1080/00949655.2016.1222391
Shah, F., Khan, Z., Aslam, M. and Kadry, S., 2021. Statistical Development of the VSQ‐Control Chart for Extreme Data with an Application to the Carbon Fiber Industry. Mathematical Problems in Engineering, 2021(1), p.9766986. doi:10.1155/2021/9766986
mxpw(n = c(5,8), alpha = 0.0027, type = "VSQ")
mxpw(n = c(5,8), alpha = 0.0027, type = "VSQ")
The mxrl
function computes key characteristics of the run length distribution for V and VSQ control charts. It calculates the Average Run Length (ARL), Standard Deviation of the Run Length (SDRL), and Median Run Length (MRL), based on the provided sample size, significance level, shift constant, and control chart type.
mxrl(n = 1, alpha = 0.0027, delta = 1, type = "V")
mxrl(n = 1, alpha = 0.0027, delta = 1, type = "V")
n |
The sample size for each subgroup (Integer). Default is |
alpha |
Probability of false alarm (type I error) for the control chart(numeric). Default is |
delta |
The shift constant representing the magnitude of the shift to detect(numeric). Default is |
type |
Specifies the type of control chart to be used. Options are |
This function computes the characteristics of the run length distribution for either the V chart or the VSQ chart. The run length distribution is an essential metric in control chart analysis as it quantifies the performance of the control chart in detecting shifts in the process. The ARL is the expected number of samples before an out-of-control signal, SDRL is the standard deviation of the run length, and MRL is the median run length.
A list with the following components:
ARL |
The Average Run Length (ARL). |
SDRL |
The Standard Deviation of the Run Length (SDRL). |
MRL |
The Median Run Length (MRL). |
Zahid Khan
Hossain, M.P., Omar, M.H. and Riaz, M., 2017. New V control chart for the Maxwell distribution. Journal of Statistical Computation and Simulation, 87(3), pp.594-606. doi:10.1080/00949655.2016.1222391
Shah, F., Khan, Z., Aslam, M. and Kadry, S., 2021. Statistical Development of the VSQ‐Control Chart for Extreme Data with an Application to the Carbon Fiber Industry. Mathematical Problems in Engineering, 2021(1), p.9766986. doi:10.1155/2021/9766986
mxrl(n = 2, alpha = 0.005, delta = 1, type = "VSQ")
mxrl(n = 2, alpha = 0.005, delta = 1, type = "VSQ")
The mxrpc
function constructs control limits for the V and VSQ control charts using real data. It allows the user to specify the value of alpha
, the type of control limit (Probability Limit Control Chart or K-Sigma Control Chart), and the type of control chart (V or VSQ). The function provides a summary of control chart parameters and plots the corresponding control chart.
mxrpc(data, alpha = 0.0027, limit = "PCL", chart = "V", summary = FALSE)
mxrpc(data, alpha = 0.0027, limit = "PCL", chart = "V", summary = FALSE)
data |
A data frame containing real-world observations for which the control charts will be constructed. |
alpha |
The false alarm probability for control limit calculation. Default is |
limit |
The type of control limit to be used: either "PCL" for Probability Limit Control Chart or "KCL" for K-Sigma Control Chart. Default is |
chart |
The type of control chart to construct: either "V" for V chart or "VSQ" for VSQ chart. Default is |
summary |
Logical value indicating whether to display a summary of control chart parameters. Default is |
This function takes a real data set and generates control charts (V or VSQ) based on the specified control limit type (PCL or KCL). The function outputs both a graphical representation of the control chart and a summary of the control chart parameters, including control limits, central line, and the values used for constructing the chart.
The function outputs a graph of the selected control chart (V or VSQ), and if summary = TRUE
, it prints a summary of the control chart parameters including the lower control limit (LCL), center line (CL), and upper control limit (UCL).
Zahid Khan
Hossain, M.P., Omar, M.H. and Riaz, M., 2017. New V control chart for the Maxwell distribution. Journal of Statistical Computation and Simulation, 87(3), pp.594-606. doi:10.1080/00949655.2016.1222391
Shah, F., Khan, Z., Aslam, M. and Kadry, S., 2021. Statistical Development of the VSQ‐Control Chart for Extreme Data with an Application to the Carbon Fiber Industry. Mathematical Problems in Engineering, 2021(1), p.9766986. doi:10.1155/2021/9766986
mxrpc(data=strength_data, alpha = 0.0027, limit = "PCL", chart = "VSQ", summary = TRUE)
mxrpc(data=strength_data, alpha = 0.0027, limit = "PCL", chart = "VSQ", summary = TRUE)
This function constructs control limits for the V and VSQ control charts based on probability limits or k-sigma limits using simulated data from the Maxwell distribution. The function allows for flexible configuration of control chart types and limit methods.
mxspc(m = 25, n = 4, alpha = 0.0027, sigma, seed = 100, limit = "PCL", chart = "V", summary = FALSE)
mxspc(m = 25, n = 4, alpha = 0.0027, sigma, seed = 100, limit = "PCL", chart = "V", summary = FALSE)
m |
The number of subgroups or samples. Default is |
n |
The size of each sample or subgroup. Default is |
alpha |
The false alarm probability for control limit calculation. Default is |
sigma |
The scale parameter of the Maxwell distribution, which must be provided by the user. |
seed |
The seed for random number generation to ensure reproducibility. Default is |
limit |
The type of control limit to be used: either "PCL" for probability limit control chart or "KLC" for k-sigma limit control chart. Default is |
chart |
The type of control chart to construct: either "V" for V chart or "VSQ" for VSQ chart. Default is |
summary |
Logical value indicating whether to display a summary of control chart parameters. Default is |
The function simulates data from the Maxwell distribution using the provided scale parameter (sigma
) and calculates control limits and statistics for the specified control chart type (V or VSQ). It allows for choosing between probability limit control charts and k-sigma control charts. A control chart is plotted based on the calculated limits and plotting statistics.
A list of control chart parameters is returned invisibly, which includes:
LCL |
The lower control limit (or probability limit). |
CL |
The center line of the control chart. |
UCL |
The upper control limit (or probability limit). |
m |
The number of subgroups. |
n |
The sample size for each subgroup. |
sigma |
The provided scale parameter for the Maxwell distribution. |
seed |
The seed used for reproducibility. |
limit |
The type of limit used ("PCL" or "KLC"). |
chart |
The type of control chart ("V" or "VSQ"). |
Zahid Khan
Hossain, M.P., Omar, M.H. and Riaz, M., 2017. New V control chart for the Maxwell distribution. Journal of Statistical Computation and Simulation, 87(3), pp.594-606. doi:10.1080/00949655.2016.1222391
Shah, F., Khan, Z., Aslam, M. and Kadry, S., 2021. Statistical Development of the VSQ‐Control Chart for Extreme Data with an Application to the Carbon Fiber Industry. Mathematical Problems in Engineering, 2021(1), p.9766986. doi:10.1155/2021/9766986
mxspc(m = 30, n = 4, alpha = 0.0027, sigma = 1777.86, seed = 100, limit="PCL", summary = TRUE, chart = "V")
mxspc(m = 30, n = 4, alpha = 0.0027, sigma = 1777.86, seed = 100, limit="PCL", summary = TRUE, chart = "V")
This dataset contains the strength measurements of carbon fiber tested under tension at various gauge levels. The data is used to construct control charts for monitoring the scale parameter of the Maxwell distribution in the carbon fiber industry.
data("strength_data")
data("strength_data")
A data frame with 12 subgroups each with 5 observations.
he dataset consists of 12 subgroups, each containing 5 measurements of carbon fiber strength. These measurements are used to compute the V-statistic, which is then applied to control charts for monitoring the Maxwell distribution's scale parameter. The data was originally reported by Badar and Priest (1982) and is slightly modified for statistical analysis.
Badar, M. G., & Priest, A. M. (1982). Statistical aspects of fiber and bundle strength in hybrid composites. In "Progress in Science and Engineering Composites ICCM-IV, Tokyo", pp. 1129-1136.
Shah, F., Khan, Z., Aslam, M. and Kadry, S., 2021. Statistical Development of the VSQ‐Control Chart for Extreme Data with an Application to the Carbon Fiber Industry. Mathematical Problems in Engineering, 2021(1), p.9766986. doi:10.1155/2021/9766986
data("strength_data") strength_data
data("strength_data") strength_data