### 9.2
# Criar uma função que retorne gráficos de análise exploratória para duas variáveis
# Dica deles, baixar pacote
install.packages("plotrix")
library("plotrix")
# e ver função multhist
help(multhist)

# Fazendo aos pouquinhos

exp.2v <- function(x,y) 
{
    if(length(x)/length(y)==1)
    {
    par(mfrow= c(2,3))
    plot(y~x, main="Plot")
    boxplot(x,y, main="Bloxplot")
    multhist(list(x,y), col=c("chartreuse4","darkblue"),xlab="Valores",ylab="Frequência",            main="Histograma")
    qqnorm(x, col="chartreuse4", xlab="Quantis teóricos", ylab="Quantis observados", 
           main="Q-Q plot")
    qqline(x)
    qqnorm(y, col="darkblue", xlab="Quantis teóricos", ylab="Quantis observados", 
           main="Q-Q plot")
    qqline(y)
    par(mfrow= c(1,1))
    summary(x) 
    summary(y) 
    correlacao <- cor(x,y) 
    resultado <- list(summary(x), summary(y), correlacao)
    return(resultado)
    } 
   
    else(length(x)/length(y)<=1 | length(x)/length(y)>=1)
    {
    stop("vetores devem conter o mesmo numero de observacoes")
    }
    
}


## Inventando x e y para testar
x1 <- seq(1,50, by=1)
y1 <- seq(1,100, by=2)
y2 <- seq(1,50, by=2)

exp.2v(x1,y1) #deve rodar
exp.2v(x1,y2) #deve dar erro