Packages

#install.packages(“AgroR”) #install.packages(‘agricolae’) #installed.packages(“devtool”) #install.packages(“readxl”) #install.packages(“dplyr”) #install.packages(“tidyverse”) #installed.packages(“ggplot2”) #install.packages(“gridExtra”) #install.packages(“ScottKnott”) #install.packages(“ggsignif”) #install.packages(“sjstats”) #install.packages(“mvtnorm”) #install.packages(‘knitr’) #install.packages(“cowplot”) #install.packages(“lme4”)

Plots

library(lme4) library(cowplot) library(agricolae) library(devtool) library(readxl) library(dplyr) library(tidyverse) library(ggplot2) library(AgroR) library(gridExtra) library(ScottKnott) library(mvtnorm) library(ggsignif) library(sjstats) library(gridExtra)

RCBD DATA SET

attach(Bio_stimulant_clean_data_1) a<-Bio_stimulant_clean_data_1 a\(Treatment<-as.character(a\)Treatment) a\(Rep<-as.factor(a\)Rep) a

attach(Biostimulant_Trial2_Data) b<-Biostimulant_Trial2_Data b\(Treatment<-as.character(a\)Treatment) b\(Rep<-as.factor(a\)Rep) b

MIXED MODEL

model <- lmer(Yield ~ Treatment + (1 | Yield), data = a) summary(model)

Diagnostic plots

plot(model)

Residuals vs Fitted values

plot(residuals(model), fitted(model))

QQ plot for residuals

qqnorm(residuals(model)) qqline(residuals(model))

##ANOVA RCBD

bb<-with(a,DBC(Treatment,Rep,Yield, mcomp = “duncan”, ylab = “Yield”, xlab = “Treatment”))# duncan

aa<-with(a,DBC(Treatment,Rep,Yield, mcomp = “sk”, ylab = “Yield”, xlab = “Treatment”))# duncan

LATIN SQUARE

#modelLatin <- lmer(Yield ~ Treatment + (1 | ID) + (1 | Rep), data = b) #summary(modelLatin)