Có cách nào khác để kiểm tra Kolmogorov-Smirnov cho dữ liệu bị ràng buộc với hiệu chỉnh không?


14

Tôi đã nhận được một loạt dữ liệu từ hai mẫu (được kiểm soát và xử lý), mỗi mẫu chứa vài nghìn giá trị phải trải qua thử nghiệm có ý nghĩa trong R. Về mặt lý thuyết, các giá trị phải liên tục, nhưng do làm tròn bằng phần mềm đo lường mà chúng phát sinh ' t và họ đã có quan hệ. Các bản phân phối chưa được biết và hình dạng của các bản phân phối được kiểm soát và xử lý có thể khác nhau, vì vậy tôi muốn sử dụng một thử nghiệm không tham số để so sánh nếu sự khác biệt giữa các mẫu có ý nghĩa đối với 10 yếu tố khác nhau.

Tôi đã nghĩ đến việc sử dụng thử nghiệm Kolmogorov-Smirnov, nhưng nó không thực sự phù hợp với các mối quan hệ. Gần đây tôi tình cờ thấy một thư viện R mới có tên là Match thực hiện phiên bản bootstrap của thử nghiệm KS và chấp nhận các mối quan hệ. Bây giờ đây thực sự là một ý tưởng tốt hay tôi nên sử dụng một thử nghiệm khác thay thế? Và tôi có cần điều chỉnh giá trị p không?


Các giấy liên kết liên quan đến phù hợp với điểm số xu hướng. Có thể là bài kiểm tra bootstrap có tính tổng quát hơn nhưng tôi không chắc chắn.
Michael R. Chernick

Tôi đã thực hiện một phiên bản ngẫu nhiên của một cái gì đó như Kolmogorov-Smirnov (thực ra, có lẽ tôi đã sử dụng Anderson-Darling hoặc Cramer-von Mises cho KS, nhưng vẫn cẩn thận với phân phối ngẫu nhiên quan hệ).
Glen_b -Reinstate Monica

Bạn đã thấy mã của Tom Waterhouse chưa?
Ray Koopman

Câu trả lời:


13

Thay vì sử dụng kiểm tra KS, bạn có thể chỉ cần sử dụng quy trình hoán vị hoặc lấy mẫu lại như được thực hiện trong oneway_testchức năng của coingói. Có một cái nhìn vào câu trả lời được chấp nhận cho câu hỏi này .

Cập nhật : Gói của tôi afexchứa hàm compare.2.vectorsthực hiện hoán vị và các thử nghiệm khác cho hai vectơ. Bạn có thể lấy nó từ CRAN:

install.packages("afex")

Đối với hai vectơ xynó (hiện tại) trả về một cái gì đó như:

> compare.2.vectors(x,y)
$parametric
   test test.statistic test.value test.df       p
1     t              t     -1.861   18.00 0.07919
2 Welch              t     -1.861   17.78 0.07939

$nonparametric
             test test.statistic test.value test.df       p
1 stats::Wilcoxon              W     25.500      NA 0.06933
2     permutation              Z     -1.751      NA 0.08154
3  coin::Wilcoxon              Z     -1.854      NA 0.06487
4          median              Z      1.744      NA 0.17867

Bất kỳ ý kiến ​​liên quan đến chức năng này đều được hoan nghênh.


3
(+1) Có thể tìm thấy mô tả về thử nghiệm này và các thử nghiệm khác trong blog này

@Henrik Cảm ơn bạn đã gợi ý và chỉ ra câu hỏi khác. Điều đó thực sự hữu ích!
AnjaM

@AnjaM Bạn được chào đón. Bạn cũng có thể muốn kiểm tra cập nhật của tôi.
Henrik
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.