Informacija

Varijansa stohastičke stope rasta za prosječnu veličinu stanovništva tijekom vremena

Varijansa stohastičke stope rasta za prosječnu veličinu stanovništva tijekom vremena


We are searching data for your request:

Forums and discussions:
Manuals and reference books:
Data from registers:
Wait the end of the search in all databases.
Upon completion, a link will appear to access the found materials.

Zanima me izračunavanje varijance stohastičke stope rasta za prosječnu veličinu populacije tijekom vremena,
$$ ln beta = lim_ {t to infty} dfrac {1} {t} ln V [N (t)] $$ što pokazuje da varijacija $ N (t) $ raste pri stopa $ beta $. Jednadžba za izračunavanje navedena je u knjizi "Matrični populacijski modeli" Hal Caswella. Varijansa se može izračunati kao, $$ beta = lambda_1^{( bf {B} _2)} $$ gdje je $ lambda_1 $ dominantna vlastita vrijednost, $$ begin {align*} bf {B_2} & = bf {F_2} [ bf {P} otimes ( bf {I} otimes bf {I})] mathrm {,} [1em] bf {F} _2 & = mathrm {diag} ({ bf A_1 otimes A_1, ldots, A_k otimes A_k}) mathrm {,} end {align*} $$ $ bf A $ matrice su projekcije populacije dimenzije $ s $ iz $ k $ neovisna okruženja, $ bf P $ je matrica prijelaza stupca-stohastike dimenzije $ k $, a $ bf I $ je matrica identiteta dimenzije $ s $. Pokušao sam to učiniti u R -u koristeći primjer goldenheather (Hudsonia) iz knjige Caswell koji se može učitati u R na sljedeći način,

knjižnica (popbio); knjižnica (Matrix); podaci (hudsonia)

Završio sam s vrijednošću $ beta $ od 0,9366351. To ne može biti točno, varijacija je prevelika. Mora da sam pogriješio ili je možda došlo do pogreške u jednoj od jednadžbi (vrlo, vrlo vjerojatno prvoj).

Evo mog koda (ako pomaže),

s = NROW (hudsonia [[1]]) k = length (hudsonia) P = Matrix (1 / k, k, k) I = Diagonal (s, 1) F2 = bdiag (lapply (hudsonia, function (i) kronecker (i, i))) I2 = kronecker (I, I) P2 = kronecker (P, I2) B2 = F2 %* % P2 lambda (B2)



Komentari:

  1. Tawil

    Po mom mišljenju, ovo - zabuna.

  2. Daim

    I je li učinkovit?

  3. Esau

    Upravo ono što trebate. Zajedno možemo doći na pravi odgovor. Siguran sam.



Napišite poruku