Abstract
In this paper we use a beta Chart to monitor fracture data. The beta Chart displays its control limits based on the beta probability distribution. This Chart was applied to a data set of contaminated peanut proportions. With toxic substances for 34 batches weighing 120 pounds and then comparing it with the traditional Shewart Chart ( p- Chart ), then a sensitivity study is performed to compare both Charts in two cases: under control and out of control. Using several values of average ratios and with different sample sizes, the evaluation is based on one of the criteria that measures the efficiency of the the evaluation is done based on one of the criteria that measures the efficiency of the Chart
which is the Average Run Length (ARL) for both cases. , and the operating average in the first case is a function of the type one error is for comparing Charts and in order to discover the shift in the proof the first type and in the second case it is a function of the type two error The second cases Sensitivity analysis using several values of the fracture rate confirmed the superior performance of the beta Chart compared to the p Chart, resulting in the proposed approximation the value of the average operating length of ARLo in the control condition slightly greater and that the value of the average operating length in an out sidesce nario ARL1 controlismuchsmaller
Keywords
Main Subjects
Highlights
When the data for the percentage of contaminated peanuts were applied to both the beta and Shewhart panels (P), it was found that the beta panel produced more accurate results than the Shewhart panel. Sensitivity analysis using several fractional values confirmed the outstanding performance of the beta panel compared to the (P) Schwart panel and the proposed approximation resulted in a slightly larger ARL0, and a much smaller ARL1.
Full Text
1 -مقدمة Introduction
تؤدی لوحات السیطرة النوعیة دوراً مهما فی عملیة المراقبة الاحصائیة ، إن عملیات السیطرة الإحصائیة )((Statistical Process Control (SPC عادة ماتستخدم فی مراقبة واکتشاف التحولات فی عملیات الإنتاج ، وتعد لوحات السیطرة للصفات ( Attribte ) أدوات مهمة موجودة فی عملیات السیطرة الإحصائیة ( (SPC لمراقبة العملیات باستخدام البیانات المتقطعة(discrete data) . وکانت لوحتا p و np أکثر شیوعًا لرصد الوحدات غیر المطابقة ، والتی طورها Shewhart فی عام 1924. ثم طورت لتکون أکثر حساسیة للکشف عن الانحرافات فی العملیة الانتاجیة . یتم حساب تقدیرات المتوسط والتباین للبیانات بافتراض ان التوزیع الاحتمالی لها هو توزیع ثنائی الحدین وبمعلمتی n و p لعدد الوحدات غیر المطابقة ، ویتم حساب حدود السیطرة بناءً على التقریب للتوزیع الطبیعی (Normal distribution approximation) للبیانات ان المتغیرات من النوع الکسریة (Fraction type) هی مشاهدات مُعَبَّر عنها فی الفترة [0،1] . وفی هذه البحث یقع الکسر ضمن الفترة [0،1] التی تمثل: (1) نسبة مئویة - النسبة بین رقمین من النوع المتقطع ( فمثلا ، عدد الأجزاء المعیبة فی دفعة إنتاج ما) ، و (2) النسبة - النسبة بین رقمین مستمرین (مثلا حجم الکحول فی النبیذ) .وهناک بعض القواعد التی تتعامل مع افتراض التماثل والتقریب للتوزیع الطبیعی حیث فی سنة (1989) أقترح کل من (Schader و Schmid) أن التقریب الطبیعی لتوزیع ذی الحدین یکون مناسباً اذا تحقق الشرطان الآتیان :
np (1- p) ≥ 9 (i)
(ii) np ≥ 5 عندما 0 < p ≤ 0.5 ≤ (1-p)
وبالمثل فی (2003) وصف کل من (Fleiss, Levin و Paik) أن التقریب یکون مناسباً عندما تکون قیمة p تقع فی هذه الفترة (0.3 ≤ p ≤ 0.7) و عندما تکون n کبیرة للغایة بالنسبة ل ≥5) (n p و n (1-p ) ≥ n بینما یبقى التباین ثابتاً np(1-p) . فی حین أشار( Montgomery) فی (2005) إلى أن التقریب الطبیعی لتوزیع ذی الحدین یکون مناسباً عندما :
(i)np ≥ 10
p ( 0.1 ≤ p ≤ 0.9 ) (ii)
وفی الکثیر من الدراسات ، تُستخدم لوحة p عندما تکون قیم المعلمة p صغیرة (الخ ; …0.01 ; 0.05 ; 0.1 p = 0.001;). فی هذه الحالات ، یکون التوزیع ذو الحدین منحرفًا تمامًا، ویکون التقریب للتوزیع الطبیعی غیرمناسب ؛ لأنه یسمح بقیم سالبة أو أکبر من واحد . وفی عام (1991) اقترح Queensberry لوحة Qذات الحدین (Binomial Q Chart ) لمراقبة البیانات الکسریة غیر المطابقة باستخدام التحویل غیر الخطی لحدود السیطرة ، وأظهر أنه یقترب من التوزیع الطبیعی أقرب إلى ذات الحدین . اما فی (1996) فقد قدّم (Heiman) تعدیلاً لحدود السیطرة فی لوحة p لأحجام العینات الکبیرة sam-ple (10000n >) مشیرًا إلى أنه فی هذه الحالة ، تکون حدود السیطرة ضیقة ، وبالتالی یزید من معدل الإنذار الخاطئ (false alarm ). کما اقترح کل من( SchwertmanوRayan) فی (1997) تعدیلا لحدود السیطرة للوحة nP لتطابق التقریب الطبیعی ( Normal approximation) عندما ( 0.03 ( p < . بینما أقترح Chen) (فی (1998) تعدیلا على حدود السیطرة للوحة P وقارنها بلوحة P التقلیدیة ولوحة Q ذات الحدین باستخدام معدل الإنذار الخاطئ . اما فی (2008) فقد قام کل من ( Sim و ( Limبتکییف لوحة السیطرة للصفات لرصد البیانات ذات القیمة الصفریة واستخدما الفاصل الزمنی (Blyth-Still ) مع 3σ لحساب حدود السیطرة على افتراض أن هذه البیانات تتبع التوزیع ذا الحدین وبواسون . وقارن (Bourke) فی عام (2008) الأداء لأربع لوحات سیطرة من خلال مراقبة تحولات الجزء غیر المطابق فی العملیات الصناعیة. ومن ثم لاحظ أوجه التشابه فی أداء لوحة السیطرة الترکیبی(Synthetic control chart) ولوحة np على مدى فترة زمنیة طویلة من عملیة السیطرة. ولأغراض مماثلة قارن (Aebtarm و Bouguila) فی (2011) الأداء لاحدى عشرة لوحة سیطرة لرصد العیوب مع توزیع بواسون . عندما یکون توزیع البیانات فی العملیة الصناعیة غیر مثماتل یزداد معدل الانذار الخاطئ مع عدم التناسق بسبب التناقض بین شکل توزیع البیانات والتوزیع الطبیعی لها . واقترح ( Nelson ) فی هذه الحالات افتراض أن توزیع البیانات معروف وانشاء لوحة سیطرة بحدود دقیقة یوفّر معدلات إندار خاطئة مرغوبة . لقد تم فی هذا البحث بتطبیق لوحة بیتا للسیطرة على متغیرات من نوع الکسر فی العملیات الصناعیة . حیث تفترض لوحة السیطرة هذه أنه یمکن تقریب بیانات الکسریة عن طریق توزیع بیتا وبناء حدود سیطرة جدیدة بناءً على هذا التوزیع . تم تطبیق لوحة بیتا على بیانات نسبة الفول السودانی الملوث على کل من لوحتی بیتا للسیطرة ولوحة شیوارت (Shewhart (p) ) التقلیدیة لأغراض المقارنة باستخدام البرامجیتین Minitab , R .
2 -هدف البحث
تم استخدام لوحة السیطرة بیتا(Beta control chart) لمراقبة العملیات التی تقیس متغیرات من النوع الکسریة. وتعتمد لوحة السیطرة هذه على توزیع احتمالیة بیتا لتقدیر حدود السیطرة للبیانات الکسریة التى لها قیم مقیدة (Restricted) فی الفترة الزمنیة [ 0 , 1 ] . ومن ثم تطبیق اللوحة على بیانات ومقارنتها مع لوحة P ، ثم حساب تحلیل الحساسیة لعدة قیم من الکسور ای النسب .
3 - لوحة بیتا للسیطرة Beta Control Chart
یستخدم توزیع بیتا لنمذجة المتغیرات المستمرة المحددة ضمن الفترة الزمنیة] 1 ، 0 [ ، مثل المعدلات ،الکسوروالنسب و کذلک توزیع بیتا ، یکون محددا أیضًا عندما یقتصر متغیر قید الاهتمام y على المدة الزمنیة المستمرة a , b] [ حیث إن a وb قیم معلومة وa < b فی هذه الحالات دون فقدان العمومیة یمکن نمذجة (y - a) / (b - a) بدلا من نمدجة y مباشرة . حیث إن y متغیر عشوائی یتوزع بیتا . وصف Johnson واخرین فی (1995) ، وان المتغیر العشوائی y یتوزع توزیع بیتا بمعلمتی شکل و بدالة کثافة احتمالیة لها تکون على النحو الآتی :
f(y; , ) = وان 0 < y < 1
حیث إن> 0 , هی معلمات الشکل للتوزیع و Γ (Ɵ) هی دالة کاما. وان متوسط توزیع بیتا و تباینه هما على التوالی:
E(y) = …. ( 1 )
Var(y) = ….. ( 2 )
وان دالة التوزیع التراکمی (CDF))) لتوزیع بیتا
F(y; ) =
حیث إنB ( , ) ، هی دالة بیتا و ) , B (y ؛ هی دالة بیتا غیر الکاملة. حیث أن الدالة الکمیة quantile function)) هی
ψ(α;θ1,θ2) = (y; )
Gupta 2004 ) و Nadarajah )
أن توزیع بیتا متعدد الاستخدامات للغایة ، وله مجموعة متنوعة من التطبیقات (Johnson وآخرون 2005 ( فی حین أن ( ten Caten Sant0Anna, ) افترض فی سنة (2012) أن المتغیرات من نوع الکسریة والنسبة تتوزع توزیع بیتا ، وأن لوحة بیتا تستوعب بشکل طبیعی عدم تماثل نوع البیانات الکسریة والنسبة ، وتکون حدود السیطرة تقتصر على الفترة [ 0 , 1] إن هذه الخصائص للوحة بیتا مفیدة مقارنة باللوحات المعتادة التی تفترض التقریب الطبیعی ، مثل لوحتی p و np ، المعروفة التی ترصد الوحدات غیر المطابقة . یتم تحدید حدی السیطرة االادنى (LCL ) والأعلى ( UCL) للوحة بیتا للسیطرة بواسطة
LCL = - W 1 ; UCL= + W 2
حیث إن و یمثلان المتوسط وتباین المتوسط المتغیر الکسری ، وان ,W12Wهما ثوابت مجدولة التی تحدد عرض منطقة السیطرة ( Control rejoin) لغرض بناء حدود السیطرة للوحة بیتا ، بالاعتماد على قیم (n و ) على الرغم من أن لوحة بیتا المقترحة یقدم أساسًا نظریًا صارمًا ، إلا أن هذا یعدّ نقطة سلبیة فی عدم تضمین هذه اللوحة فی حزم البرامج الاحصائیة المتوافرة حالیا . تعرض لوحة Beta المقترحة التطبیق العملی نفسه والبساطة التشغیلیة کما فی لوحة p التی اقترحها Shewhart حیث لا قیود على حجم العینة و اسلوب أخذ العینات. على سبیل المثال ، یمکن جمع أحجام العینات 10 أو 15 أو 20 و ستکون تحلیلات الکشف عن التحولات فی العملیة بالکفاءة نفسها عندما یکون حجم العینات من 50 أو 100 أو 200. لمراقبة العملیات طویلة الأجل ، یجب أن یتم إعادة تقییم حدود السیطرة بشکل دوری ، حتى لو لم یکن هناک دلیل على وجود اختلافات ویسمح هذا الإجراء الدینامیکی للوحة السیطرة بتوفیر مراقبة فعالة واکتشاف مصادر التباین. ( Bayer, eals , 2018 )
4 –تحلیل الحساسیة Sensitivity Analysis
هناک العدید من الأسالیب لقیاس کفاءة أداء لوحات السیطرة ، و یعّد معدل طول التشغیل ( Average Run Length (ARL) ) المقیاس الأکثر استخداماً فی عملیة المراقبة الاحصائیة لتقییم أداء حساسیة لوحات السیطرة للکشف عن الانحرفات فی مخرجات العملیة . ویعدّ العالم (Pag) اول من استخدم مفهوم معدل طول التشغیل فی عام (1954) وعرفه بأنه متوسط عدد العینات المفحوصة ، والتی تکون تحت السیطرة لحین حدوث مایشیر الى أن العملیة اصبحت خارج حدود السیطرة. (2005 , Montgomery)
5 - إستخدام معدل طول التشغیل لتقیم أداء لوحة السیطرة Using Average Run Length to determine the performance of a control chart
فی مجال اختبارالفرضیات الإحصائیة یمکننا قیاس الخطأ من النوع الاول والخطأ من النوع الثانی β . ولا یمکن استخدام هذه الطریقة للوحات السیطرة ؛ لأن لوحات السیطرة هی سلسلة من العینات بحسب اختبارات فرضیة العینة وبدلا من الخطأ التقلیدی و β المستخدم باختبار الفرضیات ، عادة مایستخدم متوسط طول التشغیل (ARL ) لقیاس أداء لوحات السیطرة . فی الشکل ( 2 - 1) تم أخد أربع عینات قبل أن تتجاوز نقطة ما حد السیطرة الأعلى ، طول التشغل هو متغیّر عشوائّی یمکن أن یأخد قیماً صحیحة فقط ، ونحن بحاجة الى استخدام القیمة المتوقعة أو متوسط طول التشغیل لقیاس الأداء . (ARL ) متوسط طول التشغیل هو القیمة المتوقعة من طول التشغیل ، والذی یمثل العدد المتوقع للمجموعات الفرعیة حتى تشیر للوحة للخروج عن السیطرة أولاً . یتم الان حساب معدل طول التشغیل فی الحالتین ، وهما :
الشکل ( 1 – 2 ) تحدید طول التشغیل للوحة
أولاً : حساب طول التشغیل فی حالة العملیة تحت السیطرة (Calculate Average Run Length In - Control )
متوسط طول التشغیل ARL0
إن متوسط طول التشغیل (ARL0 ) هو العدد المتوقع للعینات حتى تشیر لوحة السیطرة الى مؤشر خروج اللوحة عن السیطرة مع العلم أنها تحت السیطرة . و یمکننا ایضا القول بأنه قد تم ارسال إندار. خاطئ (false alarm ) ، ولتعریف متوسط طول التشغیل (ARL0 ) ریاضیاً فإنه یفسرعلى أنه دالة للخطأ من النوع الأول ( ) ویعتبر معیارا لکلفة القرارات الخاطئة ویمثل کالاتی:-
]مع العلم أن العملیة تحت السیطرة اشارات للخروج عن السیطرة Pr [=
] P0= P │ UCL < y < LCL [ P – 1 =
حیث إن P0 تمثل متوسط المتغیّر Y فی حالة استقرار العملیة وعدم حدوث تغییر فیها
ولاشتقاق متوسط طول التشغیل ARL0 طول التشغیل احتمالیة الاشارة خارج السیطرة
1
2
3
4 .......
: ........
K α K-1
وبهذا یکون معدل طول التشغیل
= E ( ) = K
= ... ( 3 )
فإذا لم یحدث تغییرفی العملیة فاحتمال وقوع نقطة خارج حدی السیطرة فان 0.0027 = α
وبذلک فإن متوسط طول التشغیل یکون370.4 = 0.0027 / 1 = وهذا یعنی أنه فی حالة استقرار العملیة یتوقع حدوث مؤشر خارج السیطرة بعد فحص 370.4 عینة .
ثانیا : حساب طول التشغیل للعملیة خارج السیطرة Calculate Average Run Length Out - Control ))
ان متوسط طول التشغیل یمثل العدد المتوقع للعینات من دون إشارة الى وجود مؤشر یشیر الى خروج اللوحة عن السیطرة مع العلم أنها خارج حدود السیطرة وضمن دالة من الخطأ من النوع الثانی . ای افتراض ان تغیرا قد حدث فی مخرجات العملیة عن مرکزها بعدد من الانحرفات الخطأ من النوع الثانی, وهو یقیس معدل عدد المواد المنتجة الردیئة ویمثل کالاتی:-
]مع العلم أن اللوحة خارج حدود السیطرة عدم الاشارة خارج السیطرة لعینة ما Pr [= β
] P1 = P │ UCL < y < LCL [ P – 1 = β
]مع العلم أن اللوحة خارج حدود السیطرة تظهر لوحة شیوارت اشارات لعینة ما Pr [= 1 – β
وحیث إن λ + P0 = P1 وان λ تمثل التغیّر الناتج فی العملیة وهو دالة للتباین و P1 تمثل متوسط المتغیرY فی العملیة الخارجة عن السیطرة . ولاشتقاق متوسط طول التشغیل
طول التشغیل احتمالیة الاشارة داخل السیطرة
1 β
2
3
: ......
K
وبهذا فإن یحسب کما یأتی :
= E ( ) = K
= .... ( 4 )
وبهذا فان معدل طول التشغیل هو عدد النقاط المرسومة ( المجموعة الجزئیة ) من بدایة تغییر العملیة وحتى حدوث اشارة خارجة عن السیطرة لمقارنة أداء لوحات سیطرة مختلفة من حیث قدرتها على کشف التغیرات . وبهذا نستطیع القول إن معدل طول التشغیل یستخدم لعرض اداء لوحة معینة للسیطرة أو لمقارنة اداء عدة لوحات للسیطرة . ( 2005 , Montgomery) على الرغم من أن لوحة بیتا المقترحة یقدم أساسًا نظریًا صارمًا ، إلا أن هذا یعدّ نقطة سلبیة فی عدم تضمین هذه اللوحة فی حزم البرامج الإحصائیة
6 – الجانب التطبیقی
ولغرض إجراء تقییم أداء لوحتی P وبیتا تم استخدام مجموعة من البیانات الحقیقیة لدراسة الفول السودانی الملوث بالمواد السامة فی 34 دفعة من 120 رطلا المتغیر المرصود هو نسبة الفول السودانی الملوث کما هو موضح فی الجدول ( 1 ) (السدحان ، 2010)
الجدول (1)نسب الفول السودانی الملوث
البیانات |
ت |
البیانات |
ت |
0.933 |
18 |
0.971 |
1 |
0.858 |
19 |
0.979 |
2 |
0.987 |
20 |
0.982 |
3 |
0.958 |
21 |
0.971 |
4 |
0.909 |
22 |
0.957 |
5 |
0.859 |
23 |
0.961 |
6 |
0.863 |
24 |
0.956 |
7 |
0.811 |
25 |
0.972 |
8 |
0.877 |
26 |
0.889 |
9 |
0.798 |
27 |
0.961 |
10 |
0.855 |
28 |
0.982 |
11 |
0.788 |
29 |
0.975 |
12 |
0.821 |
30 |
0.942 |
13 |
0.830 |
31 |
0.932 |
14 |
0.718 |
32 |
0.908 |
15 |
0.642 |
33 |
0.970 |
16 |
0.658 |
34 |
0.985 |
17 |
ولغرض دراسة تحلیل هذه البیانات قمنا اولا بعمل توصیف لهذه البیانات کما هو موضّح فی الجدول(2) ( تم احتساب التباین بالصیغة العامة له)
الجدول ( 2 ) الملخص الاحصائی لمتغیر البیانات
المتغیر |
القیمة الصغرى |
القیمة العظمى |
المتوسط |
التباین |
Y |
0.642 |
0.987 |
0.895824 |
0.00875 |
وبعد ذلک قمنا برسم کل من المدرج التکراری ودالة الکثافة الاحتمالیة للمتغیر قید الدراسة ( نسبة الفول السودانی الملوث) وکما هو موضح فی الشکلین ( 1 ) ( 2 ) .
الشکل ( 1 ) المدرج التکراری للبیانات
الشکل (2 ) دالة الکثافة الاحتمالیة للبیانات
ولتحدید التوزیع الاحتمالی للبیانات قمناباختبار الفرضیة الآتیة :
البیانات تتوزع توزیع بیتا : H0
البیانات لا تتوزع توزیع بیتا : H1
تم إجراء کل من الاختبارات التالیة اختبار Kolmogorov _ Smirnov ، اختبار حسن المطابقة Chi _ Squared واختبار Anderson_ Darling وعن طریق مقارنة قیم المختبر الاحصائی مع القیم الحرجة عند مستوى معنویة 0.05 , 0.01 = تبین أن قیم المختبر الاحصائی اصغر من القیم الحرجة ای انها تقع ضمن منطقة القبول ، وکما هو موضح فی الجدول (3) مما یؤید ان البیانات تتوزع توزیع بیتا .
الجدول ( 3 ) اختبارات حسن المطابقة
7 - تحلیلات المقارنة لوحتی بیتا وشیوارت ( p) على نسب الفول السودانی الملوث
تم حساب حدود السیطرة عن طریق الرسم البیانی لاحتمالیة الاندارات الکاذبة ( false alarm ) فی عملیة المراقبة بناءً على التوزیع الطبیعی للوحة ( p ) لشیوارت وبناءً على توزیع بیتا للوحة بیتا و الجدول ( 4 ) یبین حدود السیطرة المحسوبة باستخدام تقدیرات العینة لکلتا اللوحتین.
الجدول ( 4 ) حدود السیطرة للوحتی P وبیتا
T |
LCL |
UCL |
|
0.02364775 |
0 |
0.1087531 |
Shewhart |
0.8958235 |
0.4654641 |
0.9998359 |
Beta |
ومن ثم تم رسم کل من لوحتی P وبیتا ، لغرض مقارنتها کما فی الشکلین المتتالین ( 3 ) و ( 4 ) على التوالی وفق البیانات المستخدمة.
الشکل ( 3 ) لوحة P شیوارت |
الشکل ( 4 ) لوحة بیتا للسیطرة |
من خلال ملاحظة الشکلین ( 3 ) و( 4 )أظهرت لوحة بیتا التغیرات الحاصلة فی العملیة بشکل أکثر وضوحا من لوحة P حیث إن النقاط قبل العینة 23 کانت فوق خط الهدف ومابعد العینة 23 ظهرت معظمها تحت خط الهدف ، فهذا یدل على السلوک غیر الطبیعی لحجم السموم فی المنتج اما لوحة p فلم تتمکن من توضیح السلوک الغیر طبیعی لعملیة السیطرة .
8 - تقییم الأداء للوحتی بیتا وشیوات Evaluation performance between shewart and beta chart
نقوم الآن بإجراء دراسة لمقارنة لوحة شیوارت وبیتا فی فی حالتین هما تحت السیطرة (in – control ) وخارج نطاق السیطرة ( out of control) بحساب متوسط طول التشغیل لقیم مختلفة من و n وثم إجراء تحلیل الحساسیة باستخدام البرامجیة الحسابیة R ، بالنسبة لحالة تحت السیطرة نقوم بحساب ARL0 من المعادلة ( 3 ) ، وبالنسبة للعملیة خارج السیطرة یتم حساب ARL1من المعادلة ( 4 ) وکما أشرنا سابقا فإن ARL0 هو دالة للخطأ من النوع الأول فی حین ARL1هودالة للخطأ من النوع الثانی .
الجدول ( 5 ) مقارنة لقیم کل من ARL لکل من لوحتی شیوارت وبیتا للقیم n = 250 λ = 0.01 P0 = 0.01 و P0 =0.001 n = 2000 λ = 0.01
|
N
|
2000 LCL |
UCL |
|
β |
ARL |
n
|
250 LCL |
UCL |
|
β |
ARL |
Shewhart |
0.001 |
-0.00112 |
0.031 |
0.0032 |
|
312 |
0.001 |
-0.004996 |
0.01898 |
0.0042 |
|
238 |
|
0.002 |
-0.00096 |
0.00 4996 |
|
0.9 |
10 |
0.02 |
-0.006562 |
0.04696 |
|
0.88 |
8 |
|
0.003 |
-0.0067 |
0.006671 |
|
0.03 |
1 |
0.03 |
-0.00266 |
0.06266 |
|
0.707 |
3 |
|
0.005 |
0.0002687 |
0.0097313 |
|
0.21 |
1 |
0.05 |
0.0038 |
0.09135 |
|
0.120 |
1 |
|
0.01 |
0.003325 |
0.016675 |
|
0.02 |
1 |
0.01 |
-0.00887 |
0.02887 |
|
0.1 |
1 |
Beta |
0.001 |
0.00002 |
0.0049 |
0.0027 |
|
370 |
0.01 |
0.00028 |
0.043 |
0.0027 |
|
730 |
|
0.002 |
0.000014 |
0.0071 |
|
0.003 |
1 |
0. 02 |
0.00231 |
0.063 |
|
0.860 |
7 |
|
0.003 |
0.0004 |
0.009 |
|
0.004 |
1 |
0. 03 |
0.0059 |
0.079 |
|
0.560 |
2 |
|
0.005 |
0.0013 |
0.0123 |
|
0.006 |
1 |
0.05 |
0.0157 |
0.107 |
|
0.072 |
1 |
|
0.01 |
0.004 |
0.0194 |
|
0.001 |
1 |
0. 01 |
0.0482 |
0.173 |
|
0.001 |
|
الجدول ( 6 ) مقارنة لقیم کل من ARL لکل من لوحتی شیوارت وبیتا للقیم P0 = 0.1 n = 350 λ = 0.02 و P0= 0. 1 n = 75 λ = 0.002
|
n
|
350 LCL |
UCL |
|
β |
ARL |
n
|
75 LCL |
UCL
|
|
β |
ARL |
Shewhart |
0.1 |
-0.3563 |
0.25213 |
0.0032 |
|
333 |
0.1 |
-0.00392 |
0.20392 |
0.003 |
|
333 |
|
0.12 |
0.10329 |
0.13671 |
|
0.95 |
20 |
0.12 |
0.00743 |
0.2325 |
|
0.97 |
33 |
|
0.15 |
0.09274 |
0.20725 |
|
0.50 |
2 |
0.15 |
0.02631 |
02736 |
|
0.92 |
13 |
|
0.2 |
0.13586 |
0.2641 |
|
0.015 |
1 |
0.2 |
0.06144 |
0.33856 |
|
0.72 |
4 |
|
0.3 |
0.22652 |
0.37348 |
|
0.21 |
1 |
0.3 |
0.14126 |
0.45874 |
|
0.2 |
1 |
Beta |
0.1 |
0.054 |
0.156 |
0.0027 |
|
370 |
0.1 |
0 |
0.261 |
0.027 |
|
370 |
|
0.12 |
0.0712 |
0.172 |
|
0.966 |
31 |
0.12 |
0.0250 |
0.293 |
|
0.971 |
34 |
|
0.15 |
0.092 |
0.217 |
|
0.905 |
11 |
0.15 |
0.375 |
0.324 |
|
0.905 |
11 |
|
0.2 |
0.135 |
0.264 |
|
0.650 |
3 |
0.2 |
0.65 |
0.395 |
|
0.65 |
3 |
|
0.3 |
0.225 |
0.382 |
|
0.104 |
1 |
0.3 |
0.1312 |
0.519 |
|
0.104 |
1 |
من خلال ملاحظة الجدولین ( 5 ) و( 6 ) نجد أن حساسیة لوحات السیطرة التی تم فحصها عندما ( P = 0.001 ) وحجم العینة 2000 لاحظنا أن ARL0 = 312 للوحة شیوارت بینما کانت للوحة بیتا = 370 ARL1على الرغم من أن لوحة شیوارت توفر اکتشافات سریعة للتحولات فی قیم المعلمات :0.01) = ( 0.002 , 0.003 , 0.005 , ولکن لوحة بیتا تظهر حساسیة أکثر لاحتمال الانذار الخاطئ ( false alarm ) والاداء فی اکتشاف التغیرات ( detecting changes) مقارنة بلوحة شیوارت وبالمثل یحدث لقیمة المعلمة ( ) و n = 250 حیث تقدم لوحة بیتا أداء أفضل فی اکتشاف التغیرات عندما 0.1 P =و n = 75 فان حدود السیطرة للوحة شیوارت فی حالة تحت السیطرة تؤثر حدود السیطرة فی الاداء غیر المناسب . حیث تظهر حدود السیطرة هذه اوجه التشابه والأداء الأفضل عندما P = 0.1 و n > 250 . تفرض حدود السیطرة المقترحة فی لوحة بیتا أن متوسط طول التشغیل 370 = ARL0باحتمالیة 0.9973 ای ان مع اکتشاف اسرع للتحولات لجمیع حالات الخروج عن السیطرة ذات القیم الاقل . نظرا لکون التوزیع ذی الحدین وهو من التوزیعات المتقطعة یتم اقتطاع ( truncated) حدود السیطرة الدنیا ل 3σعند القیمة صفر. و القیمة فإن لوحة السیطرة مع التقریب بالتوزیع الطبیعی تقدّم حدوداً ضیقة ( narrow limits) تولد متوسط طول التشغیل ، ویکون اقل من 370 من العینات وانذارات خاطئة.
9 - تعریف حدود الشکل المغلق Definition of closed form limits
تم حساب حدودالسیطرة الدنیاوالعلیا لقیم النسب الصغیرة الاَتیة (P = 0.002 , 0.02 . 0.05 and0.1 ) وبأحجام عینات مختلفة حیث ان ( np = 1 , 2 , . . …. , 10 ) فعلى سبیل المثال ننظر لقیم الحدود الدنیا للوحة شیوارت وبیتا عندما n = 100 و p = 0.02 نلاحظ أن الحد الأدنى للوحة Shewhart یساوی – 0.022 فی حین أن الحد الأدنى للوحة Beta المقابلة له یساوی 0.000351 کما هو موضح فی الجدول ( 7 ) حیث إن قیم لوحة بیتا تقع ضمن المدى [ 0 , 1 ] فی حین لوحة شیوارت تکون قیم الحدود الدنیا لها موجبة عندما np ≥ 8 │p ≥ 0.2 بینما حدود السیطرة فی لوحة بیتا مناسبة لای قیم n و p ، لأن الحدود الخاصة بها لا تعتمد على حجم العینة n .
الجدول ( 7 ) حدود السیطرة للبیانات بواسطة لوحتی شیوارت وبیتا
|
Beta |
|
Shewart |
|
|
|
UCL |
LCL |
UCL |
LCL |
np |
n |
P |
0.00757 |
0.000002 |
0.00799 |
-0.00399 |
1 |
500 |
0.002 |
0.07941 |
0.000023 |
0.07939 |
-0.03939 |
1 |
50 |
0.02 |
0.4191 |
0.000062 |
0.1962 |
-0.09620 |
1 |
20 |
0.05 |
0.6541 |
0.000073 |
0.736656 |
-0.3366 |
1 |
5 |
0.2 |
0.00541 |
0.000032 |
0.006238 |
-0.00223 |
2 |
1000 |
0.002 |
0.05291 |
0.000351 |
0.062 |
-0.022 |
2 |
100 |
0.02 |
0.2234 |
0.002251 |
0.1533 |
-0.0533 |
2 |
40 |
0.05 |
0.3124 |
0.0345 |
0.5799 |
-0.1794 |
2 |
10 |
0.2 |
0.004520 |
0.000081 |
0.01238 |
-0.00146 |
3 |
1500 |
0.002 |
0.04532 |
0.000621 |
0.0542 |
-0.0142 |
3 |
150 |
0.02 |
0.17271 |
0.004541 |
0.01344 |
-0.0344 |
3 |
60 |
0.05 |
0.4672 |
0.008125 |
0.5098 |
-0.1098 |
3 |
15 |
0.2 |
0.00406 |
0.000215 |
0.00499 |
-0.000009 |
4 |
2000 |
0.002 |
0.0414 |
0.00216 |
0.0496 |
-0.0096 |
4 |
200 |
0.02 |
0.1621 |
0.00651 |
0.1231 |
-0.0231 |
4 |
80 |
0.05 |
0.3113 |
0.02173 |
0.4683 |
-0.06832 |
4 |
20 |
0.2 |
0.00375 |
0.000257 |
0.00468 |
-0.0006 |
5 |
2500 |
0.002 |
0.0376 |
0.00257 |
0.04656 |
-0.00656 |
5 |
250 |
0.02 |
0.2765 |
0.0773 |
0.11538 |
-0.01538 |
5 |
100 |
0.05 |
0.3560 |
0.074 |
0.44 |
-0.04 |
5 |
25 |
0.2 |
0.00357 |
0.000285 |
0.00444 |
-0.000447 |
6 |
3000 |
0.002 |
0.0274 |
0.00287 |
0.04428 |
-0.00424 |
6 |
300 |
0.02 |
0.0397 |
0.00495 |
0.10966 |
-0.00968 |
6 |
120 |
0.05 |
0.2571 |
0.0213 |
0.4191 |
-0.0191 |
6 |
30 |
0.2 |
0.00315 |
0.00312 |
0.004265 |
-0.000026 |
7 |
3500 |
0.002 |
0.03314 |
0.00327 |
0.04245 |
-0.00245 |
7 |
350 |
0.02 |
0.2218 |
0.0215 |
0.10525 |
-0.00525 |
7 |
140 |
0.05 |
0.0412 |
0.0326 |
0.40283 |
-0.002837 |
7 |
35 |
0.2 |
0.002935 |
0.000359 |
0.004119 |
-0.0000119 |
8 |
4000 |
0.002 |
0.02825 |
0.00359 |
0.041 |
-0.001 |
8 |
400 |
0.02 |
0.1129 |
0.02341 |
0.010169 |
-0.00169 |
8 |
160 |
0.05 |
0.0071 |
0.003567 |
0.3897 |
0.010263 |
8 |
40 |
0.2 |
0.00396 |
0.000384 |
0.00398 |
0.000003 |
9 |
4500 |
0.002 |
0.03261 |
0.00385 |
0.0397 |
0.00002 |
9 |
450 |
0.02 |
0.2063 |
0.02487 |
0.09873 |
0.00126 |
9 |
180 |
0.05 |
0.0361 |
0.03101 |
0.37888 |
0.021145 |
9 |
45 |
0.2 |
0.00322 |
0.000408 |
0.00389 |
0.0000105 |
10 |
5000 |
0.002 |
0.03215 |
0.00409 |
0.03878 |
0.001217 |
10 |
500 |
0.02 |
0.00886 |
0.01667 |
0.09623 |
0.003766 |
10 |
200 |
0.05 |
0.03166 |
0.0329 |
0.3697 |
0.0303 |
10 |
50 |
0.2 |
10 - الاستنتاجات
عندما تم تطبیق البیانات الخاصة لنسب الفول السودانی الملوث على کل من لوحتی بیتا وشیوارت (P) تبین أن لوحة بیتا تودی نتائج أکثر دقة من لوحة شیوارت . أکد تحلیل الحساسیة باستخدام العدید من القیم الکسریة الأداء المتمیز للوحة بیتا مقارنة بلوحة شیوارت (P) ونتج عن التقریب المقترح ARL0 أکبر قلیلاً ، وأن ARL1 أصغر بکثیر.