# Exercicio 9 - Construção de Função
# Camila Galheigo Coelho 

setwd("C:/Users/Camila/Documents/PESQUISA/1_DOUTORADO/Disciplinas/R")
getwd()

# Que frio! #
celsius.calculate<-function(x)
          {
        celsius=5/9*(x-32)
          return(celsius)
          }
celsius.calculate
celsius.calculate(99.9)


# analise exploratoria simultanea de duas variaveis #

?runif
?rpois
?qqline
?plot
variavel.1<-runif(100,min=10, max=100)
variavel.2<-rpois(100,25)

analise.expl<-function(x,y)
              {
               par(mfrow<-c(1,2))
               qqplot(x,y, xlab="variavel.1", ylab="variavel.2", main="analise exploratoria A")
               plot(x~y,xlab="variavel.1", ylab="variavel.2", main="analise exploratoria B")
               abline(lm(x~y))               
               descricao.variaveis<-list(summary(x), summary(y),cor(x,y))
               names(descricao.variaveis)<-c(substitute(x),substitute(y),"correlacao")
               return(descricao.variaveis)
              }

analise.expl(variavel.1,variavel.2)


## Indices de disversidade de especies ##

matx.1<-matrix( ,ncol=5,nrow=5,dimnames=list(c("sp.1","sp.2","sp.3","sp.4","sp.5"),c("area.1","area.2","area.3","area.4","area.5")))
matx.1
indice.shannon<-function(x)
         {
         if (na.rm==TRUE)
         { 
	    matriz=na.omit(x)
	    }
	     else 
	    {
		matriz=x 
	

 shannon=sum(p.i*ln(p.i))
    

           
save.image()