如何让 R 吐出基于方差分析的 F 统计量的临界值?
How do I get R to spit out the critical value for F-statistic based on ANOVA?
R 中的方差分析缺少的一件事是它不会自动显示临界值。其他一切都给定了。我可以看出我的 F 值远高于应有的值,但我想知道截止点的余量。有这个在线计算器可以根据自由度得出 F 统计量的临界值,但我希望 R 可以这样做。 http://www.danielsoper.com/statcalc/calculator.aspx?id=4
我该怎么做?
示例:
>anova(anovaModel.model1)
Df Sum Sq Mean Sq F value Pr(>F)
data$SIZE 4 0.1193 0.027926 22.056 4.55e-16 ***
Residuals 1372 1.994 0.001352
>F(variance.mod1) #??
>F(4,1372 ) #Something like this?
试试这个:
mylm <- lm(wt~mpg, data = mtcars)
myanova <- anova(mylm)
cbind(myanova, 'CriticalValue' = qf(1-.05, myanova[1,1], myanova[2,1]))
Df Sum Sq Mean Sq F value Pr(>F) CriticalValue
mpg 1 22.343135 22.3431348 91.37533 1.293959e-10 4.170877
Residuals 30 7.335613 0.2445204 NA NA 4.170877
在这种情况下,qf
函数是您的朋友。
alpha = .05
qf(1-alpha, myanova[1,1], myanova[2,1])
[1] 4.170877
R 中的方差分析缺少的一件事是它不会自动显示临界值。其他一切都给定了。我可以看出我的 F 值远高于应有的值,但我想知道截止点的余量。有这个在线计算器可以根据自由度得出 F 统计量的临界值,但我希望 R 可以这样做。 http://www.danielsoper.com/statcalc/calculator.aspx?id=4
我该怎么做?
示例:
>anova(anovaModel.model1)
Df Sum Sq Mean Sq F value Pr(>F)
data$SIZE 4 0.1193 0.027926 22.056 4.55e-16 ***
Residuals 1372 1.994 0.001352
>F(variance.mod1) #??
>F(4,1372 ) #Something like this?
试试这个:
mylm <- lm(wt~mpg, data = mtcars)
myanova <- anova(mylm)
cbind(myanova, 'CriticalValue' = qf(1-.05, myanova[1,1], myanova[2,1]))
Df Sum Sq Mean Sq F value Pr(>F) CriticalValue
mpg 1 22.343135 22.3431348 91.37533 1.293959e-10 4.170877
Residuals 30 7.335613 0.2445204 NA NA 4.170877
在这种情况下,qf
函数是您的朋友。
alpha = .05
qf(1-alpha, myanova[1,1], myanova[2,1])
[1] 4.170877