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

田间种植图的绘制有一个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.groundnuthead(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="一年多点试验")

作者介绍：

我是邓飞, 是一个读了玉米育种的研究生, 做了一个二代重测序的课题, 毕业之后进入了一家卖软件的公司, 开始了统计分析的培训, 现在做起了动物全基因组选择的数据分析工作. 

加我微信:

备注: 姓名+单位 微信(ID: yijiaobani)