ฉันพยายามทำการถดถอย logit สั่ง ฉันกำลังใช้รูปแบบเช่นนั้น (เพียงแบบจำลองโง่ ๆ ประมาณจำนวน บริษัท ในตลาดจากรายรับและมาตรการประชากร) คำถามของฉันเกี่ยวกับการทำนาย
nfirm.opr<-polr(y~pop0+inc0, Hess = TRUE)
pr_out<-predict(nfirm.opr)
เมื่อฉันใช้ทำนาย (ซึ่งฉันพยายามใช้เพื่อให้ได้ค่า y ที่คาดการณ์) ผลลัพธ์จะเป็น 0, 3, หรือ 27 ซึ่งไม่มีทางที่จะสะท้อนให้เห็นถึงสิ่งที่ควรจะเป็นคำทำนายจากค่าสัมประสิทธิ์ของฉันเองจากค่าสัมประสิทธิ์ ประมาณการและดัก ไม่มีใครรู้ว่าจะได้รับการทำนาย "ถูกต้อง" สำหรับรุ่น logit ของฉันได้อย่างไร
แก้ไข
เพื่อชี้แจงข้อกังวลของฉันข้อมูลการตอบสนองของฉันมีการสังเกตในทุกระดับ
>head(table(y))
y
0 1 2 3 4 5
29 21 19 27 15 16
ที่ซึ่งตัวแปรทำนายของฉันดูเหมือนจะพัวพัน
> head(table(pr_out))
pr_out
0 1 2 3 4 5
117 0 0 114 0 0
polr()
ฟังก์ชั่นอื่น ๆ หรือไม่? คุณอาจจะลองlrm()
จากแพคเกจ:rms
lrmFit <- lrm(y ~ pop0 + inc0); predict(lrmFit, type="fitted.ind")
ตัวเลือกหนึ่งคือvglm()
จากแพคเกจ:VGAM
vglmFit <- vglm(y ~ pop0 + inc0, family=propodds); predict(vglmFit, type="response")
ทั้งสองคืนเมทริกซ์ของความน่าจะเป็นหมวดหมู่ที่คาดการณ์ไว้ ดูคำตอบของฉันเพื่อรับหมวดหมู่ที่คาดคะเนได้จากที่นั่น
predict
ฟังก์ชันแตกต่างจากค่าที่คุณสร้างขึ้นเองอย่างไร โครงสร้างของตัวแปรตามของคุณคืออะไร? โปรดให้ตัวอย่างที่ทำซ้ำได้