Thuộc tính
Nội dung
Tiếng Anh (English)
Văn bản gốc/PDF
Lược đồ
Liên quan hiệu lực
Liên quan nội dung
Thuộc tính
Tải về

Đang tải văn bản...

Tiêu chuẩn Việt Nam TCVN 6910-5:2002 Độ chính xác phương pháp đo và kết quả đo

Số hiệu:	TCVN6910-5:2002	Loại văn bản:	Tiêu chuẩn Việt Nam
Nơi ban hành:	Bộ Khoa học Công nghệ và Môi trường	Người ký:	***
Ngày ban hành:	30/10/2002	Ngày hiệu lực:
ICS:	17.020	Tình trạng:	Đã biết

Phòng thí nghiệm	Mức
	1		2		j		q
	a	b	a	b	a	b	a	b
1
2

i

p

4.4. Mô hình thống kê

4.4.1. Mô hình cơ bản dùng trong tiêu chuẩn này được cho trong điều 5 của TCVN 6910-1:2001 (ISO 5725-1:1994). Để ước lượng độ chính xác (độ đúng và độ chụm) của phương pháp đo, thường giả thiết rằng mỗi kết quả đo là tổng của ba thành phần:

(3)

trong đó, đối với một vật liệu cụ thể được thử nghiệm:

m_j biểu thị trung bình chung ở mức cụ thể j = 1,…., q;

B_ij là biểu thị thành phần phòng thí nghiệm của độ chệch trong các điều kiện lặp lại ở một phòng thí nghiệm cụ thể i = 1,…..p, ở một mức cụ thể j = 1….., q;

e_ijk biểu thị sai số ngẫu nhiên của kết quả thử nghiệm k = 1.2….. n. nhận được trong phòng thí nghiệm i ở mức j dưới điều kiện lặp lại:

4.4.2. Đối với thí nghiệm mức tách biệt mô hình đó trở thành:

(4)

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Sự vắng mặt của chỉ số k trong B_ij có nghĩa là người ta giả thiết rằng độ chệch gắn với phòng thí nghiệm i không phụ thuộc vào vật liệu a hoặc b trong một mức. Điều này giải thích tại sao hai vật liệu phải là giống nhau.

4.4.3. Xác định trung bình ô như sau:

và độ sai khác theo ô (hiệu số ô) như sau:

4.4.4. Trung bình chung cho một mức j của một thí nghiệm mức tách biệt có thể được xác định như sau:

4.5. Phân tích thống kê các số liệu từ thí nghiệm mức tách biệt

4.5.1. Tập hợp các số liệu thành một bảng như đã chỉ ra trong bảng 1. Mỗi sự kết hợp một phòng thí nghiệm với một mức tạo ra một “ô” trong bảng đó chứa hai thánh phần của số liệu, y_ija và y_ijb.

...

a - b

và dấu của hiệu số đó cần được ghi lại.

Tính các trung bình ô y và đưa chúng vào một bảng như đã chỉ ra trong bảng 3.

4.5.2. Nếu một ô trong bảng 1 không chứa hai kết quả thử nghiệm (ví dụ: vì các mẫu đã bị hỏng hoặc các số liệu đã bị loại ra sau khi áp dụng các phép thử về giá trị bất thường sẽ được miêu tả dưới đây) thì các ô tương ứng trong bảng (2) và (3) đều là ô trống.

4.5.3. Với mỗi mức j của thí nghiệm, tính giá trị trung bình D_j và độ lệch chuẩn s_Dj của các độ sai khác trong cột j của bảng 2:

Trong đó tổng S lấy theo tất cả các phòng thí nghiệm i = 1,2,….,p.

Nếu có các ô trống trong bảng 2, thì p là số các ô có chứa số liệu trong cột j của bảng 2 và tổng được thực hiện trên tất cả các ô không trống.

4.5.4. Với mỗi mức j của thí nghiệm tính giá trị trung bình y_j và độ lệch chuẩn s_yj của giá trị trung bình ô trong cột thứ j của bảng 3:

...

Trong đó tổng S lấy theo tất cả các phòng thí nghiệm i = 1,2,…..,p.

Nếu có các ô trống trong bảng 3, thì p là số các ô có chứa số liệu trong cột j của bảng 3 và tổng được thực hiện trên tất cả các ô không trống.

4.5.5. Sử dụng bảng 2 và 3 và các thống kê được tính trong 4.5.3 và 4.5.4 để kiểm tra tính phù hợp và tính bất thường của số liệu như sẽ được mô tả trong 4.6. Nếu các số liệu bị loại bỏ, tính lại các thống kê đó.

4.5.6. Tính độ lệch chuẩn lặp lai s_rj và độ lệch chuẩn tái lặp s_R từ các công thức:

4.5.7. Nghiên cứu xem s_rj và s_RJ có phụ thuộc vào giá trị trung bình y_j hay không, nếu có phải xác định mối quan hệ hàm bằng cách sử dụng các phương pháp đã được mô tả trong 7.5 của TCVN 6910-2: 2001

Bảng 2 - Biểu mẫu nên dùng để lập bảng các độ sai khác theo ô đối với thiết kế mức tách biệt

Phòng thí nghiệm

Mức

...

Bảng 3 - Biểu mẫu nên dùng để lập bảng các giá trị trung bình ô đối với thiết kế mức tách biệt

Phòng thí nghiệm

Mức

...

4.6. Khảo sát sự phù hợp và bất thường của các số liệu

4.6.1. Để kiểm tra sự phù hợp của các số liệu người ta dùng các thống kê h đã được mô tả trong 7.3.1 của TCVN 6910-2:2001 (ISO 5725-2:1994).

Để kiểm tra sự phù hợp của các độ sai khác theo ô, tính các thống kê h như sau:

Để kiểm tra sự phù hợp của các trung bình ô, tính các thống kê h như sau

Để chỉ ra các phòng thí nghiệm không phù hợp, vẽ biểu đồ cả hai tập thống kê đó theo thứ tự của các mức, nhưng gộp lại theo phòng thí nghiệm, như đã được chỉ ra trong hình 2 và 3. Việc giải thích các đồ thị đó sẽ được đề cập đến một cách đầy đủ trong 7.3.1 của TCVN 6910-2:2001 (ISO 5725-2:1994). Nếu một phòng thí nghiệm có độ lặp lại xấu hơn các phòng thí nghiệm khác thì nó được thể hiện bởi một số lượng lớn khác thường các thống kê h có giá trị lớn trên đồ thị được vẽ từ các độ sai khác theo ô. Nếu phòng thí nghiệm đưa ra các kết quả chệch thì điều đó được thể hiện bởi sự biến thiên theo một hướng của các thống kê h trên đồ thị được dựng lên từ các trung bình ô. Trong các trường hợp như vậy phải yêu cầu các phòng thí nghiệm đó nghiên cứu và báo cáo những sai sót của những người tổ chức thí nghiệm.

4.6.2. Kiểm tra các số liệu tản mạn và bất thường bằng phép kiểm nghiệm Grubb, đã được mô tả trong 7.3.4. của TCVN 6910-2 2001 (ISO 5725-2:1994)

Để kiểm nghiệm giá trị tản mạn và bất thường của các độ sai khác theo ô, áp dụng phép kiểm nghiệm Grubb đối với các giá trị trong mỗi cột của bảng 2.

...

Các phép kiểm nghiệm đó đã được giải thích một cách đầy đủ trong 7.3.2 của TCVN 6910-2:2001 (ISO 5725-2:1994). Chúng được dùng để phát hiện các kết quả không phù hợp với các kết quả còn lại của các số liệu ghi được trong thí nghiệm, mà dựa trên các số liệu đó các kết luận rút ra từ việc tính độ lệch chuẩn lặp lại và tái lập sẽ ảnh hưởng đến các giá trị của các thống kê đó một cách thực chất. Thông thường các số liệu tỏ ra là bất thường sẽ bị loại ra khi tính toán và các số liệu tản mạn vẫn được đưa vào, trừ khi có đủ lý do để loại. Nếu các phép kiểm nghiệm chỉ ra rằng một giá trị ở trong bảng 2 hoặc bảng 3 bị loại khi tính toán độ lệch chuẩn lặp lại và tái lập thì giá trị tương ứng trong bảng kia cũng bị loại.

4.7. Báo cáo kết quả của thí nghiệm mức tách biệt

4.7.1. Khuyến nghị được đưa ra trong 7.7 của TCVN 6910-2:2001 (ISO 5725-2:1994) về:

- báo cáo kết quả của việc phân tích thống kê cho hội đồng chuyên môn;

- các quyết định của hội đồng; và

- việc chuẩn bị báo cáo đầy đủ.

4.7.2. Các khuyến nghị về hình thức công bố về độ lệch chuẩn lặp lại và tái lập của một phương pháp đo tiêu chuẩn được trình bày trong 7.1 của TCVN 6910-1:2001 (ISO 5725-1:1994).

4.8. Ví dụ 1: Thí nghiệm mức tách biệt - Xác định protein

4.8.1. Bảng 4 bao gồm các số liệu từ thí nghiệm ^[5] để xác định hàm lượng protein trong thức ăn bằng cách đốt cháy. Có chín phòng thí nghiệm tham gia và thí nghiệm bao gồm 14 mức. Trong mỗi mức sử dụng hai phần thức ăn có tỷ lệ khối lượng protein như nhau.

...

Sử dụng phương trình (8) và (9) trong 4.5.3, đối với sự sai khác trong bảng 5 thu được:

D₁₄ = 8,34 %

s_D14 = 0,4361 %

và áp dụng phương trình (10) và (11) trong 4.5.4 đối với các giá trị trung bình ô trong bảng 6 thu được:

y₁₄ = 85,46 %

s_y14 = 0,4534 %

như vậy độ lệch chuẩn lặp lại và tái lập khi sử dụng các phương trình (12) và (13) trong 4.5.6 sẽ là:

s_r14 = 0,31 %

s_R14 = 0,50 %

...

4.8.3. Hình 1 trình bày “biểu đồ Youden”' cho mẫu a tương ứng với mẫu b ở mức 14 trong bảng 4. Phòng thí nghiệm 5 cho một điểm ở góc trái phía dưới của đồ thị, phòng thí nghiệm 1 cho một điểm ở góc phải phía trên. Điều đó chứng tỏ rằng các số liệu của phòng thí nghiệm 5 có độ chệch âm và các số liệu của phòng thí nghiệm 1 có độ chệch dương trên cả hai mẫu a và b. Điều đó thường thấy khi vẽ biểu đồ các số liệu từ một thiết kế mức tách biệt như hình 1. Hình đó cũng chỉ ra rằng các kết quả của phòng thí nghiệm 4 là không bình thường, vì một điểm đối với phòng thí nghiệm này nằm cách xa đường trung bình của hai mẫu đó. Các phòng thí nghiệm khác tạo nên một nhóm ở giữa biểu đồ. Như vậy hình này nêu ra một trường hợp để nghiên cứu những nguyên nhân của các độ chệch ở ba phòng thí nghiệm đó.

Chú thích - Để có thêm thông tin khác về “biểu đồ Youden” xem các tài liệu tham khảo [7] và (8).

4.8.4. Giá trị của các thống kê h ở mức 14, được tính như trong 4.6.1, chỉ ra trong bảng 5 và 6. Các giá trị cho tất cả các mức được vẽ thành biểu đồ trong hình 2 và 3.

Trong hình 3, thống kê h đối với các trung bình theo ô chứng tỏ rằng phòng thí nghiệm 5 cho các thống kê h âm ở tất cả các mức và điều này thể hiện rõ độ chệch âm nhất quán với các số liệu của chúng. Cũng trong hình đó các phòng thí nghiệm 8 và 9 cho các thống kê h hầu như tất cả đều là dương chứng tỏ rằng có độ chệch dương nhất quán với các số liệu của chúng (nhưng nhỏ hơn độ chệch âm trong phòng thí nghiệm 5). Cũng như vậy các thống kê h đối với các phòng thí nghiệm 1, 2 và 6 thể hiện một độ chệch biến đổi theo mức trong mỗi phòng thí nghiệm đó. Sự tương tác như vậy giữa các phòng thí nghiệm và các mức có thể cung cấp các đầu mối về nguyên nhân của các độ chệch phòng thí nghiệm.

Hình 2 không thể hiện một kiểu dáng nào đáng chú ý

4.8.5. Các giá trị của thống kê Grubb được cho trong bảng 8. Các phép kiểm nghiệm đó chứng tỏ rằng các số liệu từ phòng thí nghiệm 5 là đáng nghi ngờ.

4.8.6. Ở giai đoạn này của việc phân tích, chuyên gia thống kê cần bắt đầu khảo sát các nguyên nhân gây ra số liệu nghi ngờ ở phòng thí nghiệm 5 trước khi tiến hành phân tích số liệu. Nếu không phát hiện được các nguyên nhân đó có khi phải loại tất cả các số liệu từ phòng thí nghiệm 5 khi tính độ lệch chuẩn lặp lại và tái lập. Sau đó tiếp tục sự phân tích bằng việc khảo sát các mối quan hệ hàm số có thể có giữa độ lệch chuẩn lặp lại và tái lập và trung bình chung. Trong TCVN 6910-2 (ISO 5725-2): chưa đề cập đến vấn đề nào liên quan tới điều này và như vậy nó chưa được xét ở đây.

Bảng 4 - Ví dụ 1: Xác định tỷ lệ khối lượng protein trong thức ăn, theo %

Phòng thí nghiệm

...

11,11

10,34

10,91

...

13,74

13,48

13,79

13,00

15,89

15,26

11,12

9,94

...

10,31

14,00

13,12

13,44

13,06

15,69

15,10

11,26

...

10,95

10,51

13,38

12,70

13,54

13,18

15,83

15,73

...

10,41

11,66

9,95

13,01

13,16

13,58

12,88

15,08

15,63

...

10,69

10,31

10,98

10,13

13,24

13,33

13,32

12,59

15,02

...

11,73

11,01

12,31

10,92

14,01

13,66

14,04

13,64

...

15,94

11,13

10,36

11,38

10,44

12,94

12,44

13,63

...

15,75

15,56

11,21

10,51

11,32

10,84

13,09

13,76

...

13,49

15,98

15,89

11,80

11,21

11,35

9,88

13,85

...

13,96

13,77

16,51

15,72

Phòng thí nghiệm

Mức

...

20,14

19,78

20,33

20,06

46,45

44,42

52,05

...

65,84

59,14

19,25

20,25

20,36

19,94

46,69

44,62

...

48,81

66,31

59,19

20,48

19,86

20,56

20,11

46,90

...

52,18

48,90

66,06

58,52

21,54

20,06

20,64

20,46

...

45,29

51,73

48,56

65,93

58,93

19,90

19,66

20,56

...

45,83

43,73

50,84

47,91

64,19

57,94

20,31

20,27

...

20,63

40,86

43,96

