田间种植图的绘制有一个R包叫desplot
功能很强大, 今天介绍一下这个包的使用方法
# Install the released version from CRAN:
install.packages("desplot")
# Install the development version from GitHub:
# install.packages("devtools")
devtools::install_github("kwstat/desplot")
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=("单因素试验设计田间种植图"))
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=("单因素随机区组试验田间种植图"))
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="单因素随机区组田间种植图")
require(agridat)require(desplot)
desplot(block ~ col+row, yates.oats,
col=nitro, text=gen, cex=1,main="裂区试验设计种植图",
out1=block, out2=gen)
data(nonnecke.sweetcorn.uniformity)
dat <- nonnecke.sweetcorn.uniformityrequire(desplot)
desplot(yield~col*row|loc, dat,main="一年多点试验")