--- mod1:mat_apoio:biogeo [2010/09/14 23:14] – paulo
+++ mod1:mat_apoio:biogeo [2024/01/11 15:21] (atual) – edição externa 127.0.0.1
@@ Linha 4: / Linha 4: @@
 {{:mod1:mat_apoio:arqui.jpg?300 |http://www.luontoon.fi/page.asp?Section=5353}}
-"A relação espécies-área é uma das poucas **leis** genuínas em ecologia" (Gotelli, 2007). O que você acha desta afirmação? Será que podem existir leis em biologia e especialmente em ecologia?
+"A relação espécies-área é uma das poucas **leis** genuínas em ecologia" (Gotelli, 2007). O que você acha desta afirmação? Será que podem existir leis em ecologia? Vamos conferir!
-Como toda lei que se preze, a lei da relação espécies-área é descrita por uma equação matemática:
+Que tal agora darmos um passeio pelo mundo "real"? Vamos ver alguns dados reais, coletados por pessoas reais em lugares reais. Então é tudo verdade, ok? Rode os comandos no R para ver os tais dados.
-<m>S=c A^z</m> , onde **S** é o número de espécies, **A** é a área e **c** e **z** são constantes. Esta é uma função de potência, mas podemos transformá-la em uma equação de reta:
-<m |>log(S)=log(c)+z log(A)</m> . Rode os comandos abaixo (adivinhem onde!) no **R**. Eles criarão uma função que vocês podem usar para explorar as propriedades das equações acima.
-<code>
-spp.area=function(c , z){
-	curve(expr = c*x^z , from=1, to=100, ylim=c(1,10000), xlab="Área",
-	 ylab="Número de espécies", font.lab=2, lwd=2, col=2, main=paste("c = ",c,"; z = ",z))
-	curve(expr = c*x^z , from=1, to=100, ylim=c(1,10000), xlab="Área",
-	 ylab="Número de espécies", font.lab=2, lwd=2, col=2, main=paste("c = ",c,"; z = ",z), log="xy")
-	}
-par(mfrow=c(2,2))
-spp.area(c = 1.1 , z = 2)
-</code>
-Note que os dois eixos do gráfico da direita estão em escala logarítmica. Note também que a função lhe permite variar os dois parâmetros da equação: **c** e **z**. Altere os valores desses parâmetros e veja o que acontece. Você consegue pensar em um significado biológico para esses dois parâmetros?
-Que tal agora darmos um passeio pelo mundo "real"? Vamos ver alguns dados reais, coletados por pessoas reais em lugares reais. Então é tudo verdade, ok? Rode os comandos para ver os tais dados.
 <code>
@@ Linha 40: / Linha 20: @@
 calif
-# Ajuste da reta
+# Gráfico
-mod=lm(log10(riqueza)~log10(area),data=calif)
+plot(riqueza~area,data=calif,xlab="Área",ylab="Número de espécies",main="Menina eu vou pra Califórnia...",
-mod
+ bty="l",font.lab=2)
+</code>
+A lei da relação espécies-área é empírica, ou seja, foi descrita a partir de dados como os acima. A relação acima pode ser descrita por uma função de potência:
+<m>S=c A^z</m> , onde **S** é o número de espécies, **A** é a área e **c** e **z** são constantes.
+Veja o gráfico dessa função:
+<code>
+x11()
+curve(expr = 2.39*x^.22 , from=0, to=25000, xlab="Área",
+	 ylab="Número de espécies", font.lab=2, lwd=2, col=2)
+</code>
+Agora observe o que acontece quando transformamos as escalas do gráfico em logaritmos:
+<code>
 # Gráfico
-par(mfrow=c(1,1))
+plot(riqueza~area,data=calif,log="xy",xlab="Área",ylab="Número de espécies",main=c("Menina eu vou pra Califórnia...",
-plot(riqueza~area,data=calif,log="xy",main=c("Menina eu vou pra Califórnia...",
  paste("c = ",round(mod[[1]][1],2),"; z = ",round(mod[[1]][2],2))),bty="l",font.lab=2)
+# Ajuste da reta
+mod=lm(log10(riqueza)~log10(area),data=calif)
+mod
 abline(mod,col=2)
 </code>
+Note que o que era uma função de potência, na escala logaritmica vira uma equação de reta:
+<m |>log(S)=log(c)+z log(A)</m> .
 Esses dados,  de espécies de plantas vasculares endêmicas em regiões da California foram coletados por Johnson e colaboradores e publicados em um trabalho em 1968 (ver referências no fim da página). Estão disponíveis neste [[http://math.hws.edu/~mitchell/SpeciesArea/index.html|site]] .
-Aí vai um desafio para vocês fazerem em casa: será que é possível recriar esses dados com os modelos vistos no exercício? Então se liga nos modelos que virão a seguir porque quem conseguir ganha uma jujuba! (Lembrem-se que vocês devem se empenhar pelo bem da ciência, não pelos prêmios.)
+Agora vamos brincar um pouco com a relação espécies-área variando os parâmetros **c** e **z**. Rode os comandos abaixo:
+<code>
+spp.area=function(c , z){
+	curve(expr = c*x^z , from=1, to=10^10, xlab="Área", ylim=c(1,500),
+	 ylab="Número de espécies", font.lab=2, lwd=2, col=2, main=paste("c = ",c,"; z = ",z))
+	curve(expr = c*x^z , from=1, to=10^10, xlab="Área", ylim=c(1,500),
+	 ylab="Número de espécies", font.lab=2, lwd=2, col=2, main=paste("c = ",c,"; z = ",z), log="xy")
+	}
+par(mfrow=c(2,2))
+spp.area(c = 1.5 , z = .25)
+spp.area(c = 2.1 , z = .25)
+</code>
+Note que os dois eixos do gráfico da direita estão em escala logarítmica. Note também que a função lhe permite variar os dois parâmetros da equação: **c** e **z**. Altere os valores desses parâmetros e veja o que acontece. Você consegue pensar em um significado biológico para esses dois parâmetros?
 ===== Modelo de amostragem passiva =====
 {{:mod1:mat_apoio:ilha.jpg?300 |http://www.forumfantastik.net/forum/threads/179909-Rochas-do-mundo/page2}}
@@ Linha 320: / Linha 340: @@
 ===== O equilíbrio de MacArthur & Wilson =====
 {{:mod1:mat_apoio:equilibrista.jpg?200 |http://www.1br.biz/s/papeis-de-parede/animais/ratinho-equilibrista-3171.html}}
-A idéia básica deste modelo é que o número de espécies de uma ilha é um balanço entre a //imigração// de novas espécies vindas do continente e a //extinção// das espécies presentes. Quando as taxas de imigração (<m>lambda</m>) e extinção (<m>mu</m>) se igualam, o número de espécies (**S**) atinge um equilíbrio dinâmico:
+A idéia básica deste modelo é que o número de espécies de uma ilha é um balanço entre a //imigração// de novas espécies vindas do continente e a //extinção// das espécies presentes. Quando as taxas de imigração ($$ \lambda $$) e extinção ($$\mu$$) se igualam, o número de espécies (**S**) atinge um equilíbrio dinâmico:
-$${dS}/dt=\lambda-mu$$
+$$\frac{dS}{dt} \ = \ \lambda-\mu$$
 O modelo assume que quanto mais espécies presentes na ilha, menor é a chegada de novas espécies do continente. Lembra de como vai ficando cada vez mais difícil encontrar figurinhas que não sejam repetidas? Pois é, o princípio é o mesmo. Neste modelo a função que descreve a taxa de imigração é:
@@ Linha 394: / Linha 414: @@
 <code>
-# aumentando apenas I
+# aumentando apenas E
 a=seq(from=.01,to=1,by=.01)
 for(i in 1:100){
@@ Linha 402: / Linha 422: @@
 </code>
 <code>
-# aumentando apenas E
+# aumentando apenas I
 for(i in 1:100){
 	grafeq(E=.5,I=a[i],P=100)
@@ Linha 526: / Linha 546: @@
 ===== EXTRAS =====
-<box 60% red|Atenção!> Só prossiga se você estiver muito adiantado. Se não fôr o caso vá primeiro para a seção "Teoria neutra da biodiversidade". Não fique triste, você ainda pode fazer o resto em casa, tipo na sexta ou no sábado à noite... Lembrando que o autor deste roteiro não se responsabiliza por eventuais danos causados à nota do seu questionário. ^_^</box>
+<box 60% red|Atenção!> Só prossiga se você estiver com tempo. Se não fôr o caso vá primeiro para a seção "Teoria neutra da biodiversidade". Não fique triste, você ainda pode fazer o resto em casa, tipo na sexta ou no sábado à noite... Lembrando que o autor deste roteiro não se responsabiliza por eventuais danos causados à nota do seu questionário. ^_^</box>
 ==== Complicando... ====
 Como diria o ditado, para quê simplificar se podemos complicar! Até agora consideramos o efeito da distância e do tamanho da ilha, só que a distância só tinha efeito sobre a taxa de imigração e o tamanho só tinha efeito na taxa de extinção. Mas será que distância não pode afetar também a taxa de extinção? E o tamanho da ilha não pode alterar a taxa de imigração? Veja a seguir.

mod1/mat_apoio/biogeo.1284506097.txt.gz · Última modificação: 2024/01/11 15:21 (edição externa)

Mostrar página Revisões anteriores

[unknown link type]Voltar ao topo