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:14] – [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. | + | 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 |
- | Nesse caso teríamos um //re//: o r estocástico. A realização | + | Para um crescimento discreto a construção de simulações com estocasticidade ambiental é mais intuitivo: |
< | < | ||
- | rd = rnorm(10, mean=0.3, sd=0.05) | + | npop=10 |
- | N0=rep(10,10) | + | n0=10 |
- | N1=N0*exp(re) | + | lamb.med = 1.2 |
- | re=rnorm(10, mean=0.3, sd=0.05) | + | lamb.sd= 0.4 |
- | N2=N1*exp(re) | + | lamb = rnorm(npop, mean=lamb.med, sd=lamb.sd) |
- | N3=N2*exp(rnorm(10,mean=0.3,sd=0.05)) | + | N0=rep(n0,npop) |
- | N4=N3*exp(rnorm(10, | + | N1=lamb*N0 |
- | N5=N4*exp(rnorm(10, | + | lamb=rnorm(npop, mean=lamb.med, sd=lamb.sd) |
+ | N2=N1*lamb | ||
+ | N3=N2*rnorm(npop,mean=lamb.med,sd=lamb.sd) | ||
+ | N4=N3*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. |