52,18

49,03

65,73

58,77

20,00

...

20,25

20,19

46,25

44,31

52,25

49,44

66,06

59,19

...

20,69

20,85

20,27

47,11

44,40

52,44

48,81

65,66

59,38

...

20,64

21,01

20,78

20,89

47,09

45,15

52,19

48,46

66,33

...

Phòng thí nghiệm

Mức

...

84,16

...

85,38

81,71

87,64

88,23

90,24

82,10

...

81,06

85,56

82,44

88,81

88,38

89,88

81,44

...

82,26

79,43

85,26

82,15

88,58

88,12

89,48

81,67

...

84,39

80,08

85,20

81,76

88,47

87,98

90,04

80,73

...

81,71

79,01

83,58

79,74

86,43

86,19

88,59

...

82,85

81,16

84,44

80,90

87,78

86,89

...

80,88

86,25

81,00

84,88

81,44

88,06

...

89,31

81,38

84,59

81,16

84,96

81,71

...

87,98

89,94

81,56

83,05

80,93

84,73

81,94

...

88,05

89,75

81,35

Bảng 5 - Ví dụ 1: Độ sai khác theo ô ở mức 14

Bảng 6 - Ví dụ 1: Trung bình ô ở mức 14

...

Độ sai khác trong ô

Thống kê

Phòng thí nghiệm

Độ sai khác trong ô

Thống kê

...

8,14

-0,459

86,170

1,576

8,44

...

85,660

0,451

7,81

- 1,215

...

0,263

9,31

2,224

85,385

- 0,156

...

- 0,482

84,525

- 2,052

8,52

0,413

...

85,140

- 0,696

7,93

- 0,940

85,345

- 0,244

...

8,38

0,092

85,750

0,649

8,40

0,138

...

85,550

0,208

Bảng 7 - Ví dụ 1: Giá trị của các trung bình, độ sai khác trung bình và độ lệch chuẩn tính từ các số liệu ở tất cả 14 mức trong bảng 4

Mức

Số lượng phòng thí nghiệm

Trung bình chung

...

Độ sai khác trung bình

D_j%

Độ lệch chuẩn

S_yj %

s_Dj %

s_rj %

s_Rj %

...

10,87

0,73

0,35

0,21

0,15

0,36

10,84

...

0,36

0,43

0,30

0,42

13,11

0,13

0,44

...

0,39

0,52

13,43

0,50

0,30

0,21

0,15

...

15,66

0,27

0,39

0,40

0,29

0,44

...

20,27

0,06

0,40

0,73

0,52

0,54

20,39

...

0,30

0,41

0,29

0,37

45,60

2,21

0,44

...

0,26

0,47

50,40

3,16

0,44

0,35

0,25

...

62,37

6,84

0,53

0,40

0,28

0,57

...

82,14

3,23

1,01

1,08

0,77

1,15

83,17

...

0,74

0,46

0,33

0,77

87,91

0,30

0,69

...

0,29

0,72

85,46

8,34

0,45

0,44

0,31

...

Bảng 8 - Ví dụ 1: Giá trị của các thống kê Grubb

Thống kê Grubb cho các độ sai khác

Mức

Một giá trị nhỏ nhất

Hai giá trị nhỏ nhất

Hai giá trị lớn nhất

Một giá trị lớn nhất

1,653

...

0,3139

2,125

1,418

0,3945

0,4738

1,535

1,462

...

0,5323

1,379

1,490

0,5841

0,4771

1,414

2,033

...

0,6075

1,289

1,456

0,5490

0,3210

1,947

1,185

...

0,1712

2,296* (5)

0,996

0,7571

0,1418* (6; 8)

1,876

1,458

...

0,3092

1,602

1,474

0,3360

0,4578

1,737

1,422

...

0,2943

1,865

1,418

0,6009

0,2899

1,956

2,172

...

0,6326

1,444

1,215

0,6220

0,2362

2,224* (4)

Thống kê Grubb cho các trung bình ô

Mức

...

Hai giá trị nhỏ nhất

Hai giá trị lớn nhất

Một giá trị lớn nhất

1,070

0,6607

0,1291* (6 : 9)

1,832

...

0,6288

0,2118

2,165

1,621

0,4771

0,4077

1,680

...

0,5339

0,3807

1,429

1,794

0,4018

0,5009

1,333

...

0,4947

0,4095

1,386

1,599

0,5036

0,4391

1,470

...

0,3753

0,4536

1,404

2,328* (5)

0,1317* (4 : 5)

0,7417

1,025

...

1,000

1,756

0,2469

0,5759

1,472

...

0,1063* (5 : 6)

0,7116

1,130

2,308* (5)

0,0733* (5 : 6)

0,7777

0,994

...

0,2781

0,5486

1,576

Chú thích - Số trong ngoặc đơn biểu thị các phòng thí nghiệm đưa ra các giá trị tản mạn hoặc bất thường

Giá trị tới hạn của các thống kê kiểm nghiệm Grubb cho 9 phòng thí nghiệm, áp dụng cho các độ sai khác hoặc các trung bình ô, như sau:

Giá trị tản mạn (*)

Giá trị bất thường (**)

Phép kiểm nghiệm Grubb cho một giá trị bất thường đơn lẻ

...

2.387

Phép kiểm nghiệm Grubb cho một cặp các giá trị bất thường

0,1492

0,0851

Hình 1 - Ví dụ 1: Các số liệu thu được ở mức 14

Hình 2 - Ví dụ 1: Kiểm tra sự phù hợp của các độ sai khác theo ô
(gộp lại theo phòng thí nghiệm)

...

Hình 3 - Ví dụ 1: Kiểm tra sự phù hợp trên các trung bình ô
(gộp lại theo phòng thí nghiệm)

5. Thiết kế cho các vật liệu không đồng nhất

5.1. Ứng dụng các thiết kế cho các vật liệu không đồng nhất

5.1.1. Một ví dụ về vật liệu không đồng nhất là da thuộc. Không có hai miếng da nào là giống hệt nhau và các tính chất của da thuộc thực chất thay đổi trong một miếng da. Một kiểm nghiệm thông thường đối với da thuộc là kiểm nghiệm về độ bền kéo đứt theo BS 3144^[3]. Điều đó được thực hiện trên các mẫu thử dạng cặp đôi. (BS 3144 quy định số các mẫu thử như vậy được cắt từ một miếng da cũng như vị trí và hướng của chúng trong miếng da, như vậy việc định nghĩa tự nhiên về “mẫu” để dùng khi kiểm nghiệm da thuộc là một miếng da hoàn chỉnh). Nếu một thí nghiệm về độ chụm được thực hiện khi sử dụng thiết kế đồng mức như được mô tả trong TCVN 6910-2 (ISO 5725-2), mỗi phòng thí nghiệm được gửi tới một mảnh da ở mỗi mức của thí nghiệm và thu được hai kết qua thử nghiệm trên mỗi miếng da. Sự thay đổi giữa các miếng da cộng thêm với sự thay đổi giữa các phòng thí nghiệm như vậy sẽ làm gia tăng độ lệch chuẩn tái lập. Tuy nhiên nếu mỗi phòng thí nghiệm được gửi tới hai miếng da ở mỗi mức, và hai kết quả thử nghiệm thu được trên mỗi mảnh da, thì các số liệu có thể được dùng để ước lượng sự thay đổi giữa các miếng da và để tính các giá trị cho độ lệch chuẩn tái lập của phương pháp thử nghiệm trong đó sự thay đổi giữa các miếng da đã được loại bỏ.

5.1.2. Một ví dụ khác về vật liệu không đồng nhất là cát (có thể được dùng làm bê tông). Dưới tác động của gió và nước, cát nằm theo các tầng chứa các hạt biến đổi dần theo cỡ của các hạt. Như vậy khi cát được sử dụng, phân bố cỡ hạt luôn luôn được quan tâm. Trong kỹ thuật, phân bố cỡ hạt của cát được đo bằng thử nghiệm sàng lọc (ví dụ BS 812-103^[1]). Để kiểm nghiệm cát, một mẫu chung được lấy ra từ một sản phẩm, sau đó một hoặc một số mẫu thử được tạo ra từ mẫu chung. Thông thường mẫu chung có khối lượng vào khoảng 10 kg và các mẫu thử vào khoảng 200 g. Vì sự thay đổi tự nhiên của vật liệu, luôn luôn có sự khác nhau giữa các mẫu chung của cùng một sản phẩm. Do đó cũng như với da thuộc nếu thực hiện một thí nghiệm đồng mức, trong đó mỗi phòng thí nghiệm được gửi một mẫu chung ở mỗi mức, sự thay đổi giữa các mẫu chung sẽ làm gia tăng độ lệch chuẩn tái lập tính toán của phương pháp thử nghiệm, nhưng nếu các phòng thí nghiệm được gửi tới hai mẫu chung ở mỗi mức thì giá trị của các độ lệch chuẩn tái lập có thể tính được và loại bỏ được sự thay đổi đó.

5.1.3. Các ví dụ trên cùng thể hiện rõ các đặc tính khác của các vật liệu không đồng nhất: vì sự thay đổi của vật liệu, việc chuẩn bị các mẫu thử hoặc lô thử có thể là một nguồn quan trọng của sự thay đổi. Như đối với da thuộc, quá trình cắt các mẫu thử từ một miếng da có thể ảnh hưởng lớn đến độ dài có thể kéo được và với các thử nghiệm sàng lọc trên cát quá trình chuẩn bị các mẫu thử từ mẫu chung thường là những nguồn thay đổi chủ yếu trong phương pháp thử nghiệm. Nếu mẫu hoặc các mẫu thử được chuẩn bị cho thí nghiệm về độ chụm không tương ứng đối với thực tế thông thường (cố gắng tạo ra các mẫu đồng nhất), thì các giá trị của độ lệch chuẩn lặp lại và tái lập được tạo ra bởi thí nghiệm đó sẽ không thể hiện được sự thay đổi diễn ra trong thực tế. Có những tình huống trong đó người ta mong muốn tạo ra các mẫu đồng nhất bằng một quá trình được thiết kế đặc biệt để loại bỏ càng nhiều càng tốt sự thay đổi của vật liệu ( ví dụ đối với sự thử nghiệm về độ thành thạo hoặc khi một thí nghiệm về độ chụm được sử dụng như là một phần của chương trình làm việc trong khi triển khai một phương pháp đo). Tuy nhiên, khi mục tiêu của thí nghiệm về độ chụm là phát hiện sự thay đổi sẽ diễn ra trong thực tế (ví dụ khi người bán và người mua kiểm nghiệm các mẫu của cùng một sản phẩm), sự thay đổi do tính không đồng nhất của vật liệu cần được bao gồm trong các số đo về độ chụm của phương pháp đo.

Cần phải thận trọng để đảm bảo rằng mỗi kết quả thử nghiệm trong một thí nghiệm thu được bằng cách thực hiện một quy trình thử nghiệm là độc lập với các kết quả thử nghiệm khác. Điều này sẽ không còn đúng nếu một số giai đoạn chuẩn bị mẫu thử bị chia sẻ bởi nhiều mẫu, đến nỗi độ chệch hoặc độ lệch nảy sinh ra từ sự chuẩn bị đó sẽ có ảnh hưởng chung đến các kết quả kiểm nghiệm từ các mẫu thử.

5.1.4. Việc thiết kế các vật liệu không đồng nhất được đề xuất trong điều này tạo ra các thông tin về sự thay đổi giữa các mẫu mà không thể thu được từ thiết kế đồng mức được mô tả trong TCVN 6910 -2. Có một phí tổn không thể tránh được để nhận được thông tin phụ thêm: thiết kế đề suất cần có nhiều mẫu hơn để thử nghiệm. Các thông tin phụ thêm này có thể có giá trị. Trong ví dụ về da thuộc ở 5.1.1, thông tin về sự thay đổi giữa các miếng da có thể được sử dụng để quyết định xem bao nhiêu miếng da cần phải kiểm nghiệm khi đánh giá chất lượng của mẫu gửi tới hoặc để quyết định giữa hai sự lựa chọn: kiểm nghiệm nhiều miếng da với ít mẫu thử trên một miếng hoặc kiểm nghiệm ít miếng da với nhiều mẫu thử trên một miếng. Trong ví dụ về cát ở 5.1.2, thông tin về sự thay đổi giữa các mẫu cát có thể được dùng để quyết định xem thủ tục tiến hành lấy mẫu có thỏa đáng chưa hay cần phải có sự cải tiến.

5.1.5. Thiết kế được mô tả trong điều này áp dụng được đối với các thí nghiệm bao gồm ba yếu tố được sắp xếp theo thứ bậc: yếu tố “các phòng thí nghiệm” ở mức cao nhất, yếu tố “các mẫu trong phòng thí nghiệm” ở mức sau và yếu tố “kết quả thử nghiệm trong các mẫu” ở mức thấp nhất. Một trường hợp khác có thể gặp phải trong thực tế là trường hợp xếp thứ bậc ba yếu tố: “phòng thí nghiệm” ở mức cao nhất, “các kết quả thử nghiệm trong các phòng thí nghiệm” ở mức sau. “sự xác định giá trị trong các kết quả thử nghiệm” ở mức thấp nhất. Điều này sẽ xuất hiện nếu các phòng thí nghiệm tham gia trong thí nghiệm về độ chụm được gửi tới một mẫu đơn lẻ của một vật liệu đồng nhất và yêu cầu phải thực hiện hai (hoặc có thể nhiều hơn) phép thử nghiệm trên mẫu và mỗi phép thử nghiệm gồm một số lần xác định giá trị và kết quả thử nghiệm cần tính là trung bình của các lần xác định. Các công thức cho trong 5.5, 5.6 và 5.9 có thể được áp dụng cho các số liệu thu được trong một thí nghiệm như vậy, nhưng các độ lệch chuẩn lặp lại và tái lập cần phải được tính theo cách khác hơn một chút so với cách nêu ra ở đây (xem chú thích 2 của 5.5.5). Cần định rõ số lần xác định giá trị khi tính trung bình để đưa ra một kết quả thử nghiệm, vì điều này ảnh hưởng tới giá trị của các độ lệch chuẩn lặp lại và tái lập.

...

5.2.1. Việc bố trí thiết kế các vật liệu không đồng nhất được trình bày trong bảng 9.

Mỗi phòng trong p phòng thí nghiệm tham gia được cung cấp hai mẫu ở q mức và thu được hai kết quả thử nghiệm trên mỗi mẫu. Mỗi ô như vậy trong thí nghiệm chứa bốn kết quả thử nghiệm (hai kết quả đối với mỗi mẫu).

