Kolmogorov-Smirnov hai chiều

Tôi muốn chạy một số thử nghiệm Kolmogorov - Smironov hai chiều để xác định xem phân phối hai chiều có phù hợp với tham chiếu hay không.

Có gói hoặc ứng dụng nào tôi có thể sử dụng theo cách tương đối đơn giản không? Hoặc có một thuật toán khác là thích hợp hơn? Tôi chỉ có một kiến ​​thức thống kê cơ bản.

Một phần mở rộng hai chiều của thử nghiệm Kolmogorov-Smirnov đã được Justel, Pena và Zamar mô tả trong một "Thử nghiệm đa năng Komogorov-Smirnov về sự phù hợp" . Nhận xét của @ Procrastinator cho thấy có thể có những đề xuất khác như vậy.

Tuy nhiên, tôi đã không thấy một gói với việc thực hiện đơn giản.

Tùy thuộc vào những gì bạn muốn làm, kde.test () trong gói ks của Tarn Duong cho R có thể hữu ích hơn.

Triển khai Python

Tôi đã viết một triển khai python bằng cách sử dụng numpy. Bạn có thể tìm thấy mã ở đây , bạn có thể tìm thấy nhiều thông tin hơn trong chuỗi doc trong mã.

Và đây là một cái khác (không phải bởi tôi). Notebook này cung cấp một triển khai Python cho thử nghiệm 2D KS với 2 mẫu. Các .pytập tin có thể được tải về ở đây . Mã dường như là một bản dịch Cmã thẳng , hiệu quả có thể là một vấn đề nếu kích thước mẫu lớn.

Tuy nhiên, tốt hơn hết bạn nên kiểm tra mã (bất kể là mã nào) với giấy tờ / sách gốc trước khi sử dụng. Việc triển khai python của thử nghiệm KS 2d được kiểm tra ít hơn nhiều so với các thử nghiệm trong R.

Thêm thông tin

Thuật toán được phát triển đầu tiên trong hai bài báo (như tôi thấy)

• Peacock, JA 1983, Thử nghiệm lòng tốt hai chiều trong thiên văn học
• Fasano, G. và Franceschini, A. 1987, Phiên bản đa chiều của Thử nghiệm Kolmogorov-Smirnov .

Một giới thiệu tốt đẹp và việc Cthực hiện có thể được tìm thấy trong

• Báo chí, WH và cộng sự. 1992, Công thức toán số trong C , Mục 14.7, p645.

  Bạn có thể tìm thấy C++/Fortrantriển khai trong các phiên bản khác của cuốn sách.

Đây là một bài đăng có tiêu đề Hãy xem thử nghiệm Kolmogorov-Smirnov cũng liên quan đến chủ đề này, bạn có thể muốn xem qua. Nó khuyến khích sử dụng phương pháp lấy mẫu lại để đánh giá giá trị p với khoảng cách KS đã cho.

bạn có thể thấy mã Matlab này hữu ích.

http://www.mathworks.com/matlabcentral/fileexchange/38617-two-dimensional-2d- Paired-kolmogorov-smirnov-test