R语言进行二因素随机区组方差分析及其多重比较（窗口软件GenStat也可实现）

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

【例3.3.1.2】玉米品种A有4个水平，分别是A1-A4，施肥B有两个水平，分别是B1-B2，重复三次，采用随机区组，数据如下：

数据见附件：

数据录入：

> library(reshape);library(agricolae)

>df<-read.table(file="r3.3.1.2.csv",header=T,sep=",")

> df.2<-melt(df,id=c("A","B"))

> colnames(df.2)[3:4]<-c("Blk","y")

方差分析

> duncan.model<-aov(y~A*B+Blk,data=df.2)

> summary(duncan.model) # show ANOVA table

           Df Sum Sq Mean Sq F value   Pr(>F)     

A           3  98.79   32.93  15.199 0.000111 ***

B           1  77.04   77.04  35.558 3.47e-05 ***

Blk          2 20.33   10.17   4.692 0.027567 *   

A:B          3136.46   45.49  20.994 1.90e-05 ***

Residuals   14  30.33   2.17                    

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

> anova.tab(duncan.model) #方差格式整理

           Df Sum Sq Mean Sq F value   Pr(>F)     

A           3   98.8   32.93  15.199 0.000111 ***

B           1   77.0   77.04  35.558 3.47e-05 ***

Blk          2  20.3   10.17   4.692 0.027567 *   

A:B          3 136.5   45.49  20.994 1.90e-05 ***

Residuals   14   30.3   2.17                    

Total       23  363.0                           

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

    由方差分析表可知A、B、AB互作均达到极显著水平，区组达到显著水平，下面进行多重比较。

多重比较

> duncan.test(duncan.model,"A",alpha=0.05,console=T) #对A因素进行多重比较

Study: duncan.model ~ "A"

Duncan's new multiple range test

for y

Mean Square Error:  2.166667

A,  means

          y     std r Min Max

A1 12.00000 1.264911 6  10  13

A2 17.16667 2.926887 6  12  20

A3 13.00000 5.291503 6   7  19

A4 12.66667 3.829708 6   8  17

alpha: 0.05 ; Df Error: 14

Critical Range

       2       3        4

1.822718 1.909920 1.963738

Means with the same letter are notsignificantly different.

Groups, Treatments and means

a A2 17.17

b A3 13

b A4 12.67

b A1 12

    由上表可知，A2水平和其他水平相比差异显著，处理效果最好。

> duncan.test(duncan.model,"B",alpha=0.05,console=T) #对B因素进行多重比较

Study: duncan.model ~ "B"

Duncan's new multiple range test

for y

Mean Square Error:  2.166667

B,  means

          y     std  r Min Max

B1 15.50000 2.540580 12  12  19

B2 11.91667 4.420167 12   7  20

alpha: 0.05 ; Df Error: 14

Critical Range

       2

1.288856

Means with the same letter are notsignificantly different.

Groups, Treatments and means

a B1 15.5

b B2 11.92

    由上表可知，B1和B2相比较，差异显著。

> with(df.2, duncan.test(y, A:B, DFerror = 14, MSerror =2.17,console=T)) #对AB进行多重比较

Study: y ~ A:B

Duncan's new multiple range test

for y

Mean Square Error:  2.17

A:B,  means

             y       std r Min Max

A1:B1 12.666667 0.5773503 3  12 13

A1:B2 11.333333 1.5275252 3  10 13

A2:B1 15.666667 3.5118846 3  12 19

A2:B2 18.666667 1.5275252 3  17 20

A3:B1 17.666667 1.5275252 3  16 19

A3:B2  8.333333 1.5275252 3   7 10

A4:B1 16.000000 1.0000000 3  15 17

A4:B2  9.333333 1.5275252 3   8 11

alpha: 0.05 ; Df Error: 14

Critical Range

       2       3        4        5       6        7        8

2.579695 2.703112 2.779281 2.8307672.867363 2.894153 2.914097

Means with the same letter are notsignificantly different.

Groups, Treatments and means

a A2:B2 18.67

ab A3:B1 17.67

ab A4:B1 16

b A2:B1 15.67

c A1:B1 12.67

cd A1:B2 11.33

de A4:B2 9.333

e A3:B2 8.333

    有互作的多重比较可知，A2:B2互作效果最好，显著高于A2:B1、A1:B1、A1:B2、A4:B2、A3:B2。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏9 回帖

关键词：genstat 多重比较 方差分析 R语言 significant R语言处理二因素随机区组及多重比较

非常感谢您的分享，如果能贡献一些电子版书籍就更好啦！希望您今后多为大家贡献好资料哦！

好东西，就是没有找到附件

yijiaobani 发表于 2015-7-4 19:02
【例3.3.1.2】玉米品种A有4个水平，分别是A1-A4，施肥B有两个水平，分别是B1-B2，重复三次，采用随机区组， ...

好的?

附件重新传上。论坛不支持csv格式上传，我上传的是excel格式，请把数据另存为csv格式打开。

你好，我用你的方法试了，提示并没有找到有duncan.test函数，是怎么回事呢

anova.tab(duncan.model) #方差格式整理
这句为什么提示错误呀

风味1991 发表于 2017-7-16 16:40
anova.tab(duncan.model) #方差格式整理
这句为什么提示错误呀

这里自定义了一个函数anova.tab()
anova.tab<-function(fm){

tab<-summary(fm)

k<-length(tab[[1]])-2

temp<-c(sum(tab[[1]][,1]),sum(tab[[1]][,2]),rep(NA,k))

tab[[1]]["Total",]<-temp

tab

}

好难。我在数据导入那边就出现无数问题。首先是安装了agricolae 包之后用library调用出现错误，然后说无melt 这个函数。求求大神帮帮我

没有melt 这个函数怎么办呢？