Có thể mở rộng thiết kế đơn giản đó bằng cách cho phép có nhiều hơn hai mẫu trong một phòng thí nghiệm ở mỗi mức hoặc có nhiều hơn hai kết quả thử nghiệm trên mỗi mẫu. Việc tính toán thiết kế mở rộng phức tạp hơn nhiều so với thiết kế với hai kết quả thử nghiệm trên một mẫu và hai mẫu trên một phòng thí nghiệm ở mỗi mức. Tuy nhiên các nguyên tắc thiết kế mở rộng giống hệt như các nguyên tắc thiết kế đơn giản. Vì vậy những tính toán sẽ được thiết lập chi tiết ở trong điều này cho thiết kế đơn giản. Các công thức để tính toán giá trị của độ lệch chuẩn lặp lại và tái lập với thiết kế mở rộng và ví dụ áp dụng sẽ được đưa ra trong 5.9 và 5.10 dưới đây.

5.2.2. Các số liệu từ thiết kế vật liệu không đồng nhất được biểu thị bởi

y_ijtk

trong đó

chỉ số i biểu thị phòng thí nghiệm (i = 1.... p')

chỉ số j biểu thị mức (j = 1.2,...., q).

chỉ số t biểu thị mẫu (t = 1 .2,...., g)

...

Thông thường g = 2, n = 2. trong thiết kế mở rộng cả hai số g hoặc n đều lớn hơn 2.

Chú thích - Trong TCVN 6910 -1 và TCVN 6910-2 (ISO 5725-2), p được dùng để biểu thị số các phòng thí nghiệm cũng như chỉ số trong các bảng các giá trị tới hạn đối với phép kiểm nghiệm Cochran: Với thiết kế đồng mức cả hai số là như nhau. Với thiết kế với các vật liệu không đồng nhất chỉ số đối với phép kiểm nghiệm Cochran có thể là bội số của số các phòng thí nghiệm. Vì vậy ở đây p' dùng để chỉ số lượng các phòng thí nghiệm còn p biểu thị chỉ số đối với phép kiểm nghiệm Cochran.

5.3. Tổ chức thí nghiệm với vật liệu không đồng nhất

5.3.1. Phải tuân theo hướng dẫn đã nêu trong điều 6 của TCVN 6910 -1: 2001 khi thiết kế thí nghiệm với vật liệu không đồng nhất, cần xét thêm các vấn đề sau:

Cần phải chuẩn bị bao nhiêu mẫu cho mỗi phòng thí nghiệm ở mỗi mức?

Thông thường, vì phải tính đến chi phí nên số mẫu cần chuẩn bị là hai.

Có thể sử dụng các công thức, bảng, hình vẽ trong điều 6 và phụ lục B của TCVN 6910 -1: 2001 để chọn số các phòng thí nghiệm, số mẫu và số lần lặp lại, nhưng với sự cải biên sẽ được thiết lập trong 5.3.2 và 5.3.5.

5.3.2. Có thể đánh giá độ không đảm bảo của ước lượng độ lệch chuẩn lặp lại và tái lập từ một thí nghiệm trên vật liệu không đồng nhất bằng cách tính đại lượng A_r như đã trình bầy trong 6.3 của TCVN 6910 -1: 2001 (ISO 5725-1:1994):

...

Chú thích - Công thức cho A_r ở trên ( và A_R dưới đây) đã được dẫn ra theo phương pháp đã mô tả trong chú thích 24 của 6910 - 1: 2001.

5.3.3. Có thể đánh giá độ không đảm bảo của ước lượng độ lệch chuẩn lặp lại và tái lập từ một thí nghiệm trên một vật liệu không đồng nhất bằng cách tính đại lượng A_R như đã được trình bày trong 6.3 của TCVN 6910-1: 2001 (ISO 5725-1: 1994) như sau:

thay cho đại lượng được xác định bởi phương trình (10) trong TCVN 6910-1: 2001 (ISO 5725-1:1994).

Trong đó:

Có thể dẫn ra các giá trị của F và g từ các ước lượng sơ bộ của các đô lệch chuẩn s_H, s_R và s_r thu được trong quá trình chuẩn hóa phương pháp đo.

5.3.4. Chú ý tới hướng dẫn được đưa ra trong điều 5 và 6 của TCVN 6910-2:2001 (ISO 5725-2:1994) với việc xem xét các chi tiết về tổ chức thí nghiệm với vật liệu không đồng nhất.

Điều 5.1.2 của TCVN 6910-2:2001 (ISO 5725-2:1994) bao gồm các yêu cầu “một nhóm gồm n thử nghiệm” hoặc “một nhóm gồm n lần đo” ( ví dụ một nhóm gồm n thử nghiệm cần phải được tiến hành với các điều kiện lặp lại). Trong một thí nghiệm với vật liệu không đồng nhất, các yêu cầu đó có liên quan tới nhóm gồm g x n thử nghiệm trong một ô, tức là liên quan tới tất cả các thử nghiệm trong một phòng thí nghiệm ở một mức.

...

5.4. Mô hình thống kê cho một thí nghiệm với vật liệu không đồng nhất

5.4.1. Mô hình cơ bản được dùng trong tiêu chuẩn này chính là phương trình (3) trong 4.4.1 đã được trình bày lại. Đối với một thí nghiệm với vật liệu không đồng nhất mô hình đó được mở rộng thành

Các số hạng m, B và e có cùng ý nghĩa như trong phương trình (3) trong 4.4.1. Nhưng phương trình (19) còn có thêm số hạng H_ijt biểu thị sự thay đổi giữa các mẫu và chỉ số t biểu thị các mẫu trong các phòng thí nghiệm (ý nghĩa của các chỉ số khác đã được nêu ra trong 5.2.2)

Giả thiết sau đây là hợp lý sự thay đổi giữa các mẫu là ngẫu nhiên và không phụ thuộc vào phòng thí nghiệm, nhưng có thể phụ thuộc vào mức thí nghiệm, như vậy số hạng H_ijt trung bình bằng không và phương sai được cho bởi:

5.4.2. Đối với trường hợp thông thường với hai mẫu trong một phòng thí nghiệm và hai kết quả thử nghiệm trên một mẫu (g = n = 2), ta xác định:

a) Trung bình mẫu, phạm vi giữa các kết quả thử nghiệm cho phòng thí nghiệm i, mức j và mẫu t (t= 1 hoặc 2)

...

c) Trung bình chung, độ lệch chuẩn của các trung bình ô cho mức j

trong đó tổng được lấy theo các phòng thí nghiệm i = 1,2,.....,p'.

5.5. Phân tích thống kê các số liệu của một thí nghiệm với vật liệu không đồng nhất

5.5.1. Ở đây khảo sát chi tiết trường hợp thông thường khi hai mẫu được chuẩn bị cho mỗi phòng thí nghiệm ở mỗi mức và trên mỗi mẫu thu được hai kết quả thử nghiệm. (Trường hợp tổng quát sẽ được khảo sát ở 5.9 và 5.10).

Tập hợp các số liệu thành bảng như đã chỉ ra trong bảng 9. Mỗi sự kết hợp của một phòng thí nghiệm và một mức sẽ cho một “ô” trong bảng này, mỗi ô chứa bốn kết quả thử nghiệm.

Sử dụng các phương trình (21) đến (26) trong 5.4.2:

a) tính các phạm vi giữa kết quả thử nghiệm và đưa chúng vào một bảng như bảng 10

...

c) tính các giá trị trung bình ô và đưa chúng vào một bảng như bảng 12

Ghi lại các phạm vi như là các giá trị dương (không kể dấu).

Bảng 9 - Biểu mẫu nên dùng để tập hợp các số liệu của thiết kế cho vật liệu không đồng nhất

Phòng thí nghiệm

Mẫu

Mức 1

Mức 2

Mức j

...

Mức q

Số các kết quả thử nghiệm

...

Bảng 10 - Biểu mẫu nên dùng để lập bảng các phạm vi giữa các kết quả thử nghiệm cho thiết kế với vật liệu không đồng nhất

Phòng thí nghiệm

Mẫu

...

Mức 2

Mức j

Mức q

...

Phòng thí nghiệm

Mức 1

Mức 2

Mức j

Mức q

...

Phòng thí nghiệm

Mức 1

Mức 2

Mức j

Mức q

...

a) hoặc sử dụng công thức cho trường hợp tổng quát được nêu ra về sau...

b) hoặc bỏ qua tất cả các số liệu trong ô.

Lựa chọn a) được ưa chuộng hơn. Lựa chọn b) lãng phí số liệu nhưng cho phép sử dụng công thức đơn giản.

5.5.3. Với mỗi mức j của thí nghiệm tính các đại lượng sau

a) Tổng các bình phương của các phạm vi giữa các kết quả thử nghiệm trong cột j của bảng 10 (lấy tổng trên p' phòng thí nghiệm và trên 2 mẫu)

b) Tổng các bình phương các phạm vi giữa các mẫu trong cột j của bảng 11 (lấy tổng trên p' phòng thí nghiệm)

c) Trung bình và độ lệch chuẩn của các trung bình ô trong cột j của bảng 12, khi sử dụng phương trình (25) và (26) trong 5.4.2.

...

5.5.5. Tính toán các độ lệch chuẩn lặp lại s_r và độ lệch chuẩn tái lập s_Rj từ:

nếu xẩy ra

thì xác định

Tính một ước lượng s_Hj của độ lệch chuẩn mà nó đo sự biến đổi giữa các mẫu như sau:

...

1) Nên thực hiện việc kiểm nghiệm về ý nghĩa để thấy rõ sự thay đổi giữa các mẫu là có ý nghĩa về mặt thống kê hay không. Tuy nhiên điều đó không phải là một phần cần thiết của việc phân tích. Việc sử dụng phép kiểm nghiệm đó để quyết định xem sự thay đổi giữa các mẫu có thể bỏ qua được hay không (tức là quyết định xem các kết quả kiểm nghiệm trong mỗi ô có được xử lý như khi tất cả chúng nhận được trên cùng một mẫu hay không trong việc pnân tích là điều không đúng. Điều này có thể dẫn đến độ chệch trong việc ước lượng độ lệch chuẩn lặp lại, bởi vì việc phát hiện sự biến đổi giữa các mẫu là không có ý nghĩa về mặt thống kê không chứng tỏ rằng sự thay đổi giữa các mẫu là không đáng kể.

2) Trong trường hợp đã mô tả ở 5.1.5 (khi ba yếu tố là “phòng thí nghiệm”, “các thử nghiệm trong phòng thí nghiệm”, “xác định giá trị trong các thử nghiệm”) độ lệch chuẩn lặp lại và tái lập cần được tính như sau:

Các công thức này được áp dụng khi các kết quả thử nghiệm là trung bình của hai lần xác định giá trị.

5.5.6. Nghiên cứu xem s_rj, và s_Rj có phụ thuộc vào trung bình chung y_j hay không, nếu phụ thuộc, hãy xác định quan hệ hàm số đó bằng cách sử dụng phương pháp đã nêu trong 7.5 của TCVN 6910-2:2001 (ISO 5725-2:1994)

5.6. Khảo sát sự phù hợp và sự bất thường của các số liệu

5.6.1. Sử dụng các thông kê h và k như đã mô tả trong 7.3.1 của TCVN 6910-2:2001 (ISO 5725- 2:1994)

Để kiểm tra sự phù hợp của các trung bình ô, tính thống kê h như sau:

...

Để kiểm tra sự phù hợp của các phạm vi giữa các mẫu, tính các giá trị thống kê k như sau:

Vẽ biểu đồ các giá trị thống kê đó theo thứ tự các mức, nhưng gộp lại theo phòng thí nghiệm để chỉ ra phòng thí nghiệm không phù hợp

Để kiểm tra sự phù hợp của các phạm vi giữa các kết quả thử nghiệm, tính các giá trị thống kê k như sau:

Vẽ biểu đồ các giá trị thống kê đó theo thứ tự các mức, nhưng gộp lại theo phòng thí nghiệm để chỉ ra phòng thí nghiệm không phù hợp.

Việc giải thích các đồ thị đó được bàn đến một cách đầy đủ trong 7.3.1 của TCVN 6910-2:2001 (ISO 5725-2:1994) Nếu một phòng thí nghiệm báo cáo các kết quả nói chung bị chệch, thì hầu hết các thống kê h đối với phòng thí nghiệm đó sẽ lớn và theo cùng một hướng được thể hiện trong đồ thị lập nên từ các trung bình ô. Nếu một phòng thí nghiệm không thực hiện việc thử nghiệm trong các mức với các điều kiện lặp lại (và cho phép các nhân tố lạ làm tăng thêm sự thay đổi giữa các mẫu) thì các thống kê k với các giá trị lớn khác thường sẽ hiện rõ trên đồ thị được dựng lên từ các phạm vi giữa các mẫu. Nếu một phòng thí nghiệm có độ lặp lại xấu thì nó sẽ đưa ra các thống kê k với các giá trị lớn khác thường được hiện rõ trên đồ thị được dựng lên từ các phạm vi giữa các kết quả thử nghiệm.

5.6.2. Sử dụng phép kiểm nghiệm Cochran như được mô tả trong 7.3.3 và 7.3.4 của TCVN 6910- 2:2001 (ISO 5725-2:1994) để kiểm tra sự tản mạn và bất thường của các số liệu.

Để kiểm tra sự tản mạn và bất thường của các phạm vi giữa các kết quả thử nghiệm, tính giá trị thống kê Cochran cho mỗi mức j:

...

trong đó w_max là giá trị lớn nhất của các phạm vi giữa các kết quả thử nghiệm w_ijt cho mức j.

Để sử dụng bảng các giá trị tới hạn trong 8.1 của TCVN 6910-2:2001 (ISO 5725-2:1994) tra bảng đó theo hàng ứng với p = 2p' theo rìa biên bên trái và theo cột có đầu đề n = 2.

Để kiểm nghiệm các giá trị tản mạn và bất thường theo các phạm vi giữa các mẫu, tính giá trị thống kê Cochran cho mỗi mức j:

trong đó w_max là giá trị lớn nhất trong số các phạm vi giữa các mẫu w_ij đối với mức j.

Để sử dụng các giá trị tới hạn trong TCVN 6910-2:2001 (ISO 5725-2:1994) tra bảng đó theo hàng ứng với p = p' theo rìa biên bên trái và theo cột có đầu đề n = 2.

Để kiểm nghiệm các giá trị tản mạn và bất thường theo các trung bình ô, ta tính giá trị thống kê Grubb từ các trung bình ô như đã chỉ ra trong 7.3.4 của TCVN 6910-2:2001 (ISO 5725-2:1994) cho mỗi mức j (trong đó s trong TCVN 6910-2: 2001 (ISO 5725-2: 1994) chính là s_yj được xác định bởi phương trình (26) trong 5.4.2)

