มีความดีของแอนเดอร์สันที่รักการทดสอบพอดีสำหรับชุดข้อมูลสองชุดหรือไม่?


10

ฉันรู้ว่า ad.test () สามารถใช้สำหรับทดสอบมาตรฐานได้

เป็นไปได้ไหมที่จะรับ ad.test เพื่อเปรียบเทียบการแจกแจงจากตัวอย่างข้อมูลสองชุด?

x <- rnorm(1000)
y <- rgev(2000)
ad.test(x,y)

ฉันจะทำการทดสอบ Anderson-Darling กับ 2 ตัวอย่างได้อย่างไร


2
บทความวิกิพีเดียในการทดสอบการโฆษณากล่าวถึงนี้ภายใต้หัวข้อ "Non-พาราทดสอบ K-ตัวอย่าง." เอกสารอ้างอิงปี 1987 ของ JASA โดย Sholz และ Stephens สามารถดูได้ที่cithep.caltech.edu/~fcp/statistics/hypothesisTest/ฟรี
whuber

หากคำถามคือ: ฉันจะทำอย่างไรใน R (ตามที่แท็กแนะนำ): คำถามที่ดี (+1) (และคำตอบน่าจะเป็น: rig ด้วยตัวคุณเอง) แม้ว่าจะถูกวางผิดที่ค่อนข้างที่นี่ ( StackOverflowเป็นสถานที่ที่ดีกว่าสำหรับประเภทนี้ ของคำถาม)
Nick Sabbe

@Nick การค้นหาหรือใช้การทดสอบ GoF ไม่ว่าจะในภาษา R หรือภาษาอื่น ๆ นั้นเหมาะสมกับความสนใจของเราในทุกสิ่งที่เกี่ยวกับสถิติ
whuber

1
@whuber: ฉันยืนแก้ไข: ฉันเพิ่งอ่านส่วนที่เกี่ยวข้องของคำถามที่พบบ่อย ถึงกระนั้นก็เป็นเส้นบาง ๆ ระหว่างความรักและความเกลียดชัง แต่ฉันไม่ได้ลงคะแนนให้โยกย้าย :-)
Nick Sabbe

2
@Nick ฉันเห็นด้วยกับเส้นบาง ๆ เมื่อคำถามมุ่งเน้นไปที่กลไกการเขียนโปรแกรมล้วนๆความเหมาะสมของมันจะกลายเป็นที่น่าสงสัย คุณสามารถค้นหาการสนทนาเป็นระยะเกี่ยวกับเรื่องนี้ในเมตา
whuber

คำตอบ:


7

แพคเกจadkถูกแทนที่ด้วยแพคเกจkSamples:

ลอง:

install.packages("kSamples")  
library(kSamples)
ad.test(runif(50), rnorm(30))

kSamples::ad.testฟังก์ชั่นค่อนข้างช้า มีทางเลือกที่มีประสิทธิภาพมากกว่านี้ไหม?
Nemesi

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.