日韩黑丝制服一区视频播放|日韩欧美人妻丝袜视频在线观看|九九影院一级蜜桃|亚洲中文在线导航|青草草视频在线观看|婷婷五月色伊人网站|日本一区二区在线|国产AV一二三四区毛片|正在播放久草视频|亚洲色图精品一区

分享

科學(xué)網(wǎng)-單因素方差分析之R篇

 楓林秋2016 2019-08-19

2013-11-4 10:12


程序代碼如下:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
######### 代碼清單   #########
#### 建立數(shù)據(jù)集df ####
yield <- scan()
 24 30 28 26
 27 24 21 26
 31 28 25 30
 32 33 33 28
 21 22 16 21
Treat <- rep( paste("A", 1:5, sep = ""), rep(4, 5) )
df <- data.frame( Treat, yield)
######## 方差分析  #####
fit <- aov(yield ~ Treat, data = df)
summary(fit)

運(yùn)行結(jié)果如下:

1
2
3
4
5
6
> summary(fit)
           Df Sum Sq Mean Sq F value   Pr(>F)    
Treat        4  301.2   75.30   11.18 0.000209 ***
Residuals   15  101.0    6.73                    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

當(dāng)F檢驗(yàn)顯著或極顯著時(shí),就有必要進(jìn)行處理平均數(shù)兩兩之間的比較,以判斷不同處理平均數(shù)之間的差異顯著性。這種多個(gè)平均數(shù)兩兩間的相互比較稱為多重比較(multiple comparison)。

運(yùn)行結(jié)果如下:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
##### 多重比較之ducan 方法  ####
library( agricolae )
duncan.test( fit, "Treat", alpha = 0.05 )
##### 運(yùn)行結(jié)果 #####
> duncan.test( fit, "Treat", alpha = 0.05 )
Study:
Duncan's new multiple range test for yield
Mean Square Error:    6.733333
Treat,  means  
  yield  std.err replication
A1  27.0 1.290994           4
A2  24.5 1.322876           4
A3  28.5 1.322876           4
A4  31.5 1.190238           4
A5  20.0 1.354006           4
alpha: 0.05 ; Df Error: 15
Critical Range
      2        3        4        5
3.910886 4.099664 4.216980 4.296902
Means with the same letter are not significantly different.
Groups, Treatments and means
a    A4    31.5
ab   A3    28.5
b    A1    27
b    A2    24.5
c    A5    20

由結(jié)果可知,A4顯著高于A1、A2和A5,但與A3之間無顯著差異。



http://blog.sciencenet.cn/blog-1114360-738841.html 

    本站是提供個(gè)人知識管理的網(wǎng)絡(luò)存儲空間,所有內(nèi)容均由用戶發(fā)布,不代表本站觀點(diǎn)。請注意甄別內(nèi)容中的聯(lián)系方式、誘導(dǎo)購買等信息,謹(jǐn)防詐騙。如發(fā)現(xiàn)有害或侵權(quán)內(nèi)容,請點(diǎn)擊一鍵舉報(bào)。
    轉(zhuǎn)藏 分享 獻(xiàn)花(0

    0條評論

    發(fā)表

    請遵守用戶 評論公約

    類似文章 更多