如何在 R 中找到单向方差分析的 F 的临界值？

要找到 R 中单向方差分析的 F 临界值，我们可以按照以下步骤操作 -

• 首先，创建一个包含一个分类列和一个数字列的数据框。

• 然后，使用 aov 函数查找方差分析表。

• 之后，使用 qf 函数为单向方差分析找到 F 的临界值。

创建数据框

让我们创建一个数据框，如下所示 -

Grp<-sample(LETTERS[1:4],20,replace=TRUE)
Score<-rnorm(20)
df<-data.frame(Grp,Score)
df

执行时，上述脚本生成以下内容output(this output will vary on your system due to randomization)-

   Grp   Score
1  B    1.75508031
2  D   -1.43867197
3  B   -0.18409783
4  C   -1.72435769
5  C   -0.95996448
6  B    0.19776077
7  B    0.52247374
8  C   -1.28337249
9  C   -0.63236147
10 B    0.93966870
11 B    0.23925102
12 B    0.86565239
13 B    0.07353123
14 A   -0.62096596
15 A    -1.76680335
16 A    0.12203536
17 B    0.66276852
18 D    -0.50199349
19 A    -1.30960082
20 D    -1.32413279

创建方差分析表

使用 aov 函数和汇总函数查找方差分析表 -

Grp<-sample(LETTERS[1:4],20,replace=TRUE)
Score<-rnorm(20)
df<-data.frame(Grp,Score)
ANOVA<-aov(Score~Grp,data=df)
ANOVA<-summary(ANOVA)
ANOVA

输出

           Df Sum    Sq   Mean Sq F value Pr(>F)
Grp        3  12.870 4.29 11.6 0.000274 ***
Residuals 16  5.917 0.37
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

求单向方差分析的 F 临界值

使用 qf 函数为单向方差分析找到 F 的临界值 -

Grp<-sample(LETTERS[1:4],20,replace=TRUE)
Score<-rnorm(20)
df<-data.frame(Grp,Score)
ANOVA<-aov(Score~Grp,data=df)
ANOVA<-summary(ANOVA)
qf(1-0.05,ANOVA[[1]][1,1],ANOVA[[1]][2,1])

输出

[1] 3.238872