Dane

Nagrania video alczyków - zachowań i obecności w kolonii (Hornsund, Spitsbergen), sezon 2019; wczesny okres pisklęcy (pisklęta w wieku 2-7 dni) i środkowy okres pisklęcy (pisklęta w wieku 8-14 dni).

Ptaki znakowane indywidualnie za pomocą kombinacji kolorowych obrączek i znaków barwnych malowanych na piórach piersi. Kamery ustawione na wejściu do gniazda, w trybie poklatkowych, rejestracja ciągła przez 72h. Płeć ptaków oznaczana molekularnie.

Materiał video przetwarzany z pomocą programu VLC przez Antoine Grissot, Emilia Zalewska, Katarzyna Wojczulanis-Jakubas i Alan Szydłowski. Dokładność obserwacji sprawdzona za pomocą dwóch obserwatorów rejestrujących ten sam materiał dla kilku nagrań wideo - wysoka powtarzalność odczytu (>95%).

Liczba ptaków w kolonii

Ogólna obecność (wszystkie rodzaje aktywności połączone)

Obencość ptaków w kolonii z wyłączeniem aktywności związanej z żerowaniem/karmieniem

Obencość ptaków w kolonii związana z karmieniem (przyloty z pokarmem)

Dobowy rozkład lotów żerowiskowych - modalność

Modalność testowana za pomocą testu Silverman’a. Przeprowadzony osobno dla grup płciowych oraz etapów okresu pisklęcego. Na podstawie średniej liczby znakowanych ptaków, które wyleciały z kolonii, w kazdej godzinie doby (dwie/trzy doby obserwacji połączone).

(1) Wczesny okres pisklęcy_f (1) Wczesny okres pisklęcy_m (2) Środkowy okres pisklęcy_f (2) Środkowy okres pisklęcy_m
k=1 0.9619620 0.3543544 0.4514515 0.6736737
k=2 0.6596597 0.3003003 0.6006006 0.5965966
k=3 0.2992993 0.0420420 0.2832833 0.3663664

Interpretacja - wartości w tabeli zawierają wartość p dla testowanej hipotezy zerowej k modalnych (\(k \in <1,2,3>\)). Testy nie wykazują żadnej istotnej modalności dla żadnej z grup płciowych oraz etapów okresu pisklęcego.

Dobowy rozkład lotów żerowiskowych - różnice między rozkładami.

Różnice między rozkładami testowane za pomocą testów Kołmogorowa-Sminrova. Przeprowadzone parami dla płci i etapów okresu pisklęcego. Rozkłady utworzone na podstawie średniej liczby znakowanych ptaków, które wyleciały z kolonii, w kazdej godzinie doby (dwie/trzy doby obserwacji połączone).

loop_vecA loop_vecB ks_res
(1) Wczesny okres pisklęcy_f (1) Wczesny okres pisklęcy_m 0.99798
(1) Wczesny okres pisklęcy_f (2) Środkowy okres pisklęcy_f 0.44131
(1) Wczesny okres pisklęcy_f (2) Środkowy okres pisklęcy_m 0.25906
(1) Wczesny okres pisklęcy_m (2) Środkowy okres pisklęcy_f 0.17377
(1) Wczesny okres pisklęcy_m (2) Środkowy okres pisklęcy_m 0.04180
(2) Środkowy okres pisklęcy_f (2) Środkowy okres pisklęcy_m 0.89278

Interpretacja - wartości w tabeli zawierają wartość p dla testowanej hipotezy zerowej (brak różnic między rozkładami). Istotne różnice dla p <0.05

Actogramy

Osobne wykresy dla płci i etapów sezonu.

Interpretacja - zastosowanie raczej analityczne, nie do raportowania (co najwyżej w aneksie)

Modelowanie

Liczba znakowanych ptaków zaobserwowanych w poszczególnych godzinach doby (zmienna wyjaśniana; ogólna liczba bez uzwględniana poszczegolnych aktywności) w odniesieniu do płci, etapu sezonu i poszczególnych godzin - zmienne wyjaśniające; (godzina traktowana jako zmienna ciągła; nie do końca poprawne, ale z odpowiednimi zastrzeżeniami przy interpretacji dopuszczalne; czas/godziny to zmienne kołowe!). Wartości wyliczone na podstawie średniej liczby ptaków dla kilku (2/3) dni nagrań w każdej sesji.