Việc diễn giải các phép kiểm nghiệm đó được thảo luận một cách đầy đủ trong 7.3.2 của TCVN 6910- 2:2001 (ISO 5725-2:1994). Trong thí nghiệm với vật liệu không đồng nhất, các kết quả của việc áp dụng các phép kiểm nghiệm đó cần phải tiến hành theo thứ tự sau đây. Trước tiên cần áp dụng phép kiểm nghiệm Cochran đối với các phạm vi giữa các kết quả trung bình. Dựa trên phép kiểm nghiệm đó nếu quyết định rằng một phạm vi giữa các kết quả thử nghiệm là bất thường và bị loại bỏ, thì hai kết quả thử nghiệm xác định phạm vi bất thường đó cũng cần phải được loại bỏ khi tính độ lệch chuẩn lặp lại và tái lập (nhưng các kết quả thử nghiệm khác trong ô vẫn cần phải được giữ lại). Tiếp theo hãy áp dụng phép kiểm nghiệm Cochran đối với các phạm vi giữa các mẫu và cuối cùng áp dụng phép kiểm nghiệm Grubb đối với các trung bình ô. Nếu quyết định rằng một phạm vi giữa các mẫu là bất thường hoặc một trung bình ô là bất thường và rằng các kết quả gây ra những bất thường như vậy cần loại bỏ, thì cùng loại bỏ tất cả các kết quả đối với các ô tương ứng khi tính độ lệch chuẩn lặp lại và tái lập.

5.7. Báo cáo kết quả của một thí nghiệm trên vật liệu không đồng nhất

...

5.8. Ví dụ 2: Thí nghiệm trên vật liệu không đồng nhất

5.8.1. Tập hợp vật chất được sử dụng làm vật liệu (lớp bê tông hoặc bitum) để tạo nên bề mặt của các phi trường và các xa lộ cần phải có khả năng chịu thấm nước và chống đóng băng. Phương pháp thường dùng để đo khả năng đó là phương pháp thử nghiệm về sự không bị hư hỏng bởi sunphat manhê ^[2] trong đó một mẫu thử của tập hợp chất bị nhúng vào dung dịch sunphát manhê bão hòa và sau đó được làm khô một số lần. Ban đầu mẫu thử được chuẩn bị sao cho tất cả các hạt của tập hợp vật chất đều được giữ lại trên sàng với lưới cỡ 10,0 mm. Việc xử lý trong thử nghiệm đó sẽ làm cho các hạt bị nhỏ dần và kết quả thử nghiệm là tỷ lệ khối lượng hạt của lô thử lọt qua sàng đó ở cuối cuộc thử nghiệm. Một tỷ lệ cao trong thử nghiệm (vượt quá 10 % đến 25 %) chứng tỏ rằng tập hợp các chất đó có độ không hư hại thấp.

5.8.2. Các số liệu chỉ ra trong bảng 13 là các số liệu thu được trong một thí nghiệm, trong đó các cặp mẫu của tám tập hợp chất được gửi tới mười một phòng thí nghiệm và hai kết quả thử nghiệm về độ không bị hư hại bởi sunphat manhê thu được trên mỗi mẫu. Các mẫu có khối lượng xấp xỉ 100 kg (chúng được dùng cho một số lần thử nghiệm khác nhau) và các mẫu thử có khối lượng xấp xỉ 350 g.

5.8.3. Bảng 14, 15 và 16 chỉ ra phạm vi giữa các kết quả thử nghiệm, phạm vi giữa các mẫu và các trung bình ô được tính theo các công thức (21) đến (24) trong 5.4 2 chỉ đối với mức 6 của thí nghiệm.

Nếu sử dụng các phương trình (27) và (28) trong 5.5.3, phạm vi giữa các kết quả thử nghiệm trong bảng 14 và phạm vi giữa các mẫu trong bảng 15 sẽ cho:

Nếu áp dụng các phương trình (25) và (26) trong 5.4.2 đối với các trung bình ô trong bảng 16 ta được:

y₆ = 19.0% (trung bình chung)

s_r6 = 5.03%

...

Bảng 17 đưa ra các kết quả tính toán cho các mức khác

5.8.4. Hình 4 chỉ rõ biểu đồ của các phạm vi giữa các kết quả thử nghiệm, các phạm vi giữa các mẫu, và các trung bình ô ở mức 6. Các đồ thị loại này đưa ra một bức tranh dễ hiểu về lượng biến đổi nẩy sinh từ các nguồn khác nhau (giữa các kết quả thử nghiệm, giữa các mẫu và giữa các phòng thí nghiệm). Hình 4 chỉ ra rằng, trong thí nghiệm này, ở mức 6 có sự thay đổi lớn giữa các trung bình ô mà nếu phương pháp thử nghiệm đó được sử dụng trong một quy định sẽ dẫn đến sự bất đồng do sự khác nhau giữa các kết quả của chúng giống như sự bất đồng giữa người bán và người mua, phạm vi giữa các mẫu nhỏ hơn phạm vi giữa các kết quả thử nghiệm, điều này gợi ý rằng sự biến đổi giữa các mẫu ở mức 6 là không đáng kể.

5.8.5. Giá trị của các thống kê h và k được tính như ở 5.6.1 cũng được chỉ ra trong bảng 14, 15 và 16 cho mức 6. Các giá trị đó cho tất cả các mức được vẽ thánh biểu đồ trong hình 5 đến 7. (Trong các hình đó các mức đã được sắp xếp lại theo thứ tự tăng dần của các trung bình chung như được chỉ ra trong bảng 17). Hình 5 chỉ ra rằng phòng thí nghiệm 6 đã thu được nhiều thống kê k với giá trị lớn đối với phạm vi giữa các kết quả thử nghiệm, điều này chứng tỏ rằng phòng thí nghiệm 6 có độ lặp lại xấu hơn so với các phòng thí nghiệm khác. Hình 6 chỉ ra rằng 3 phòng thí nghiệm (1, 6 và 10) thu được các thống kê k với các giá trị lớn đối với phạm vi giữa các mẫu: Điều này có lẽ là do các phòng thí nghiệm đó đã không tuân theo các thủ tục đã được khuyến cáo để chuẩn bị các mẫu thử từ các mẫu chung. Hình 7 chỉ ra rằng các thống kê h dương và âm trong hầu hết các phòng thí nghiệm (với các phòng thí nghiệm 1, 6,10, các thống kê đó có giá trị lớn nhất). Rõ ràng rằng có các độ chệch nhất quán theo hầu hết các phòng thí nghiệm, điều đó chứng tỏ rằng phương pháp thử nghiệm đó chưa được quy định một cách thích hợp.

5.8.6. Việc áp dụng các phép kiểm nghiệm Cochran và Grubb như trong 5.6.2 đối với các số liệu sẽ thu được các kết quả ở bảng 18. Hai giá trị bất thường được phát hiện. Do thiếu các thông tin khác, nên số liệu ứng với các giá trị bất thường đó cần phải được loại ra và cần phải tính lại. Sau đó cần tiếp tục việc phân tích với sự khảo sát mối quan hệ hàm số theo cách giống như với thiết kế đồng mức trong TCVN 6910-2 (ISO 5725-2).

Các trung bình ô

...

Phạm vi giữa các kết quả thử nghiệm

Hình 4 - Ví dụ 2 : Biểu đồ của các phạm vi và trung bình trích từ các bảng 14, 15 và 16 cho mức 6

Hình 5 - Ví dụ 2: Kiểm tra sự phù hợp đối với phạm vi giữa các kết quả thử nghiệm

(được ghép nhóm theo phòng thí nghiệm)

Hình 6 - Ví dụ 2: Kiểm tra sự phù hợp đối với các phạm vi giữa các mẫu

(được ghép nhóm theo phòng thí nghiệm)

...

(được ghép nhóm theo phòng thí nghiệm)

Bảng 13 - Ví dụ 2: Xác định sự không bị hư hại bởi sunphat manhê (%)

Phòng thí nghiệm

Mẫu

Mức1

Mức 2

Mức 3

Mức 4

Số kết quả thử nghiệm

...

67,0

7,4

8,0

4,1

3,5

10,4

10,1

69,7

...

6,6

5,7

10,5

13,1

13,9

13,8

66,5

...

1,9

2,1

3,0

3,2

8,7

6,7

65,7

65,8

...

3,3

1,9

1,1

8,3

4,8

68,7

69,5

...

5,8

2,4

2,9

11,7

7,0

67,7

77,7

9,7

...

2,1

3,3

7,9

12,0

77,5

75,3

2,0

...

2,4

1,4

9,4

7,1

76,3

77,2

4,7

3,8

...

2,3

10,7

7,7

55,4

63,2

3,8

4,1

...

0,8

3,7

6,3

65,9

54,7

2,1

3,1

0,7

...

3,3

3,7

64,8

70,9

8,4

6,1

6,0

...

16,5

12,3

78,2

73,4

8,3

10,6

12,4

9,8

...

16,8

64,8

63,4

4,3

5,7

2,9

3,0

...

9,3

67,0

63,4

7,7

3,9

4,3

6,4

11,1

...

64,9

68,4

4,4

2,8

1,3

2,8

5,7

...

65,4

65,5

5,4

6,7

2,7

2,8

4,8

5,5

...

1,1

0,0

6,6

7,0

...

0,7

3,7

4,9

6,3

...

57,0

57,7

3,3

0,4

2,1

2,4

5,5

5,8

...

52,7

4,2

2,3

3,6

3,5

3,9

5,7

...

75,2

5,3

6,4

5,7

1,9

9,5

7,2

77,9

...

3,5

7,1

1,4

3,0

8,1

7,4

Phòng thí nghiệm

Mẫu

Mức 5

...

Mức 7

Mức 8

Số kết quả thử nghiệm

...

8,9

7,4

31,1

28,5

38,7

41,7

...

4,1

7,6

9,1

23,0

23,1

44,2

41,1

7,3

...

3,2

3,5

16,5

15,4

36,6

45,2

3,2

...

2,8

4,0

10,3

12,8

43,2

40,5

1,7

2,5

...

4,4

6,1

24,3

16,7

38,9

43,1

3,7

7,7

...

6,0

20,8

22,2

46,1

47,4

3,5

5,6

...

2,7

3,1

20,2

16,2

32,0

35,5

2,9

2,2

...

2,9

20,0

11,9

26,5

35,7

3,2

2,3

...

1,4

13,8

15,1

36,7

39,5

1,1

1,2

1,5

...

11,5

13,3

37,6

34,1

0,6

1,7

8,2

...

20,3

24,7

49,4

50,6

11,9

18,5

3,7

4,6

...

18,9

48,2

52,4

14,9

8,1

3,1

5,5

...

23,3

38,9

29,9

1,7

5,6

5,5

21,5

...

34,4

38,3

2,2

5,0

1,8

2,2

13,6

...

27,0

37,0

0,3

2,2

4,0

15,6

16,7

...

34,6

3,6

3,7

3,8

17,7

17,1

...

33,1

1,8

2,0

3,5

2,8

21,4

16,8

26,5

...

2,5

1,6

3,5

3,0

21,7

23,9

35,3

...

0,5

4,3

3,2

3,5

27,0

32,5

18,0

18,2

...

2,1

3,5

2,5

11,0

18,4

27,0

33,5

...

3,9

2,0

2,8

16,4

8,1

35,4

29,3

2,1

...

Bảng 14 - Ví dụ 2: Phạm vi giữa các kết quả thử nghiệm cho mức 6

Phòng thí nghiệm

Mẫu

Phạm vi giữa các kết quả thử nghiệm

Thống kê k

2,6

...

0,1

0,024

1,1

0,264

2,5

...

7,6

1,825

1,4

0,336

...

0,960

8,1

1,945

1,3

0,312

...

0,432

4,4

1,056

2,1

0,504

...

3,9

0,936

1,2

0,288

1,6

0,384

...

1,1

0,264

0,6

0,144

4,6

1,104

...

2,2

0,528

5,5

1,320

7,4

...

8,1

1,945

Bảng 15 - Ví dụ 2: Phạm vi giữa các mẫu cho mức 6

Bảng 16 - Ví dụ 2: Các trung bình ô cho mức 6

Phòng thí nghiệm

Phạm vi giữa các mẫu

...

Thống kê k

Phòng thí nghiệm

Trung bình ô

Thống kê h

6,75

1,776

...

26,425

1,475

4,40

1,152

13,750

...

1,00

0,262

21,000

0,397

2,25

...

17,075

-0,382

2,05

0,537

...

-1,108

2,55

0,668

21,225

0,442

...

0,825

23,675

0,929

3,35

0,877

...

14,475

-0,899

1,70

0,445

18,250

-0,149

...

6,95

1,819

26,275

1,445

2,55

0,668

...

13,425

1,108

Bảng 17 - Ví dụ 2: Các giá trị của các trung bình, tổng bình phương của các phạm vi, và độ lệch chuẩn được tính từ các số liệu cho tất cả 8 mức trong bảng 13

(loại ra các ô vài các số liệu bị thiếu)

Mức

Số các phòng thí nghiệm

...

y_j%

Tổng bình phương các phạm vi

Các độ lệch chuẩn

SS_rj(%)²

SS_HJ(%)²

s_yj%

s_rj%

s_Rj%

s_Hj%

...

3,7

82,99

96,3725

2,62

1,37

2,56

1,85

...

4,0

34,70

11,2550

1,88

0,89

2,01

0,34

...

155,39

29,4225

3,49

1,97

3,92

0,00

5,0

...

25,2375

1,95

1,44

2,29

0,47

8,2

131,07

...

3,10

1,73

3,47

0,00

19,0

381,66

160,5300

...

2,95

5,51

1,72

36,5

636,19

305,4775

7,28

...

7,78

2,58

67,4

529,71

92,9225

6,23

3,64

...

0,00

Bảng 18 - Ví dụ 2: Các giá trị của thống kê Grubb và Cochran

Mức

Số các phòng thí nghiệm

Thống kê Cochran cho các phạm vi giữa các kết quả thử nghiệm

Thống kê Cochran cho các phạm vi giữa các mẫu

...

0,203

0,664 * (1)

0,461 ** (6)

0,374

0,298

...

0,232

0,238

0,169

0,550

...

0,172

0,301

0,157

0,536

0,237

...

Các thống kê Grubb cho các trung bình ô

Mức

Số các phòng thí nghiệm

Một giá trị nhỏ nhất

Hai giá trị nhỏ nhất

Hai giá trị lớn nhất

Một giá trị lớn nhất

...

0,970

0,791

0,098** (1;6)

2,219

1,396

0,709

...

2,266

0,849

2,643**(6)

...

0,614

0,466

1,713

1,290

0,681

0,294

2,082

...

1,108

0,700

0,479

1,475

1,649

0,562

...

1,875

1,808

0,345

0,590

1,476

Chú thích - Các số trong ngoặc biểu thị các phòng thí nghiệm đưa ra các giá trị tản mạn tiêu chuẩn là như sau:

Phép thử nghiệm thống kê

...

Số các phòng thí nghiệm

Chỉ số trong bảng TCVN 6910-2 (ISO 5725-2)

Giá trị tản mạn

Giá trị bất thường

p’

(*)

...

Phép kiểm nghiệm Cochran

