Import dat: výsledkem je, že v CANOCO budete mít tři matice, dvě typu compositional (molluscs
1) a vascular
) a jednu typu general (matice env
). Všechny matice musí obsahovat 43 řádků (ploch).
Analýza 1: sestává ze dvou kroků, a není na ni možné použít přednastavenou šablonu. Při vytvoření této analýzy postupujte takto:
vascular
), None
a Center and standardize by samples Hellinger standardization
.moluscs
), Explanatory data - sem je třeba vložit skóre na ordinačních osách z předchozí analýzy, a to takto: klikněte na tlačítko Advanced role v sekci Explanatory data, v novém okně na tlačítko New, v okně Use score results from: vyberte Current analysis
a Step 1
, vpravo v boxu Select score type: vyberte Case scores
a pod ním Use these scores from vyberte all axes
. Klikněte na OK, a v políčku pro Explanatory variables se objeví Case scores from step 1
, což je přesně to, co potřebujeme. Log
a Center and standardize by samples None
(na plžích není třeba dělat Hellingerovu standardizaci, je ale potřeba je log ztransformovat, protože jde o počty).Forward selection
, v možnosti Permutation Test Parameters vyberte Unrestricted permutations
.CaseR.1
a CaseR.2
), je ale možné, že vám druhá osa už nevyjde signifikantní (s Bonferroniho korekcí).Analysis 'Analysis 1', step 'Step 2' Forward Selection Results: Name Explains % Contribution % pseudo-F P P(adj) CaseR.1 25.4 56.4 14.0 0.002 0.016 CaseR.2 4.5 10.0 2.6 0.006 0.048
env
, a to takto: Cases (1)
2) označte sloupce s osami, které byly vybrány postupným výběrem v analýze Step 2 (v mém případě CaseR.1
a CaseR.2
), a nakopírujte je na clipboard (CTRL+C)env
, posuňte se na poslední sloupeček, pravým tlačítkem myši klikněte na záhlaví prvního prázdného sloupce ([C15]) a vyberte Append ... from Clipboard
- vloží se vám sloupečky se skóre na vegetačních ordinačních osách.
Analýza 2 - Analysis > Add new analysis > Canoco adviser, vyberte matice molluscs
a env
, analýzu Interactive forward selection ...
.
molluscs
a jako Explanatory data vyberte matici env
a v ní všechny proměnné kromě skóre na vegetačních osách z předchozí analýzy (chceme dělat výběr jen mezi měřenými proměnnými prostředí). Ca
) a marginálně signifikantní i konduktivita (cond
) - kterou jsem nakonec zahrnul do modelu, i když její P(adj) byl mírně vyšší než 0.05.Analysis 'Interactive-forward-selection', step 'Forward Selection' Forward Selection Results: Name Explains % Contribution % pseudo-F P P(adj) Ca 20.2 41.5 10.4 0.002 0.028 conduct 5.3 10.9 2.8 0.004 0.056
Analýza 3 - vlastní rozklad variance: Analysis > Add new analysis > Canoco adviser, vyberte matice molluscs
a env
a analýzu v sekci Variation Partitioning Analyses
, konkrétně s testováním parciálních (conditional efektů (Var-part-2groups-Conditional-effects-tested...
).
molluscs
, Explanatory data tabulka env
(která nyní obsahuje jak vlastní měřené proměnné, tak i ordinační skóre z vegetačních dat)First group
bude CaseR1
a CaseR2
(tedy vliv vegetace), Second group
bude Ca
a cond
(vliv prostředí, proměnné které vyšli jako důležité při postupném výběru).Variation Partitioning Results for Two Groups in Analysis 'Var-part-2groups-Conditional-effects-tested' Explained variation Fraction Variation(adj) % of Explained % of All DF Mean Square a 0.062197 22.2 6.2 2 0.04677 b 0.015341 5.5 1.5 2 0.02446 c 0.20208 72.3 20.2 -- -- Total Explained 0.27962 100.0 28.0 4 0.08706 All Variation 1 -- 100.0 42 --
Variation Partitioning Results for Two Groups in Analysis 'Var-part-2groups-Conditional-effects-tested' Significance tests Tested Fraction F P a+b+c 5.1 0.002 a 2.7 0.002 b 1.4 0.072
Obě skupiny proměnných vysvětlí celkem 28% variability (adjustovaný R2), z toho čistý efekt vegetace je 6.2% (a ten je průkazný na p < 0.01), zatímco čistý efekt proměnných prostředí je jen 1.5% (a není průkazný, p = 0.072). Sdílená složka je největší - 20.2%, protože obě sady proměnných jsou výrazně korelované (prostředí ovlivňuje vegetaci).
Interpretace by byla asi následující: vegetace je lepší prediktor druhového složení než měřené proměnné prostředí, a obsahuje i informaci, kterou proměnné prostředí neobsahují. Pokud bych měl jistotu, že v proměnných prostředí jsou všechny faktory, které mohou být relevantní pro složení plžích společenstev, pak by se dalo říci, že variabilita kterou vegetace vysvětlí navíc odráží biotické interakce mezi plži a rostlinami. Protože ale s jistotou nevím, jestli jsou měřeny všechny proměnné prostředí, je pravděpodobné, že variabilita kterou vegetace vysvětlí navíc je způsobena některým z faktorů prostředí který nebyl změřen (proto není ve frakci [b] nebo [c]) a kytky na něj reagují stejně jako plži. Celý model a parciální efekt vlivu vegetace jsou průkazné, parciální vliv proměnných prostředí už průkazný není.
Show brief version of notebooks with analysis results
. Když analýzu provedete znovu, objeví se vám i další záložky se skóre na osách a dalšími informacemi.