Tiêu chuẩn Việt Nam TCVN9603 2013 Công ty luật Minh Khuê www luatminhkhue vn TIÊU CHUẨN QUỐC GIA TCVN 9603 2013 ISO 5479 1997 GIẢI THÍCH CÁC DỮ LIỆU THỐNG KÊ KIỂM NGHIỆM SAI LỆCH SO VỚI PHÂN BỐ CHUẨN[.]
Công ty luật Minh Khuê www.luatminhkhue.vn TIÊU CHUẨN QUỐC GIA TCVN 9603:2013 ISO 5479:1997 GIẢI THÍCH CÁC DỮ LIỆU THỐNG KÊ - KIỂM NGHIỆM SAI LỆCH SO VỚI PHÂN BỐ CHUẨN Statistical interpretation of data - Tests for departure from the normal distribution Lời nói đầu TCVN 9603:2013 hồn tồn tương đương với ISO 5479:1997; TCVN 9603:2013 Ban kỹ thuật tiêu chuẩn quốc gia TCVN/TC 69 Ứng dụng phương pháp thống kê 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ệ công bố Lời giới thiệu Nhiều phương pháp thống kê khuyến nghị tiêu chuẩn này, phương pháp mô tả ISO 2854 [1], dựa giả định (các) biến ngẫu nhiên áp dụng với phương pháp độc lập có phân bố chuẩn với hai tham số chưa biết Do nảy sinh câu hỏi sau Phân bố thể mẫu đủ gần với phân bố chuẩn để sử dụng cách tin cậy phương pháp tiêu chuẩn không? Không có câu trả lời đơn giản có khơng cho câu hỏi có hiệu lực trường hợp Vì lý này, nhiều “kiểm nghiệm tính chuẩn” xây dựng, phép kiểm nghiệm nhiều nhạy với đặc trưng cụ thể phân bố xem xét; ví dụ độ bất đối xứng hay độ nhọn Nói chung, phép kiểm nghiệm sử dụng thiết kế để tương ứng với rủi ro tiên nghiệm xác định trước giả thuyết tính chuẩn bị bác bỏ (sai lầm loại một) Mặt khác, xác định xác suất giả thuyết khơng bị bác bỏ khơng (sai lầm loại hai) đối giả thuyết (nghĩa ngược với giả thuyết tính chuẩn) xác định xác Điều nhìn chung khơng thể và, nữa, địi hỏi nỗ lực tính toán Đối với phép kiểm nghiệm riêng rẽ, rủi ro đặc biệt lớn cỡ mẫu nhỏ GIẢI THÍCH CÁC DỮ LIỆU THỐNG KÊ - KIỂM NGHIỆM SAI LỆCH SO VỚI PHÂN BỐ CHUẨN Statistical interpretation of data - Tests for departure from the normal distribution Phạm vi áp dụng 1.1 Tiêu chuẩn đưa hướng dẫn phương pháp phép kiểm nghiệm để sử dụng việc xác định có nên bác bỏ giả thuyết phân bố chuẩn hay không, giả định quan trắc độc lập 1.2 Bất có nghi ngờ việc quan trắc có phân bố chuẩn hay không, việc sử dụng phép kiểm nghiệm sai lệch so với phân bố chuẩn hữu ích hay chí cần thiết Tuy nhiên, trường hợp phương pháp ổn định (nghĩa kết thay đổi phân bố xác suất thực tế quan trắc phân bố chuẩn), phép kiểm nghiệm sai lệch so với phân bố chuẩn khơng hữu ích Đây trường hợp, ví dụ, trung bình mẫu ngẫu nhiên đơn quan trắc kiểm tra dựa giá trị lý thuyết cho trước cách sử dụng phép kiểm nghiệm t 1.3 Không thiết phải sử dụng phép kiểm nghiệm đề cập đến phương pháp thống kê dựa giả thuyết tính chuẩn Có khả khơng nghi ngờ phân bố chuẩn quan trắc cho dù có lý lý thuyết (ví dụ vật lý) khẳng định giả thuyết giả thuyết coi chấp nhận theo thơng tin trước 1.4 Các phép kiểm nghiệm sai lệch so với phân bố chuẩn lựa chọn tiêu chuẩn chủ yếu dùng cho liệu đầy đủ, khơng phải liệu phân nhóm Chúng khơng thích hợp với liệu bị theo dõi 1.5 Các phép kiểm nghiệm sai lệch so với phân bố chuẩn lựa chọn tiêu chuẩn áp dụng cho giá trị quan trắc hàm chúng, logarit bậc hai 1.6 Phép kiểm nghiệm sai lệch so với phân bố chuẩn hiệu cỡ mẫu nhỏ tám Vì vậy, tiêu chuẩn giới hạn cỡ mẫu từ tám trở lên LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê www.luatminhkhue.vn Tài liệu viện dẫn Các tài liệu viện dẫn tiêu chuẩn cần thiết cho việc áp dụng tiêu chuẩn Đối với tài liệu có ghi năm cơng bố áp dụng nêu Đối với tài liệu khơng ghi năm cơng bố áp dụng phiên nhất, bao gồm sửa đổi TCVN 8244-1 (ISO 3534-1), Thống kê - Từ vựng ký hiệu - Phần 1: Thuật ngữ chung xác suất thống kê Thuật ngữ, định nghĩa ký hiệu 3.1 Thuật ngữ định nghĩa Tiêu chuẩn áp dụng thuật ngữ định nghĩa TCVN 8244-1 (ISO 3534-1) 3.2 Ký hiệu Tiêu chuẩn sử dụng ký hiệu ak hệ số phép kiểm nghiệm Shapiro-Wilk A đại lượng phụ trợ cho phép kiểm nghiệm Epps-Pulley b2 độ nhọn thực nghiệm độ bất đối xứng thực nghiệm B đại lượng phụ trợ cho phép kiểm nghiệm Epps-Pulley E kỳ vọng Gj đại lượng phụ trợ cho phép kiểm nghiệm kết hợp sử dụng nhiều mẫu độc lập h số lượng mẫu liên tiếp H0 giả thuyết không H1 đối giả thuyết k số giá trị quan trắc x mẫu, xếp theo thứ tự không giảm mj mômen trung tâm bậc j mẫu n cỡ mẫu p xác suất kèm với phân vị p phân bố P xác suất Pk xác suất kèm với phân vị X(k) S đại lượng phụ trợ cho phép kiểm nghiệm Shapiro-Wilk T thống kê kiểm nghiệm TEP thống kê kiểm nghiệm phép kiểm nghiệm Epps-Pulley up p-phân vị phân bố chuẩn chuẩn hóa vj đại lượng phụ trợ cho phép kiểm nghiệm kết hợp sử dụng nhiều mẫu độc lập W thống kê kiểm nghiệm phép kiểm nghiệm Shapiro-Wilk Wj đại lượng phụ trợ cho phép kiểm nghiệm kết hợp sử dụng nhiều mẫu độc lập x giá trị X X biến ngẫu nhiên x(j) giá trị thứ j mẫu, xếp theo thứ tự không giảm x(k) giá trị thứ k mẫu, xếp theo thứ tự khơng giảm trung bình số học α mức ý nghĩa β xác suất sai lầm loại hai β2 độ nhọn tổng thể LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê www.luatminhkhue.vn β2-3 độ tù tổng thể độ bất đối xứng tổng thể γ đại lượng phụ trợ cho phép kiểm nghiệm kết hợp sử dụng nhiều mẫu độc lập γ (n) hệ số phép kiểm nghiệm kết hợp sử dụng nhiều mẫu độc lập δ đại lượng phụ trợ cho phép kiểm nghiệm kết hợp sử dụng nhiều mẫu độc lập δ(n) hệ số phép kiểm nghiệm kết hợp sử dụng nhiều mẫu độc lập ε đại lượng phụ trợ cho phép kiểm nghiệm kết hợp sử dụng nhiều mẫu độc lập ε(n) hệ số phép kiểm nghiệm kết hợp sử dụng nhiều mẫu độc lập µ kỳ vọng µ2 phương sai tổng thể µ3 mơmen trung tâm bậc ba mẫu µ4 mômen trung tâm bậc bốn mẫu σ độ lệch chuẩn tổng thể ( ) Khái quát 4.1 Có nhiều loại phép kiểm nghiệm sai lệch so với tính chuẩn Trong tiêu chuẩn xem xét phương pháp đồ thị, kiểm nghiệm mômen, kiểm nghiệm hồi quy kiểm nghiệm hàm đặc trưng Phép kiểm nghiệm khi-bình phương thích hợp với liệu nhóm việc phân nhóm dẫn đến thơng tin nên chúng không xem xét tiêu chuẩn 4.2 Nếu khơng có sẵn thơng tin bổ sung mẫu khuyến nghị trước tiên vẽ đồ thị xác suất chuẩn; nghĩa vẽ hàm phân bố tích lũy giá trị quan trắc biểu đồ xác suất chuẩn gồm hệ trục tọa độ hàm phân bố tích lũy phân bố chuẩn thể đường thẳng Phương pháp này, mô tả điều 5, cho phép “thấy” phân bố quan trắc có gần với phân bố chuẩn hay khơng Với thơng tin bổ sung định tiến hành phép kiểm nghiệm định hướng hay tiến hành kiểm nghiệm hồi quy kiểm nghiệm hàm đặc trưng, hay khơng thực kiểm nghiệm Ngồi ra, cách trình bày đồ thị khơng coi kiểm nghiệm chặt chẽ thông tin tổng hợp mà đưa bổ trợ thiết yếu cho phép kiểm nghiệm độ chệch so với phân bố chuẩn Trong trường hợp bác bỏ giả thuyết khơng cách thường có khả ước đốn loại đối giả thuyết áp dụng 4.3 Phép kiểm nghiệm sai lệch so với phân bố chuẩn kiểm nghiệm giả thuyết không mẫu gồm n quan trắc độc lập từ phân bố chuẩn Phép kiểm nghiệm gồm việc tính tốn hàm T quan trắc gọi thống kê kiểm nghiệm Giả thuyết không phân bố chuẩn sau khơng bị bác bỏ bác bỏ tùy thuộc vào việc giá trị T có nằm phạm vi tập hợp giá trị gần giá trị dự kiến tương ứng với phân bố chuẩn hay không 4.4 Miền tới hạn phép kiểm nghiệm tập hợp giá trị T dẫn đến việc bác bỏ giả thuyết không Mức ý nghĩa kiểm nghiệm xác suất P thu giá trị T nằm miền tới hạn giả thuyết không Mức cho xác suất bác bỏ sai giả thuyết không (sai lầm loại một) Ranh giới miền tới hạn (hoặc trường hợp kiểm nghiệm hai phía, ranh giới miền tới hạn) (các) giá trị tới hạn thống kê kiểm nghiệm 4.5 Hiệu lực phép kiểm nghiệm xác suất bác bỏ giả thuyết không khơng Hiệu lực cao ứng với xác suất thấp việc không bác bỏ giả thuyết không cách sai lầm (sai lầm loại hai) Cần nhấn mạnh hiệu lực phép kiểm nghiệm (nghĩa trường hợp định, xác suất giả thuyết không phân bố chuẩn bị bác bỏ điều sai) tăng số quan trắc tăng Ví dụ, độ lệch khỏi phân bố chuẩn trở nên rõ ràng sử dụng phép kiểm nghiệm sai lệch so với phân bố chuẩn mẫu lớn khơng phát phép kiểm nghiệm có quan trắc 4.6 Có khác biệt hai loại kiểm nghiệm sai lệch so với phân bố chuẩn Khi dạng sai lệch so với phân bố chuẩn quy định đối giả thuyết kiểm nghiệm có hướng Tuy nhiên, dạng sai lệch so với phân bố chuẩn không quy định đối giả thuyết kiểm nghiệm kiểm nghiệm vô hướng LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê www.luatminhkhue.vn Trong kiểm nghiệm có hướng, miền tới hạn xác định cho hiệu lực phép kiểm nghiệm đạt giá trị cực đại Trong phép kiểm nghiệm vô hướng, cần chia miền tới hạn cho miền tới hạn bao gồm giá trị thống kê kiểm nghiệm nằm cách xa giá trị dự kiến Nếu giả định thể loại sai lệch so với phân bố chuẩn, nghĩa phân bố nghiên cứu có độ bất đối xứng hay độ nhọn khác so với phân bố chuẩn, cần áp dụng kiểm nghiệm có hướng hiệu lực lớn hiệu lực kiểm nghiệm vô hướng 4.7 Chú ý kiểm nghiệm có hướng thiết phía Trong trường hợp bất đối xứng, ví dụ, tập trung vào bất đối xứng dương bất đối xứng âm Tuy nhiên, xem xét nhiều lựa chọn kiểm nghiệm đa hướng Đây trường hợp đặc biệt độ bất đối xứng khác không độ nhọn khác so với phân bố chuẩn xem xét 4.8 Các Bảng đến Bảng 14 Hình cho phép thực kiểm nghiệm cho hầu hết mức α thông thường; nghĩa α = 0,05 α = 0,01 Mức ý nghĩa phải quy định trước thực kiểm nghiệm Chú ý kiểm nghiệm dẫn đến bác bỏ giả thuyết khơng mức 0,05 khơng bác bỏ giả thuyết mức 0,01 4.9 Trong q trình tính tốn thống kê kiểm nghiệm, cần sử dụng sáu chữ số có nghĩa Các tổng phụ, kết trung gian đại lượng phụ trợ không làm trịn đến sáu chữ số có nghĩa Phương pháp đồ thị 5.1 Hàm phân bố tích lũy giá trị quan trắc vẽ giấy đồ thị xác suất chuẩn Trên đồ thị này, trục (trong tiêu chuẩn trục tung) có thang phi tuyến tính theo vùng nằm hàm phân bố chuẩn chuẩn hóa ghi giá trị tương ứng tần suất tích lũy Trục cịn lại có thang tuyến tính cho giá trị X theo thứ tự Hàm phân bố tích lũy biến X gần đường thẳng Đơi khi, hai trục thay đổi cho Ngoài ra, có chuyển đổi chuẩn hóa biến X, thang tuyến tính thay thang logarit, bậc hai, nghịch đảo thang đo khác Hình đưa ví dụ giấy đồ thị xác suất chuẩn Trên trục tung, giá trị tần suất tích lũy cho theo phần trăm, trục hồnh có thang đo tuyến tính tùy ý Giấy đồ thị xác suất chuẩn để trống cho Phụ lục A Nếu đồ thị bảng đưa tập hợp điểm xuất rải rác quanh đường thẳng điều cung cấp hỗ trợ sơ cho giả định mẫu xem xét cách hợp lý xuất phát từ phân bố chuẩn Tuy nhiên, có độ lệch hệ thống so với đường thẳng đồ thị thường gợi ý loại phân bố cần đưa vào xem xét Tầm quan trọng cách tiếp cận cung cấp thơng tin rõ ràng loại sai lệch so với phân bố chuẩn Nếu đồ thị liệu đến từ phân bố định dạng (ví dụ đồ thị hàm phân bố tích lũy thể Hình 6) việc chuyển đổi liệu dẫn đến phân bố chuẩn Nếu đồ thị liệu không đến từ phân bố đơn giản mà từ pha trộn hai nhiều tổng thể (ví dụ đồ thị hàm phân bố tích lũy thể Hình 7) khuyến nghị cần nhận biết tổng thể tiếp tục phân tích tổng thể cách riêng rẽ Cần lưu ý đồ thị kiểm nghiệm sai lệch so với phân bố chuẩn cách chặt chẽ Trong trường hợp mẫu nhỏ, đường cong rõ rệt xuất phân bố chuẩn, mẫu lớn, đường cong khơng rõ nét phân bố phân bố chuẩn 5.2 Quy trình vẽ đồ thị bao gồm xếp giá trị quan trắc (x(1), x(2), …,x(n)) theo thứ tự khơng giảm sau vẽ đồ thị Pk = (k - 3/8)/(n + 1/4) theo x(k) giấy đồ thị xác suất chuẩn CHÚ THÍCH 1: Các thay cho công thức (1) thường sử dụng Pk = (k - 1/2)/n Pk = k/(n + 1) LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 (1) Công ty luật Minh Khuê www.luatminhkhue.vn Đây phép gần cho hàm phân bố chuẩn thống kê thứ tự dự kiến, F [E(X(k))] việc sử dụng chúng không khuyến nghị 5.3 Ví dụ cách sử dụng giấy đồ thị xác suất chuẩn cho Hình Hình - Giấy đồ thị xác suất chuẩn giải LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Cơng ty luật Minh Kh www.luatminhkhue.vn Hình - Bảng đồ thị xác suất chuẩn giải LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê www.luatminhkhue.vn Bảng thể giá trị x(k) theo thứ tự không giảm kết loạt 15 phép thử chịu uốn luân phiên độc lập Bảng - Kết quả, x(k) loạt 15 phép thử chịu uốn luân phiên giá trị tương ứng lg (10x(k)) k x(k) lg(10x(k)) 0,041 0,200 0,301 0,107 0,330 0,519 0,172 0,445 0,648 0,238 0,490 0,690 0,303 0,780 0,892 0,369 0,920 0,964 0,434 0,950 0,978 0,500 0,970 0,987 0,566 1,040 1,017 10 0,631 1,710 0,233 11 0,697 2,220 1,346 12 0,762 2,275 1,357 13 0,828 3,650 1,562 14 0,893 7,000 1,845 15 0,959 8,800 1,944 CHÚ THÍCH 2: Trong Bảng ví dụ đây, đơn vị quan trắc bỏ qua chúng khơng thích hợp cho phép thử tiêu chuẩn Bằng cách kết hợp xác suất Pk = (k - 3/8)/(n + 1/4) với giá trị x(k) nhỏ thứ k, thu loạt điểm thể Hình 2a) Có thể thấy từ đồ thị điểm không tạo thành đường thẳng Tuy nhiên, x(k) thay lg(10x(k)) đồ thị [Hình 2b)] dẫn đến loạt điểm lúc nằm gần với đường thẳng Do đó, giả thuyết phân bố chuẩn logarit quan trắc thích hợp 5.4 Cần ý giá trị quan trắc cực trị có phương sai lớn giá trị Do thang đo tần suất tích lũy mở rộng phía cực trị, nên số giá trị nằm hai đầu phân bố tích lũy tách biệt khỏi đường thẳng xác định giá trị không coi thị sai lệch so với phân bố chuẩn Cỡ mẫu lớn kết luận rút từ hình dạng đồ thị đáng tin cậy Nếu đồ thị hàm phân bố tích lũy giá trị quan trắc giá trị lớn có xu hướng nằm phía đường thẳng xác định giá trị khác việc chuyển đổi y = log x hay thường dẫn đến đồ thị phù hợp với đường thẳng [xem Hình 2b) Hình 5] Phần Hình đến cho thấy hàm phân bố tích lũy so sánh với hàm mật độ tương ứng thể phần hình Nếu đồ thị hàm phân bố tích lũy giá trị quan trắc thể Hình phân bố LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê www.luatminhkhue.vn tần suất tương ứng độ nhọn bè (platykurtic) độ nhọn tù (leptokurtic) Đồ thị hàm phân bố tích lũy thể Hình ứng với hàm mật độ có độ bất đối xứng dương độ bất đối xứng âm Hình thể hàm phân bố tích lũy hàm mật độ xếp chồng hai hàm mật độ khác Hình - Hàm mật độ với độ nhọn bè Hình - Hàm mật độ với độ nhọn tù Hình - Hàm mật độ với độ bất đối xứng dương Hình - Hàm mật độ với độ bất đối xứng âm LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê www.luatminhkhue.vn Hình - Xếp chồng hai hàm mật độ khác Kiểm nghiệm có hướng 6.1 Quy định chung 6.1.1 Kiểm nghiệm có hướng xem xét liên quan đến đặc trưng độ nhọn độ bất đối xứng phân bố quan trắc Chúng dựa thực tế trường hợp biến ngẫu nhiên chuẩn X có trung bình µ = E(X), mơmen trung tâm bậc ba µ3 = E [(X - µ)3] = …(2) mơmen trung tâm chuẩn hóa bậc ba …(3) mơmen trung tâm chuẩn hóa bậc bốn …(4) µ = E[(X - µ)2] µ4 = E[X - µ) ] .(5) (6) bất đối xứng tổng thể lớn hơn, nhỏ khơng; β2 độ nhọn tổng thể dương; β2 - độ tù tổng thể; LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê bất đẳng thức www.luatminhkhue.vn 6.1.2 Trong kiểm nghiệm độ bất đối xứng, đối giả thuyết H1:µ3 > hoặc, tương ứng >0 có nghĩa bất đối xứng dương (xem Hình 5) H1:µ3 < hoặc, tương ứng có nghĩa độ nhọn tù (hàm mật độ leptokurtic) (xem Hình 4) H1:β2 < có nghĩa độ nhọn bè (hàm mật độ platykurtic) (xem Hình 3) So với phân bố chuẩn, phân bố có độ nhọn tù có xu hướng có nhiều giá trị biến gần với trung bình hướng tới hai phía cực trị Ngược lại trường hợp độ nhọn bè 6.1.4 Việc sử dụng kiểm nghiệm có hướng hợp lý có thơng tin cụ thể cách thức phân bố thực khác biệt so với phân bố chuẩn Thơng tin có từ tính chất tự nhiên liệu loại nhiễu ảnh hưởng tới trình tạo liệu Ví dụ, thực tế biến khơng âm, có trung bình gần với khơng so với giá trị độ lệch chuẩn, lý bất đối xứng dương phân bố thực Tương tự, nhiễu trình tạo liệu gây pha trộn tổng thể chuẩn trung bình khác phương sai dẫn đến phân bố khơng chuẩn có β2 > 6.1.5 Trong trường hợp, việc lựa chọn kiểm nghiệm hướng cần dựa xem xét chung liên quan đến tính chất quan trắc q trình tạo chúng dựa dạng phân bố cụ thể giá trị quan trắc Trong trường hợp đề cập sau, kết kiểm nghiệm vô hướng coi khách quan 6.1.6 Nếu x1, x2,… xn loạt quan trắc …(7) …(8) j = 2, 3, thống kê kiểm nghiệm độ bất đối xứng độ nhọn tương ứng đại lượng …(9) …(10) LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê www.luatminhkhue.vn xj zj xj zj đo chuyển đổi đo chuyển đổi 141 1,799 176 1,447 183 1,322 160 1,643 190 1,146 174 1,477 123 1,908 153 1,708 155 1,690 162 1,623 164 1,602 167 1,568 183 1,322 179 1,398 150 1,732 78 2,100 134 1,845 173 1,491 170 1,531 168 1,556 144 1,778 Từ Bảng tìm TEP(x) = 0,612 sử dụng chương trình máy tính ngắn đơn giản Đối với n = 25, cách nội suy Bảng 12 tìm p phân vị p = - α = 0,99 0,567 Giá trị TEP(x) tính vượt giá trị tới hạn Vì giả thuyết không bị bác bỏ mức ý nghĩa 0,01 giá trị xj Ngoài ra, từ Bảng tìm TEP(z) = 0,006 sử dụng chương trình máy tính Vì giá trị nhỏ giá trị tới hạn n = 25 nội suy từ Bảng 12 nên không bác bỏ giả thuyết khơng giá trị zj Ví dụ minh họa thực tế biết rõ độ bền kéo đứt sợi tơ nhân tạo phân bố theo phân bố chuẩn logarit VÍ DỤ 6: Ví dụ sau minh họa chi tiết cách tính thống kê kiểm nghiệm TEP theo công thức (15) Cột thứ hai Bảng thể n = 10 giá trị xj sử dụng để thực kiểm nghiệm Epps-Pulley Theo công thức (13) (14), = 10,4 m2 = 11,858 tính Tổng kép số hạng thứ ba công thức (15) chuỗi hữu hạn (n - 1) chuỗi con, chuỗi số có số hạng chuỗi cuối có (n - 1) số hạng Đối với chuỗi đầu tiên, số cố định k = số hạng chuỗi thu j = Trong chuỗi thứ hai, số cố định k = 3; chuỗi có hai số hạng thu j = j = Trong chuỗi cuối cùng, số cố định k = 10 chín số hạng , …, thu j = 1, 2, 3, , Các số hạng chuỗi n - = liệt kê cột thứ ba đến mười Bảng LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê www.luatminhkhue.vn Cột thứ mười hai thể n = 10 số hạng tổng số hạng thứ tư công thức (15) Bảng - Độ bền kéo đứt sợi tơ nhân tạo - Tính thống kê kiểm nghiệm TEP k=2 k=3 k=4 k=5 k=6 k=7 k=8 k=9 k = 10 j=1 j = 1,2 j = j = j = j = j = j = j =1 j = 10 0,9996 0,8977 0,2192 0,2083 0,1684 0,0769 0,0587 0,0304 0,0205 0,5285 j xj 4,9 5,0 - 6,5 - - 10,9 - - - 11,0 - - - - 11,4 - - - - - 12,7 - - - - - - 13,1 - - - - - - - 14,0 - - - - - - - - 0,9895 0,7609 10 14,5 - - - - - - - - - 0,7016 0,9095 0,2304 0,2192 0,1778 0,0821 0,0629 0,0329 0,0222 0,4421 0,4258 0,3633 0,1977 0,1593 0,0933 0,0673 0,9996 0,9895 0,8723 0,8154 0,6668 0,5790 0,9933 0,8853 0,8303 0,6842 0,5966 0,9312 0,8853 0,7520 0,6668 0,9933 0,9312 0,8723 0,9664 0,9207 Tổng 104,0 0,9996 1,8072 0,8916 1,8528 2,6923 3,0455 3,8052 4,1573 4,7350 Tổng cộng 0,5407 0,7257 0,9947 0,9924 0,9791 0,8945 0,8575 7,9757 23,9865 Mỗi số mười cột cuối Bảng 6, tổng chúng tính nhập cuối cột Tất 45 số hạng thuộc tổng số hạng thứ ba công thức (15) cộng lại thành giá trị tổng cộng = 23,9865 Cuối cơng thức (15) tính = 0,2914 Đối với n = 10 Bảng 12 cho thấy p phân vị p = - α = 0,95 0,357 Giá trị TEP = 0,2914 khơng vượt q giá trị tới hạn Vì không bác bỏ giả thuyết không mức ý nghĩa 0,05 ví dụ Kiểm nghiệm sử dụng đồng thời nhiều mẫu độc lập Kiểm nghiệm áp dụng cho nhiều mẫu, mẫu cỡ n với n ≥ 8, nhiên, lý thực tế, Bảng 13 giới hạn n ≤ 50 Cơ sở giả định mẫu độc lập lấy từ tổng thể Trong nhiều trường hợp, cần kiểm nghiệm sai lệch so với phân bố chuẩn cách sử dụng nhiều mẫu độc lập mẫu độc lập nhỏ để phát sai lệch đáng kể so với phân bố chuẩn Trong trường hợp này, kiểm nghiệm Shapiro-Wilk áp dụng Đối với h mẫu liên tiếp lấy từ tổng thể mẫu có cỡ n, giá trị Wj (j = 1, 2, …, h) tính theo cơng thức (12) Đối với kiểm nghiệm kết hợp giá trị tương ứng Gj tính từ quan hệ sau đây: … (16) LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê www.luatminhkhue.vn … (17) Các hệ số γ (n), δ(n) ε(n) dùng để chuyển đổi Wj thành biến Gj lấy từ Bảng 13 Trong trường hợp phân bố xét chuẩn biến Gj gần tuân theo phân bố chuẩn chuẩn hóa Giá trị trung bình biến Gj … (18) thống kê kiểm nghiệm Giả thuyết không bị bác bỏ mức ý nghĩa α (19) up = u1-α p phân vị phân bố chuẩn chuẩn hóa VÍ DỤ 7: Ví dụ việc sử dụng kiểm nghiệm kết hợp sử dụng nhiều mẫu độc lập h = 22 mẫu ngẫu nhiên, mẫu cỡ n = 20 lấy từ tổng thể đặc trưng X 20 cá thể đo Đặc trưng không giả định có phân bố chuẩn Đối với số mẫu này, giá trị tương ứng Wj (j = 1, 2, … 22) tính theo cơng thức (12) Trong Bảng liệt kê 22 giá trị Wj Từ Bảng 13, hệ số sau lấy ra: γ (20) = -5,153; δ(20) = 1,802; ε (20) = 0,2359 Sử dụng số này, 22 giá trị tương ứng Gj tính theo cơng thức (16) (17), đồng thời liệt kê Bảng Theo Bảng 11 giá trị tới hạn thống kê W 0,868 n = 20 mức ý nghĩa α = 0,01 Từ Bảng 14, giá trị tới hạn cho = - 2,326 mức ý nghĩa α = 0,01 Bảng - Giá trị Wj Gj 22 mẫu cỡ n = 20 lấy từ tổng thể Mẫu số Wj Gj 0,9543 -0,189 0,9645 +0,292 0,9148 -1,413 0,8864 -2,008 0,9573 -0,059 0,9158 -1,389 0,9462 -0,503 0,9277 -1,083 0,9639 +0,260 10 0,9363 -0,833 11 0,9067 -1,598 12 0,9218 -1,240 13 0,9551 -0,155 14 0,9338 -0,909 15 0,9584 -0,009 j LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 Công ty luật Minh Khuê www.luatminhkhue.vn Mẫu số Wj Gj 16 0,9088 -1,552 17 0,9028 -1,683 18 0,8947 -1,849 19 0,9488 -0,407 20 0,9445 -0,563 21 0,9471 -0,470 22 0,9451 -0,542 j Tổng -17,902 Nếu số 22 mẫu xử lý riêng khơng mẫu cho thấy sai lệch so với phân bố chuẩn mức ý nghĩa α = 0,01 khơng giá trị Wj nhỏ giá trị tới hạn 0,868 không giá trị Gj nhỏ giá trị tới hạn - 2,326 Tuy nhiên, đánh giá kết hợp tất 22 mẫu = - 17,902 / 22 = - 0,814 = - 3,82 Giá trị so với giá trị tới hạn -up = - 2,326 mức ý nghĩa α = 0,01 cho Vì giá trị tính - 3,82 nằm giá trị tới hạn nên giả thuyết không bị bác bỏ mức ý nghĩa α = 0,01 10 Bảng thống kê Bảng - Kiểm nghiệm độ bất đối xứng, (p phân vị n p p = - α = 0,95 0,99) n 0,95 0,99 0,99 1,42 0,97 10 p 0,95 0,99 400 0,20 0,28 1,41 450 0,19 0,27 0,95 1,39 500 0,18 0,26 12 0,91 1,34 550 0,17 0,24 15 0,85 1,26 600 0,16 0,23 20 0,77 1,15 650 0,16 0,22 25 0,71 1,06 700 0,15 0,22 30 0,66 0,98 750 0,15 0,21 35 0,62 0,92 800 0,14 0,20 40 0,59 0,87 850 0,14 0,20 45 0,56 0,82 900 0,13 0,19 50 0,53 0,79 950 0,13 0,18 60 0,49 0,72 1000 0,13 0,18 70 0,46 0,67 1200 0,12 0,16 80 0,43 0,63 1400 0,11 0,15 LUẬT SƯ TƯ VẤN PHÁP LUẬT 24/7 GỌI 1900 6162 ... giá trị quan trắc vẽ giấy đồ thị xác suất chuẩn Trên đồ thị này, trục (trong tiêu chuẩn trục tung) có thang phi tuyến tính theo vùng nằm hàm phân bố chuẩn chuẩn hóa ghi giá trị tương ứng tần suất... xét chuẩn biến Gj gần tuân theo phân bố chuẩn chuẩn hóa Giá trị trung bình biến Gj … (18) thống kê kiểm nghiệm Giả thuyết không bị bác bỏ mức ý nghĩa α (19) up = u1-α p phân vị phân bố chuẩn chuẩn... mẫu µ4 mơmen trung tâm bậc bốn mẫu σ độ lệch chuẩn tổng thể ( ) Khái quát 4.1 Có nhiều loại phép kiểm nghiệm sai lệch so với tính chuẩn Trong tiêu chuẩn xem xét phương pháp đồ thị, kiểm nghiệm