Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anteriorRevisão anteriorPróxima revisão | Revisão anterior | ||
exercicios:exerc2 [2012/05/16 14:18] – [Estocasticidade Ambiental] adalardo | exercicios:exerc2 [2024/01/09 18:18] (atual) – edição externa 127.0.0.1 | ||
---|---|---|---|
Linha 284: | Linha 284: | ||
===== Estocasticidade Ambiental ===== | ===== Estocasticidade Ambiental ===== | ||
- | Flutuações ambientais podem exercer efeito na taxa de crescimento instantâneo da população. De uma forma simples, podemos imaginar que essa variação funcione como um ruído no //r//, como se a população em média tivesse uma taxa, mas a cada realização ela pudesse ser um tanto diferente devido a condições externar a ela própria. A realização | + | Flutuações ambientais podem exercer efeito na taxa de crescimento instantâneo da população. De uma forma simples, podemos imaginar que essa variação funcione como um ruído no //r//, como se a população em média tivesse uma taxa, mas a cada realização ela pudesse ser um tanto diferente devido a condições externar a ela própria. A implementação |
+ | Para um crescimento discreto a construção de simulações com estocasticidade ambiental é mais intuitivo: | ||
< | < | ||
- | lamb = rnorm(10, mean=1.2, sd=0.8) | + | npop=10 |
- | N0=rep(10,10) | + | n0=10 |
+ | lamb.med = 1.2 | ||
+ | lamb.sd= 0.4 | ||
+ | lamb = rnorm(npop, mean=lamb.med, sd=lamb.sd) | ||
+ | N0=rep(n0,npop) | ||
N1=lamb*N0 | N1=lamb*N0 | ||
- | lamb=rnorm(10, mean=1.2, sd=0.8) | + | lamb=rnorm(npop, mean=lamb.med, sd=lamb.sd) |
N2=N1*lamb | N2=N1*lamb | ||
- | N3=N2*rnorm(10,mean=0.3,sd=0.05) | + | N3=N2*rnorm(npop,mean=lamb.med,sd=lamb.sd) |
- | N4=N3*exp(rnorm(10, | + | N4=N3*rnorm(10, |
- | N5=N4*exp(rnorm(10, | + | N5=N4*rnorm(10, |
Nt< | Nt< | ||
- | matplot(0: | + | matplot(0: |
- | lines(1:6, exp(1: | + | |
</ | </ | ||
==== Desafio ==== | ==== Desafio ==== | ||
- | É possível adaptar a nossas | + | É possível adaptar a nossas |
<box 70% green |Dicas> | <box 70% green |Dicas> | ||
O primeiro passo sempre e pensar quais argumentos vamos precisar | O primeiro passo sempre e pensar quais argumentos vamos precisar | ||
- | Nesse caso, temos apenas mais um argumento o **//dpr//** : o desvio padrão de //r//. O resto continua o mesmo, lembre-se que se o **//dpr//** for 0, nosso população é determinística! Ou seja, a mesma função pode se prestar para simular ambos cenários. | + | Nesse caso, temos apenas mais um argumento o **//lamb.dp//** : o desvio padrão de //lambda//. O resto continua o mesmo, lembre-se que se o **//lamb.dp//** for 0, nosso população é determinística! Ou seja, a mesma função pode se prestar para simular ambos cenários. |