descritiva <- function(dados,tit.grafico){ limites <- c(1.04,1.49) dad = dados par(mfrow=c(2,3)) temp1 <- hist(dad$vari, plot=FALSE)[[3]] temp2 <- density(dad$vari, na.rm=T) hist(dad$vari, prob=TRUE, ylab="Frequencia", xlab=tit.grafico, main="A", ylim=c(0, max(temp1, temp2[[2]]))) lines(temp2) ## Boxplot comparativo por ponto de coleta boxplot(dad$vari ~ dad$LOCAL,xlab="Pontos de coleta",ylab=tit.grafico,ylim=c(-0.5,2.6),main="B",notch=TRUE,outline=FALSE) abline(h=limites) boxplot(dad$vari ~ dad$TRIM,xlab="Trimestre",ylab=tit.grafico,ylim=c(-0.5,2.6),main="C",notch=TRUE,outline=FALSE) abline(h=limites) if(as.numeric(dim(dad)[1]) > 900){ boxplot(dad$vari ~ dad$GRUPO,xlab="Grupos",ylab=tit.grafico,ylim=c(-0.5,2.6)) abline(h=limites)} else { boxplot(dad$vari ~ dad$Ano, xlab= "Ano", ylab=tit.grafico,ylim=c(-0.5,2.6),main="D",notch=TRUE,outline=FALSE) abline(h=limites)} boxplot(dad$vari ~ dad$UHE,xlab="",ylab=tit.grafico,xaxt="n",ylim=c(-0.5,2.6),main="E",notch=TRUE,outline=FALSE) labels <- levels(dad$UHE) abline(h=limites) text(1:length(labels), par("usr")[3] - 0.05, srt = 90, adj = 1, labels = labels, xpd = TRUE, cex = 0.9) teste <- data.frame(TRIM = dad$TRIM,LOCAL = dad$LOCAL,vari = dad$vari) teste <- na.exclude(teste) interaction.plot(teste$TRIM, teste$LOCAL, teste$vari, ylab = "medias", xlab = "Trimestre", legend=F, lty=1:3, ylim=c(1, 1.6),main="F") legend("topright", c("M","R","J"), lty=1:3,cex=0.7) abline(h=limites) return(invisible())}