如何在R中找到單因素方差分析的臨界值F?


要在R中找到單因素方差分析的臨界值F,我們可以按照以下步驟進行:

  • 首先,建立一個包含一個分類列和一個數值列的資料框。
  • 然後,使用aov函式查詢方差分析表。
  • 之後,使用qf函式查詢單因素方差分析的臨界值F。

建立資料框

讓我們建立一個如下所示的資料框:

 線上演示

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

執行上述指令碼後,將生成以下輸出(由於隨機化,此輸出在您的系統上會有所不同):

   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函式和summary函式查詢方差分析表:

 線上演示

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

更新於:2021年8月14日

2K+ 次檢視

啟動您的職業生涯

完成課程後獲得認證

開始學習
廣告
© . All rights reserved.