Xây dựng phần mềm sử dụng biểu đồ shewhart để kiểm soát khối lượng trung bình viên trong quá trình dập viên

Trước tình hình cạnh tranh gay gắt của thị trường dược phẩm trong và ngoài nước, cùng với nhu cầu của người sử dụng thuốc đòi hỏi về chất lượng và hiệu quả ngày càng cao. Hầu hết các công ty dược phẩm trong nước xem chất lượng sản phẩm là yếu tố sống còn của đơn vị. Thời gian qua, Bộ Y tế ban hành quy định về áp dụng nguyên tắc tiêu chuẩn GMPWHO để các doanh nghiệp nâng cấp nhà xưởng, đầu tư máy móc thiết bị sản xuất, nghiên cứu cải tiến kỹ thuật bào chế thuốc...... bảo đảm thuốc được sản xuất có chất lượng và hiệu quả điều trị như mong muốn của người sử dụng, đồng thời chất lượng thuốc cũng có thể được so sánh ngang hàng với các nước trong khu vực cũng như thế giới. Và một trong số các quy định trong GMP là kiểm soát được quá trình sản xuất, có nghĩa là phải kiểm soát được các thông số vật lý của quá trình như: Hình thức cảm quan, đồng đều khối lượng, độ cứng, bề dầy viên, độ mài mòn, thời gian rã viên….Trong đó, đồng đều khối lượng viên là một trong những yếu tố đầu tiên cần phải được kiểm soát trong quá trình nhằm thỏa mãn yêu cầu về chất lượng sản phẩm. Trong thời gian qua, việc kiểm soát các thông số vật lý của sản phẩm trong quá trình được thực hiện bằng thủ công, phụ thuộc chủ yếu vào người vận hành thiết bị. Các dữ liệu thu được có thể mang tính khách quan lẫn chủ quan, nên việc thu thập và lưu trữ dữ liệu đôi khi không chính xác và mang tính đại diện, mà đặc biệt chưa tính được khối dữ liệu có phân bố chuẩn hay không, độ rộng của từng cỡ mẫu, khuynh hướng dịch chuyển của các giá trị trung bình là như thế nào. Trước nhu cầu phải nâng cao hệ thống kiểm soát quá trình để giải quyết các nhu cầu trên, chúng tôi tiến hành nghiên cứu: “Xây dựng phần mềm sử dụng biểu đồ Shewhart để kiểm soát khối lượng trung bình viên trong quá trình dập viên” , giúp cho nhà sản xuất nâng tầm quản lý lên một bước.

