Obrigado Walmes.
Esqueci de anexar o arquivo com a base de dados para rodar o script. Acabo de anexar nessa mensagem.
Tentei reproduzir o passo para salvar o arquivo em formato png, mas em todos os testes o arquivo saia em branco, com a informação de que é impossível de ser reproduzida a imagem.
Mais uma vez, agradeço pela ajuda!
rm(list = ls())
setwd('C:\\Users\\Yury\\Desktop\\Mestrado\\1_TESE\\Model_Data\\Maize\\ENSEMBLE\\Simulacoes\\PA\\50')
a = read.table('sdesiderio.txt', header = T, sep = '\t')
media = matrix(NA, nrow = 36, ncol = 4)
colnames(media) = c('Desv.Pad', 'Media', 'M.mais', 'M.menos')
for(i in 1:36){
d = sd(a[i,1:(ncol(a)-1)])
media[i,1] = d
m = rowMeans(a, na.rm = T, dims = 1)
media[i,2] = m[i]
media[i,3] = media[i,2]+media[i,1]
media[i,4] = media[i,2]-media[i,1]
ifelse(media[i,]<0, 0, media[i,])
}
a.box = t(a)
colnames(a.box) = c('J1','J2','J3','F1','F2','F3','M1','M2','M3','A1','A2','A3','M1','M2','M3',
'J1','J2','J3','J1','J2','J3','A1','A2','A3','S1','S2','S3','O1','O2','O3',
'N1','N2','N3','D1','D2','D3')
png("grafico-b1.png", width = 1200, height = 900, res = 196)
b1 = boxplot(a.box, las = 1, xlab = 'Semeadura', ylab = 'kg/ha', cex=0.5, pch=16)
title('Variação de Ya de acordo com a data de semeadura')
dev.off()