str(mrbp1)
par(mfrow = c(2,1),mar = c(4,4,2,2))
plot(rep(mrbp1$trtmt,times = 3)+rep(c(-.3,0,.3),each = 18),c(mrbp1$spp1,mrbp1$spp2,
mrbp1$spp3),col = rep(c("blue","red","black"),each = nrow(mrbp1)),
pch = rep(c(19,16,17),each = nrow(mrbp1)),
xlab = "Treatment Group Indicator",cex = 1.5,ylab = "Dried Biomass",
main = "Dried Biomass for each species and treatment group")
abline(v = seq(from = 1.5,to = 5.5,by = 1),lty = "dashed",col = "grey")
legend(x = .65, y = 23, legend = c("Species 1","Species 2","Species 3"),
col = c("blue","red","black"),
pch = c(19,16,17),bg = "white")
par(mar = c(5,4,3,2))
plot(rep(mrbp1$block,times = 3)+rep(c(-.3,0,.3),each = 18),
c(mrbp1$spp1,mrbp1$spp2,mrbp1$spp3),col = rep(c("blue","red","black"),
each = nrow(mrbp1)),pch = rep(c(19,16,17),each = nrow(mrbp1)),
xlab = "Block Indicator",cex = 1.5,ylab = "Dried Biomass",
main = "Dried Biomass for each species and block",xaxp = c(1,3,2))
abline(v = seq(from = 1.5,to = 5.5,by = 1),lty = "dashed",col = "grey")
legend(x = .8, y = 24, legend = c("Species 1","Species 2","Species 3"),
col = c("blue","red","black"),
pch = c(19,16,17),bg = "white")
Run the code above in your browser using DataLab