#Exercício Funções simples
conversor <- function(dados)
	{
	rest = dados - 32
	celsius = 5/9 * rest	
	resulta <- data.frame(dados,celsius)
	return(resulta)
	}
#Exercício Análise exploratória simultânea de duas variáveis
x <- round(runif(50,0,100))
y <- round(runif(50,0,100))
Plot <- function(x,y)
	{
	par(mfrow=c(3,2))
	boxplot(x,y, ylab= "Valores dos Grupos", xlab="Grupos", ylim=c(0,100), main="Valores das Amostras")
	hist(y, main="Histograma", xlab="Valores do Y", ylab= "Valores", prob=T)
	qqnorm(x,xlab="Quantis Teóricos", ylab="Quantis Amostrais", main="Valores do X")
	qqline(x, col="red")
	qqnorm(y,xlab="Quantis Teóricos", ylab="Quantis Amostrais", main="Valores do Y")
	qqline(y, col="green")
	plot(y~x, main= "Regresão", ylab="Dados y", xlab="Dados x")
  	lm.correlacao<-lm(y~x)
  	abline(lm.correlacao, col="blue")
  	Coef<-coef(lm.correlacao)
  	summary<-list (summary(x),summary(y),Coef)
 	return(summary)
	}
#Índices de Diversidade de Espécies
shannon <- function(abundancia)
	{
	abun1 = na.omit(abundancia)
	abun2 = abun1[abun1>0]
	pi = abun2/sum(abun2)
	ln.pi = log(pi)
	shannon = -(sum((pi)*ln.pi))
	return(shannon)
	}

simpson <- function(abundancia)
	{
	abund1 = na.omit(abundancia)
	abund2 = abund1[abund1>0]
	Pi = abund2/sum(abund2)
	simpson = sum(Pi^2)
	return(simpson)
	}

abundancia <- matrix(round(runif(30,0,100)), 5,6)
rownames(abundancia)=paste("especie", c(1:5))
colnames(abundancia)=paste("amostra", c(1:6))

diversidade <- function(abundancia,indices=c("shannon", "simpson"))
	{
	if(indices=="shannon")
  	{
    	shannon <- apply(abundancia,MARGIN=2,FUN=shannon)
    	indices=shannon
	return(shannon)
	} 
	else(indices=="simpson")
	{ 
	simpson <- apply(abundancia, MARGIN=2,FUN=simpson)
    	indices=simpson
  	return(simpson)
	}
  	}