Ambos lados da revisão anteriorRevisão anteriorPróxima revisão | Revisão anterior |
mod1:mat_apoio:biogeo [2010/09/14 23:17] – paulo | mod1:mat_apoio:biogeo [2024/01/11 15:21] (atual) – edição externa 127.0.0.1 |
---|
{{:mod1:mat_apoio:arqui.jpg?300 |http://www.luontoon.fi/page.asp?Section=5353}} | {{: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> | <code> |
calif | 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 | # 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) | 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) | abline(mod,col=2) |
</code> | </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]] . | 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 ===== | ===== Modelo de amostragem passiva ===== |
{{:mod1:mat_apoio:ilha.jpg?300 |http://www.forumfantastik.net/forum/threads/179909-Rochas-do-mundo/page2}} | {{:mod1:mat_apoio:ilha.jpg?300 |http://www.forumfantastik.net/forum/threads/179909-Rochas-do-mundo/page2}} |
| |
===== EXTRAS ===== | ===== 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... ==== | ==== 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. | 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. |