install.packages("plotrix") 
library(plotrix)
x<-c(14,34,40,60,80)
y<-c(150,100,120,130,200)

ana.exp=function(x,y,rmna=TRUE)
{
  if(length(x)!=length((y)))
  { cat("\t", " os vetores não possuem mesmo tamanho")}
  if(rmna==TRUE)
  {x.1<-na.omit(x)##ESTUDAR NA.OMIT
  x.1
  y.1<-na.omit(y)
  y.1}
 
  par(mfrow = c(2,2)) ##passo1
  lista<-list(x,y) ##so para multilist
  ##histograma
  multhist(lista)
  ##qqnorm para duas variaveis em um grafico só ##PASSO2
  (a<-min(x.1,y.1))
  (b<-max(x.1,y.1))
  qqnorm(x, col="red",pch=16,xlab="quantis", ylab="quantis",ylim=c(a,b))
  par(new=TRUE)
  qqnorm(y, xlab="",ylab="")
  plot(y~x, xlab= "valores de x", ylab="valores do y",pch=17, col="blue") 
  boxplot(x,y,col=c("red","blue"))
  sumx<-summary(x.1)
  sumy<-summary(y.1)
  coef.co<-cor(x.1,y.1)
  taboa=list(sumx,sumy,coef.co)
  names(taboa)=c("Sumário de x", "Sumário de y", "Coef.correlação")
  taboa
  
  return(taboa)
 
}
ana.exp(x,y,rmna=T)

##tive dificuldades para finalizar a tarefa.