欢迎您访问 最编程 本站为您分享编程语言代码,编程技术文章!
您现在的位置是: 首页

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

最编程 2024-04-27 14:30:47
...

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

1. 安装desplot包

# 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 单因素试验设计

library(agricolae)
cul = paste0("cul",1:10)
r= 3

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

book = design.crd(cul,r)$book
book$row = rep(1:5,each=6)
book$col = rep(1:6,5)
library(desplot)
desplot(cul~row+col,data=book,main=("单因素试验设计田间种植图"))

2.2 单因素随机区组试验

library(agricolae)
cul = paste0("cul",1:10)
r= 3

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

book = design.rcbd(cul,r)$book
book$row = rep(1:5,each=6)
book$col = rep(1:6,5)
library(desplot)

desplot(cul~row+col,data=book,out1=block,main=("单因素随机区组试验田间种植图"))

2.2 单地点多重复试验

library(agridat)library(desplot)
data(ryder.groundnut)
gnut <- ryder.groundnut
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

desplot(wet~row+col,text=gen,cex=1,data=gnut,main="单因素随机区组田间种植图")

2.4 裂区试验

require(agridat)require(desplot)
desplot(block ~ col+row, yates.oats,
        col=nitro, text=gen, cex=1,main="裂区试验设计种植图",
        out1=block, out2=gen)

2.5 一年多点试验

data(nonnecke.sweetcorn.uniformity)
dat <- nonnecke.sweetcorn.uniformityrequire(desplot)
desplot(yield~col*row|loc, dat,main="一年多点试验")

推荐阅读