มีใครช่วยอธิบายโมเดล Cox ของฉันให้ฉันฟังเป็นภาษาอังกฤษธรรมดาได้ไหม
ฉันติดตั้งโมเดลการถดถอยของ Cox ต่อไปนี้กับข้อมูลทั้งหมดของฉันโดยใช้cph
ฟังก์ชั่น Data
ข้อมูลของฉันจะถูกบันทึกไว้ในวัตถุที่เรียกว่า ตัวแปรw
, x
และy
มีความต่อเนื่อง z
เป็นปัจจัยสองระดับ เวลามีหน่วยวัดเป็นเดือน ผู้ป่วยบางรายของฉันขาดข้อมูลสำหรับตัวแปรz
( หมายเหตุ : ฉันได้ระบุไว้อย่างชัดเจนว่าคำแนะนำของดร. ฮาร์เรลด้านล่างนี้ว่าฉันใส่ค่าเหล่านี้เพื่อหลีกเลี่ยงการทำให้ลำเอียงแบบของฉันและจะทำในอนาคต)
> fit <- cph(formula = Surv(time, event) ~ w + x + y + z, data = Data, x = T, y = T, surv = T, time.inc = 12)
Cox Proportional Hazards Model
Frequencies of Missing Values Due to Each Variable
Surv(time, event) w x y z
0 0 0 0 14
Model Tests Discrimination
Indexes
Obs 152 LR chi2 8.33 R2 0.054
Events 64 d.f. 4 g 0.437
Center 0.7261 Pr(> chi2) 0.0803 gr 1.548
Score chi2 8.07
Pr(> chi2) 0.0891
Coef S.E. Wald Z Pr(>|Z|)
w -0.0133 0.0503 -0.26 0.7914
x -0.0388 0.0351 -1.11 0.2679
y -0.0363 0.0491 -0.74 0.4600
z=1 0.3208 0.2540 1.26 0.2067
ฉันพยายามทดสอบสมมติฐานเกี่ยวกับอันตรายตามสัดส่วนโดยใช้cox.zph
คำสั่งด้านล่าง แต่ไม่ทราบวิธีตีความผลลัพธ์ การใส่plot()
คำสั่งให้ข้อความแสดงข้อผิดพลาด
cox.zph(fit, transform="km", global=TRUE)
rho chisq p
w -0.1125 1.312 0.2520
x 0.0402 0.179 0.6725
y 0.2349 4.527 0.0334
z=1 0.0906 0.512 0.4742
GLOBAL NA 5.558 0.2347
ปัญหาแรก
- ใครสามารถอธิบายผลลัพธ์ของผลลัพธ์ข้างต้นให้ฉันเป็นภาษาอังกฤษธรรมดาได้หรือไม่ ฉันมีพื้นฐานทางการแพทย์และไม่มีการฝึกอบรมอย่างเป็นทางการในสถิติ
ปัญหาที่สอง
ตามที่ดร. ฮาร์เรลล์แนะนำฉันต้องการตรวจสอบแบบจำลองของฉันภายในด้วยการตรวจสอบความถูกต้องแบบข้าม 10 ครั้ง 100 ครั้งโดยใช้
rms
แพ็คเกจ (จากสิ่งที่ฉันเข้าใจสิ่งนี้จะนำ100 * 10 = 1000
ไปสู่ ของผู้ป่วยที่พวกเขาไม่เคยเห็น)ฉันลองใช้
validate
ฟังก์ชั่นดังที่แสดง> v1 <- validate(fit, method="crossvalidation", B = 10, dxy=T) > v1 index.orig training test optimism index.corrected n Dxy -0.2542 -0.2578 -0.1356 -0.1223 -0.1320 10 R2 0.0543 0.0565 0.1372 -0.0806 0.1350 10 Slope 1.0000 1.0000 0.9107 0.0893 0.9107 10 D 0.0122 0.0128 0.0404 -0.0276 0.0397 10 U -0.0033 -0.0038 0.0873 -0.0911 0.0878 10 Q 0.0155 0.0166 -0.0470 0.0636 -0.0481 10 g 0.4369 0.4424 0.6754 -0.2331 0.6700 10
คุณจะทำการสุ่มตัวอย่าง 100x ใหม่ได้อย่างไร ฉันคิดว่าโค้ดด้านบนของฉันทำการตรวจสอบไขว้เพียงครั้งเดียวเท่านั้น
จากนั้นฉันก็อยากรู้ว่าแบบจำลองของฉันดีแค่ไหนที่ทำนายไว้ ฉันพยายามต่อไปนี้:
> c_index <- abs(v1[1,5])/2 + 0.5 > c_index [1] 0.565984
นี่หมายความว่าแบบจำลองของฉันดีกว่าการพลิกเหรียญเพียงเล็กน้อยเท่านั้นหรือ
ปัญหาที่สาม
ดร. ฮาร์เรลล์ชี้ให้เห็นว่าฉันสันนิษฐานว่าเป็นเส้นตรงสำหรับเอฟเฟกต์โควาเรียตและจำนวนของเหตุการณ์ในตัวอย่างของฉันมีขนาดใหญ่พอที่จะพอดีกับแบบจำลองที่เชื่อถือได้หากเอฟเฟกต์ covariate ทั้งหมดเป็นเส้นตรง
- นี่หมายความว่าฉันควรจะรวมคำศัพท์บางอย่างไว้ในแบบจำลองของฉันหรือไม่? ถ้าเป็นเช่นนั้นคำแนะนำใด ๆ เกี่ยวกับสิ่งที่จะนำไป?
cph
ผลลัพธ์ข้างต้นให้ฉันเป็นภาษาอังกฤษธรรมดาหรือชี้ให้ฉันอ้างอิงที่จะทำ ดร. ฮาร์เรลขอบคุณมากสำหรับความช่วยเหลือของคุณจนถึงตอนนี้!