getwd() setwd("C:/Users/Claire/Documents/RUSP") #---------------- 9.2 Exercício Análise exploratória simultânea de duas variáveis #Tarefa: Crie uma função para saída gráfica de análises exploratórias de duas variáveis. #Essa saída pode ter por exemplo boxplot, histograma, qq norm e y por x. #A função deve permitir a entrada de dois objetos vetores de mesmo tamanho (x e y). #Caso os valores de x e y forem desenhados em um mesmo gráfico, defina cores diferentes. #As legendas devem ser em português e a saída (return) deve apresentar o sumário das duas #variáveis e o coeficiente de correlação. ?rnorm graficos=function(x,y) { if(class(x)!= "numeric" & class(y)!= "numeric"){ # premissa1- se a classe dos vetores for diferente de numérico, a função não será executada stop("as classes dos vetores não são numéricos") } if(length(x)!= length(y)){ # premissa2- se o comprimento do vetor x for diferente do vetor y, a função não será executada. stop("os vetores tem diferententes tamanhos") } par(mfrow=c(3,2),las=1) #criando a áre de plotagem com duas linhas e colunas boxplot(x,y,border=c("blue","green"), outline = TRUE, names=c("Temperatura", "Atividade")) qqnorm(y, xlab = "Quantis teóricos", ylab = "Quantis amostrais") #plotando um gráfico qqnorm com a variável y qqline(y,col="blue") #inserindo a linha de normalidade em axul no qqnorm de y qqnorm(x,xlab="Quantis teóricos", ylab = "Quantis amostrais") #plotando um gráfico qqnorm com a variável x qqline(x,col="green") #inserindo a linha de normalidade em axul no qqnorm de x install.packages("plotrix") #instalando pacote para plotar dois histogramas em um mesmo gráfico library("plotrix") #carregando o conteúdo do pacote multhist(list(x,y),col=c("blue","green"),ylab="Frequência",main="", ylim=c(0,40)) #criando o gráfico com os histogramas das duas variáveis text(18,35,"Temperatura", cex=1) #inserindo legenda do histograma text(18,32,"Atividade",cex=1) #inserindo legenda do histograma symbols(15,35, squares = 0.1, inches =0.1, bg = "blue", add = TRUE) #inserindo símbolo para fazer referência entre cor e legenda symbols(15,32,squares=0.1,inches=0.1,bg="green",add=TRUE) #inserindo símbolo para fazer referência entre cor e legenda mxy=lm(y~x) #regressão linear entra as duas variáveis plot(y~x, xlab="Temperatura",ylab="Atividade",las=1) #gráfico y por x abline(lmxy, col="red", lty=2,lwd=2) #linha de tendência entre as duas variáveis par(mfrow=c(1,1)) #salvando as imagem como uma só. sumario.x=summary(x) #sumário da variável x sumario.y=summary(y) #sumário da variável y coeficientes=coef(lmxy) #coeficientes da linha da regressão linear return(list(sumario.x,sumario.y,coeficientes)) #dizendo ao R para retornar os resultados em formato de lista, já que ele não retorna mais de uma info por vez. } graficos(x,y) #TESTE #set.seed(21) #x=rnorm(100) #y=rnorm(100) #TESTE2 #x=rnorm(100) #y=c("a","b","c","d","e")