Ambos lados da revisão anteriorRevisão anteriorPróxima revisão | Revisão anterior |
mod1:mat_apoio:biogeo [2010/09/29 12:11] – marcel | mod1:mat_apoio:biogeo [2024/01/11 15:21] (atual) – edição externa 127.0.0.1 |
---|
plot(riqueza~area,data=calif,xlab="Área",ylab="Número de espécies",main="Menina eu vou pra Califórnia...", | plot(riqueza~area,data=calif,xlab="Área",ylab="Número de espécies",main="Menina eu vou pra Califórnia...", |
bty="l",font.lab=2) | 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> | </code> |
| |
Como toda lei que se preze, a lei da relação espécies-área é descrita por uma equação matemática: | |
<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> . | |
| |
Vamos ver então o que acontece quando transformamos as escalas do gráfico de espécies-área em escalas logaritmicas: | Agora observe o que acontece quando transformamos as escalas do gráfico em logaritmos: |
<code> | <code> |
# Ajuste da reta | |
mod=lm(log10(riqueza)~log10(area),data=calif) | |
mod | |
| |
# Gráfico | # Gráfico |
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",xlab="Área",ylab="Número de espécies",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]] . |
| |
Rode os comandos abaixo. Eles criarão uma função que vocês podem usar para explorar as propriedades das equações acima. | Agora vamos brincar um pouco com a relação espécies-área variando os parâmetros **c** e **z**. Rode os comandos abaixo: |
| |
<code> | <code> |