Monise Terra Cerezini Exercícios 9 - Construção de Funções Simples #Exercício Programar 1.: QUE FRIO! celsius<- function(x) { frio<- 5/9*(x-32) cat ("\t","Converte Fahrenheit em Celsius\n") return(frio) } ****Testando celsius(10) temperatura=c(20,30,40,50) celsius(temperatura) #Exercício Programar 2. : Análise exploratória simultânea de duas variáveis explor<-function(x,y) { x11() par(mfrow=c(1,2)) plot(y~x) modelo<-lm(y~x) par(new=TRUE) abline(modelo,col="red") boxplot(y~x) par(mfrow=c(1,1)) cat("\t", "sumario x","\n") cat("\t", "sumario y","\n") cat("\t", "coeficiente","\n") resultados<-list(summary(x),summary(y), coef(modelo)) return(resultados) } dados.1<-c(15,18,22,24,25,30,31,34,37,39,41,45) dados.2<-c(2,4.5,6,10,30,34,50,56,60,77.5,80,85) explor(dados.1,dados.2) #Programar 3: Índices de Diversidade de Espécies Diversidade<-function(x) { Pi=(x[,2:ncol(x)]) / apply(x[,2:ncol(x)],2,sum,na.rm=T) H=-apply(pi*log(pi,base=exp(1)),2,sum,na.rm=T) D=apply(pi^2,2,sum,na.rm=T) return(list(colnames(x),H,D)) } Parc.1=abs(round(rnorm(15,4,3))) Parc.2=abs(round(rnorm(15,3,4))) Parc.3=abs(round(rnorm(15,6,3))) Parc.4=abs(round(rnorm(15,4.3,2.1))) Sp=rep("especie",75) tabela=data.frame(Sp,Parc.1,Parc.2,Parc.3,Parc.4) indices(tabela) ###Acabooooouuuuu!!!!!!!!!!!!!