วิธีรับค่าของ Mean squared error ในการถดถอยเชิงเส้นใน R


20

ปล่อยให้โมเดลการถดถอยเชิงเส้นที่ได้รับจากฟังก์ชัน R ต้องการทราบว่าเป็นไปได้ที่จะได้รับจากคำสั่ง Mean Squared Error

ฉันมีตัวอย่างผลลัพธ์ต่อไปนี้

> lm <- lm(MuscleMAss~Age,data)
> sm<-summary(lm)
> sm

Call:
lm(formula = MuscleMAss ~ Age, data = data)

Residuals:
     Min       1Q   Median       3Q      Max 
-16.1368  -6.1968  -0.5969   6.7607  23.4731 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) 156.3466     5.5123   28.36   <2e-16 ***
Age          -1.1900     0.0902  -13.19   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1   1

Residual standard error: 8.173 on 58 degrees of freedom
Multiple R-squared:  0.7501,    Adjusted R-squared:  0.7458 
F-statistic: 174.1 on 1 and 58 DF,  p-value: < 2.2e-16

หลาย R-squared เป็นข้อผิดพลาดกำลังสอง? ถ้าคำตอบคือไม่สามารถอธิบายความหมายของ Multiple R-squared และ Multiple R-squared

คำตอบ:


25

ค่า R หลายค่าที่รายงาน R คือสัมประสิทธิ์การตัดสินใจซึ่งกำหนดโดยสูตร

R2=1SSresSStot.

ผลรวมของ squared ข้อผิดพลาดจะได้รับ (ขอบคุณก่อนหน้านี้คำตอบ ) sum(sm$residuals^2)โดย

mean(sm$residuals^2)ข้อผิดพลาดยกกำลังสองเฉลี่ยจะได้รับจาก คุณสามารถเขียนฟังก์ชันเพื่อคำนวณสิ่งนี้ได้เช่น:

mse <- function(sm) 
    mean(sm$residuals^2)

5
sm$sigma^2 * sm$fstatistic[3]/(1+sum(sm$fstatistic[2:3]))8.173581+18.1732×58/(1+1+58)=64.57

6
Antepenultimate เป็นคำที่ยอดเยี่ยม
fbt

1
ตามปกติวิธี neater มากคือการเขียนmse <- function(sm) mean(sm$residuals^2)
Marius Hofert

4
> นั่นหมายความว่าข้อผิดพลาดกำลังสองถูกกำหนดโดยค่าส่วนที่เหลือ ^ 2 / ข้อผิดพลาด df จากตาราง ANOVA แทนที่จะเป็นค่าเฉลี่ย (ค่าส่วนที่เหลือ ^ 2) หลังคือค่าเฉลี่ยข้อผิดพลาดการทำนายกำลังสอง ไม่แน่ใจว่าฉันขาดความเข้าใจไปหรือเปล่า ขออภัยฉันไม่มีคะแนนชื่อเสียงเพียงพอที่จะโพสต์ความคิดเห็น
SrikanthRaja

1
แต่ฉันเห็นด้วยกับสิ่งที่ผู้ใช้ 45409 พูด ตามสิ่งที่ฉันเรียนรู้ MSE = ส่วนที่เหลือ ^ 2 / ข้อผิดพลาด df
vtshen
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.