农业试验设计中田间种植图的绘制方法

2019-06-13 20:29:30 浏览数 (2)

田间种植图的绘制有一个R包叫desplot功能很强大, 今天介绍一下这个包的使用方法

1. 安装desplot包

代码语言：javascript复制
# Install the released version from CRAN:
install.packages("desplot")

# Install the development version from GitHub:
# install.packages("devtools")
devtools::install_github("kwstat/desplot")

2. 功能介绍

  • 单因素方试验设计
  • 单因素随机区组试验设计
  • 单点多重复试验(RCBD)
  • 裂区试验设计
  • 一年多点试验设计
2.1 单因素试验设计
代码语言：javascript复制
library(agricolae)
cul = paste0("cul",1:10)
r= 3

试验描述: 10个品种, 重复3次, 共30个小区, 5行6列

代码语言：javascript复制
book = design.crd(cul,r)$book
book$row = rep(1:5,each=6)
book$col = rep(1:6,5)
代码语言：javascript复制
library(desplot)
desplot(cul~row col,data=book,main=("单因素试验设计田间种植图"))
2.2 单因素随机区组试验
代码语言：javascript复制
library(agricolae)
cul = paste0("cul",1:10)
r= 3

试验描述: 10个品种, 3个区组, 共30个小区, 5行6列. 随机区组和完全随机的区别在于一个区组地块是连在一起的.

代码语言：javascript复制
book = design.rcbd(cul,r)$book
book$row = rep(1:5,each=6)
book$col = rep(1:6,5)
代码语言：javascript复制
library(desplot)

desplot(cul~row col,data=book,out1=block,main=("单因素随机区组试验田间种植图"))
2.2 单地点多重复试验
代码语言：javascript复制
library(agridat)library(desplot)
data(ryder.groundnut)
gnut <- ryder.groundnut
代码语言：javascript复制
head(gnut)

block

row

col

gen

wet

dry

B1

4

1

F

3.8

2.3

B1

4

2

A

5.2

3.3

B1

4

3

D

3.0

1.8

B1

4

4

B

2.6

1.8

B1

4

5

E

4.2

2.8

B1

4

6

C

2.4

1.4

代码语言：javascript复制
desplot(wet~row col,text=gen,cex=1,data=gnut,main="单因素随机区组田间种植图")
2.4 裂区试验
代码语言：javascript复制
require(agridat)require(desplot)
desplot(block ~ col row, yates.oats,
        col=nitro, text=gen, cex=1,main="裂区试验设计种植图",
        out1=block, out2=gen)
2.5 一年多点试验
代码语言：javascript复制
data(nonnecke.sweetcorn.uniformity)
dat <- nonnecke.sweetcorn.uniformityrequire(desplot)
desplot(yield~col*row|loc, dat,main="一年多点试验")
block dry row

0 人点赞