1.请根据表1数据,以灌溉模式为主区,生物炭用量为子区,利用R语言agricolae包生成裂区试验设计试验处理布置图;针对裂区试验设计进行方差分析和Turkey HSD进行均值对比,并利用ggplot2工具包绘制主效应和交互效应图。请将布置图、R代码(包括注释)、方差分析表、均值对比结果和主效应图以及交互效应图粘贴到答题册上。(25分)
#### 导入相关程序包
library(broom)
library(ggpmisc)
library(readxl)
library(minpack.lm