Le presenti analisi tengono in considerazione solo la vegetazione erbacea. I codice per tipi di vegetazione sono:

brcp: vegetazione dominata da Brachypodium genuense

nard: vegetazione dominata da Nardus stricta

fest: vegetazione dominata da Festuca circummediterranea

glob: vegetazione dominata da Globularia meridionalis

sesl: vegetazione dominata da Sesleria juncifolia

- analisi #1

relazione tra le vegetazioni studiate e la profondità di suolo

boxplot(df_noNA$prof_suolo~df_noNA$vegetazione, main="vegetazione e profondità di suolo",
  xlab="tipi di vegetazione", ylab="profondità di suolo (cm)")

testo le differenze tra livelli (pairwaise Wilcox test)

kruskal.test(prof_suolo ~ vegetazione, data = df_noNA)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  prof_suolo by vegetazione
## Kruskal-Wallis chi-squared = 25.354, df = 4, p-value = 4.269e-05
pairwise.wilcox.test(df_noNA$prof_suolo, df_noNA$vegetazione,
                     exact=F, p.adjust.method = "bonferroni")
## 
##  Pairwise comparisons using Wilcoxon rank sum test with continuity correction 
## 
## data:  df_noNA$prof_suolo and df_noNA$vegetazione 
## 
##      brcp   fest   glob   nard  
## fest 0.1334 -      -      -     
## glob 0.0055 1.0000 -      -     
## nard 1.0000 0.0504 0.0094 -     
## sesl 0.0195 0.8228 1.0000 0.0102
## 
## P value adjustment method: bonferroni

risultato: esistono relazioni significative tra la profondità di suolo e i diversi tipi di vegetazione.

Le vegetazioni a dominanza di Brachypodium genuense e Nardus stricta preferiscono suoli profondi, mentre le vegetazioni a dominanza di Sesleria juncifolia e Globularia meridionalis preferiscono suoli sottili. Una posizione intermedia è occupata dalle vegetazioni a dominanza di Festuca circummediterranea

- analisi #2

relazione tra le vegetazioni studiate e la matrice (rocciosità percentuale all’interno del plot)

testo le differenze tra livelli (pairwaise Wilcox test)

kruskal.test(petr ~ vegetazione, data = df_noNA)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  petr by vegetazione
## Kruskal-Wallis chi-squared = 30.02, df = 4, p-value = 4.848e-06
pairwise.wilcox.test(df_noNA$petr, df_noNA$vegetazione,
                     exact=F, p.adjust.method = "bonferroni")
## 
##  Pairwise comparisons using Wilcoxon rank sum test with continuity correction 
## 
## data:  df_noNA$petr and df_noNA$vegetazione 
## 
##      brcp   fest   glob   nard  
## fest 1.0000 -      -      -     
## glob 0.0032 0.0329 -      -     
## nard 1.0000 0.2847 0.0020 -     
## sesl 0.0144 0.0530 1.0000 0.0075
## 
## P value adjustment method: bonferroni

risultato: esistono relazioni significative tra la quantità di matrice (petrosità) e i diversi tipi di vegetazione.

Le vegetazioni a dominanza di Brachypodium genuense e Nardus stricta preferiscono suoli con poca matrice, mentre le vegetazioni a dominanza di Sesleria juncifolia e Globularia meridionalis preferiscono suoli ricchi di matrice. Una posizione intermedia è occupata dalle vegetazioni a dominanza di Festuca circummediterranea.

- analisi #3

associazione tra le categorie di tessitura di suolo e la vegetazione

mosaicplot che mostra l’associazione tra variabili categoriche. I valori sulle ascisse indicano: A: argilloso, S: sabbioso, L: limo.

Le categorie di suolo sono stati attribuite secondo il rango più alto della tassonomia dei suoli di USDA (partendo dalle attribuzioni in classi fatte dalle analisi). la fonte è shorturl.at/qCSW2 (documento su Microsoft teams)

testo l’associazione delle mie categorie usando il GTest perchè il campione è molto piccolo rispetto alle categorie prese in considerazione. I test Chi-sqr e Fisher potrebbero dare risultati alterati

GTest(df_noNA$vegetazione, df_noNA$tessitura_class)
## 
##  Log likelihood ratio (G-test) test of independence without correction
## 
## data:  df_noNA$vegetazione and df_noNA$tessitura_class
## G = 19.697, X-squared df = 8, p-value = 0.01155

risultato: esistono relazioni significative tra la categorie di tessitura e la vegetazione. Nardeti e Brachypodieti tendono ad associarsi a suoli con tessitura argillosa, Glbularieti, Festuceti e Seslerieti tendono ad associarsi inferiormente a suoli argillosi ma maggiormente a suoli limosi.

I suoli sabbiosi sono principalmente associati a Seslerieti e Festuceti.

- analisi #4

relazione tra le vegetazioni studiate e la pressione di pascolo (misurata tramite proxy)

testo le differenze tra livelli (pairwaise Wilcox test)

kruskal.test(morso ~ vegetazione, data = df_noNA)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  morso by vegetazione
## Kruskal-Wallis chi-squared = 26.727, df = 4, p-value = 2.257e-05
pairwise.wilcox.test(df_noNA$petr, df_noNA$vegetazione,
                     exact=F, p.adjust.method = "bonferroni")
## 
##  Pairwise comparisons using Wilcoxon rank sum test with continuity correction 
## 
## data:  df_noNA$petr and df_noNA$vegetazione 
## 
##      brcp   fest   glob   nard  
## fest 1.0000 -      -      -     
## glob 0.0032 0.0329 -      -     
## nard 1.0000 0.2847 0.0020 -     
## sesl 0.0144 0.0530 1.0000 0.0075
## 
## P value adjustment method: bonferroni

risultato: esistono relazioni significative tra la pressione di pascolo e tutti i tipi di vegetazione.

Le vegetazioni a dominanza di Nardus stricta si trovano in aree più pascolate. Un’alta pressione di pascolo è presente anche nelle vegetazioni a dominanza di Festuca circummediterranea. Una posizione intermedia è occupata dalle vegetazioni a dominanza di Globularia meridionalis. Infine le vegetazioni a dominanza di Sesleria juncifolia e Brachypodium genuense si trovano in aree meno o per nulla pascolate.