การทดสอบตัวผู้ใน R [ปิด]


14

ในการทดสอบสมมติฐานการถดถอยแบบขนานในการถดถอยโลจิสติกอันดับฉันพบว่ามีหลายวิธี ฉันใช้ทั้งวิธีกราฟิก (ตามรายละเอียดในหนังสือของ Harrell) และวิธีรายละเอียดโดยใช้ แพ็คเกจลำดับใน R

อย่างไรก็ตามฉันต้องการรันการทดสอบ Brant (จาก Stata) สำหรับทั้งตัวแปรเดี่ยวและสำหรับโมเดลทั้งหมด ฉันได้ดูไปรอบ ๆ แต่ไม่พบการใช้งานใน R

มีการนำการทดสอบ Brant ไปใช้ใน R หรือไม่?


ไม่ทราบเกี่ยวกับการนำไปใช้ใด ๆ แต่ในส่วนนี้จากหนังสือของ J. Scott Long ให้คำอธิบายโดยละเอียดเกี่ยวกับวิธีการคำนวณการทดสอบที่ฉันเชื่อว่าคุณอ้างถึง
NRH

ขอบคุณ - ฉันเคยดูกระดาษต้นฉบับและไฟล์. ada จาก stata ด้วย อย่างไรก็ตามการเขียนโปรแกรมที่จำเป็นอยู่นอกเหนือระดับของฉัน
Misha

คำตอบ:


5

ฉันใช้การทดสอบตัวผู้ในอาร์แพคเกจและฟังก์ชั่นเรียกว่าแบรนต์และตอนนี้ใช้ได้ใน CRAN แล้ว

การทดสอบตัวผู้ถูกกำหนดโดย Rollin Brant เพื่อทดสอบสมมติฐานการถดถอยแบบขนาน (Brant, R. (1990) การประเมินสัดส่วนในรูปแบบอัตราต่อรองแบบสัดส่วนสำหรับการถดถอยโลจิสติกตามลำดับชีวภาพ , 46 , 1171–1178)

นี่คือตัวอย่างรหัส:

data = MASS::survey
data$Smoke = ordered(MASS::survey$Smoke, levels=c("Never","Occas","Regul","Heavy"))
model1 = MASS::polr(Smoke ~ Sex + Height, data=data, Hess=TRUE)
brant(model1)

ในตัวอย่างสมมติฐานการถดถอยแบบขนานจะเก็บไว้เนื่องจากค่า p ทั้งหมดมีค่ามากกว่า 0.05 รถโดยสารสำหรับทั้งรุ่นส่วนที่เหลือสำหรับค่าสัมประสิทธิ์ indvidual


2

ใช่ - อันที่จริงแพ็คเกจลำดับที่คุณเชื่อมโยงสามารถทำได้ (แม้ว่าพวกเขาจะไม่เรียกมันว่าการทดสอบตัวผู้) ลองดูที่หน้า 6 และ 7 ของลิงค์ของคุณซึ่งแสดงให้เห็นว่า "การทดสอบอัตราส่วนความน่าจะเป็นของความลาดชันที่เท่ากันหรือการคาดคะเนอัตราต่อรองแบบสัดส่วน" ซึ่งเป็นสิ่งที่คุณกำลังมองหา


ฉันได้เปรียบเทียบผลลัพธ์ระหว่างสองแนวทาง แต่ไม่เหมือนกัน ฉันเชื่อว่าการทดสอบตัวผู้เป็นแบบทดสอบคะแนนมากกว่า
Misha

5
ไม่ในตัวอย่างที่ จำกัด วิธีการทั้งหมดเหล่านี้แตกต่างกัน แต่ไม่มีอาการพวกเขาควรจะเหมือนกัน การทดสอบตัวผู้ประเมินการประมาณค่าของโมเดลที่ไม่มีข้อ จำกัด โดยใช้การถดถอยแบบแยกส่วนและทำการทดสอบแบบ Wald การเปรียบเทียบวิธีการต่าง ๆ สามารถพบได้ที่นี่
Maarten Buis

2

บางบันทึกในหัวข้อ

แพคเกจ R VGAMในCumulativeคำสั่ง (ลำดับการถดถอยที่มีความน่าจะเป็นที่สะสม) parallel=FALSEช่วยให้การเปลี่ยนแปลงอัตราต่อรองสัดส่วนสมมติฐานที่มีตัวเลือก

เป็นที่ทราบกันว่าเป็นปัญหาที่พบบ่อย (จากหนังสือ: แบบจำลองการถดถอยสำหรับตัวแปรตามหมวดหมู่โดยใช้ Stata, Second Edition, โดย J. Scott Long, Jeremy Freese)

"ข้อแม้เกี่ยวกับสมมติฐานการถดถอยแบบขนาน: เราพบว่ามีการละเมิดสมมติฐานการถดถอยแบบขนาน (PRA) บ่อยครั้งเมื่อมีการปฏิเสธแบบจำลองทางเลือกอื่น ๆ ที่ไม่ได้กำหนดข้อ จำกัด ของการถดถอยแบบขนานควรได้รับการพิจารณา การถดถอย usig OLS เนื่องจากข้อสันนิษฐานโดยนัยของการใช้ LRM กับข้อมูลลำดับนั้นแข็งแกร่งยิ่งขึ้นรูปแบบทางเลือกที่สามารถพิจารณาได้รวมถึงแบบจำลองสำหรับผลลัพธ์ที่ระบุ [... ] Stereotype Logistic model หรือ Stereotype Model model; Logized Ordered General General Model; แบบจำลองอัตราส่วนความต่อเนื่องเป็นทางเลือก "(หน้า 221)

บทความนี้มีความลึกในหัวข้อนี้ชัดเจนและเขียนได้ดี แต่ไม่พิจารณาแพ็คเกจ VGAM หรือคำสั่ง "สะสม": การถดถอยโลจิสติกอันดับในการศึกษาทางระบาดวิทยา


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