ความแตกต่างระหว่าง AIC () และ extractAIC () ใน R คืออะไร?


11

เอกสาร R สำหรับทั้งไม่ส่องแสงมาก สิ่งที่ฉันได้รับจากลิงค์นี้คือการใช้อย่างใดอย่างหนึ่งควรจะดี สิ่งที่ฉันไม่ได้รับคือทำไมพวกเขาไม่เท่ากัน

ความจริง: ฟังก์ชั่นแบบขั้นตอนการถดถอยในการวิจัย, การใช้งานstep()extractAIC()

ที่น่าสนใจการทำงานlm()รูปแบบและglm()'null' รุ่น (เฉพาะตัด) ในชุดข้อมูล 'mtcars ของ R ให้ผลลัพธ์ที่แตกต่างกันสำหรับและAICextractAIC()

> null.glm = glm(mtcars$mpg~1)
> null.lm = lm(mtcars$mpg~1)

> AIC(null.glm)
[1] 208.7555
> AIC(null.lm)
[1] 208.7555
> extractAIC(null.glm)
[1]   1.0000 208.7555
> extractAIC(null.lm)
[1]   1.0000 115.9434

เป็นเรื่องแปลกเพราะทั้งสองรุ่นข้างต้นเหมือนกันและAIC()ให้ผลลัพธ์ที่เหมือนกันสำหรับทั้งคู่

ทุกคนสามารถโยนแสงในปัญหาหรือไม่

คำตอบ:


6

ตามความช่วยเหลือสำหรับทั้งสองฟังก์ชั่น (ใช้? AIC และ? extractAIC) นี้คาดว่า

โปรดทราบว่า AIC นั้นถูกกำหนดให้มีค่าคงที่ที่เพิ่มขึ้นเนื่องจากนี่เป็นกรณีสำหรับโอกาสในการบันทึก ซึ่งหมายความว่าคุณควรตรวจสอบว่า

extractAIC(full.modell) - extractAIC(null.modell)

และ

AIC(full.modell) - AIC(null.modell)

ให้ผลเหมือนกัน ตราบใดที่พวกเขาทำฟังก์ชั่นทั้งสองจะเทียบเท่าสำหรับวัตถุประสงค์ในทางปฏิบัติ


2
ฉันอาจหายไปบางอย่าง แต่ผมก็ยังไม่เข้าใจว่าทำไมextractAIC(null.lm) != AIC(null.lm)ในขณะที่extractAIC(null.glm) == AIC(null.glm)แม้ว่าจะเป็นรูปแบบเดียวกับnull.lm null.glmคุณช่วยขยายคำตอบหน่อยได้ไหม?
smillig

2
@smillig extractAICใช้วิธีการที่แตกต่างกันสำหรับlmพอดีและglmเหมาะคือและextractAIC.lm extractAIC.glmคุณสามารถใช้getAnywhereเพื่อศึกษารหัสของพวกเขา AICใช้วิธีการเดียวกันสำหรับทั้งคู่
Roland

ฉันมีโมเดลหลายคู่ (พร้อมตัวทำนายหลายตัว) ซึ่งทั้งสองฟังก์ชั่นให้ผลลัพธ์ที่แตกต่างกัน รุ่น 1: y = x1 + x2, รุ่น 2: y = z + x1 + x2 * z extractAIC()ให้ค่า (ลบ) ที่ลดลงสำหรับรุ่น 1 ในขณะที่ AIC ให้ค่า (บวก) ที่ลดลงสำหรับรุ่น 2
Maxim.K

1
@ Maxim.K คุณให้ข้อมูลเล็กน้อยเกี่ยวกับประเภทของตัวแปรและโมเดลที่ใช้ หากคุณทำและมีความแตกต่างบางอย่างกับคำถามนี้มันอาจคุ้มค่าที่จะโพสต์คำถามนี้เป็นคำถามใหม่ ยากที่จะพูดโดยไม่ทราบรายละเอียด
Erik

@ เอริคฉันสงสัยว่ามันจะคุ้มค่ามากถ้าฉันบอกว่า z ต่อเนื่องและ x2 เป็นหมวดหมู่ (dummified) จะต้องใช้ข้อมูลในการทำซ้ำและฉันไม่สามารถเผยแพร่พวกเขาฉันกลัว
Maxim.K
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.