Phạm vi giữa các kết quả thử nghiệm

0,389

0,365

0,480

...

Phép kiểm nghiệm Cochran

Phạm vi giữa các mẫu

0,602

0,570

0,718

...

Phép kiểm nghiệm Grubb cho một giá trị bất thường

Trung bình ô

2,290

2,355

2,482

...

Phép kiểm nghiệm Grubb cho cặp giá trị bất thường

Trung bình ô

0,186 4

0,2213

0,1150

...

5.9. Các công thức tổng quát để tính toán với thiết kế cho vật liệu không đồng nhất

Tính các thống kê sau đây ở mỗi mức j:

a) Trung bình chung (với tổng lấy theo i, t, k)

trong đó n_j là số các kết quả thử nghiệm bao gồm trong tổng

b) Các hiệu quả phòng thí nghiệm với mỗi i (với tổng lấy theo t và k):

= trung bình phòng thí nghiệm - trung bình chung.

trong đó n_ij là số các kết quả thử nghiệm bao gồm trong tổng:

...

= trung bình mẫu - trung bình phòng thí nghiệm.

trong đó n_ijt là số các kết quả thử nghiệm bao gồm trong tổng.

d) Các phần dư đối với i, t, k:

= kết quả thử nghiệm - trung bình mẫu.

e) Tổng bình phương đối với các mẫu (với tổng lấy theo i):

f) Tổng bình phương đối với các mẫu (với tổng lấy theo i, t):

...

g) Tổng bình phương đối với sự lặp lại (với tổng lấy theo i, t và k):

h) Các bậc tự do:

trong đó

p’_j là số phòng thí nghiệm có báo cáo ít nhất một kết quả thử nghiệm.

g_i là số các mẫu có ít nhất một kết quả thử nghiệm được báo cáo.

n là tổng số các kết quả thử nghiệm.

i) Các thừa số đối với mỗi i (với tổng lấy theo t):

...

j) Các thừa số (với tổng lấy theo i):

k) Độ lệch chuẩn lặp lại s_r, độ lệch chuẩn giữa các mẫu s_Hj, độ lệch chuẩn giữa các phòng thí nghiệm s_Lj và độ lệch chuẩn tái lập s_RJ (khi sử dụng ) được cho bởi:

Chú thích - Các công thức trên được dẫn ra bằng cách sử dụng lý thuyết thống kê của Scheffé [4]

5.10 Ví dụ 3: Áp dụng các công thức tổng quát

5.10.1. Các số liệu của ví dụ 2, mức 4, được sử dụng để làm ví dụ áp dụng các công thức tổng quát bằng cách bỏ qua một số các kết quả thử nghiệm (xem bảng 19). Công thức trong 5.9 cho trung bình chung tính theo công thức được cho ở 5.9 phía dưới bảng 19; tổng các bình phương, các bậc tự do, các thừa số được chỉ ra ở phía dưới bảng 20, 21 và 22.

5.10.2. Khi áp dụng các phương trình (52) -(55) ở bước k của 5.9 ta thu được:

...

và

như vậy

và

như vậy

...

Bảng 19 - Ví dụ 3: Xác định độ không bị hư hại bởi sunphát manhê ở mức 4

Phòng thí nghiệm

Mẫu

Kết quả thử nghiệm k = 1

Kết quả thử nghiệm k = 2

...

10,1

13,9

13,8

...

8,3

4,8

7,0

...

12,0

9,4

...

3,7

6,3

3,3

3,7

16,5

12,3

...

13,2

16,8

7,5

9,3

11,1

8,3

...

5,7

6,8

4,8

5,5

6,6

...

4,9

6,3

5,5

5,8

3,9

...

9,5

7,2

8,1

7,4

Trung bình chung

m_j= 8,1111 %

...

Số các kết quả thử nghiệm

n_j = 36

Bảng 20 - Ví dụ 3: Tính tổng các bình phương cho các phòng thí nghiệm

Phòng thí nghiệm

Trung bình, phòng thí nghiệm

Số các kết quả thử nghiệm

...

Hiệu quả phòng thí nghiệm

B_ij %

Thừa số

K_ij

12,600

4,488 9

...

6,550

- 1,561 1

9,500

1,388 9

...

9,400

1,288 9

4,250

- 3,861 1

...

14,700

- 6,588 9

9,050

0,938 9

...

5,700

- 2,411 1

6,200

- 1,911 1

...

5,225

- 2,886 1

8,050

- 0,061 1

...

SS_Lj = 378,853 1 (%)²,

Bậc to do cho các phòng thí nghiệm,

y_Lj = 11 - 1 = 10

Các thừa số

K_j = 130 K’_j = 68 K”_j = 19,667

Bảng 21 - Ví dụ 3: Tính tổng bình phương cho các mẫu

Phòng thí nghiệm

Mẫu

...

Trung bình mẫu

Số các kết quả thử nghiệm

n_ijt

Hiệu quả mẫu

H_ijt

...

- 2,500

13,85

1,250

...

6,55

0,000

7,00

...

12,00

2,500

9,40

0,000

...

5,00

0,750

...

- 0,750

14,40

- 0,300

15,00

...

0,300

8,40

- 0,650

9,70

...

6,25

0,550

5,15

- 0,550

...

6,80

0,600

5,60

- 0,600

...

5,65

0,425

4,80

-0,425

...

0,300

7,55

- 0,300

Tổng bình phương cho các mẫu

SS_HJ = 29,907 5 (%)²,

...

V_HJ = 20 - 11 = 9

Bảng 22 - Ví dụ 3: Tính tổng bình phương cho độ lặp lại

Phòng thí nghiệm

Mẫu

Kết quả thử nghiệm k = 1

...

Kết quả thử nghiệm k = 2

0,00

0,05

-0.05

...

1,75

...

0,00

...

-1,30

1,30

-0,20

0,20

...

-2,10

-1,80

1,80

-0,90

0,90

...

-1,40

-0,55

0,55

-0,35

0,35

...

-0,20

0,20

-0,70

0,70

-0,15

0,15

...

-0,90

0,90

1,15

-1,15

0,35

-0,35

...

SS_rj = 36,895 (%)².

Bậc tự do cho độ lặp lại

v_rj = 36 - 20 = 16

6. Phương pháp ổn định để phân tích số liệu

6.1. Áp dụng các phương pháp ổn định để phân tích số liệu

6.1.1. TCVN 6910-2 (ISO 5725-2), khuyến nghị cần áp dụng hai phép kiểm nghiệm (Cochran và Grubb) về giá trị bất thường đối với các số liệu thu được trong thí nghiệm về độ chụm. Bất kỳ số liệu nào tạo ra một thống kê kiểm nghiệm trong hai phép kiểm nghiệm trên có giá trị vượt quá giá trị tới hạn ở mức ý nghĩa 1 % cần phải được loại bỏ (trừ khi nhà thống kê có lý do xác đáng để giữ lại số liệu đó). Trên thực tế thường không dễ áp dụng cách này. Xét kết quả của các phép kiểm nghiệm giá trị bất thường trong ví dụ 1 của 4.8. Các kết quả đó được đưa ra trong bảng 8. Phòng thí nghiệm 5 đưa ra một trung bình ở vượt trội duy nhất (ở mức 10) và nó cần được phân loại là giá trị bất thường, theo phép kiểm nghiệm Grubb, nhưng nó đưa ra ba giá trị tản mạn khác và thể hiện rõ trong hình 3 rằng có điều gì đó sai sót trong phòng thí nghiệm này. Trong tình huống này nhà thống kê cần phải quyết định lựa chọn giữa:

a) giữ lại tất cả số liệu của phòng thí nghiệm 5

...

c) loại bỏ tất cả các số liệu của phòng thí nghiệm 5.

Quyết định này sẽ ảnh hưởng thực sự đến các giá trị tính được cho độ lệch chuẩn lặp lại và tái lập. Theo kinh nghiệm chung, khi phân tích số liệu từ các thí nghiệm về độ chụm để phát hiện các số liệu nằm trên biên giới giữa các giá trị tản mạn và bất thường, cần có các phán quyết và điều này ảnh hưởng đến kết quả tính toán. Nhưng điều đó có thể không thỏa đáng. Các phương pháp ổn định được mô tả trong điều này cho phép phân tích các số liệu theo cách sao cho không đòi hỏi thực hiện các quyết định có ảnh hưởng đến kết quả tính toán. Như vậy nếu có lý do để cho rằng các kết quả của thí nghiệm về độ chụm sẽ chứa các số liệu bất thường, phương pháp ổn định có thể được ưa chuộng hơn.

6.1.2. Mô hình cơ bản trình bày trong điều 5 của TCVN 6910-1: 2001 (ISO 5725-1:1994) dựa trên giả thiết rằng đã thiết lập được đúng đắn một giá trị chung cho độ lệch chuẩn lặp lại và tái lập đối với tất cả các phòng thí nghiệm sử dụng phương pháp đo đó. Trong thực tế thường thấy một số phòng thí nghiệm có độ lặp lại xấu hơn các phòng thí nghiệm khác. Ví dụ, hãy xem hình 5 của ví dụ 2 trong 5.8, phòng thí nghiệm 6 rõ ràng có độ lặp lại xấu hơn nhiều so với phòng thí nghiệm 9 trong thí nghiệm đó. Như vậy giả thiết các phòng thí nghiệm đạt được độ lặp lại tương tự tỏ ra không còn đúng trong trường hợp này. Một số người tham gia trong thí nghiệm về độ chụm có thể đạt độ lặp lại xấu trong thời gian đầu, khi một phương pháp đo được quy định cho một thí nghiệm về độ chụm hoặc khi họ có ít kinh nghiệm về phương pháp đo đó, và trong các tình huống đó việc sử dụng phương pháp ổn định sẽ đặc biệt thích hợp.

6.1.3. Mục đích của việc sử dụng phương pháp ổn định khi phân tích số liệu của một thí nghiệm độ chụm là tính các giá trị của độ lệch chuẩn lặp lại và tái lập theo một cách sao cho chúng không bị ảnh hưởng bởi các số liệu bất thường. Nếu những người tham gia thí nghiệm đó coi như được phân thành hai nhóm, một nhóm tạo ra các số liệu với chất lượng tốt, nhóm kia tạo ra các số liệu với chất lượng xấu, khi đó phương pháp ổn định sẽ tạo ra các các giá trị của độ lệch chuẩn lặp lại và tái lập phù hợp với lớp số liệu có chất lượng tốt và không bị ảnh hưởng bởi các số liệu có chất lượng xấu (nếu như lớp các số liệu có chất lượng xấu không quá lớn).

6.1.4. Việc sử dụng phương pháp ổn định để phân tích số liệu không ảnh hưởng đến việc lập kế hoạch, tổ chức hoặc thực thi thí nghiệm về độ chụm. Việc quyết định nên sử dụng phương pháp ổn định hay phương pháp đòi hỏi phải loại bỏ các số liệu bất thường cần được tiến hành bởi các chuyên gia thống kê và phải báo cáo lên hội đồng chuyên môn. Khi sử dụng phương pháp ổn định cần phải áp dụng các phép kiểm nghiệm giá trị bất thường và kiểm tra sự phù hợp được mô tả trong TCVN 6910-2 (ISO 5725-2) hoặc TCVN 6910 - 5 cho các số liệu và các nguyên nhân của giá trị bất thường hoặc các kiểu dáng của các thống kê h và k cần phải được nghiên cứu. Tuy nhiên không cần phải loại bỏ các số liệu do các kết quả của việc kiểm nghiệm và kiểm tra.

6.1.5. Theo phương pháp tính toán các thống kê h, k mô tả trong TCVN 6910-2 (ISO 5725-2) mẫu số trong các thống kê đó là độ lệch chuẩn được tính từ các số liệu đã được báo cáo. Nếu có số liệu bất thường trong các số liệu đó thì chúng sẽ làm cho mẫu số lớn lên và tạo ra sai lệch trong đồ thị của các thống kê đó. Ví dụ, nếu ở một mức của thí nghiệm, một phòng thí nghiệm nào đó đưa ra giá trị trung bình ô bất thường vượt trội hẳn so với các giá trị bất thường khác ở cùng mức đó thì đồ thị của thống kê h sẽ biểu hiện một giá trị h lớn khác thường. Tuy nhiên các thống kê h đối với tất cả các phòng thí nghiệm khác ở mức đó sẽ là nhỏ, thậm chí một số trong các phòng thí nghiệm đã đưa ra các giá trị bất thường. Việc sử dụng trung bình chung để tính thống kê h có thể dẫn tới một hiệu quả tương tự. Việc sử dụng các ước lượng ổn định của độ lệch chuẩn cho các mẫu số trong thống kê h và k và việc sử dụng ước lượng ổn định của các giá trị trung bình chung để tính các thống kê h sẽ tránh được sai lệch đó. Vì vậy nên dùng các ước lượng đó cho mục đích này.

6.1.6. Các số liệu của thí nghiệm độ chụm cho phép tính hai loại thống kê.

a) các trung bình ô, từ đó độ lệch chuẩn có thể tính được, độ lệch chuẩn này là thước đo sự biến đổi giữa các phòng thí nghiệm.

b) các độ lệch chuẩn hoặc các phạm vi (hoặc độ sai khác trong thiết kế mức tách biệt) trong các ô và chúng được kết hợp lại để cho một thước đo về sự biến đổi trong các phòng thí nghiệm.

...

Ví dụ, với các số liệu ở một mức của thiết kế đồng mức đã được trình bày trong TCVN 6910-2 (ISO 5725-2), giai đoạn đầu tiên của việc phân tích là tính giá trị trung bình, độ lệch chuẩn của các kết quả đo trong mỗi ô. Sau đó sử dụng trung bình ô để tính độ lệch chuẩn, dùng làm thước đo sự biến đổi giữa các phòng thí nghiệm. Khi phương pháp ổn định của điều này được sử dụng, việc tính toán này được thực hiện theo “thuật toán A” và các trung bình ô không bị loại ra khi tính toán dù cho kết quả của việc áp dụng phép kiểm nghiệm Grubb là như thế nào. Cũng như vậy trong thiết kế này, các độ lệch chuẩn ô được kết hợp lại với nhau để đưa ra một ước lượng của độ lệch chuẩn lặp lại. Với việc phân tích ổn định, điều đó được thực hiện theo “thuật toán S” và các độ lệch chuẩn ô không bị loại ra khi áp dụng phép kiểm nghiệm Cochran. Với mỗi cách tiếp cận (được mô tả trong TCVN 6910-2 (ISO 5725-2) hoặc ở đây) hai thước đo đó được sử dụng để tính ước lượng của độ lệch chuẩn lặp lại và tái lập theo cùng một cách.