Uogólniony model liniowy (GLM) z rozkładem Poisson’a dla błędów (nie uwzględnia czynnika losowego w postaci osobnika/gniazda, które występują wielokrotnie w zestawie danych - dane niekoniecznie niezależne - potencjalny wpływ pseudoreplikacji; ale można przyjąć, że efekt tej pseudoreplikacji jest niewielki bo liczba powtórzeń dla każdego osobnika/gniazda jest podobna). W modelu uwzględniona interakcja między etapem sezonu i płcią.


Call:
glm(formula = Nind ~ Session * Sx + Hr, family = "poisson", data = td_model)

Deviance Residuals: 
     Min        1Q    Median        3Q       Max  
-1.92355  -0.43503  -0.05721   0.52416   1.18936  

Coefficients:
                                         Estimate Std. Error z value Pr(>|z|)
(Intercept)                             2.126e+00  9.028e-02  23.550   <2e-16
Session(2) Środkowy okres pisklęcy      9.032e-02  9.757e-02   0.926    0.355
Sxm                                     1.128e-01  9.705e-02   1.162    0.245
Hr                                     -6.004e-05  4.901e-03  -0.012    0.990
Session(2) Środkowy okres pisklęcy:Sxm -9.925e-02  1.358e-01  -0.731    0.465
                                          
(Intercept)                            ***
Session(2) Środkowy okres pisklęcy        
Sxm                                       
Hr                                        
Session(2) Środkowy okres pisklęcy:Sxm    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for poisson family taken to be 1)

    Null deviance: 43.265  on 95  degrees of freedom
Residual deviance: 41.559  on 91  degrees of freedom
AIC: 438.81

Number of Fisher Scoring iterations: 4

Interpretacja - żaden z testowanych paramterów nie wydaje się mieć większego wpływu na obserowaną liczbę ptaków. Prawdopodobnie wielkość próby jest za mała aby ujawnić statystycznie istotne różnice, które obserwuje się na wykresie.

Statystyki kołowe - jednorodność rozkładu

Zmienna rozważana jest teraz godziną, w której ptak odleciał na żerowisko. Modalność analizowana za pomocą testu Rayleigh’a. Analiza wykonana osobno dla płci i etapów okresu pisklęcego.

c_gr ray_p_val
(1) Wczesny okres pisklęcy_f 0.0397453
(1) Wczesny okres pisklęcy_m 0.0402743
(2) Środkowy okres pisklęcy_f 0.1054384
(2) Środkowy okres pisklęcy_m 0.0786227

Interpretacja: Wartości w tabeli zawierają wartość p dla testowanej hipotezy zerowej o jednorodności rozkładu. Wczesny okres pisklęcy jest bardziej niejednorodny.

Statystyki kołowe - porównanie rozkładów

Zmienne zgodne z opisem powyżej dla jednorodności rozkładu. Analiza porównawcza rozkładów wykonana za pomocą testów Watson’a.

loop_vecA loop_vecB wats_stat Sigf
(1) Wczesny okres pisklęcy_f (1) Wczesny okres pisklęcy_m 0.3416375 0.05
(1) Wczesny okres pisklęcy_f (2) Środkowy okres pisklęcy_f 0.0780274 ns.
(1) Wczesny okres pisklęcy_f (2) Środkowy okres pisklęcy_m 0.0629210 ns.
(1) Wczesny okres pisklęcy_m (2) Środkowy okres pisklęcy_f 0.2561469 0.05
(1) Wczesny okres pisklęcy_m (2) Środkowy okres pisklęcy_m 0.2874610 0.05
(2) Środkowy okres pisklęcy_f (2) Środkowy okres pisklęcy_m 0.0389027 ns.

Interpretacja: Istotne różnice w rozkładzie dla płci we wczesnym okresie pisklecym i dla samców miedzy okresami. Co jest dość zgodne z tym co widać na wykresie.

W sumie analizowanie godziny jako zmiennej, daje mocniejsze wnioskowanie niż analizowanie liczby osobników w danej godzine (pierwsza część raportu).

```