农业试验设计中田间种植图的绘制方法
田间种植图的绘制有一个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="一年多点试验")
本文参与 腾讯云自媒体分享计划，分享自微信公众号。
原始发表：2019-04-01，如有侵权请联系 cloudcommunity@tencent.com 删除
评论
登录后参与评论
推荐阅读
目录