Một ví dụ phức tạp hơn là thiết kế xếp chồng so le 6 yếu tố được nêu trong phụ lục C của TCVN 6910 - 3: 2001. Với thiết kế này giai đoạn đầu tiên của việc phân tích là tính giá trị trung bình của các số liệu cho mỗi phòng thí nghiệm (ở mỗi mức) và chúng được ký hiệu bởi y_i(1),….. y_i(5), và dẫy các phạm vi được ký hiệu bởi w_i(1),…..w_i(5), chúng chứa các thông tin về sự biến đổi gắn với các yếu tố khác nhau được kiểm tra trong thí nghiệm. Để phân tích các số liệu khi sử dụng phương pháp ổn định được trình bày ở đây, “thuật toán A” được áp dụng các các trung bình ô và “thuật toán S” được áp dụng cho dẫy các phạm vi. Sử dụng các thống kê thu được bởi các thao tác đó được sử dụng để thu được các ước lượng của độ chụm lặp lại, độ chụm trung gian và độ lệch chuẩn tái lập, theo cùng một cách như khi áp dụng phương pháp phân tích được nêu trong TCVN 6910 - 3.

6.1.7. Các phương pháp ổn định được chọn trong tiêu chuẩn này vì chúng có thể áp dụng cho tất cả các thiết kế thí nghiệm được nêu ra trong các phần 2, 3, 4 và 5 của TCVN 6910 và vì việc tính toán tương đối đơn giản. Tuy nhiên cần lưu ý rằng các phương pháp này đưa ra những cách kết hợp một cách ổn định các trung bình ô, độ lệch tiêu chuẩn ở và các phạm vi ô. Chúng không kết hợp các kết quả thử nghiệm đơn lẻ theo cách ổn định, tức là chúng bắt đầu với các trung bình số học của các trung bình ô và của các độ lệch chuẩn ô. Có những phương pháp ổn định kết hợp các kết quả thử nghiệm trong các ô theo cách ổn định và việc áp dụng chúng trở nên phức tạp hơn trong thực hành.

6.2. Phân tích ổn định: Thuật toán A

6.2.1. Thuật toán này tạo ra các giá trị ổn định của giá trị trung bình và độ lệch chuẩn của các số liệu được dùng và nó được áp dụng cho:

a) các trung bình ô, theo một thiết kế bất kỳ;

b) độ sai khác ô, theo thiết kế mức tách biệt.

6.2.2. Ký hiệu p thành phần của số liệu được sắp xếp theo thứ tự tăng bởi:

x₁, x₂,……, x_i,……, x_p

...

6.2.3. Tính các giá trị ban đầu của x* và s* như sau:

x* = trung vị của của các x_i (i = 1, 2,…, p) (56)

s*= 1,483 x trung vị của các |x₁ - x*|, (i = 1, 2,…, p) (57)

6.2.4. Cập nhật các giá trị của x* và s* như sau:

Tính

j = 1,5s* (58)

với mỗi giá trị x₁ (i = 1, 2,…, p), tính:

Tính các giá trị mới của x* và s* từ:

...

6.2.5. Các ước lượng ổn định x* và s* có thể được dẫn ra bằng cách tính lặp, tức là bằng cách lặp lại việc tính toán trong 6.2.4 nhiều lần, cho đến khi sự biến đổi của các ước lượng của x* và s* giữa hai lần tính liên tiếp là nhỏ. Đó là một phương pháp đơn giản cho việc lập trình trên máy vi tính.

6.2.6. Một phương pháp khác không cần tính lặp và như vậy rất tiện dụng khi tính toán bằng tay, được đưa ra bởi nhận xét rằng các phương trình (60) và (61) trong 6.2.4 có thể viết dưới dạng:

trong đó

u_L là số các thành phần x_i của số liệu mà x_i < x* - j

u_U là số các thành phần x_i của số liệu mà x_i > x* + j

x’ và s' là trung bình và độ lệch chuẩn của (p - u_L - u_U) các thành phần x_i của số liệu đối với nó |x_i - x*| ≤ j

Các đại lượng này được dùng để tính x* và s* một cách trực tiếp nếu u_L và u_U đã biết. Một cách tiếp cận nữa là hãy thử các khả năng khác nhau theo một thứ tự có hệ thống (tức là hãy thử u_L = 0, u_U = 0; tiếp đó u_L = 0, u_U = 1; u_L = 1, u_U = 0; u_L = 1, u_U = 1 và tiếp tục) cho đến khi tìm được một lời giải hợp lý trong đó số lượng thực sự của các thành phần của số liệu mà chúng sai khác với x* một lượng lớn hơn 1,5s* bằng các giá trị u_L và u_U được dùng để tính x* và s*. Trong thực tế, người phân tích có thể sử dụng các biểu đồ trong hình 4 để nhận ra các giá trị tỏ ra khác với x* một lượng lớn hơn 1,5s* và như vậy sẽ tìm được lời giải bằng cách đánh giá với một số ít các trường hợp.

Một khả năng nữa là sử dụng phương pháp lặp để tìm lời giải xấp xỉ, sau đó giải các phương trình (62) và (63) để tìm lời giải chính xác. Đó là một cách tiếp cận được sử dụng trong ví dụ dưới đây.

...

6.3.1. Thuật toán này được áp dụng cho các độ lệch chuẩn trong phòng thí nghiệm (hoặc các phạm vi trong phòng thí nghiệm) trong một thiết kế bất kỳ. Nó tạo ra một giá trị kết hợp ổn định của các độ lệch chuẩn và các phạm vi mà thuật toán áp dụng với chúng.

6.3.2. Ký hiệu p thành phần của dữ liệu được sắp xếp theo thứ tự tăng dần bởi:

w₁, w₂ …… w_p….. w_p’

(Các đại lượng này có thể là phạm vi hoặc độ lệch chuẩn).

Ký hiệu giá trị kết hợp ổn định bởi w* và bậc tự do ứng với mỗi w_i, bởi v. (Khi w_i là phạm vi, v = 1. Khi w_i là độ lệch chuẩn của n kết quả, v = n -1). Các giá trị x và h cần thiết cho thuật toán này được cho trong bảng 23.

6.3.3. Tính giá trị ban đầu cho w* như sau:

w* = trung vị của các w_i (i = 1,2,…, p) (64)

6.3.4. Cập nhật giá trị của w* như sau

Tính

...

với mỗi giá trị w_i (i = 1,2,…, p), tính:

Tính giá trị mới của w* từ:

6.3.5. Ước lượng ổn định w* có thể được dẫn ra bảng cách tính lặp lại các tính toán trong 6.3.4 nhiều lần cho đến khi sự biến đổi theo ước lượng của w* giữa hai lần tính toán liên tiếp là nhỏ. Đó là một phương pháp đơn giản cho việc lập trình trên máy tính.

6.3.6. Một phương pháp khác không cần tính lặp và như vậy rất tiện dụng khi tính toán bằng tay, tương tự với phương pháp đã nêu trong 6.2.6. Phương trình (67) trong 6.3.4 có thể được viết dưới dạng:

trong đó

S’ là tổng lấy theo các w_i sao cho w_i ≤ Y

...

Điều này có thể giải quyết bằng cách thử u_U = 0, u_U =1, u_U = 2, v.v... cho đến khi có được một lời giải hợp lý, trong đó số lượng thực sự các w_i vượt quá h x w* bằng u_U. Trong thực tế, nhà phân tích sẽ có thể sử dụng biểu đồ trong hình 4 để nhận ra các phạm vi tỏ ra vượt h x w* và tìm được lời giải bảng cách đánh giá với một số ít trường hợp.

Cách tiếp cận dùng trong ví dụ sau là sử dụng phương pháp lặp để tìm một lời giải xấp xỉ, sau đó giải phương trình (68) để tìm lời giải chính xác.

Bảng 23 - Các yếu tố cần thiết cho phân tích ổn định: Thuật toán S

Bậc tự do

Thừa số giới hạn

Thừa số hiệu chỉnh

...

1,645

1,097

1,517

1,054

1,444

0,039

...

1,032

1,359

1,027

1,332

1,024

1,310

...

1,292

1,019

1,277

10,18

1.264

1,017

...

6.4. Các công thức: Phân tích ổn định ở một mức cụ thể của thiết kế đồng mức

6.4.1. Với thiết kế đồng mức, một ước lượng của độ lệch chuẩn lặp lại s_r ở một mức có thể thu được bằng cách áp dụng thuật toán S đối với các phạm vi ở hoặc các độ lệch tiêu chuẩn ô đề tính giá trị ổn định của w từ phương trình (67) trong 6.3.4. Nếu thuật toán S được áp dụng đối với các độ lệch chuẩn ô thì

s_r = w* (69)

Nếu có hai kết quả đo trên một ô và thuật toán S được áp dụng cho các phạm vi ô, thì

(70)

6.4.2. Ước lượng ổn định của độ lệch chuẩn của các trung bình ô s_d cho một mức có thể thu được bằng cách áp dụng thuật toán A đối với các trung bình ô cho mức đó để thu được một giá trị ổn định s* từ phương trình (61) trong 6.2.4 và sau đó sử dụng

s_d = s (71)

6.4.3. Tiếp theo độ lệch chuẩn giữa phòng thí nghiệm s_L có thể được dẫn ra khi sử dụng

...

Nếu biểu thức dưới dấu căn bậc hai là âm thì đặt

s_D = 0 (73)

Tính độ lệch chuẩn tái lập ở mức đó

6.5. Ví dụ 4: Phân tích ổn định ở một mức cụ thể của thiết kế đồng mức

6.5.1. Ví dụ 3 trong TCVN 6910-2:2001 (ISO 5725-2:1994) là một ví dụ của thiết kế đồng mức trong đó số liệu chứa các giá trị tản mạn và bất thường. Mức 5 trong ví dụ đó là một mức đáng lưu ý vì phòng thí nghiệm 1 đưa ra một trung bình ô tỏ ra là gần tản mạn theo phép kiểm nghiệm Grubb và phòng thí nghiệm 6 đưa ra một trung bình ô tỏ ra là gần tản mạn theo phép kiểm nghiệm Cochran. Các số liệu đó được tái tạo lại ở đây trong bảng 24.

6.5.2. Nếu số liệu từ tất cả các phòng thí nghiệm đều được giữ lại, các độ lệch chuẩn lặp lại và tái lập có thể được ước lượng bằng cách sử dụng các công thức trong 7.4 của TCVN 6910-2:2001 (ISO 5725- 2:1994) và cho ta

p = 9

m = 20,511

...

s_d = 1,727

s_L = 1,677

s_R = 1,776

6.5.3. Tuy nhiên, theo TCVN 6910-2 (ISO 5725-2): nhà phân tích số liệu đã sử dụng thông tin từ các mức khác trong thí nghiệm và các nghi ngờ liên quan tới việc nhận dạng các mẫu được thử nghiệm bởi phòng thí nghiệm 6 để lý giải việc loại bỏ các số liệu của cả hai phòng thí nghiệm 1 và 6 khi tính toán và tính được

p = 7

m = 20,412

s_r = 0,393

s_d = 0,573

s_L = 0,501

...

Rõ ràng rằng việc quyết định loại các số liệu từ hai phòng thí nghiệm đó đã có ảnh hưởng thực sự đến các ước lượng của độ Iêch chuẩn lặp lại và tái lập.

6.5.4. Bước đầu tiên trong việc phân tích là thu được một ước lượng ổn định của độ lệch chuẩn lặp lại. Việc tính toán có thể được thực hiện một cách tiện lợi như được chỉ ra trong bảng 25, trong đó các phạm vi ô đã được sắp xếp theo thứ tự tăng. Việc sử dụng thuật toán S và bằng cách lặp lại đưa đến các kết quả được chỉ ra trong bảng 25. Trong ví dụ này bậc tự do của mỗi phạm vi ô là v = 1, như vậy x = 1,097 và h =1,645. Sau bốn bước lặp được chỉ ra trong bảng, giá trị ổn định thu được w* » 0,7 và chỉ có một phạm vi ô (w₉ = 1,98) vượt quá y. Nếu các tính toán được thực hiện trên máy tính, quá trình này có thể được tiếp tục cho đến khi sự biến đổi theo giá trị của w* giữa hai lần lặp liên tiếp là nhỏ.

Lời giải cũng có thể được dẫn ra một cách trực tiếp như sau: việc sử dụng phương trình (68) trong 6.3.6 với

u_U = 1

cho

và đưa đến lời giải (nếu giả thiết u_U = 1 là đúng)

w* = 0,69 % dấu criôsốt (dầu lấy từ cây giẻ gai)

...

Như vậy ước lượng của độ lệch chuẩn lặp lại là:

Giá trị này nằm giữa hai ước lượng được cho trong 6.5.2 và 6.5.3.

6.5.5. Bước tiếp theo trong việc phân tích là xác định ước lượng ổn định của độ lệch chuẩn của các trung bình ô. Việc áp dụng thuật toán A đối với các trung bình ô cho ta các kết quả được chỉ ra trong bảng 26. trong đó các trung bình ô được sắp xếp theo thứ tự tăng. Sau bốn bước lặp được chỉ ra trong bảng đó các giá trị ổn định sẽ là x* = 20,412 và s* = 1,1. và chỉ có hai trung bình ô vượt trội (x₁* = 17,570: x*₉ = 24,140) sai khác với x* một lượng lớn hơn j. Nếu việc tính toán được thực hiện trên máy tính, quá trình đó có thể được tiếp tục cho đến khi sự biến đổi theo x* và s* giữa hai lần lặp liên tiếp là nhỏ.

Nếu việc tính toán được thực hiện bằng tay, người phân tích số liệu cần sử dụng phương pháp trực tiếp đã được mô tả trong 6.2.6 và cần thử với u_L = u_U = 1, bằng cách đó thu được

x' = 20,412 và s' = 0,573 % dầu criôsốt.

Do đó, từ phương trình (62). (63) trong 6.2.6

(s*)² = 6 x (0,573)²/ [8/1,134² - 1,5² (9+9-4)/7]

từ đó thu được:

...

và

x* = x' = 20,412 % dấu criôsốt.

Có thể khẳng định rằng, giá trị đó của s* sẽ cho j = 1,605 (như đã giả định, chỉ có x₁* và x₉* là khác với x* = 20,412 một lượng lớn hơn j), và nếu thay x₁* bởi 18,807 và x₉* bởi 22,017 lại thu được các giá trị mới x* = 20,412 và s* = 0,944 x 1,134 = 1,070. Như vậy lời giải đó tỏ ra là hợp lý.

Từ phương trình (72) trong 6.4.3, thu được ước lượng của độ lệch chuẩn giữa các phòng thí nghiệm là

và ước lượng độ lệch chuẩn tái lập là

Một lần nữa, giá trị này nằm giữa hai ước lượng được cho trong 6.5.2 và 6.5.3.

Bảng 24 - Ví dụ 4: Chuẩn độ nhiệt của dầu criôsốt

...

Số liệu
% dấu criôsốt

Trung bình ô
% dấu

Phạm vi ô
% dầu