1 ĐẶT VẤN ĐỀ Trước tình hình cạnh tranh gay gắt thị trường dược phẩm nước, với nhu cầu người sử dụng thuốc đòi hỏi chất lượng hiệu ngày cao Hầu hết công ty dược phẩm nước xem chất lượng sản phẩm yếu tố sống đơn vị Thời gian qua, Bộ Y tế ban hành quy định áp dụng nguyên tắc tiêu chuẩn GMPWHO để doanh nghiệp nâng cấp nhà xưởng, đầu tư máy móc thiết bị sản xuất, nghiên cứu cải tiến kỹ thuật bào chế thuốc bảo đảm thuốc sản xuất có chất lượng hiệu điều trị mong muốn người sử dụng, đồng thời chất lượng thuốc so sánh ngang hàng với nước khu vực giới Và số quy định GMP kiểm sốt q trình sản xuất, có nghĩa phải kiểm sốt thơng số vật lý q trình như: Hình thức cảm quan, đồng khối lượng, độ cứng, bề dầy viên, độ mài mòn, thời gian rã viên….Trong đó, đồng khối lượng viên yếu tố cần phải kiểm sốt q trình nhằm thỏa mãn u cầu chất lượng sản phẩm Trong thời gian qua, việc kiểm sốt thơng số vật lý sản phẩm q trình thực thủ cơng, phụ thuộc chủ yếu vào người vận hành thiết bị Các liệu thu mang tính khách quan lẫn chủ quan, nên việc thu thập lưu trữ liệu đơi khơng xác mang tính đại diện, mà đặc biệt chưa tính khối liệu có phân bố chuẩn hay khơng, độ rộng cỡ mẫu, khuynh hướng dịch chuyển giá trị trung bình Trước nhu cầu phải nâng cao hệ thống kiểm sốt q trình để giải nhu cầu trên, tiến hành nghiên cứu: “Xây dựng phần mềm sử dụng biểu đồ Shewhart để kiểm sốt khối lượng trung bình viên q trình dập viên” , giúp cho nhà sản xuất nâng tầm quản lý lên bước Trang TỔNG QUAN TÀI LIỆU Lợi ích việc kiểm sốt biểu đồ Shewhart cho phép nhà sản xuất theo dõi trình đo lường độ tin cậy phân tích xu hướng quy trình, trở thành công cụ hữu hiệu để xác định loại bỏ nguyên nhân ảnh hưởng xấu đến biến thiên quy trình Từ đó, giúp ích nhiều cho việc cải tiến quy trình quản lý chất lượng sau Biểu đồ Shewhart loại biểu đồ thống kê dùng phổ biến kiểm sốt quy trình sản xuất, khơng để phát dịch chuyển giá trị trung bình mà giúp theo dõi biến thiên giá trị trung bình phạm vi 2-3sigma Biểu đồ Shewhart giúp phân biệt biến thiên quy trình nguyên nhân phổ biến hay nguyên nhân đặc biệt Hiệu việc phòng ngừa sai lỗi Các giới hạn kiểm soát cảnh báo biểu đồ giúp phòng ngừa sai lỗi, ví dụ: sản phẩm vượt ngồi giới hạn kiểm sốt Cung cấp thơng tin hiệu trình Dữ liệu sử dụng cho biểu đồ kiểm sốt hữu ích cho việc tính tốn số hiệu q trình Việc phân tích liệu với biểu đồ Shewhart có ý nghĩa liệu thu thập đạt phân phối chuẩn, phải xem xét phân phối chuẩn khối liệu trước xem xét vẽ biện luận biểu đồ Shewhart Trang 2.1 KHẢO SÁT PHÂN PHỐI CHUẨN CỦA DỮ LIỆU BẰNG TRẮC NGHIỆM CHI BÌNH PHƯƠNG Trong thực tế, thường gặp loại trắc nghiệm χ2: - Trắc nghiệm tính tương thích (test of goodness-of-fit): Để khảo sát liệu quan sát có thuộc phân phối chuẩn? - Trắc nghiệm tính độc lập (test of independence): Để khảo sát liệu quan sát có độc lập với nhau? - Trắc nghiệm tính đồng (test of homogeneity): Để khảo sát liệu quan sát có nguồn gốc xuất phát từ dân số? Như vậy, để khảo sát liệu có thuộc phân phối chuẩn hay không ta dùng trắc nghiệm tính tương thích - Để kiểm tra tính phân bố chuẩn quần thể mẫu, phương pháp phổ biến cho kết có tính định lượng dùng thử nghiệm chi bình phương để so sánh số cá thể lớp biểu đồ mật độ phân bố (histogram) với số cá thể lý thuyết tương ứng phải có - Nếu χ2 thực tế < χ2 lý thuyết  Có phân bố chuẩn quần thể mẫu - Thực thử nghiệm χ2 cần đến liệu số lớp, khoảng rộng lớp, giới hạn để vẽ biểu đồ mật độ phân bố Để khơng dẫn tới kết sai lệch, việc tính tốn số liệu phải tuân thủ theo quy định 2.1 VẼ BIỂU ĐỒ SHEWHART X VÀ R - Các biểu đồ Shewhart X R thường sử dụng theo cặp đôi - Biểu đồ X nhằm phát thay đổi giá trị trung bình quy trình Trang - Biểu đồ R nhằm giám sát độ phân tán quy trình cách theo dõi biến thiên khoảng (ranges) mẫu Những thay đổi khoảng báo hiệu thay đổi biến thiên q trình - Trình tự thích hợp vẽ biện luận biểu đồ kiểm soát R trước Nếu biện luận cho thấy độ phân tán quy trình kiểm sốt, vẽ biện luận biểu đồ kiểm soát X - Thơng thường người ta chọn giới hạn kiểm sốt vị trí σ σ cách đường trung tâm - Nếu biết giá trị trung bình µ độ lệch chuẩn σ quy trình việc tính tốn vẽ biểu đồ Shewhart đơn giản Bảng Dạng tổng quát biều đồ Shewhart - µ, σ giá trị trung bình độ lệch chuẩn quy trình, n cỡ mẫu Trang MỤC TIÊU NGHIÊN CỨU 3.1 Mục tiêu tổng qt Xây dựng phần mềm kiểm sốt khối lượng viên/gói ba sản phẩm thuộc ba dạng bào chế: viên nang Cephalexin 500 mg, viên nén cefuroxim 500 mg, thuốc gói Cefaclor 125 mg Cơng ty CPDP IMEXPHARM biểu đồ Shewhart qua hệ thống cân online 3.2 Mục tiêu cụ thể - Xây dựng sở liệu - Khảo sát khối lượng trung bình viên tính tốn thủ cơng - Xây dựng phần mềm - Kiểm tra hoạt động phần mềm - Nối mạng cân Online Trang ĐỐI TƯỢNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 3.1 ĐỐI TƯỢNG 3.1.1 Ba sản phẩm thuộc dạng bào chế khác sản xuất Công ty CPDP IMEXPHARM: - Viên nang Cephalexin 500 mg - Viên nén Cefuroxim 500 mg - Thuốc gói Cefaclor 125 mg 3.1.2 Phần mềm theo Shewhart PHƯƠNG PHÁP NGHIÊN CỨU 3.2.1 Xây dựng sở liệu Tiến hành giai đoạn xây dựng áp dụng biểu đồ: - Trong thực tế giá trị trung bình µ độ lệch chuẩn σ quy trình sản xuất thường khơng biết nên chúng ước tính từ cỡ mẫu thực nghiệm Việc thiết lập biểu đồ theo bước tiến hành sau: + Thu thập 20 mẫu có cỡ mẫu n + Kiểm tra loại giá trị bất thường thử nghiệm DIXON + Tính giá trị trung bình Xtb khoảng quan sát R mẫu + Tính giá trị trung bình Xtb giá trị trung bình khoảng quan sát R mẫu - Vẽ đường trung tâm giới hạn 1σ, 2σ, 3σ Trang Chuẩn bò liệu Khơng đạt Loại giá trò khác thường Đánh giá phân phối chuẩn Đạt Khảo sát biểu ñoà Shewhart Đánh giá phần mềm Nối mạng cân online 3.2.2 Kiểm tra loại bỏ giá trị bất thường thử nghiệm Dixon - Giá trị bất thường mẫu giá trị lớn hay nhỏ Nếu khơng bị loại, giá trị bất thường làm sai lệch kết phân tích thống kê - Cơ sở tốn thơng kê: - Đặt giả thuyết: H0: X1 không bất thường HA: X1 bất thường - Phân tích thống kê: - Chuỗi liệu: Các giá trị thống kê xếp thành chuỗi liệu thừ nhỏ đến lớn hay từ lớn đến nhỏ X1 X2 X3 X4 ………… X5 XN - XN - XN Giá trị thống kê Q tính tùy theo cỡ mẫu N: Q = ( X2 - X1)/ ( XN - X1) N=3–7 ( X2 - X1)/ ( XN-1 - X1) N = – 12 ( X3 - X1)/ ( XN-2 - X1) N = 13 – 40 Biện luận kết Nếu Q > Q α (tra bảng) => Bác bỏ H0 => Giá trị bất thường Trang 3.2.3 Phân phối chuẩn Một quần thể gọi có tính phân phối chuẩn phân phối giá trị chúng tuân theo đường congcủa quy luật chuẩn (dạng hình chng) Nếu gọi µ σ giá trị trung bình độ lệch chuẩn dân số gồm N giá trị quan sát Xi Giá trị Z biến số ngẫu nhiên cho biểu thức: Z= ( X i − µ) σ Phương trình phân phối chuẩn có dạng hình chng: −z2 / e 2π f (Z ) = Diện tích đường cong phân phối chuẩn từ - ∞ đến Z0 (bên trái): ∏ Z0 = ∫ −∞ −Z / e dZ 2π Kính nhờ thầy kiệt vẽ dùm hình chng (khơng biết vẽ nào, mong thầy giúp) Hình Diện tích đường cong ∏ o phân phối chuẩn từ - ∞ đến Z0 Diện tích đường cong phân phối chuẩn từ Z1 đến Z2 : ∏ (Z < Z < Z ) = ∏ (−∞ < Z < Z ) - ∏ (−∞ < Z < Z ) ∏ (−∞ < Z < Zi ) tìm bảng diện tích bên trái ( ∏ i ) phân phối chuẩn (tra 2 bảng) Lý luận để tính tóan đường giới hạn kiểm sóat biểu đồ Shewhart tính số hiệu quy trình thíêt lập giả thuyết dân số khảo sát có phân phối chuẩn Trong thực tế, phải sử dụng liệu từ dân số mẫu để đến kết luận cho quy trình mà phân bố tính chất đầu quy trình chưa biết trước Để không bị sai lệch việc biện luận kết quả, yêu cầu phải kiểm tra tính phân phối chuẩn dân số mẫu Khi khảo sát dân số hay mẫu lớn, việc quan sát giá trị cá thể thời gian, vậy: Các giá trị gộp lại thành nhóm gọi lớp Số giá trị quan sát lớp gọi tần số thực nghiệm lớp Tòan quan sát xếp bảng tần số thực nghiệm Trang Nhật đồ tần số biểu diễn phân bố tần số theo lớp Vẽ biểu đồ mật độ phân bố cách tính: - Số lớp biểu đồ k = (1 + 10logN)/3, với N: Số giá trị quần thể mẫu - Khoảng rộng lớp D = (Xmax – Xmin)/k - Giới hạn lớp L1 = Xmin - 1/2 độ phân giải phép đo - Giới hạn lớp U1 = L1+ D - Giới hạn lớp 2: L2 = L1+ D - Giới hạn lớp 2: U2 = L2 + D = U1+ D - Tiếp tục tính giới hạn dưới/ lớp 3,……và k tương tự - Tính giới hạn lớp tương tự - Từ giới hạn lớp, công cụ Excel vẽ biểu đồ mật độ phân bố cung cấp tần suất giá trị theo lớp (giá trị Ni thử nghiệm χ2) Hình Minh họa nhật đồ tần số Cơ sở tóan thống kê Đặt giả thuýêt: Ho: P1 = ∏ ; P2= ∏ ;… ; Pk = HA : trường hợp Pi ≠ ∏k “thuộc phân phối chuẩn” ∏ i ” không thuộc phân phối chuẩn” Giả thuýêt Ho tức tần số thực nghiệm không khác với tần số lý thuýêt Phân tích thống kê: χ2=  (Oi − Ei )   Ei i =1   k ∑ Trang Oi: tần số thực nghiệm Ei: tần số lý thuyết = ∏ i N ∏ i : tỷ lệ lý thuyết (tra bảng diện tích bên trái ( ∏ i ) phân phố chuẩn) N: tổng số giá trị Xi χ α2 = CHIINV (α , γ ) γ = (K-1) γ lý thuýêt tìm bảng γ ứng với bậc tự γ = (k-1) theo ngưỡng α Biện luận kết quả: Nếu γ > χα bác bỏ Ho ngược lại 3.2.4 Khảo sát biểu đồ Shewhart - Giá trị trung bình quy trình (mức trung tâm) xem ổn định cho dù giá trị trung bình mẫu biến thiên xung quanh - Đối với biểu đồ Shewhart X việc biện luận dựa quy tắc nêu bên Trong ý điểm nằm ngồi hay áp sát gới hạn trên/dưới, dịch chuyển vùng giá trị cao/thấp, hay lên/xuống - Đối với biểu đồ Shewhart R, việc biện luận dựa nguyên tắc từ 1-4 song ý điểm nằm hay áp sát giới hạn trên, dịch chuyển vùng giá trị cao hay lên Bảng Biểu thức tính tốn đường trung tâm đường giới hạn Các đường quy ước Giới hạn +3 σ Giới hạn +2 σ Giá trị bình X + A2 R trung Khoảng quan sát D4 R X + 2/3 A2 R X + 1/3 A2 R (1+2D4)/3* R Đường trung tâm X R Giới hạn -1 σ X - 1/3 A2 R X - 2/3 A2 R (2+D3)/3* R Giới hạn +1 σ Giới hạn -2 σ (2+D4)/3* R (1+2D3)/3* R Trang 10 Ghi X : ước lượng µ A2 R : ước lượng 3σ* n D4* R : ước lượng R + 3σ D3* R : ước lượng R - 3σ A2, D3,D4: Tra bảng Giới hạn -3 σ X - A2 R D3* R - X + A2 R - X + 2/3 A2 R - X + 1/3 A2 R - X X - 1/3 A2 R X - 2/3 A2 R X - A2 R Hình Minh họa biểu đồ Shewhart X (n ≥ 2) với mức giới hạn ước tính -D4 R -(1+2D4)/3* R -(2+D4)/3* R - R (2+D3)/3* R (1+2D3)/3* R D3* R 3.2.3.1 Khảo sát khối lượng trung bình viên tính tốn thủ cơng Trang 11 Tiến hành thu thập liệu khối lượng trung bình lô liên tiếp dạng bào chế - Thời điểm lấy mẫu: 20 phút trình - Số mẫu lấy: mẫu/lần - Cỡ mẫu: 10 viên - Cách lấy mẫu: ngẫu nhiên Phương pháp: - Sản xuất theo QTSX số… - Cân ghi lại liệu khối lượng viên - Phân tích liệu, vẽ biểu đồ Shewhart X biểu đồ R - Đánh giá theo bảng ( riêng biểu đồ R đánh giá giới hạn trên) Bảng 3: Các quy tắc biện luận Quy tắc 1a Giới hạn +3σ 1b Giới hạn - 3σ 2a Vùng giá trị cao 3a Đi lên 3b Đi xuống 4a Giới hạn +3σ u cầu (khơng được) điểm ngồi giới hạn +3σ điểm giới hạn - 3σ Liên tiếp 6,7,8,9….điểm nằm X Liên tiếp 6,7,8,9….điểm nằm X Liên tiếp 6,7,8,9….điểm lên Liên tiếp 6,7,8,9….điểm xuống vùng +A; vùng +A; 4b Giới hạn - 3σ 10 vùng +A; vùng +B; vùng - A; vùng - A; 2b Vùng giá trị thấp 10 vùng - A; vùng - B; 2.2 Xây dựng phần mềm Yêu cầu: - Lập bảng liệu (số mẫu tối thiểu 20 lần, cỡ mẫu) theo yêu cầu thời gian thực tế quy trình Trang 12 - Nối mạng truyền liệu khối lượng mẫu viên 10 viên từ cân máy tính - Dữ liệu mẫu liệt kê Max, Min, Average, RSD - Tính tích lũy liệu đánh giá dựa quy tắc bảng Phần mềm dừng lại vi phạm quy tắc Quy trình phép vận hành điều chỉnh lại - Hiện khởi đầu quy trình tần số kiểm tra khối lượng viên, có chế độ hẹn chuông báo đến cân kiểm tra 2.3 Đánh giá kết biều đồ R, Xtb: Đạt, không đạt (theo quy tắc nào) Đánh giá phần mềm TRIỂN VỌNG ĐỀ TÀI Kết nghiên cứu làm sở khoa học để nghiên cứu triển khai áp dụng thực tế cho nhu cầu quản lý trình cơng ty CPDP IMEXPHARM biểu đồ Shewhart thông qua hệ thống cân Online NƠI THỰC HIỆN ĐỀ TÀI − Bộ môn Bào chế - Khoa Dược – ĐH Y Dược TP HCM − Bộ môn vi tính dược - Khoa Dược – ĐH Y Dược TP HCM − Công ty CPDP IMEXPHARM DỰ TRÙ KINH PHÍ – TRANG THIẾT BỊ - NGUYÊN LIỆU 4.1 Dự trù kinh phí − Kinh phí dự kiến: 20.000.000đ Trang 13 − Nguồn kinh phí: Tự túc 4.2 Trang thiết bị STT Tên thiết bị Cân điện tử số lẻ Nhà sản xuất Model Ohaus Corp Pine OHAUS Brook KJ, USA Ghi Item: AR 2130 Max cap: 210 g Readability: 0,001g SN: 1226250276 Bộ điều khiển kết nối Máy tính Phần mềm Shewhart Máy đóng nang Máy ép gói Bosch China Volpak GFK 700S SL 110 packaging Máy dập viên machines, Spain CILT JEMKAY ENGINEERS PVT.LTD Hệ thống mạng KẾ HOẠCH THỰC HIỆN STT CÁC PP SỬ THỜI GIAN KẾT DỤNG THỰC HIỆN QUẢ CÁC HOẠT ĐỘNG NGHIÊN CỨU Khảo sát biểu đồ kiểm soát khối lượng Tham khảo 12/2009 trung bình viên theo Shewhart tài liệu liên phương pháp thủ công, bao gồm nội quan Trang 14 dung: Thu thập - xử lý liệu: loại sai số thơ, khảo sát tính phân phối chuẩn, khảo sát vẽ biểu đồ Shewhart R , khảo sát vẽ biểu đồ Shewhart X Viết phần mềm theo Shewhart, nối mạng 02/2010 nghiên cứu truyền liệu máy tính – link xử lý nhờ phần mềm Áp dụng thử liệu thực tế sản xuất 04/2009 dạng bào chế khác nhau, dạng bào chế khảo sát tối thiểu lô Đánh giá phần mềm Viết hồn chỉnh luận văn 05/2009 06/2010 PHỤ ĐÍNH Hình Minh họa biểu đồ Shewhart R (n ≥ 2) với mức giới hạn ước tính A2, D3, D4 (tra bảng) Các quy tắc biện luận Quy tắc 1a Có điểm giới hạn +3σ Trang 15 685 680 675 670 665 660 655 650 645 10 11 12 13 14 15 16 17 18 19 20 Quy tắc 1b Có điểm ngồi giới hạn - 3σ 685 680 675 670 665 660 655 650 645 10 11 12 13 14 15 16 17 18 19 20 Quy tắc 2a Tiến vùng giá trị cao: Liên tiếp 6,7,8,9….điểm nằm X Trang 16 Quy tắc 2b Tiến vùng giá trị thấp: Liên tiếp 6,7,8,9….điểm nằm X Quy tắc 3a Liên tiếp 6,7,8,9….điểm lên Trang 17 Quy tắc 3b Liên tiếp 6,7,8,9….điểm xuống Quy tắc 4a: vùng +A; vùng +A; 10 vùng +A; vùng +B Trang 18 Quy tắc 4b: vùng - A; vùng - A; 10 vùng - A; vùng – B Ý KIẾN THẦY HƯỚNG DẪN Trang 19 TÀI LIỆU THAM KHẢO Trang 20 Tiếng Việt Bộ Y tế (2005), Hướng dẫn thực hành tốt sản xuất thuốc theo khuyến cáo tổ chức y tế giới, Nhà xuất y học, tr 309 Đặng Văn Giáp (2003), Phân tích thống kê kiểm nghiệm dược phẩm Đặng Văn Giáp (1997), Phân tích liệu khoa học chương trình Ms-Excel, NXB Giáo dục, tr 63-66 Hiệp hội Tu nghiệp Kỹ thuật Hải ngoại Nhật ( AOTS ), Bảy công cụ QC, tr.5 – Huỳnh Văn Hóa, Đặng Văn Giáp, Hồng Minh Châu (2002), Thẩm định lui quy trình sản xuất viên nén Paracetamol, Tạp chí dược học số 6 Huỳnh Văn Hóa, Chung Khang Kiệt, Đỗ Quang Dương, Đặng Văn Giáp, Hoàng Minh Châu (2003), Thành lập phần mềm thống kê Pharmasoft TĐ phục vụ cho việc kiểm sốt quy trình sản xuất thuốc, Tạp chí Dược học, tr 22-25 Huỳnh Văn Hóa (2005), Xây dựng số mơ hình thẩm định quy trình sản xuất thuốc Việt nam, luận án tiến sĩ dược học Tiếng Anh Pramote Cholayudth, Concepts and practices of Pharmaceutical process validation, Statistical process control for QA, Chapter Trang 21 Trang 22 ... giá trị trung bình độ lệch chuẩn quy trình, n cỡ mẫu Trang MỤC TIÊU NGHIÊN CỨU 3.1 Mục tiêu tổng quát Xây dựng phần mềm kiểm soát khối lượng viên/ gói ba sản phẩm thuộc ba dạng bào chế: viên nang... lỗi Các giới hạn kiểm soát cảnh báo biểu đồ giúp phòng ngừa sai lỗi, ví dụ: sản phẩm vượt ngồi giới hạn kiểm sốt Cung cấp thơng tin hiệu q trình Dữ liệu sử dụng cho biểu đồ kiểm soát hữu ích cho... giới hạn để vẽ biểu đồ mật độ phân bố Để không dẫn tới kết sai lệch, việc tính tốn số liệu phải tn thủ theo quy định 2.1 VẼ BIỂU ĐỒ SHEWHART X VÀ R - Các biểu đồ Shewhart X R thường sử dụng theo