Tiêu chuẩn Việt Nam TCVN 6910-3:2001 quy định bốn thước đo trung gian độ chụm dựa theo những thay đổi về điều kiện quan trắc trong một phòng thí nghiệm (thời gian, hiệu chuẩn, người thao tác và thiết bị). Có thể xác định những thước đo trung gian này bằng thí nghiệm trong một phòng thí nghiệm cụ thể hoặc bằng thí nghiệm so sánh liên phòng.
TIÊU CHUẨN VIỆT NAM TCVN 6910-3 : 2001 ISO 5725-3 : 1994 ĐỘ CHÍNH XÁC (ĐỘ ĐÚNG VÀ ĐỘ CHỤM) CỦA PHƯƠNG PHÁP ĐO VÀ KẾT QUẢ ĐO – PHẦN 3: CÁC THƯỚC ĐO TRUNG GIAN ĐỘ CHỤM CỦA PHƯƠNG PHÁP ĐO TIÊU CHUẨN Accuracy (trueness and precision) of measurement methods and results – Part 3: Intermediate measures of the precision of a standard measurement method Lời nói đầu TCVN 6910-3 : 2001 hoàn toàn tương đương với ISO 5725-3 : 1994 Các phụ lục A, B, C tiêu chuẩn quy định, phụ lục D, E để tham khảo TCVN 6910-3 : 2001 Tiểu ban Kỹ thuật Tiêu chuẩn TCVN/ TC69/SC6 Phương pháp Kết đo biên soạn, Tổng cục Tiêu chuẩn Đo lường Chất lượng đề nghị, Bộ Khoa học Công nghệ Môi trường ban hành Lời giới thiệu 0.0 TCVN 6910-3 : 2001 phần TCVN 6910, tiêu chuẩn gồm phần tên chung “Độ xác (độ độ chụm) phương pháp đo kết đo”: - Phần 1: Nguyên tắc định nghĩa chung - Phần 2: Phương pháp xác định độ lặp lại độ tái lập phương pháp đo tiêu chuẩn - Phần 3: Các thước đo trung gian độ chụm phương pháp đo tiêu chuẩn - Phần 4: Các phương pháp xác định độ phương pháp đo tiêu chuẩn - Phần 5: Các phương pháp khác xác định độ chụm phương pháp đo tiêu chuẩn - Phần 6: Sử dụng giá trị độ xác thực tế 0.1 TCVN 6910 sử dụng hai thuật ngữ “độ đúng” “độ chụm” để diễn tả độ xác phương pháp đo “Độ đúng” gần trung bình số học số lớn kết thử nghiệm giá trị thực giá trị quy chiếu chấp nhận “Độ chụm” gần kết thử nghiệm 0.2 Việc xem xét chung đại lượng trình bày TCVN 6910-1 nên khơng nhắc lại Cần nhấn mạnh phải đọc TCVN 6910-1 với phần khác TCVN 6910, định nghĩa nguyên tắc chung quan trọng trình bày tiêu chuẩn 0.3 Nhiều yếu tố khác (không kể thay đổi mẫu thử xem giống hệt nhau) đóng góp vào thay đổi kết phương pháp đo, bao gồm: a) người thao tác; b) thiết bị sử dụng; c) việc hiệu chuẩn thiết bị; d) môi trường (nhiệt độ, độ ẩm, ô nhiễm khơng khí…) e) mẻ thuốc thử f) khoảng thời gian phép đo Sự thay đổi phép đo người thao tác khác thực và/ thiết bị khác thường lớn thay đổi phép đo thực khoảng thời gian ngắn với người thao tác, sử dụng thiết bị 0.4 Trong nhiều trường hợp thực tế, để thuận tiện cho việc mô tả thay đổi phương pháp đo, cần đưa hai điều kiện độ chụm gọi điều kiện lặp lại điều kiện tái lập Theo điều kiện lặp lại yếu tố a) đến f) thuộc mục 0.3 xem không đổi khơng đóng góp vào thay đổi theo điều kiện tái lập yếu tố thay đổi đóng góp vào thay đổi kết thử nghiệm Như điều kiện lặp lại tái lập hai cực độ chụm Điều kiện lặp lại mô tả khả thay đổi tối thiểu điều kiện tái lập mô tả thay đổi tối đa kết Cũng hình dung điều kiện trung gian hai điều kiện có nhiều yếu tố a) đến f) biến đổi sử dụng tình định Độ chụm thường biểu thị qua độ lệch chuẩn 0.5 Tiêu chuẩn tập trung vào thước đo trung gian độ chụm phương pháp đo Gọi thước đo trung gian độ lớn chúng nằm hai cực độ chụm phương pháp đo: độ lệch chuẩn lặp lại độ lệch chuẩn tái lập Để minh họa cần thiết thước đo trung gian độ chụm xem xét hoạt động phòng thí nghiệm gắn liền với q trình sản xuất, ví dụ hệ thống làm việc ba ca, người thao tác khác thực phép đo thiết bị khác Người thao tác, thiết bị nhiều yếu tố khác tạo thay đổi kết thử nghiệm Những yếu tố cần phải tính đến đánh giá độ chụm phương pháp đo 0.6 Thước đo trung gian độ chụm định nghĩa phần tiêu chuẩn có ích chủ yếu việc ước lượng chúng thành phần quy trình nhằm phát triển, tiêu chuẩn hóa kiểm sốt phương pháp đo phạm vi phòng thí nghiệm Cũng ước lượng phép đo nghiên cứu so sánh liên phòng thí nghiệm thiết kế đặc biệt, việc diễn giải áp dụng chúng đòi hỏi thận trọng lý giải thích 1.3 9.1 0.7 Bốn yếu tố gần ảnh hưởng nhiều đến độ chụm phương pháp đo là: a) Thời gian: khoảng thời gian phép đo ngắn dài b) Hiệu chuẩn: thiết bị chưa hiệu chuẩn lại nhóm phép đo c) Người thao tác: người thao tác nhiều người thao tác khác thực phép đo d) Thiết bị: thiết bị nhiều thiết bị khác (mẻ thuốc thử nhiều mẻ thuốc thử khác nhau) sử dụng phép đo 0.8 Để tiện lợi cho việc xem xét thay đổi điều kiện đo phòng thí nghiệm (thời gian, hiệu chuẩn, người thao tác thiết bị), người ta đưa điều kiện chụm gian M yếu tố khác sau: a) M = 1: có bốn yếu tố khác nhau; b) M = 2: hai bốn yếu tố khác nhau; c) M = 3: ba bốn yếu tố khác nhau; d) M = 4: bốn yếu tố khác Điều kiện chụm trung gian khác dẫn đến việc độ lệch chuẩn trung gian khác độ chụm ký hiệu SI( ); điều kiện cụ thể liệt kê dấu ngoặc Ví dụ S I (TO) độ lệch chuẩn trung gian độ chụm với thời gian (T) người thao tác (O) khác 0.9 Đối với phép đo thực điều kiện chụm trung gian nhiều yếu tố liệt kê mục 0.7 khác Với điều kiện lặp lại, yếu tố giả thiết không đổi Độ lệch chuẩn kết thử nghiệm nhận điều kiện lặp lại thường nhỏ độ lệch chuẩn kết thử nghiệm nhận điều kiện chụm trung gian Nói chung phân tích hóa học độ lệch chuẩn điều kiện chụm trung gian lớn hai đến ba lần độ lệch chuẩn điều kiện lặp lại Tất nhiên độ lệch khơng thể vượt q độ lệch chuẩn tái lập Ví dụ xác định lượng đồng quặng đồng, thí nghiệm phối hợp 35 phòng thí nghiệm phát độ lệch chuẩn điều kiện chụm trung gian với yếu tố khác (cùng người thao tác thiết bị thời gian đo khác nhau) lớn 1,5 lần so với độ lệch chuẩn điều kiện lặp lại Điều hai phương pháp: đo trọng lực điện phân chuẩn độ Na2S2O3 ĐỘ CHÍNH XÁC (ĐỘ ĐÚNG VÀ ĐỘ CHỤM) CỦA PHƯƠNG PHÁP ĐO VÀ KẾT QUẢ ĐO – PHẦN 3: CÁC THƯỚC ĐO TRUNG GIAN ĐỘ CHỤM CỦA PHƯƠNG PHÁP ĐO TIÊU CHUẨN Accuracy (trueness and precision) of measurement methods and results – Part 3: Intermediate measures of the precision of a standard measurement method Phạm vi áp dụng 1.1 Tiêu chuẩn quy định bốn thước đo trung gian độ chụm dựa theo thay đổi điều kiện quan trắc phòng thí nghiệm (thời gian, hiệu chuẩn, người thao tác thiết bị) Có thể xác định thước đo trung gian thí nghiệm phòng thí nghiệm cụ thể thí nghiệm so sánh liên phòng Ngồi ra, tiêu chuẩn còn: a) Thảo luận hàm ý định nghĩa thước đo trung gian độ chụm; b) Diễn giải ứng dụng việc ước lượng thước đo trung gian độ chụm tình thực tế; c) Không đưa thước đo sai số việc ước lượng thước đo trung gian độ chụm; d) Không đề cập đến việc xác định độ thân phương pháp đo có thảo luận mối liên hệ độ điều kiện đo 1.2 Tiêu chuẩn liên quan tới phương pháp đo tạo phép đo thang đo liên tục cho giá trị đơn kết thử nghiệm giá trị đơn kết tính tốn từ tập hợp quan trắc 1.3 Thực chất việc xác định thước đo trung gian độ chụm đánh giá khả phương pháp đo lặp lại kết phép thử điều kiện xác định 1.4 Các phương pháp thống kê thực tiêu chuẩn dựa vào giả thuyết tập hợp thông tin từ điều kiện đo “tương tự” để nhận thơng tin xác thước đo trung gian độ chụm Giả thuyết có hiệu lực yêu cầu “tương tự” phải “tương tự” Nhưng khó khăn để giữ vững giả thuyết thước đo trung gian độ chụm ước lượng từ nghiên cứu so sánh liên phòng Ví dụ, khó khăn để kiểm soát ảnh hưởng “thời gian” “người thao tác” qua nhiều phòng nghiệm cách làm cho chúng “tương tự” để việc tập hợp thơng tin từ phòng thí nghiệm khác có ý nghĩa Như yêu cầu phải thận trọng sử dụng kết nghiên cứu so sánh liên phòng thước đo trung gian độ chụm Những nghiên cứu phạm vi phòng thí nghiệm phải dựa vào giả thuyết này, nghiên cứu thực tế nhà phân tích kiểm sốt hiểu biết rõ tác động thực yếu tố 1.5 Ngồi kỹ thuật mơ tả tiêu chuẩn để ước lượng xác nhận thước đo trung gian độ chụm phòng thí nghiệm có kỹ thuật khác, ví dụ biểu đồ kiểm sốt (xem TCVN 6910-6) Tiêu chuẩn khơng đòi hỏi mơ tả cách tiếp cận để ước lượng thước đo trung gian độ chụm phạm vi phòng thí nghiệm cụ thể Chú thích – Tiêu chuẩn đề cập tới thiết kế thí nghiệm thiết kế xếp Một số thông tin cho phụ lục B C Những tài liệu tham khảo khác thuộc lĩnh vực cho phụ lục E Tiêu chuẩn viện dẫn ISO 3534-1 : 1993, Thống kê học – Từ vựng ký hiệu – Phần 1: Thuật ngữ xác suất thống kê đại cương TCVN 6910-1 : 2001, Độ xác (độ độ chụm) phương pháp đo kết đo – Phần 1: Những nguyên tắc định nghĩa chung TCVN 6910-2 : 2001, Độ xác (độ độ chụm) phương pháp đo kết đo – Phần 2: Phương pháp xác định độ lặp lại độ tái lập phương pháp đo tiêu chuẩn ISO Guide 33 : 1989, Sử dụng mẫu chuẩn chứng nhận ISO Guide 35 : 1989, Chứng nhận mẫu chuẩn Những nguyên tắc chung thống kê Định nghĩa Trong tiêu chuẩn áp dụng định nghĩa cho ISO 3534-1 TCVN 6910-1 Các ký hiệu sử dụng TCVN 6910 cho Phụ lục A Yêu cầu chung Để phép đo thực nhau, phương pháp đo phải tiêu chuẩn hóa Tất phép đo phần hợp thành thí nghiệm phòng thí nghiệm cụ thể thí nghiệm so sánh liên phòng phải tiến hành theo tiêu chuẩn Các yếu tố quan trọng 5.1 Bốn yếu tố (thời gian, hiệu chuẩn, người thao tác thiết bị) thuộc điều kiện đo phòng thí nghiệm coi thành phần tạo khả thay đổi phép đo (xem bảng 1) 5.2 “Những phép đo tiến hành thời gian” phép đo thực thời gian ngắn nhằm giảm thiểu thay đổi điều kiện đo, điều kiện môi trường, điều kiện đảm bảo luôn không đổi “Những phép đo tiến hành thời gian khác nhau” phép đo thực khoảng thời gian dài, phép đo chứa ảnh hưởng thay đổi điều kiện môi trường Bảng – Bốn yếu tố quan trọng trạng thái chúng Yếu tố Điều kiện phòng thí nghiệm Trạng thái (Giống nhau) Trạng thái (Khác nhau) Thời gian Các phép đo tiến hành thời gian Các phép đo tiến hành thời gian khác Hiệu chuẩn Giữa phép đo không tiến hành hiệu chuẩn Giữa phép đo có tiến hành hiệu chuẩn Người thao tác Cùng người thao tác Người thao tác khác Thiết bị Cùng thiết bị chưa hiệu chuẩn lại Thiết bị khác 5.3 “Hiệu chuẩn” không đề cập tới phép hiệu chuẩn mà phương pháp đo yêu cầu phần cần có phương pháp để nhận kết thử nghiệm Ở nói đến q trình hiệu chuẩn thực vào khoảng thời gian đặn nhóm phép đo phòng thí nghiệm 5.4 Trong vài thao tác “người thao tác” thực tế nhóm người thao tác, thành viên nhóm thực phần việc cụ thể q trình Trong trường hợp vậy, nhóm phải xem người thao tác thay đổi thành viên phân cơng nhiệm vụ nhóm phải xem có người thao tác khác 5.5 Trong thực tế “thiết bị” thường thiết bị thay đổi thành phần quan trọng thiết bị phải coi thiết bị khác Phải định rõ thành phần quan trọng thiết bị Sự thay đổi nhiệt kế coi thành phần quan trọng đáng kể, việc sử dụng bình khác đơi chút khả chứa nước coi bình thường Sự thay đổi mẻ thuốc thử phải xem thành phần quan trọng đáng kể Sự thay đổi dẫn đến “thiết bị” khác hiệu chuẩn lại thay đổi kéo theo việc hiệu chuẩn 5.6 Trong điều kiện lặp lại tất bốn yếu tố trạng thái bảng Đối với điều kiện chụm trung gian có nhiều yếu tố trạng thái bảng chúng xác định “điều kiện chụm với M yếu tố khác nhau”, M số yếu tố trạng thái Trong điều kiện tái lập kết thu từ phòng thí nghiệm khác không tất bốn yếu tố trạng thái mà có ảnh hưởng phụ khác phòng thí nghiệm quản lý trì phòng thí nghiệm, trình độ đào tạo chung người thao tác, ổn định kiểm tra kết thử nghiệm v.v… 5.7 Trong điều kiện độ chụm trung gian với M yếu tố khác cần phải xác định yếu tố trạng thái bảng số Ví dụ: - Độ lệch chuẩn độ chụm trung gian khác thời gian, S I(T); - Độ lệch chuẩn độ chụm trung gian khác hiệu chuẩn, S I(C); - Độ lệch chuẩn độ chụm trung gian khác người thao tác, S I(O); - Độ lệch chuẩn độ chụm trung gian khác thời gian + người thao tác, S I(TO); - Độ lệch chuẩn độ chụm trung gian khác thời gian + người thao tác + thiết bị, SI(TOE); - Và nhiều trường hợp khác theo cách đánh số tương tự Mơ hình thống kê 6.1 Mơ hình sở Để đánh giá độ xác (độ độ chụm) phương pháp đo, thường giả thiết kết thử nghiệm y tổng thành phần: y=m+B+e … (1) đó, với vật liệu cụ thể thử nghiệm: m trung bình chung (kỳ vọng) B thành phần phòng thí nghiệm độ chệch điều kiện lặp lại e sai số ngẫu nhiên xuất phép đo điều kiện lặp lại Dưới đề cập tới thành phần việc mở rộng mơ hình sở 6.2 Trung bình chung, m 6.2.1 Trung bình chung, m, trung bình lấy tồn kết thí nghiệm Giá trị m thu nghiên cứu phối hợp (xem TCVN 6910-2) phụ thuộc vào “giá trị thực” phương pháp đo Nó khơng phụ thuộc vào phòng thí nghiệm, thiết bị, người thao tác thời gian kết thử nghiệm thu Trung bình chung vật liệu cụ thể đo gọi “mức phép thử” Ví dụ: mẫu với độ tinh khiết khác chất hóa học chất liệu khác (ví dụ loại thép khác nhau) tương ứng với mức thử khác Trong nhiều trường hợp, khái niệm giá trị thực µ áp dụng, chẳng hạn nồng độ thực dung dịch chuẩn độ Mức thử m thường không giá trị thực (µ) Hiệu (m - µ) gọi “độ chệch” phương pháp đo Trong số trường hợp, mức thử định nghĩa riêng phương pháp đo khái niệm giá trị thực độc lập khơng áp dụng Ví dụ độ cứng Vickers thép số Micum than cốc thuộc phạm trù Song cách tổng quát, độ chệch ký hiệu δ (δ = không tồn giá trị thực) trung bình chung m là: m=µ+δ Chú thích – Thảo luận độ chệch TCVN 6910-4 …(2) mơ tả thí nghiệm độ cho 6.2.2 Khi khảo sát khác kết thử nghiệm thu phương pháp đo, độ chệch phương pháp khơng có ảnh hưởng bỏ qua trừ độ chệch hàm mức thử Độ chệch phương pháp đo cần tính đến so sánh kết thử nghiệm với giá trị quy định hợp đồng với giá trị tiêu chuẩn hợp đồng quy định kỹ thuật đề cập tới giá trị thực µ, không đề cập tới mức thử m, so sánh kết thử nghiệm thu cách sử dụng nhiều phương pháp đo khác 6.3 Số hạng B 6.3.1 B hạng biểu thị độ lệch phòng thí nghiệm so với m nhiều nguyên nhân gây ra, không kể sai số ngẫu nhiên xuất kết thử nghiệm Trong điều kiện lặp lại phòng thí nghiệm, B xem không đổi gọi “thành phần phòng thí nghiệm độ chệch” 6.3.2 Tuy nhiên sử dụng thường xuyên phương pháp đo, điều thấy rõ số lượng lớn ảnh hưởng thể giá trị chung B Những ảnh hưởng do, ví dụ thay đổi người thao tác, thiết bị sử dụng việc hiệu chuẩn thiết bị môi trường (nhiệt độ, độ ẩm, nhiễm khơng khí v.v…) gây Mơ hình thống kê (phương trình 1) viết lại dạng: y = m + B0 + B(1) + B(2) + … + e … (3) Hoặc y=µ+ + B0 + B(1) + B(2) + … + e (4) B hợp thành biến B0 , B(1) , B(2)… tính theo số yếu tố độ chụm trung gian Trong thực tế, mục tiêu việc nghiên cứu cân nhắc độ nhậy phương pháp đo định phạm vi sử dụng mơ hình Trong nhiều trường hợp sử dụng dạng rút gọn 6.4 Các số hạng B0 , B(1) , B(2)… 6.4.1 Trong điều kiện lặp lại, tất thành phần không đổi cộng thêm vào độ chệch kết thử nghiệm Trong điều kiện chụm trung gian, B0 ảnh hưởng cố định yếu tố lại giống (trạng thái bảng 1) B(1) , B(2)… ảnh hưởng ngẫu nhiên yếu tố biến đổi (trạng thái bảng 1) Các thành phẩn khơng góp phần vào độ chệch làm tăng độ lệch chuẩn chụm trung gian tới mức trở thành lớn so với độ lệch chuẩn lặp lại 6.4.2 Những ảnh hưởng khác biệt người thao tác bao gồm thói quen cá nhân việc thực phương pháp đo (ví dụ việc đọc giá trị độ chia thang đo v.v…) Một vài khác biệt phải loại trừ việc tiêu chuẩn hóa phương pháp đo, đặc biệt phải mơ tả xác, rõ ràng kỹ thuật có phương pháp Thậm chí độ chệch có kết phép thử thu từ cá nhân người thao tác khơng phải ln số (ví dụ độ lớn độ chệch thay đổi tùy theo trạng thái sinh lý người ngày hơm đó) khơng thể hiệu định cỡ sai lệch cách xác Phải giảm độ lớn sai lệch cách đào tạo sử dụng tài liệu hướng dẫn thao tác rõ ràng Trong hoàn cảnh ảnh hưởng thay đổi người thao tác xem mang chất ngẫu nhiên 6.4.3 Những ảnh hưởng khác biệt thiết bị gồm ảnh hưởng vị trí lắp đặt khác nhau, đặc biệt dao động thị v.v… Một vài ảnh hưởng khác biệt thiết bị hiệu việc hiệu chuẩn xác Những sai khác nguyên nhân hệ thống thiết bị phải hiệu hiệu chuẩn quy trình phải bao gồm phương pháp tiêu chuẩn Ví dụ thay đổi nhóm thuốc thử xử lý cách Cần có giá trị quy chiếu chấp nhận cần tham khảo ISO Guide 33 ISO Guide 35 cho công việc Ảnh hưởng lại thiết bị sau hiệu chuẩn mẫu chuẩn xem ảnh hưởng ngẫu nhiên 6.4.4 Những ảnh hưởng thời gian gây khác mơi trường, thay đổi nhiệt độ, độ ẩm phòng v.v… Phải cố gắng tiêu chuẩn hóa điều kiện môi trường nhằm giảm thiểu ảnh hưởng 6.4.5 Ảnh hưởng kỹ mệt mỏi người thao tác xem tác động qua lại người thao tác với thời gian Tính tập hợp thiết bị khác thời điểm bắt đầu sử dụng thời điểm sau nhiều ví dụ tác động qua lại thiết bị thời gian, Khi số người thao tác chí tập hợp thiết bị nhỏ ảnh hưởng gây nên yếu tố đánh ảnh hưởng cố định (không phải ngẫu nhiên) 6.4.6 Những quy trình cho TCVN 6910-2 xây dựng với giả thiết phân bố thành phần phòng thí nghiệm độ chệch gần chuẩn Nhưng thực tế, chúng thực với hầu hết phân bố phân bố đỉnh Phương sai B gọi “phương sai thành phần phòng thí nghiệm độ chệch” biểu thị bằng: var (B) = σ2L … (5) Đương nhiên phương sai chứa ảnh hưởng thay đổi người thao tác, thiết bị, thời gian môi trường Phương sai độ chụm trung gian tính từ thí nghiệm độ chụm sử dụng người thao tác, lần đo, mơi trường… khác thiết kế thí nghiệm xếp Phương sai var (B) xem bao gồm thành phần độc lập phòng thí nghiệm, người thao tác, ngày thí nghiệm, môi trường v.v… var (B) = var (B0) + var (B(1)) + var (B(2)) + … … (6) Các phương sai ký hiệu bởi: var (B0) = σ2(0) var (B(1)) = σ2(1) var (B(2)) = σ2(2) ,… … (7) Phương sai var (B) ước lượng nhờ thành phần thực tế S 2L ước lượng độ chụm trung gian tương tự thu từ thí nghiệm thiết kế thích hợp 6.5 Số hạng sai số, e 6.5.1 Thuật ngữ biểu thị sai số ngẫu nhiên xuất kết thử nghiệm quy trình đề cập tới toàn tiêu chuẩn xây dựng với giả thiết phân bố biến số sai số gần phân bố chuẩn thực tế quy trình thực cho hầu hết phân bố chúng phân bố đỉnh 6.5.2 Trong phòng thí nghiệm riêng rẽ, phương sai gọi phương sai phòng thí nghiệm biểu thị là: var (e) = σ2w .(8) 6.5.3 Có thể σ2w có giá trị khác phòng thí nghiệm khác có khác kỹ người thao tác, tiêu chuẩn điều giả thiết phương pháp đo tiêu chuẩn hóa cách hồn chỉnh khác phòng thí nghiệm phải nhỏ cho phép thiết lập giá trị chung phương sai phòng thí nghiệm cho tất các phòng thí nghiệm sử dụng phương pháp đo Giá trị chung ước lượng trung bình phương sai phòng thí nghiệm gọi “phương sai lặp lại”, ký hiệu bằng: Giá trị trung bình lấy tất phòng thí nghiệm tham gia thí nghiệm độ xác lại sau loại trừ giá trị bất thường Chọn điều kiện đo 7.1 Khi áp dụng phương pháp đo, nhiều điều kiện đo phòng thí nghiệm trùng như: a) điều kiện lặp lại (4 yếu tố không đổi); b) số điều kiện chụm trung gian với yếu tố khác nhau; c) số điều kiện chụm trung gian với hai yếu tố khác nhau; d) số điều kiện chụm trung gian với ba yếu tố khác nhau; e) điều kiện chụm trung gian với bốn yếu tố khác Trong tiêu chuẩn phương pháp đo không thiết phải cơng bố tất thước đo độ chụm có, độ lệch chuẩn lặp lại ln ln phải xác định Trong thực tế thương mại chung phải điều kiện thường gặp thước đo trung gian độ chụm chúng phải đủ để xác định thước đo trung gian độ chụm thích hợp với quy định chi tiết điều kiện đo cụ thể gắn liền với Yếu tố điều kiện đo thay đổi phải xác định cẩn thận Đặc biệt độ chụm trung gian thời gian khác nhau, khoảng thời gian thực tế trung bình lần đo liên tiếp phải xác định 7.2 Giả thiết phương pháp đo tiêu chuẩn sai lệch đơi chút sai lệch vốn có phương pháp phải xem xét phương tiện kỹ thuật Vì tiêu chuẩn nghiên cứu độ chệch phát sinh từ điều kiện đo 7.3 Sự thay đổi yếu tố điều kiện đo (thời gian, hiệu chuẩn, người thao tác thiết bị) so với điều kiện lặp lại (nghĩa từ trạng thái đến trạng thái bảng 1) làm tăng khả thay đổi kết thử nghiệm Tuy nhiên kỳ vọng giá trị trung bình số lượng kết thử nghiệm sai lệch điều kiện lặp lại Sự tăng lên độ lệch chuẩn điều kiện độ chụm trung gian khắc phục cách lấy kết trích dẫn cuối khơng phải kết thử nghiệm đơn mà trung bình số kết thử nghiệm 7.4 Trong đa số phòng thí nghiệm việc khảo sát thực tế khảo sát độ chụm (độ lệch chuẩn) mong muốn kết trích dẫn cuối chi phí thực phép đo định số yếu tố việc lựa chọn yếu tố mà thay đổi nghiên cứu tiêu chuẩn hóa phương pháp đo Nghiên cứu phân tích thước đo trung gian độ chụm phòng thí nghiệm 8.1 Cách tiếp cận đơn giản Phương pháp đơn giản để xác định độ lệch chuẩn độ chụm trung gian phòng thí nghiệm gồm việc lấy mẫu (hoặc thử nghiệm phá hủy tập hợp mẫu giả thiết giống nhau) thực loạt n phép đo với thay đổi yếu tố lần đo, n khuyến nghị phải 15 Điều khơng thỏa mãn phòng thí nghiệm phương pháp ước lượng thước đo trung gian độ chụm phòng thí nghiệm coi hiệu so với phương pháp khác Tuy nhiên việc phân tích đơn giản nên dùng phương pháp để nghiên cứu độ chụm trung gian với thời gian khác cách thực phép đo liên tiếp mẫu vào ngày liên tiếp nghiên cứu ảnh hưởng việc hiệu chuẩn lần đo Nên dùng đồ thị (yk – ) theo số k phép đo để nhận biết giá trị bất thường, yk kết thử nghiệm thứ k giá trị trung gian n lần kết thử nghiệm lặp lại Một phép thử thức giá trị bất thường bao gồm việc áp dụng phép thử Grubb cho 7.3.4 TCVN 6910-2 : 2001 Ước lượng độ lệch chuẩn trung gian độ chụm với M yếu tố khác cho bởi: Sl ( ) n n 1k (yk y) (10) Trong ký hiệu biểu thị điều kiện chụm trung gian phải để dấu ngoặc đơn 8.2 Phương pháp luân chuyển 8.2.1 Phương pháp luân chuyển khảo sát t nhóm phép đo, nhóm chứa n kết thử nghiệm lặp lại Ví dụ phạm vi phòng thí nghiệm, vật liệu thuộc tập hợp t vật liệu đo sau yếu tố chụm trung gian thay đổi t vật liệu đo lại, quy trình lặp lại có n kết thử nghiệm vật liệu thuộc t vật liệu Mỗi nhóm n kết thử nghiệm phải thu mẫu (hoặc tập hợp mẫu giả thiết giống hệt phép thử phá hủy) không thiết vật liệu phải giống hệt Yêu cầu tất t vật liệu nằm phạm vi mức thử xem xét để sử dụng giá trị độ lệch chuẩn trung gian độ chụm với M yếu tố khác Giá trị t (n – 1) khuyến nghị phải 15 Ví dụ: Một người thao tác thực phép đo đơn vật liệu thuộc t vật liệu, sau phép đo lặp lại người thao tác thứ hai lặp lại người thao tác thứ ba, thứ tư v.v… cho phép tính ước lượng sI(O) 8.2.2 Nên dùng đồ thị (yjk – j) theo số j vật liệu để nhận biết giá trị bất thường yjk kết thử nghiệm thứ k vật liệu thứ j ȳ j trung bình n kết vật liệu thứ j Một phép thử thức giá trị bất thường bao gồm việc áp dụng phép thử Grubb cho 7.3.4 TCVN 6910 – : 2001 tổ hợp lại cho nhóm riêng biệt cho tất tn kết thử nghiệm Ước lượng độ lệch chuẩn trung gian độ chụm với M yếu tố khác nhau, s I( ) cho công thức: Sl ( ) t (n 1) t m j 1k ( y jk y j )2 (11) Với n = (nghĩa có kết thử nghiệm vật liệu), công thức đơn giản là: Sl ( ) 2t t j ( y j1 y j2 ) (12) 8.3 Ảnh hưởng điều kiện đo kết trích dẫn cuối 8.3.1 Kỳ vọng y khác tổ hợp tổ hợp khác thời gian, hiệu chuẩn, người thao tác thiết bị, chí có bốn yếu tố thay đổi Đây điều hạn chế việc sử dụng giá trị trung bình Trong phân tích hóa học thí nghiệm vật lý, y thơng báo kết trích dẫn cuối Trong mua bán vật liệu thơ, kết trích dẫn cuối thường sử dụng để đánh giá chất lượng vật liệu thô ảnh hưởng đáng kể đến giá thành sản phẩm Ví dụ: Trong mua bán than thị trường quốc tế, sức chứa kho thường 70 000 hàm lượng tro xác định cuối mẫu thử g Theo hợp đồng quy định phần % tro tương ứng với giá 1,5 USD cho than Lượng mg khối lượng tro cân tro cân hóa học tương ứng với 0,1 % tro thành phần, nghĩa 0,15 USD cho than, cho kho thay đổi tương đương với 10 500 USD (0,1 x 1,5 x 70 000) 8.3.2 Do đó, kết trích dẫn cuối phân tích hóa học thí nghiệm vật lý phải có độ xác cần thiết, độ tin cậy cao, đặc biệt khả tái lập tính phổ qt Kết trích dẫn cuối đảm bảo điều kiện người thao tác xác định, thiết bị hay thời gian không đủ tốt xem xét thương mại Nghiên cứu phân tích liên phòng thí nghiệm thước đo trung gian độ chụm 9.1 Những giả thiết Việc ước lượng thước đo trung gian độ chụm từ nghiên cứu liên phòng dựa vào giả thiết ảnh hưởng yếu tố cụ thể tất phòng thí nghiệm Ví dụ, thay đổi người thao tác phòng thí nghiệm có ảnh hưởng thay đổi người thao tác phòng thí nghiệm khác, thay đổi thời gian tất phòng thí nghiệm Nếu giả thiết bị vi phạm khái niệm thước đo trung gian độ chụm biện pháp kỹ thuật đề xuất mục để ước lượng thước đo trung gian độ chụm khơng ý nghĩa Cần phải lưu ý thích đáng đến giá trị bất thường (không cần thiết phải loại bỏ giá trị này) điều giúp ích cho việc phát khác biệt với giả thiết để tập hợp thơng tin từ tất phòng thí nghiệm Một kỹ thuật tốt phát giá trị bất thường mơ tả đồ thị phép đo hàm mức khác yếu tố hàm nhiều phòng thí nghiệm khác tham gia vào nghiên cứu 9.2 Cách tiếp cận đơn giản Nếu vật liệu q mức gửi đến p phòng thí nghiệm, phòng thí nghiệm thực phép đo mức q mức thử có thay đổi yếu tố chụm trung gian phép đo n phép đo, việc phân tích thực phương pháp tính tốn trình bày TCVN 6910-2, phải ước lượng độ lệch chuẩn trung gian độ chụm thay cho độ lệch chuẩn lặp lại 9.3 Thí nghiệm xếp chồng Nghiên cứu sâu phương pháp ước lượng thước đo trung gian độ chụm nhằm thực thí nghiệm phức tạp Những thí nghiệm thí nghiệm xếp chồng so le thí nghiệm xếp chồng đầy đủ (xem định nghĩa thuật ngữ ISO 3534 - 3) Lợi ích việc tham gia thiết kế thí nghiệm xếp chồng thời điểm, thí nghiệm liên phòng, khơng ước lượng độ lệch chuẩn lặp lại độ lệch chuẩn tái lập mà ước lượng vài độ lệch chuẩn trung gian độ chụm Tuy nhiên, số vấn đề cần phải cân nhắc giải thích 9.8 9.4 Thí nghiệm xếp chồng đầy đủ Sơ đồ bố trí thí nghiệm xếp chồng đầy đủ mức thử cụ thể phép thử cho hình Bằng thực thí nghiệm xếp chồng đầy đủ ba yếu tố số phòng thí nghiệm nhận thước đo trung gian độ chụm đồng thời với độ lệch chuẩn lặp lại độ lệch chuẩn tái lập, nghĩa ước lượng σ (0), σ(1) σr Tương tự, sử dụng thí nghiệm xếp chồng đầy đủ yếu tố để thu hai thước đo độ chụm trung gian nghĩa ước lượng σ(0), σ(1), σ(2) σr Hình – Sơ đồ bố trí thí nghiệm xếp chồng đầy đủ ba yếu tố bốn yếu tố Các ký hiệu sử dụng số C Sự khác hiệu chuẩn E Sự khác thiết bị i Chỉ số phòng thí nghiệm cụ thể I( ) gian Chỉ số thước đo trung gian độ chụm, dấu ngoặc loại tình trung j Chỉ số mức cụ thể (TCVN 6910-2) Chỉ số nhóm phép thử nghiệm yếu tố (TCVN 6910-3) k Chỉ số kết thử nghiệm cụ thể phòng thí nghiệm i mức j L Liên phòng thí nghiệm (liên phòng) m Chỉ số độ chệch biết M Mẫu thử liên phòng O Sự khác người thao tác P Xác suất r Độ lặp lại R Độ tái lập T Sự khác thời gian W Phòng thí nghiệm thành viên 1, 2, Đối với kết thử nghiệm, đánh số theo thứ tự thu nhận chúng (1), (2), (3) Đối với kết thử nghiệm, đánh số theo thứ tự tăng độ lớn Phụ lục B (quy định) Phân tích phương sai thí nghiệm xếp chồng đầy đủ Sự phân tích phương sai mơ tả phụ lục thực riêng biệt với mức thử thí nghiệm liên phòng Để đơn giản, khơng cần thêm số mức thử sau liệu Cần ý số j sử dụng tiêu chuẩn cho yếu tố thứ (yếu tố phòng thí nghiệm) phần khác TCVN 6910, sử dụng cho mức thử Phải áp dụng phương pháp mô tả 7.3 TCVN 6910-2 : 2001 để kiểm tra số liệu quán giá trị bất thường Với thiết kế đươc mô tả phụ lục này, phân tích xác số liệu phức tạp có số kết thử nghiệm phòng thí nghiệm bị thiếu Nếu định rằng, số kết thử nghiệm phòng thí nghiệm giá trị tản mạn bất thường cần phải loại trừ khỏi việc phân tích, khuyến nghị tất số liệu thuộc phòng thí nghiệm (tại mức bị ảnh hưởng) phải loại trừ khỏi việc phân tích B1 Thí nghiệm xếp chồng đầy đủ ba yếu tố Các số liệu thu thí nghiệm biểu thị y ijk, giá trị trung bình độ rộng chúng là: ij = (yij1 + yij2) wij(1) = |yij1 – yij2| i = ( i1 + ) i2 wi(2) = | y i1 – y i2| p số phòng thí nghiệm tham gia thí nghiệm liên phòng Tổng tồn bình phương, SST, chia nhỏ sau: y p i yi Trong Khi bậc tự tổng bình phương SS0, SS1, SSe tương ứng p-1, p 2p, bảng ANOVA cho bảng B1 Bảng B1 – Bảng ANOVA cho thí nghiệm xếp chồng đầy đủ ba yếu tố Nguồn Tổng bình phương Bậc tự Bình phương trung bình Kỳ vọng bình phương trung bình SS0 p–1 MS0 = SS0 / (p – 1) σ2r + σ2(1) + σ2(0) SS1 p MS1 = SS1 / p σ2r + σ2(1) Phần dư SSe 2p MSe = SSe / 2(p) σ2r Tổng SST 4p – Các ước lượng không chệch σ2(0), σ2(1) σ2r s2(0), s2(1) s2r nhận từ bình phương trung bình MS0, MS1 MSe sau: s2(0) = (MS0 – MS1) s2(1) = (MS1 – MSe) s2r = MSe Ước lượng phương sai lặp lại, phương sai trung gian yếu tố khác độ chụm phương sai tái lập tương ứng là: s2r s2i(1) = s2r + s2(1) s2R = s2r - s2(1) + s2(0) B2 Thí nghiệm xếp chồng đầy đủ bốn yếu tố Các số liệu thu thí nghiệm biểu thị y ijkl, giá trị trung bình độ rộng chúng là: = (yijk1 + yijk2) wijk(1) = |yijk1 – yijk2| ij = ( ij1 wij(2) = |yij1 – yij2| i = ( i1 ijk + + ) ij2 ) i2 wi(3) = | i1 - | i2 = i i p số phòng thí nghiệm tham gia thí nghiệm liên phòng Tổng tồn bình phương, SST, chia nhỏ sau: Trong Với p – 1, p, 2p 4p bậc tự tổng bình phương SS0, SS1, SS2 SSe, bảng ANOVA cho bảng B2 Bảng B2 – Bảng ANOVA cho thí nghiệm xếp chồng đầy đủ bốn yếu tố Nguồn Tổng bình phương Bậc tự Bình phương trung bình Bình phương trung bình kỳ vọng SS0 p–1 MS0 = SS0 / (p - 1) σ2r + σ2(2) + σ2(1) + σ2(0) SS1 p MS1 = SS1 / p σ2r + σ2(2) + σ2(1) SS2 2p MS2 = SS2 / (2p) σ2r + σ2(2) Phần dư SSe 4p MSe = SSe / (4p) σ2r Tổng SST 8p – Các ước lượng không chệch σ2(0), σ2(1), σ2(2) σ2r s2(0), s2(1), s2(2) s2r nhận từ bình phương trung bình MS0, MS1 MSe là: s2(0) = (MS0 – MS1) s2(1) = (MS2 – MS1) s2(2) = (MS1 – MSe) s2r = MSe Ước lượng phương sai lặp lại, phương sai trung gian yếu tố khác độ chụm, phương sai trung gian hai yếu tố khác độ chụm phương sai tái lập tương ứng là: s2r s2i(1) = s2r + s2(2) s2i(2) = s2r + s2(2) +s2(1) s2R = s2r + s2(2) + s2(1) + s2(0) Phụ lục C (quy định) Phân tích phương sai thí nghiệm xếp chồng so le Phân tích phương sai mơ tả phụ lục thực riêng biệt mức thử đưa vào thí nghiệm liên phòng Để đơn giản, không cần thêm số mức thử sau số liệu Cần ý số j sử dụng tiêu chuẩn cho yếu tố lặp lại thí nghiệm Còn phần khác TCVN 6910 số j sử dụng để mức thử Phải áp dụng phương pháp mô tả 7.3 TCVN 6910-2 : 2001 để kiểm tra số liệu quán giá trị bất thường Với thiết kế đươc mô tả phụ lục này, phân tích xác số liệu phức tạp có số kết thử nghiệm phòng thí nghiệm bị thiếu Nếu định rằng, số kết thử nghiệm phòng thí nghiệm giá trị tản mạn bất thường cần phải loại trừ khỏi việc phân tích, khuyến nghị loại trừ tất số liệu thuộc phòng thí nghiệm (tại mức bị ảnh hưởng) khỏi việc phân tích C.1 Thí nghiệm xếp chồng so le ba yếu tố Số liệu thí nghiệm thu phòng thí nghiệm i biểu thị y ij (j = 1,2,3), giá trị trung bình độ rộng chúng là: i(1) = (yi1 + yi2) wi(1) = |yi1 – yi2| i(2) = (yi1 + yi2 + yi3) wi(2) = | = i i1 – yi3| 1(2) p số phòng thí nghiệm tham gia thí nghiệm liên phòng Tổng tồn bình phương, SST, chia nhỏ sau: Trong Với p – 1, p p bậc tự tổng bình phương SS0, SS1, SSe, bảng ANOVA cho bảng C.1 Bảng C.1 – Bảng ANOVA cho thí nghiệm xếp chồng so le ba yếu tố Nguồn Tổng bình phương Bậc tự Bình phương trung bình Bình phương trung bình kỳ vọng SS0 p–1 SS0 / (p – 1) SS1 p SS1 / p Phần dư SSe p SSe / p Tổng SST 3p – σ2r + σ2(1) + σ2(0) σ2r + σ2(1) σ2r Các ước lượng không chệch σ2(0), σ2(1), σ2(2) σ2r s2(0), s2(1) … s2r nhận từ bình phương trung bình MS0, MS1 MSe: s2(0) = MS0 - s2(1) = MS1 - MS1 + MSe MSe s2r = MSe Ước lượng phương sai lặp lại, phương sai trung gian yếu tố khác độ chụm phương sai tái lập tính sau: s2r s2i(1) = s2r - s2(2) s2R = s2r - s2(1) + s2(0) C.2 Thí nghiệm xếp chồng so le bốn yếu tố Số liệu thí nghiệm thu phòng thí nghiệm i biểu thị y ij (j = 1,2,3,4), giá trị trung bình độ rộng chúng là: i(1) = (yi1 + yi2) wi(1) = |yi1 – yi2| i(2) = (yi1 + yi2 + yi3) wi(2) = | i(3) = (yi1 - yi2 + yi3 + yi4) wi(3) = | y p i i(1) i(2) – yi3| – yi4| y i ( 3) p số phòng thí nghiệm tham gia thí nghiệm liên phòng Bảng ANOVA cho bảng C.2 Bảng C.2 – Bảng ANOVA cho thí nghiệm xếp chồng so le bốn yếu tố Nguồn Tổng bình phương Bậc tự Bình phương trung bình Bình phương trung bình kỳ vọng ( y i ( ) ) 4p( y) p–1 SS0 / (p - 1) σ2r + σ2(2) + σ2(1) + σ2(0) p SS1 / p σ2r + σ2(2) + σ2(1) p SS2 / p σ2r + σ2(2) p SSe / p σ2r i Phần dư w i2( 3) i w i2( ) i w i2(1) i Tổng i j ( yij y) 4p – C.3 Thí nghiệm xếp chồng so le năm yếu tố Số liệu thí nghiệm thu phòng thí nghiệm i biểu thị y ij (j = 1,2,3,4,5), giá trị trung bình độ rộng chúng là: i(1) = (yi1 + yi2) wi(1) = |yi1 – yi2| i(2) = (yi1 + yi2 + yi3) wi(2) = | i(3) = (yi1 + yi2 + yi3 + yi4) i(4) = (i1 + yi2 + yi3 + yi4 + yi5) y p i wi(3) = | wi(4) = | i(1) i(2) i(3) – yi3| – yi4| – yi5| yi (4) p số phòng thí nghiệm tham gia thí nghiệm liên phòng Bảng ANOVA cho bảng C.3 Bảng C.3 – Bảng ANOVA cho thí nghiệm xếp chồng so le năm yếu tố Nguồn Tổng bình phương Bậc tự Bình phương trung bình p–1 SS0 / (p – 1) Bình phương trung bình kỳ vọng σ2r+σ2(3) + σ2(2) + σ2(1) +5σ2(0) σ2r + σ2(3) + σ2(2) + p SS1 / p p SS2 / p σ2r + σ2(3) + p SS3 / p σ2r + σ2(3) σ2(1) σ2(2) Phần dư σ2r p Tổng i j ( yij SSe / p y) 5p – C.4 Thí nghiệm xếp chồng so le sáu yếu tố Số liệu thí nghiệm thu phòng thí nghiệm i biểu thị y ij (j = 1,2,3,4,5,6), giá trị trung bình độ rộng chúng là: i(1) = (yi1 + yi2) wi(1) = |yi1 – yi2| i(2) = (yi1 + yi2 + yi3) wi(2) = | i(3) = (yi1 + yi2 + yi3 + yi4) i(4) = (i1 + yi2 + yi3 + yi4 + yi5) wi(4) = | i(3) – yi5| i(5) = (i1 + yi2 + yi3 + yi4 + yi5 + yi6) wi(5) = | i(4) – yi6| y p i wi(3) = | i(1) i(2) – yi3| – yi4| yi (4) p số phòng thí nghiệm tham gia thí nghiệm liên phòng Bảng ANOVA cho bảng C.4 Bảng C.4 – Bảng ANOVA cho thí nghiệm xếp chồng so le sáu yếu tố Nguồn Tổng bình phương Bậc tự Bình phương trung bình Bình phương trung bình kỳ vọng p–1 SS0 / (p – 1) σ2 +σ2 +2σ2 +3σ2 + r (4) (3) (2) p SS1 / p σ2r + σ2(4) + σ2(3) + σ2(2) + p SS2 / p σ2r + σ2(2) p SS3 / p σ2r + σ2(4) + p SS4 / p σ2r + σ2(4) Phần dư p SSe / p Tổng i j ( yij σ2(4) + σ2(3) + σ2(1) + σ2(0) σ2(1) σ2(3) σ2r 6p – y) Phụ lục D (tham khảo) Những ví dụ phân tích thống kê thí nghiệm độ chụm trung gian D.1 Ví dụ – Về độ lệch chuẩn trung gian [thời gian + người thao tác] khác độ chụm, Si (TO), phòng thí nghiệm cụ thể mức cụ thể phép thử D.1.1 Cơ sở a) Phương pháp đo: xác định lượng carbon thép quang phổ xạ chân không với kết thử nghiệm biểu thị dạng % khối lượng b) Nguồn: báo cáo hàng ngày phân xưởng thép (tháng 11, 1984) c) Thiết kế thí nghiệm: phòng thí nghiệm cụ thể, lấy cách ngẫu nhiên mẫu vật liệu đem phân tích, mẫu hơm sau lại phân tích thực phân tích viên khác Trong vòng tháng, thu 29 cặp số liệu (xem bảng D.1) Bảng D.1 – Số liệu gốc – Hàm lượng carbon, % (m/ m) Chỉ số mẫu Ngày j yj1 Ngày Độ rộng wj Chỉ số mẫu Ngày j yj1 yj2 Ngày Độ rộng wj yj2 0.130 0.127 0.003 16 0.149 0.144 0.005 0.140 0.132 0.008 17 0.044 0.044 0.000 0.078 0.080 0.002 18 0.127 0.122 0.005 0.110 0.113 0.003 19 0.050 0.048 0.002 0.126 0.128 0.002 20 0.042 0.146 0.104 0.036 0.032 0.004 21 0.150 0.145 0.005 0.050 0.047 0.003 22 0.135 0.133 0.002 0.143 0.140 0.003 23 0.044 0.045 0.001 0.091 0.089 0.002 24 0.100 0.161 0.061 10 0.040 0.030 0.010 25 0.132 0.131 0.001 11 0.110 0.113 0.003 26 0.047 0.045 0.002 12 0.142 0.145 0.003 27 0.168 0.165 0.003 13 0.143 0.150 0.007 28 0.092 0.088 0.004 14 0.169 0.165 0.004 29 0.041 0.043 0.002 15 0.169 0.173 0.004 D.1.2 Phân tích Số liệu yj1, yj2 wj = | yj1 - yj2 | cho bảng D.1 Phân tích theo quy trình cho 8.2 Đồ thị số liệu [độ lệch khỏi giá trị trung bình phép đo ngày (y jk - ȳj) theo số lượng mẫu số thứ tự mẫu j] cho hình D.1 Đồ thị với việc áp dụng phép kiểm nghiệm Cochran cho ta thấy phạm vi mẫu thử số 20 24 giá trị bất thường Có khác nhiều phép đo hàng ngày mẫu chủ yếu sai số việc ghi số liệu Các giá trị hai mẫu 20, 24 loại bỏ tính độ lệch chuẩn trung gian [thời gian + người thao tác khác độ chụm], S i(TO) Độ lệch chuẩn tính cơng thức 12) sau: Sl ( TO ) x 27 27 j w 2j 2,87 x 10 Hình D.1 – Hàm lượng carbon thép – Độ lệch khỏi giá trị trung bình phép đo hai ngày theo số lượng mẫu thử D.2 Ví dụ – Về độ lệch chuẩn trung gian với thời gian khác độ chụm thí nghiệm liên phòng D.2.1 Cơ sở: a) Phương pháp đo: xác định lượng vanadi thép phương pháp quang phổ hấp thụ miêu tả tài liệu hướng dẫn thí nghiệm Kết thí nghiệm biểu thị dạng % khối lượng b) Nguồn: ISO/TC 17, Steel//SC 1, Phương pháp xác định thành phần hóa học Thí nghiệm thực tháng 5/1985 c) Thiết kế thí nghiệm: thí nghiệm xếp chồng so le yếu tố thực 20 phòng thí nghiệm, phòng thí nghiệm báo cáo kết thử nghiệm nhận điều kiện lặp lại ngày thứ ngày thứ hai mức mức thử đưa vào thí nghiệm Tất phép đo phòng thí nghiệm thực người thao tác với thiết bị đo D.2.2 Phân tích Số liệu tất mức thử cho bảng D.2 Bảng D.2 – Hàm lượng vanadi thép Đồ thị kết thử nghiệm ngày ngày mức thử theo số thứ tự phòng thí nghiệm Bảng D.2 – Số liệu gốc – Hàm lượng vanadi, % (m/m) Chỉ số phòng thí nghiệm Mức (0,01 %) Mức (0,04 %) Mức (0,1 %) Mức (0,2 %) Mức (0,5 %) Mức (0,75 %) Ngày Ngày Ngày Ngày Ngày Ngày yi1 yi2 Ngày yi3 yi1 yi2 Ngày yi3 yi1 yi2 Ngày yi3 yi1 yi2 Ngày yi3 yi1 yi2 Ngày yi3 yi1 yi2 Ngày yi3 0,0091 0,0102 0,0098 0,0382 0,0388 0,038 0,101 0,103 0,102 0,214 0,211 0,210 0,514 0,510 0,513 0,755 0,753 0,751 0,0100 0,0100 0,0090 0,0410 0,0410 0,0390 0,111 0,111 0,108 0,220 0,220 0,215 0,512 0,540 0,540 0,800 0,755 0,750 0,0095 0,0090 0,0094 0,0390 0,0380 0,0370 0,108 0,110 0,107 0,213 0,215 0,215 0,500 0,514 0,504 0,738 0,730 0,724 0,0080 0,0083 0,0077 0,0374 0,0361 0,0382 0,109 0,106 0,104 0,214 0,222 0,201 0,519 0,518 0,518 0,744 0,742 0,732 0,0100 0,0100 0,0100 0,0350 0,0370 0,0370 0,103 0,103 0,110 0,210 0,210 0,205 0,495 0,500 0,512 0,743 0,753 0,750 0,0089 0,0094 0,0094 0,0368 0,0368 0,0377 0,106 0,106 0,108 0,232 0,240 0,221 0,526 0,532 0,513 0,733 0,740 0,746 0,0098 0,0099 0,0101 0,0376 0,0380 0,0384 0,107 0,105 0,108 0,215 0,215 0,216 0,521 0,519 0,526 0,754 0,756 0,756 0,0096 0,0094 0,0099 0,0379 0,0366 0,0379 0,108 0,107 0,108 0,193 0,195 0,210 0,507 0,493 0,511 0,732 0,729 0,732 0,0104 0,0094 0,0100 0,0365 0,0370 0,0367 0,104 0,106 0,105 0,211 0,205 0,213 0,509 0,515 0,515 0,734 0,738 0,747 10 0,0107 0,0118 0,0097 0,0370 0,0375 0,0380 0,105 0,110 0,105 0,210 0,220 0,225 0,520 0,520 0,525 0,760 0,760 0,765 11 0,0100 0,0100 0,0090 0,0380 0,0380 0,0375 0,102 0,102 0,102 0,213 0,211 0,214 0,513 0,516 0,514 0,746 0,748 0,746 12 0,0109 0,0115 0,0097 0,0390 0,0390 0,0390 0,101 0,108 0,105 0,208 0,215 0,210 0,509 0,528 0,510 0,758 0,748 0,750 13 0,0100 0,0095 0,0095 0,0375 0,0375 0,0375 0,103 0,104 0,108 0,212 0,222 0,215 0,510 0,520 0,505 0,735 0,755 0,750 14 0,0096 0,0096 0,0100 0,0374 0,0374 0,0389 0,104 0,106 0,110 0,218 0,218 0,212 0,520 0,528 0,522 0,740 0,735 0,742 15 0,0099 0,0091 0,0082 0,0381 0,0375 0,0392 0,109 0,106 0,107 0,214 0,210 0,211 0,510 0,510 0,515 0,749 0,729 0,744 16 0,0098 0,0100 0,0095 0,0373 0,0377 0,0397 0,105 0,105 0,104 0,215 0,212 0,218 0,519 0,517 0,531 0,754 0,751 0,759 17 0,0105 0,0102 0,0112 0,0389 0,0382 0,0373 0,107 0,108 0,104 0,214 0,210 0,209 0,517 0,515 0,514 0,735 0,728 0,741 18 0,0103 0,0105 0,0118 0,0382 0,0380 0,0374 0,103 0,104 0,103 0,224 0,218 0,217 0,515 0,514 0,517 0,788 0,798 0,787 19 0,0098 0,0096 0,0104 0,0383 0,0375 0,0366 0,110 0,109 0,104 0,217 0,215 0,215 0,530 0,525 0,520 0,755 0,745 0,740 20 0,0140 0,0140 0,0100 0,0370 0,0408 0,0369 0,104 0,106 0,107 0,214 0,214 0,203 0,518 0,518 0,481 0,730 0,737 0,658 Việc phân tích phương sai trình bày cho mức thử mức Đồ thị số liệu (sự phụ thuộc kết thử ngày ngày vào số thứ tự phòng thí nghiệm) trình bày hình D.2 Đồ thị phòng thí nghiệm số 20 phòng thí nghiệm bất thường Có khác nhiều kết thử nghiệm ngày giá trị trung bình ngày 1, sai lệch lớn so với kết thử nghiệm phòng thí nghiệm khác Phòng thí nghiệm loại bỏ tính tốn thước đo độ chụm Các đại lượng wi(1), wi(2) yi(2) tính tốn theo C.1 phụ lục C kết cho bảng D.3 Bảng D.3 – Giá trị wi(1), wi(2) Chỉ số phòng thí nghiệm i(2) wi(1) wi(2) 0,001 0,000 15 0,009 700 0,000 0,001 00 0,009 667 0,000 0,000 15 0,009 300 0,000 0,000 45 0,008 000 0,000 0,000 00 0,010 000 0,000 0,000 25 0,009 233 0,000 0,000 25 0,009 933 0,000 0,000 40 0,009 633 0,001 0,000 10 0,009 933 10 0,001 0,001 55 0,010 733 11 0,000 0,001 00 0,009 667 12 0,000 0,001 50 0,010 700 13 0,000 0,000 25 0,009 667 14 0,000 0,000 40 0,009 733 15 0,000 0,001 30 0,009 067 16 0,000 0,000 40 0,009 767 17 0,000 0,000 85 0,010 633 18 0,000 0,001 40 0,010 867 19 0,000 0,000 70 0,009 933 i(2) i Tổng bình phương wi(1), wi(2) i i i y i(2) giá trị trung bình i(2), tính sau: w i2(1) = 5.52 x 10-6 w i2( ) = 12.44 x 10-6 ( y1( ) ) = 832.16 x 10-6 16 i y i ( ) = 0.009 798 25 Từ giá trị thu tổng bình phương SS0, SS1 SSe bảng ANOVA cho bảng D.4 Bảng D.4 – Bảng ANOVA – Hàm lượng vanadi Nguồn Tổng bình phương Bậc tự Bình phương trung bình (phòng thí nghiệm) 24,16 x 10-6 18 1,342 x 10-6 8,29 x 10-6 (ngày) 0,436 x 10-6 19 Bình phương trung bình kỳ vọng σ2r + σ2(1) + σ2(0) σ2r + σ2(1) 2,76 x 10-6 0,145 x 10-6 Phần lại 19 35,21 x 10-6 Tổng σ2r 56 Ước lượng khơng chệch phương sai phòng thí nghiệm, s 2(0), ngày phòng thí nghiệm s2(1) phương sai lặp lại ước lượng s2r thu là: s2(0) = 0,278 x 10-6 s2(1) = 0,218 x 10-6 s2r = 0,145 x 10-6 Độ lệch chuẩn tái lập sR, độ lệch chuẩn trung gian thời gian khác độ chụm s I(T) độ lệch chuẩn lặp lại s thu sau: S2r S(21) S(20 ) SR S2r Sl ( T ) S2r SR S(21) 0.801 x 10 0.603 x 10 0.381 x 10 3 Giá trị độ lệch chuẩn mức thử hàm lượng vanadium tóm tắt bảng trình bày hình D.3 Bảng D.5 – Giá trị sr, sI(T) sR sáu mức hàm lượng vanadi thép Mức Chỉ số phòng thí nghiệm bất thường Trung bình (%) sr (%) sI(T) (%) sR (%) 20 0,009 0,381 x 10-3 0,603 x 10-3 0,801 x 10-3 2 0,037 0,820 x 10-3 0,902 x 10-3 0,954 x 10-3 - 0,105 1,739 x 10-3 2,305 x 10-3 2,650 x 10-3 0,213 3,524 x 10-3 4,710 x 10-3 4,826 x 10-3 20 0,516 6,237 x 10-3 6,436 x 10-3 9,412 x 10-3 20 0,748 9,545 x 10-3 9,545 x 10-3 15,862 x 10-3 Hình D.3 – Hàm lượng vanadi thép – Độ lệch chuẩn lặp lại s r, độ lệch chuẩn trung gian thời gian khác độ chụm sI(T) độ lệch chuẩn tái lập sR hàm mức nồng độ Phụ lục E (tham khảo) Tài liệu tham khảo [1] ISO 3534-2 : 1993, Statistic – Vocabulary and symbols – Part 2: Statistical quality control [2] ISO 3534-3 : 1993, Statistic – Vocabulary and symbols – Part 3: Design of experiments [3] TCVN 6910-4, Độ xác (độ độ chụm) phương pháp đo kết đo – Phần 4: Phương pháp xác định độ phương pháp đo tiêu chuẩn [4] TCVN 6910-5, Độ xác (độ độ chụm) phương pháp đo kết đo – Phần 5: Các phương pháp khác để xác định độ chụm phương pháp đo tiêu chuẩn [5] TCVN 6910-6, Độ xác (độ độ chụm) phương pháp đo kết đo – Phần 6: Sử dụng giá trị độ xác thực tế [6] Winer, B.J Statistical principles in experimental design, - McGraw - Hill 1962 [7] Snedecor, G.W Cochran, W.G, Statistical methods, lowa University Press, 1967 ... 1 0-3 0,603 x 1 0-3 0,801 x 1 0-3 2 0,037 0,820 x 1 0-3 0,902 x 1 0-3 0,954 x 1 0-3 - 0,105 1,739 x 1 0-3 2,305 x 1 0-3 2,650 x 1 0-3 0,213 3,524 x 1 0-3 4,710 x 1 0-3 4,826 x 1 0-3 20 0,516 6,237 x 1 0-3 ... đo tiêu chuẩn ISO Guide 33 : 1989, Sử dụng mẫu chuẩn chứng nhận ISO Guide 35 : 1989, Chứng nhận mẫu chuẩn Những nguyên tắc chung thống kê Định nghĩa Trong tiêu chuẩn áp dụng định nghĩa cho ISO. .. dụ: - Độ lệch chuẩn độ chụm trung gian khác thời gian, S I(T); - Độ lệch chuẩn độ chụm trung gian khác hiệu chuẩn, S I(C); - Độ lệch chuẩn độ chụm trung gian khác người thao tác, S I(O); - Độ