24,28

24,00

24,140

0,28

...

20,40

19,91

20,155

0,49

19,30

19,70

19,500

0,40

...

20,30

20,300

0,00

20,53

20,88

20,705

0,35

...

18,56

16,58

17,570

1,98

19,70

20,50

20,100

0,80

...

21,10

20,78

20,940

0,32

20,71

21,66

21,185

0,95

...

(v = 1; x = 1,097; h = 1,645)

Lần lặp

0¹⁾

...

0,86

1,00

1,09

w₁*

0,00

...

0,28

w₃*

0,32

...

0,32

w₄*

0,35

w₅*

0,40

...

0,40

w₆*

0,49

...

0,80

0,66

0,80

w₈*

0,95

0,66

0,86

...

0,95

w₉*

1,98

0,66

0,86

1,00

1,09

w kết hợp

0,83

...

0,56

0,60

0,62

w* mới

0,40 ²⁾

0,52

0,61

0,66

0,68

...

2) 0,40 là phạm vi trung vị [ xem phương trình (64) trong 6.3.3]

Bảng 26 - Ví dụ 4 : Áp dụng thuật toán A đối với các trung bình ô ( % dầu criôsốt)

Lần lặp

0¹⁾

...

1,424

1,478

1,514

1,539

x* - j

18,876

18,909

18,893

...

x* + j

21,724

21,865

21,921

21,950

x₁*

17,570

18,876

...

18,893

18,872

x₂*

19,500

x₃*

...

20,100

x₄*

20,155

...

x₅*

20,300

x₆*

20,705

...

20,705

x₇*

20,940

x₈*

...

21,185

x₉*

24,140

21,724

21,865

21,921

...

Trung bình

20,511

20,387

20,407

20,411

20,412

Đô lệch chuẩn

1,727

0,869

...

0,905

0,916

x* mới

20,300 ²⁾

20,387

20,407

20,411

20,412

s* mới

...

0,985

1,009

1,026

1,039

1) Cột với lần lặp 0 được trích từ bảng 24 sau khi sắp xếp theo thứ tự tăng

2) 20,300 và 0,949 thu được khi sử dụng phương trình (56) và (57) trong 6.2.3

6.6. Các công thức : Phân tích ổn định ở một mức cụ thể của thiết kế mức tách biệt

6.6.1. Với thiết kế mức tách biệt, ước lượng ổn định của độ lệch chuẩn lặp lại s_r ở một mức có thể thu được bằng cách áp dụng thuật toán A đối với độ sai khác ô ở mức đó để suy ra giá trị ổn định s*, từ phương trình (61): và sau đó tính s_r theo công thức

...

Có thể sử dụng các công thức trong 4.5.6 để tính ước lượng của độ lệch chuẩn tái lập ở mức đó.

6.7. Ví dụ 5: Phân tích ổn định ở một mức cụ thể của thiết kế mức tách biệt

6.7.1. Các số liệu của ví dụ 1 trong 4.8 chứa một số giá trị tản mạn và bất thường (xem bảng 8). Cũng như vậy, hình 3 thể hiện một độ chênh âm nhất quán theo các kết quả của phòng thí nghiệm 5. Nếu không phát hiện ra các nguyên nhân của sự bất bình thường đó, người phân tích số liệu sẽ bị đặt vào tình huống khó xử khi cần quyết định xem số liệu nào phải bị loại khi tính toán độ lệch chuẩn lặp lại và tái lập. Các số liệu ở mức 14 (xem bảng 4) được dùng ở đây để minh họa các kết quả thu được bằng việc phân tích ổn định.

6.7.2. Để nhận được ước lượng ổn định của độ lệch chuẩn lặp lại, ta hãy áp dụng thuật toán A đối với các độ sai khác ô (từ bảng 5). Các kết quả thu được chỉ ra trong bảng 27, trong đó các độ sai khác ô đã được sắp xếp theo thứ tự tăng. Sau bốn lần lặp được chỉ ra trong bảng đó, rõ ràng rằng các giá trị ổn định là x* » 8,29 và s* » 0,36 và chỉ có x₉* khác với x* một lượng lớn hơn j.

Việc áp dụng phương pháp đã mô tả trong 6.2.6 với u_L = 0, u_U = 1 sẽ cho

x’= 8,219 và s'= 0,257 % protein.

như vậy các phương trình (62) và (63) có thể được viết :

x* = 8,219 + 1,5 x s*/8

...

(s*)² = 7 x (0,257)²/ [8/1,134² - 1,5²(0 + 9 - 0)/ 8]

do đó

s* = 0,354 % protein.

và khi sử dụng phương trình (75) trong 6.6.1 ta được

Trung bình ổn định của độ sai khác ô là

Với các giá trị đó của x* và s* thì

j = 1,5 x 0,354 = 0,531

...

x* - j = 7,754 và x* - j = 8,816 % protein.

Như đã giả định, để tính x* va s*, rằng chỉ có x₉* nằm ngoài các giới hạn đó. Có thể khẳng định rằng điều đó là đúng và như vậy lời giải thu được là có hiệu lực (chấp nhận được).

6.7.3. Việc áp dụng thuật toán A đối với các trung bình ô (từ bảng 6) cho ta các kết quả trong bảng 28, trong đó các trung bình ô đã được sắp xếp theo thứ tự tăng. Tình huống ở đây tương tự với tình huống được phát hiện ở bảng 26, trong đó x₁* và x₉* sai khác với x* một lượng lớn hơn j, và x* nằm gần trung bình của các giá trị từ x₂* đến x₈* tức là giá trị 85,486. Nếu áp dụng một lần nữa phương pháp của 6.2.6 nhưng với u_L = u_U = 1, giá trị trung bình và độ lệch chuẩn của các giá trị từ x₂* đến x₈* là :

x' = 85,486 và s' = 0,209

Do đó có thể dùng phương trình (63) trong 6.2.6 để suy ra s* từ

(s*)² = 6 x (0,209)²/ [8/1,134² - 1,5²(9 + 9 - 4 )/ 7]

như vậy

s* = 0,390 % protein.

Bây giờ có thể sử dụng phương trình (62) trong 6.2.6 để nhận được

...

Để kiểm tra xem lời giải đó có hiệu lực hay không, hãy tính

j = 1,5 x 0,390 = 0,585

x* - j = 84,901 và x* + j = 86,071 % protein.

Ta thấy rằng, như đã giả thiết, chỉ có x₁* và x₉* nằm ngoài các giới hạn đó.

Để nhận được độ lệch chuẩn tái lập, sử dụng phương trình (76) trong 6.6.2 để có:

s_r = 0,390 % protein.

và sử dụng phương trình (13) trong 4.5.6 để có:

s_R = 0,410 % protein.

Do đó, trong ví dụ này, phương pháp ổn định sẽ cho các ước lượng s_r và s_R nhỏ hơn một chút so với các giá trị thu được khi sử dụng tất cả các số liệu đã báo cáo (cho trong bảng 7).

...

Lần lặp

0¹⁾

0,53

...

0,55

0,54

x* - j

7,85

7,74

7,75

x* + j

...

8,91

8,86

8,84

8,83

x₁*

7,81

7,85

7,81

...

x₂*

7,93

x₃*

8,13

...

8,13

x₄*

8,14

x₅*

...

8,38

x₆*

8,40

...

x₇*

8,44

x₈*

8,52

...

8,52

x₉*

9,31

8,91

8,86

8,84

8,83

Trung bình

...

8,300

8,290

3,288

8,287

Độ lệch chuẩn

0,436

0,326

0,322

0,317

...

x* mới

8,380 ¹⁾

8,300

8,290

8,288

8,287

s* mới

0,356 ¹⁾

0,370

...

0,359

0,357

1) Thu được khi sử dụng phương trình (56), (57) trong 6.2.3.

Bảng 28 - Ví dụ 5: Áp dụng thuật toán A đối với các trung bình ô (% protein)

Lần lặp

0¹⁾

...

0,446

0,492

0,519

0,537

x* - j

85,104

...

84,971

84,950

x* + j

85,996

85,993

86,009

86,024

x₁*

...

85,104

85,009

84,971

84,950

x₂*

85,140

...

x₃*

85,345

x₄*

85,385

...

85,385

x₅*

85,550

x₆*

...

85,575

x₇*

85,660

...

x₈*

85,750

x₉*

86,170

85,996

...

86,009

86,024

Trung bình

85,456

85,501

85,490

85,487

Độ lệch chuẩn

...

0,289

0,305

0,316

0,324

x* mới

85,550 ¹⁾

85,501

85,490

85,487

...

s* mới

0,297 ¹⁾

0,328

0,346

0,358

0,367

1) Thu được khi áp dụng phương trình (56), (57) trong 6.2.3.

6.8. Các công thức : Phân tích ổn định ở một mức cụ thể của một thí nghiệm trên vật liệu không đồng nhất

6.8.1. Với thiết kế dùng vật liệu không đồng nhất, trong trường hợp bình thường khi hai mẫu được chuẩn bị cho một trong p' phòng thí nghiệm ở mỗi mức và hai kết quả thử nghiệm thu được ở mỗi mẫu, các ước lượng ổn định của độ lệch chuẩn lặp lại và tái lập có thể thu được bằng cách áp dụng ba lần các thuật toán A và S như sau

...

SS_r = 2p(w*)² (77)

b) Áp dụng thuật toán S đối với phạm vi giữa các mẫu để dẫn ra giá trị ổn định khác w* từ phương trình (67) và tính:

SS_H=p'(w*)² (78)

c) Áp dụng thuật toán A đối với trung bình ô để suy ra một giá trị ổn định s* từ phương trình (61) trong 6.2.4 và đặt

s_y = s* (79)

Các tính toán đó có thể được sắp đặt một cách tiện lợi dưới dạng bảng, với cột thứ nhất gồm các phạm vi hoặc các trung bình được đưa vào theo thứ tự tăng, như được chỉ ra trong ví dụ dưới đây.

6.8.2. Các công thức nêu ra trong 5.5 có thể được sử dụng để tính các ước lượng của độ lệch chuẩn lặp lại và tái lập và của độ lệch chuẩn s_H đo độ thay đổi giữa các mẫu.

6.9. Ví dụ 6 : Phân tích ổn định ở một mức cụ thể của thí nghiệm trên vật liệu không đồng nhất

6.9.1. Các số liệu ở mức 6 của Ví dụ 2 trong 5.8 không chứa bất kỳ giá trị bất thường hoặc tản mạn nào và vì vậy chúng được sử dụng ở đây để minh họa cho các kết quả thu được bằng phương pháp ổn định trong trường hợp đó.

...

u_U = 4

å'(w₁*)²/p = 137,92 / 22 = 6,2691.

Như vậy phương trình (68) trong 6.3.6 trở thành

(w*)² = 1,097² x 6,2691 + 4(1,097 x 1,645w*)²/22

và do đó

w* = 4,30 %

Có thể khẳng định rằng với giá trị này của w*, y =7,1 và bốn giá trị từ w₁₉* đến w₂₂* đều vượt quá y, như vậy giá trị đó của w* là một lời giải hợp lý.

Việc sử dụng phương trình (77) trong 6.8.1 cho ta

SS_r = 22 x 4,30² = 406,78 (%)²

...

u_v = 2

å'(w₁*)²/p' = 66,665/11 = 6,0605

như vậy phương trình (68) trong 6.3.6 trở thành:

(w*)² = 1,097² x 6,0605 + 2(1,097 x 1.645w*)²/11

và do đó

w* = 4,23 %

Đáng tiếc rằng, giá trị này ứng với y = 1,645 x 4,23 = 6,96, do đó lời giải này không hợp lý vì w₁₀* và w₁₁* không vượt quá giá trị đó. Điều này gợi ý rằng lời giải đòi hỏi phải có u_U = 1 hoặc u_U = 0.

Bằng cách thử cho u_U =1 ta có :

å'(w_i*)²/p' = 112,2275/11 = 10,2025

...

(w*)² = 1,097² x 10,2025 + (1,097 x 1,645w*)²/11

và từ đó

w* = 4,18 %

Bây giờ y = 1,645 x 4,18 = 6,88 và có thể thấy đó là lời giải hợp lý vì w₁₁* vượt quá 6,88.

Việc sử dụng phương trình (78) trong 6.8.1 cho ta

SS_H=11 x 4,18² = 192,20 (%)²

6.9.4. Áp dụng thuật toán A đối với các trung bình ô (từ bảng 16) cho ta các kết quả được chỉ ra trong bảng 31. Việc tính toán sẽ hội tụ sau hai lần lặp và cho ta s* = 5,70 (không có x_i* nào khác với x* một lượng lớn hơn j).

Nếu sử dụng phương trình (79) trong 6.8.1 thu được:

s_y = 5,70 %

...

s_r² = 406,78 / 44 (%) ²

S_R² = 5,70² + (406,78 - 192,20)/ 44 (%)²

s_H² = 192,20 / 22 - 406,78 / 88 (%)².

do đó

s_y = 3,04 %

s_R = 6,11 %

s_H = 2,03%

Do đó trong ví dụ này, phương pháp ổn định cho các ước lượng s_r, s_R, s_H lớn hơn các giá trị thu được khi sử dụng tất cả các số liệu đã được báo cáo (cho trong 5.8.3 và bảng 17).

Bảng 29 - Ví dụ 6: Áp dụng thuật toán S đối với phạm vi giữa các kết quả thử nghiệm

...

Lần lặp

3,9

...

5,9

6,4

w₁*

0,1

w₂*

...

0,6

w₃*

1,1

...

w₄*

1,1

w₅*

1,2

...

1,2

w₆*

1,3

w₇*

...

1,4

w₈*

1,6

...

w₉*

1,8

w₁₀*

2,1

...

2,1

w₁₁*

2,2

w₁₂*

...

2,5

w₁₃*

2,6

...

w₁₄*

3,9

w₁₅*

4,0

3,9

...

4,0

w₁₆*

4,4

3,9

4,4

w₁₇*

...

3,9

4,6

w₁₈*

5,5

3,9

5,1

5,5

...

w₁₉*

7,4

3,9

5,1

5,9

6,4

w₂₀*

7,6

3,9

...

5,9

6,4

w₂₁*

8,1

3,9

5,1

5,9

6,4

w₂₂*

...

3,9

5,1

5,9

6,4

w mới

4,17

2,80

3,29

3,55

...

w* mới

2,35 ¹⁾

3,07

3,61

3,89

4,06

1) Thu được khi sử dụng phương trình (64) trong 6.3.3

Bảng 30 - Ví dụ 6: Áp dụng thuật toán S đối với các phạm vi giữa các mẫu (%)

(v = 1; x = 1,097; h = 1,645)

...

4,19

5,43

...

6,45

w₁*

1,00

w₂*

1,70

...

1,70

w₃*

2,05

...

2,25

w₅*

2,55

...

2,55

w₆*

2,55

w₇*

