写在前面的话,学习Python时,发现notebook很好用,又发现里面可以运行R代码。本身也正在学习markdown语法,就以此为契机,将cookbook里面的东西翻译一下。希望对后来者有用。

安装agricolae

如果没有安装的话,需要键入下面代码安装agricolae软件包:

install.packages("agricolae")

如果已经安装了,直接键入下面代码载入agricolae包:

library(agricolae)

随机区组试验设计

假设20个品种A1~A20,种植4个区组,怎么设计,代码如下

library(agricolae)
cul <- paste("A",1:20,"")
rcbd <- design.rcbd(trt = cul,r = 4)
rcbd$book
plots block cul
101 1 A 10
102 1 A 8
103 1 A 18
104 1 A 1
105 1 A 4
106 1 A 2
107 1 A 17
108 1 A 12
109 1 A 7
110 1 A 13
111 1 A 11
112 1 A 16
113 1 A 19
114 1 A 5
115 1 A 6
116 1 A 15
117 1 A 3
118 1 A 14
119 1 A 20
120 1 A 9
201 2 A 6
202 2 A 4
203 2 A 17
204 2 A 16
205 2 A 8
206 2 A 18
207 2 A 15
208 2 A 20
209 2 A 14
210 2 A 10
311 3 A 12
312 3 A 6
313 3 A 11
314 3 A 19
315 3 A 2
316 3 A 1
317 3 A 14
318 3 A 10
319 3 A 5
320 3 A 8
401 4 A 5
402 4 A 7
403 4 A 8
404 4 A 11
405 4 A 17
406 4 A 18
407 4 A 19
408 4 A 4
409 4 A 2
410 4 A 6
411 4 A 12
412 4 A 3
413 4 A 1
414 4 A 16
415 4 A 20
416 4 A 14
417 4 A 15
418 4 A 13
419 4 A 9
420 4 A 10

田间种植图如下:
其中,每一行表示一个区组,一个区组有A1~A20个品种

rcbd$sketch
A 10 A 8 A 18 A 1 A 4 A 2 A 17 A 12 A 7 A 13 A 11 A 16 A 19 A 5 A 6 A 15 A 3 A 14 A 20 A 9
A 6 A 4 A 17 A 16 A 8 A 18 A 15 A 20 A 14 A 10 A 3 A 1 A 7 A 19 A 9 A 12 A 5 A 13 A 2 A 11
A 15 A 9 A 4 A 3 A 7 A 17 A 16 A 13 A 18 A 20 A 12 A 6 A 11 A 19 A 2 A 1 A 14 A 10 A 5 A 8
A 5 A 7 A 8 A 11 A 17 A 18 A 19 A 4 A 2 A 6 A 12 A 3 A 1 A 16 A 20 A 14 A 15 A 13 A 9 A 10