Päivä: 25
Pyydettyjä smoltteja: 38
Merkattu: 37
Takaisinpyydettyjä merkittyjä smoltteja: 0
Seuraavassa esitellään tuoreella havaintoaineistolla päivitetyt päivän 25 (3.5.) tulokset käytössä olevista kolmesta eri mallista missä muuttointensiteettiä on kuvattua 1) log-normaali-, 2) normaali- ja 3) tasajakaumilla. Taustat N:n prioriin löytyvät aikaisemmasta tulosbloggauksesta ja lähtokohdat malleihin tästä.
Merkittävin muutos seuranta-asetelmissa, ja tarkemmin kokonaispyyntitehon määrittelemisessä, on Smolttirysän poistaminen käytöstä päivästä 21 eteenpäin. Eli tällä hetkellä pyynnissä on taas ainoastaan ruuvi. Tämän osalta malleja on päivitetty seuraavasti:
lambda[i] <- q_s*n[i]*step(15-i) + q*n[i]*step(i-16) + q_s*n[i]*step(i-22) – q*n[i]*step(i-22)
Malli 1 – “log-normaali“
Log-normaalimalli lähtemisintensiteetille tuottaa taimenten kokonaismäärän odotusarvon piste-estimaatiksi 1411 kappaletta, jolle 95% luottamusväli on 402–3722 (Taulukko 1). Todennäköisin yksittäinen arvio kokonaismäärästä on (posteriorijakauman moodi) on jossain 500 ja 1000 välissä (Kuva 1).
Maanantaina 4.5. vaeltavien taimensmolttien lukumäärän (n[26]) odotusarvon piste-estimaatti on 33 ja 95 % luottamusväli 7–90 (Taulukko 1). Todennäköisin yksittäinen arvio (posteriorijakauman moodi) vaeltavien smolttien määrästä on 16 (Kuva 2).
Maanantaina 4.5. ruuviin jäävien taimensmolttien lukumäärän (x[26]) odotusarvon piste-estimaatti on 1 ja 95 % luottamusväli 0–4 (Taulukko 1). Todennäköisin yksittäinen arvio (posteriorijakauman moodi) on, että ruuviin ei jää yhtään taimensmolttia (Kuva 3).
Kokonaispyydystettävyyden (q) estimaatti on mallissa 0,073, eli 7,3 % pyydystettävissä olevista kaloista joutuu pyydykseen (Taulukko 1). Toistaiseksi pyydykset eivät ole onnistuneet uudelleenpyydystämään ainuttakaan 37 merkatusta kalasta.
Kuva 1. Malli 1:n posteriorijakauma vaeltavien taimensmolttien kokonaismäärästä 2015.
Kuva 2. Malli 1:n posteriorijakauma vaeltavien taimensmolttien määrästä maanantaina 4.5.2015
Kuva 3. Malli 1:n posteriorijakauma ruuviin jäävien taimensmolttien määrästä maanantaina 4.5.2015
Taulukko 1. Parametrien N, n[26], x[26] ja q posteriorijakaumien tunnusluvut malli 1:ssa.
mean | sd | MC_error | val2.5pc | median | val97.5pc | start | sample | |
N | 1411.0 | 873.2 | 19.46 | 402.0 | 1177.0 | 3722.0 | 51 | 161900 |
n[26] | 32.53 | 21.87 | 0.4414 | 7.0 | 27.0 | 90.0 | 51 | 161900 |
x[26] | 0.9503 | 1.053 | 0.00358 | 0.0 | 1.0 | 4.0 | 51 | 161900 |
q | 0.07326 | 0.03892 | 6.637E-4 | 0.02079 | 0.06546 | 0.1682 | 51 | 161900 |
Malli 2 – “normaali”
Normaalijakauma lähtemisintensiteetille tuottaa taimenten kokonaismäärän odotusarvon piste-estimaatiksi 2359 kappaletta, jolle 95% luottamusväli on 707–5075 (Taulukko 2). Tulos on selvästi korkeampi kuin log-normaalin mallin saman päivän tulos (1411) ja kuin saman mallin edellinen tulos (1199). Todennäköisin yksittäinen arvio kokonaismäärästä on (posteriorijakauman moodi) on jossain 1500 ja 2000 välissä (Kuva 4).
Maanantaina 4.5. vaeltavien taimensmolttien lukumäärän (n[26]) odotusarvon piste-estimaatti on 68 (95 % luottamusväli 20–150, Taulukko 2), joka on yli kaksinkertainen log-normaalin mallin vastaavaan estimaattiin (33) verrattuna. Todennäköisin yksittäinen arvio (posteriorijakauman moodi) vaeltavien smolttien määrästä on 31–32 (Kuva 5).
Maanantaina 4.5. ruuviin jäävien taimensmolttien lukumäärän (x[26]) odotusarvon piste-estimaatti on 2 ja 95 % luottamusväli 0–6 (Taulukko 2). Todennäköisin yksittäinen arvio (posteriorijakauman moodi) ruuviin jäävien smolttien määrästä on 1 (Kuva 6).
Kokonaispyydystettävyyden (q) estimaatti on mallissa 0,074, eli 7,4 % pyydystettävissä olevista kaloista joutuu pyydykseen (Taulukko 2).
Kuva 4. Malli 2:n posteriorijakauma vaeltavien taimensmolttien kokonaismäärästä 2015.
Kuva 5. Malli 2:n posteriorijakauma vaeltavien taimensmolttien määrästä maanantaina 4.5.2015
Kuva 6. Malli 2:n posteriorijakauma ruuviin jäävien taimensmolttien määrästä maanantaina 4.5.2015
Taulukko 2. Parametrien N, n[26], x[26] ja q posteriorijakaumien tunnusluvut malli 2:ssa.
mean | sd | MC_error | val2.5pc | median | val97.5pc | start | sample | |
N | 2359.0 | 1154.0 | 19.21 | 707.0 | 2163.0 | 5075.0 | 51 | 125850 |
n[26] | 67.57 | 33.88 | 0.5517 | 20.0 | 61.0 | 150.0 | 51 | 125850 |
x[26] | 2.292 | 1.716 | 0.007054 | 0.0 | 2.0 | 6.0 | 51 | 125850 |
q | 0.07382 | 0.03459 | 4.894E-4 | 0.02592 | 0.06734 | 0.1582 | 51 | 125850 |
Malli 3 – “tasa”
Tasajakauma lähtemisintensiteetille tuottaa taimenten kokonaismäärän odotusarvon piste-estimaatiksi 2506 kappaletta, jolle 95% luottamusväli on 707–5075 (Taulukko 3). Tulos on korkeampi kuin normaali-mallin saman päivän tulos (2359) ja kuin saman mallin edellinen tulos (1578). Todennäköisin yksittäinen arvio kokonaismäärästä on (posteriorijakauman moodi) on jossain 1500 ja 2000 välissä (Kuva 7).
Maanantaina 4.5. vaeltavien taimensmolttien lukumäärän (n[26]) odotusarvon piste-estimaatti on 42 (95 % luottamusväli 13–90, Taulukko 3), joka on pienempi kuin normaali-mallin tulos (68) mutta suurempi kuin log-normaalin mallin vastaava estimaatti (33). Todennäköisin yksittäinen arvio (posteriorijakauman moodi) vaeltavien smolttien määrästä on 31–32 (Kuva 8).
Maanantaina 4.5. ruuviin jäävien taimensmolttien lukumäärän (x[26]) odotusarvon piste-estimaatti on 1 ja 95 % luottamusväli 0–4 (Taulukko 2). Todennäköisin yksittäinen arvio (posteriorijakauman moodi) ruuviin jäävien smolttien määrästä on 1 (Kuva 9).
Kokonaispyydystettävyyden (q) estimaatti on mallissa 0,082, eli 8,2 % pyydystettävissä olevista kaloista joutuu pyydykseen (Taulukko 2). Tämän estimaatti on korkeampi kuin mallin 1 ja 2 vastaavat tulokset.
Kuva 7. Malli 3:n posteriorijakauma vaeltavien taimensmolttien kokonaismäärästä 2015.
Kuva 8. Malli 3:n posteriorijakauma vaeltavien taimensmolttien määrästä maanantaina 4.5.2015
Kuva 9. Malli 3:n posteriorijakauma ruuviin jäävien taimensmolttien määrästä maanantaina 4.5.2015.
Taulukko 3. Parametrien N, n[26], x[26] ja q posteriorijakaumien tunnusluvut malli 3:ssa.
mean | sd | MC_error | val2.5pc | median | val97.5pc | start | sample | |
N | 2506.0 | 1134.0 | 11.63 | 902.0 | 2302.0 | 5232.0 | 10000 | 3420903 |
n[26] | 41.78 | 19.97 | 0.1914 | 13.0 | 38.0 | 90.0 | 10000 | 3420903 |
x[26] | 1.109 | 1.096 | 6.281E-4 | 0.0 | 1.0 | 4.0 | 10001 | 3150900 |
q | 0.08204 | 0.03619 | 2.966E-4 | 0.03185 | 0.07518 | 0.1702 | 10000 | 3420903 |
Seuraavissa malliajoissa tulisi tarkemmin kiinnittää huomiota simulaatioketjujen konvergoitumiseen, autokorrelaatioon ja MC-error:iin, ja näiden seikkojen tarkastelemiseen. Sen pohjalta tulisi tehdä tarvittavat toimenpiteet (burn-in, thin, iteraatioiden määrän lisääminen) tulosten laadun parantamiseksi.