3,15

...

3,15

w₈*

3,35

...

4,40

4,19

4,40

w₁₀*

6,75

4,19

5,43

...

6,45

w₁₁*

6,95

4,19

5,43

6,10

6,45

w mới

3,82

...

3,38

3,58

3,69

w* mới

2,55 ¹⁾

3,30

3,71

3,92

4,05

...

Bảng 31 - Ví dụ 4 : Áp dụng thuật toán A đối với các trung bình ô ( %)

Lần lặp

...

8,550

x* - j

8,245

10,450

...

28,255

27,550

x₁*

13,425

...

x₂*

13,425

x₃*

13,750

...

13,750

x₄*

14,475

...

17,075

x₆*

18,250

...

x₇*

21,000

x₈*

21,225

...

21,225

x₉*

23,675

...

26,275

x₁₁*

26,425

...

Trung bình

19,00

Độ lệch chuẩn

5,03

...

5,03

x* mới

18,25 ¹⁾

19,00

...

6,67 ¹⁾

5,70

1) Thu được khi sử dụng phương trình (56) và (57) trong 6.2.3

PHỤ LỤC A

(quy định)

...

Phần bị chắn trong mối quan hệ

s = a + bm

Yếu tố dùng để tính độ không đảm bảo của ước lượng

Độ dốc trong mối quan hệ

s = a + bm

...

B_o

Thành phần của B biểu thị tất cả các yếu tố không thay đổi trong điều kiện chụm trung gian

B₍₁₎, B₍₂₎....

Các thành phần của B biểu thị những yếu tố thay đổi trong điều kiện chụm trung gian

Phần bị chắn trong mối quan hệ

lg s = c + d lg m

C, C', C''

Các thống kê kiểm nghiệm

...

Các giá trị tới hạn đối với những phép kiểm nghiệm thống kê

CD_P

Độ sai khác tới hạn với xác suất P

CR_P

Phạm vi tới hạn với xác suất P

Độ dốc trong mối liên hệ

lg s = c + d lg m

...

Yếu tố phạm vi tới hạn

F_p(v₁, v₂)

Phân vị mức p của phân bố F với các bậc tự do v₁ và v₂

Thống kê kiểm nghiệm Grubb

Thống kê kiểm nghiệm nhất quán giữa các phòng thí nghiệm của Mandel

...

LCL

Giới hạn kiểm soát dưới (hoặc giới hạn hành động hoặc giới hạn cảnh báo)

Trung bình chung của đặc tính thử; mức

Yếu tố được xem xét trong điều kiện chụm trung gian

Số phép lặp

...

Số phòng thí nghiệm tham gia thí nghiệm liên phòng

Xác suất

Số lượng các mức của đặc tính thử nghiệm trong thí nghiệm liên phòng

Giới hạn lặp lại

...

Mẫu chuẩn

Ước lượng của độ lệch chuẩn

Độ lệch chuẩn dự đoán

Tổng thể hoặc tổng của biểu thức nào đó

...

UCL

Giới hạn kiểm soát trên (hoặc giới hạn hành động hoặc giới hạn cảnh báo)

Yếu tố trọng số sử dụng trong tính toán hồi quy trọng số

Độ rộng của tập hợp các kết quả thử nghiệm

Dữ liệu sử dụng cho thử nghiệm Grubb

...

Trung bình số học của kết quả thử nghiệm

Trung bình chung của kết quả thử nghiệm

Mức ý nghĩa

Xác suất sai lầm loại II

...

Độ chệch phòng thí nghiệm

Ước lượng của D

Độ chệch của phương pháp đo

Ước lượng của d

...

Giá trị thực hoặc giá trị quy chiếu được chấp nhận của đặc tính thử nghiệm

Số bậc tự do

Tỷ số phát hiện được giữa độ lệch chuẩn lặp lại của phuơng pháp B và phương pháp A

Giá trị thực của độ lệch chuẩn

...

Tỷ số phát hiện được giữa căn bậc hai của bình phương trung bình giữa các phòng thí nghiệm của phương pháp B và phương pháp A

c²_r(n)

Phân vị mức r của phân bố c²với bậc tự do n

Các ký hiệu được sử dụng như chỉ số

Sự khác nhau về hiệu chuẩn

Sự khác nhau về thiết bị

...

Chỉ số của một phòng thí nghiệm cụ thể

I( )

Chỉ số của thước đo trung gian của độ chụm, trong dấu ngoặc chỉ loại tình huống trung gian

Chỉ số của một mức cụ thể (TCVN 6910-2 (ISO 5725-2))

Chỉ số của một nhóm phép thử nghiệm hoặc một yếu tố (TCVN 6910-3)

Chỉ số của một kết quả thử nghiệm cụ thể trong phòng thí nghiệm / ở mức j

...

Liên phòng thí nghiệm (liên phòng)

Chỉ số của độ chệch có thể biết được

Mẫu thử liên phòng

Sự khác nhau về người thao tác

Xác suất

...

Độ lặp lại

Độ tái lập

Sự khác nhau về thời gian

Phòng thí nghiệm thành viên

1,2,3 ...

Đối với các kết quả thử nghiệm, đánh số theo thứ tự thu nhận chúng

...

Đối với các kết quả thử nghiệm, đánh số theo thứ tự tăng độ lớn

Các ký hiệu và chữ viết tắt dùng trong TCVN 6910-5 : 2002

Độ sai khác trong ô của thiết kế mức tách biệt

Số mẫu thử nghiệm trong một phòng thí nghiệm ở một mức.

Thành phần trong kết quả đo biểu thị sai số ngẫu nhiên của mẫu.

...

Số các phòng thí nghiệm tham gia trong thí nghiệm liên phòng

Tổng các bình phương

u_L

Số các giá trị của dữ liệu nhỏ hơn giới hạn dưới trong phân tích ổn định

u_U

Số các giá trị của dữ liệu lớn hơn giới hạn trên trong phân tích ổn định

...

Tỷ số của các độ lệch chuẩn.

Giới hạn được dùng trong phân tích ổn định (thuật toán A)

Thừa số giới hạn được dùng trong phân tích ổn định (thuật toán S)

Giới hạn được dùng trong phân tích ổn định (thuật toán S)

...

Các ký hiệu phụ thêm dùng làm chỉ số trong TCVN 6910 - 5:2002

a,b

Chỉ số phân biệt các mẫu trong thí nghiệm mức tách biệt

Chỉ số phân biệt một mẫu trong phòng thí nghiệm i ở mức j

Giữa các mẫu

Các ký hiệu phụ thêm được sử dụng như chỉ số trong TCVN 6910-5 : 2002

...

PHỤ LỤC B

(tham khảo)

DẪN RA CÁC THỪA SỐ SỬ DỤNG TRONG THUẬT TOÁN A VÀ S

B.1. Mở đầu

Việc sử dụng phương pháp ổn định để phân tích số liệu của các thí nghiệm độ chụm đã được đề suất bởi Hội đồng về các phương pháp phân tích của Hội Hóa học Hoàng gia Anh (UK) xem [6]. Thuật toán A trong tiêu chuẩn này được dẫn ra từ bài báo của họ và thừa số 1,134 dùng để tính s* trong thuật toán A cũng được lấy ra từ bài báo của họ (đó là giá trị 1/ cho trường hợp c = 1,5 theo ký hiệu của bài báo đó )

“Thuật toán S” là thủ tục tương tự với thủ tục được đưa ra trong [6] cho trường hợp đặc biệt, trong đó mỗi phòng thí nghiệm báo cáo n = 2 lần đo ở một mức. Nó đưa ra một phương pháp tiện lợi cho việc áp dụng việc phân tích ổn định đối với các thí nghiệm độ chụm với nhiều hơn hai thừa số (như thiết kế với các vật liệu không đồng nhất trong điều 5 của tiêu chuẩn này và thiết kế xếp chồng so le của TCVN 6910 -3.). Việc dẫn ra các thừa số sử dụng trong thuật toán S được thực hiện dưới đây.

B.2. Các ký hiệu sử dụng trong phụ lục này

...

độ lệch chuẩn ước lượng của s

bậc tự do của s

n + 2

Thừa số hiệu chính của thuật toán S

...

c_n²

Biến khi - bình phương với n bậc tự do

B.3. Dẫn ra các thừa số giới hạn h và thừa số hiệu chính x

Thừa số hiệu chính x được định nghĩa như là giá trị cần thiết để hiệu chính s* sao cho (s*)² cho ta một ước lượng không chệch của s² tức là:

Yêu cầu đó có thể được viết dưới dạng

trong đó các đại lượng trong ngoặc { } có liên quan chặt chẽ với đại lượng n (s/s)², biến khi - bình phương, c_n²

...

như vậy

vì s £ hs tương đương với n(s /s)² £ nh²

Số hạng thứ hai trong về phải của B4 là:

Đối với thuật toán S, thừa số giới hạn h được chọn sao cho hs là điểm ở trên mức 10 % của phân số s, tức là điểm sao cho:

Các bảng Biometrica (sinh trắc) của phân bố khi - bình phương cho ta các giá trị của h ở trong bảng 23 tiêu chuẩn này. Từ các phương trình (B5), (86) suy ra rằng thành phần thứ hai của vế phải của (B4) là 0,1 nh². Chú ý rằng h phụ thuộc vào bậc tự do của s

...

Theo tính chất của hàm gama với w = n + 2 ta có

như vậy số hạng đầu đó trở thành:

Như vậy, với số bậc tự do n đã cho, h có thể tính như trên và sau đó cũng tính được z khi sử dụng bảng phân bố khi - bình phương cho trong các bảng Biometrika. Như vậy cả hai ở vế phải của (B4) có thể xác định được.

Nếu thay các phương trình (B2) đến (B7) vào (B4) ta được.

hoặc

...

Hệ thức này dùng để thu được các giá trị của thừa số hiệu chỉnh x cho trong bảng 23 của tiêu chuẩn này.

PHỤ LỤC C

(tham khảo)

DẪN RA CÁC PHƯƠNG TRÌNH DÙNG ĐỂ PHÂN TÍCH ỔN ĐỊNH

Các phương trình (62) và (63) dùng để tính các giá trị ổn định của trung bình và độ lệch chuẩn bằng phương pháp đã mô tả trong 6.2.6 có thể được suy ra từ phương trình (60) và (61) trong 6.2.4 (Thuật toán A) như sau:

Với các ký hiệu như trong 6.2.4 và 6.2.6

và

...

trong đó S' biểu thị tổng lấy trên (p - u_L - u_U) các số hạng x₁ của số liệu sao cho |x₁ - x₁*| £ j

Do đó có thể viết phương trình (C1) như sau:

như vậy

và từ đó:

đó chính là phương trình (62) trong 6.2.6.

Để suy ra phương trình (63) từ phương trình (61), chú ý rằng có thể khai triển tổng trong phương trình (61) như sau.

...

Nếu thay x* trong số hạng tổng ở đây bởi vế phải của (C4), sau một số tính toán đại số ta thu được:

Nếu sử dụng định nghĩa của s' trong (C3), tổng trên có thể viết thành:

Sau đó nếu thay phương trình (C7) vào phương trình (61) ta thu được phương trình (63).

PHỤ LỤC D

(tham khảo)

Tài liệu tham khảo

...

[2] BS 812-121:1989, Testing aggregates - Part 121: Methods for determination of soundness. Britisl Standards Institution.

[3] BS 3144:1968, Methods of sampling and physical testing of leather. British Standards Institution.

[4] Scheffé, H. The analysis of variance. Wiley, New York, 1959.

[5] Sweeney, An inter-laboratory study of the determination of protein by combustion in feeds. Journal of the Association of Official Analytical Chemists, 72, 1989, pp. 770-774.

[6] Analytical Methods Committee. Robust statistics - How not to reject outliers. Part 1: Basic concepts Part 2: Inter-laboratory trials. The Analyst, 114, 1989, pp. 1653-1697 (part 1), pp. 1699-1702 (part 2) Royal Society of Chemistry, London.

[7] Youden, W.J. The Youden plot. Industrial Quality Control, 15, 1959, pp. 133-137.

[8] Mandel, J and Lashof, T.W. Interpretation and Generalization of Youden's Two-Sample Diagram. Journal of Quanlity Technology, 6, 1974, pp. 22-36.

MỤC LỤC

...

1. Phạm vi áp dụng

2. Tiêu chuẩn viện dẫn

3. Định nghĩa

4. Thiết kế mức tách biệt

4.1. Các ứng dụng của thiết kế mức tách biệt

4.2. Bố trí thiết kế mức tách biệt

4.3. Tổ chức thí nghiệm mức tách biệt

4.4. Mô hình thống kê

4.5. Phân tích thống kê các số liệu từ thí nghiệm mức tách biệt

...

4.7. Báo cáo kết quả thí nghiệm mức tách biệt

4.8. Ví dụ 1: Thí nghiệm mức tách biệt - Xác định protein

5. Thiết kế cho các vật liệu không đồng nhất

5.1. Ứng dụng các thiết kế cho các vật liệu không đồng nhất

5.2. Bố trí thiết kế các vật liệu không đồng nhất

5.3. Tổ chức thí nghiệm với vật liệu không đồng nhất

5.4. Mô hình thống kê cho một thí nghiệm với vật liệu không đồng nhất

5.5. Phân tích thống kê các số liệu từ một thí nghiệm với vật liệu không đồng nhất

5.6. Khảo sát sự phù hợp và sự bất thường của các số liệu

...

5.8. Ví dụ 2: Thí nghiệm trên vật liệu không đồng nhất

5.9. Các công thức tổng quát để tính toán với thiết kế cho vật liệu không đồng nhất

5.10. Ví dụ 3: Áp dụng các công thức tổng quát

6. Phương pháp ổn định để phân tích số liệu

6.1. Áp dụng các phương pháp ổn định để phân tích số liệu

6.2. Phân tích ổn định: Thuật toán A

6.3. Phân tích ổn định: Thuật toán S

6.4. Các công thức: Phân tích ổn định ở một mức cụ thể của thiết kế đồng mức

6.5. Ví dụ 4: Phân tích ổn định ở một mức cụ thể của thiết kế đồng mức

...

6.7. Ví dụ 5: Phân tích ổn định ở một mức cụ thể của thiết kế mức tách biệt

6.8. Các công thức: Phân tích ổn định ở một mức cụ thể của thí nghiệm trên vật liệu không đồng nhất

6.9. Ví dụ 6: Phân tích ổn định ở một mức cụ thể của thí nghiệm trên vật liệu không đồng nhất

Phụ lục A (quy định) Các kí hiệu và chữ viết tắt dùng trong TCVN 6910

Phụ lục B (tham khảo) Dẫn ra các thừa số sử dụng trong thuật toán A và S

Phụ lục C (tham khảo) Dẫn ra các phương trình dùng để phân tích ổn định

Phụ lục D (tham khảo) Tài liệu tham khảo