/Tortorelli-etal-2022-Ecology-Invasion-response-to-community-traits-and-environment

This repository contains the R code and data to support the Manuscript "Community invasion resistance is influenced by interactions between plant traits and site productivity", Tortorelli et al. 2022, Ecology

Invasion-response-to-community-traits-and-environment

Metadata for sript and data files

R Script files

  • Calculating community weighted means: Code walks through calculating community weighted means for each subplot

  • Calculating nearest species metrics: Code walks through calculating nearest species distance metrics for each subplot

  • Modeling community trait metrics: All models are present here (all traits by all community metrics). Plots were initiated using this code and cleaned up in inkscape as svg files

  • Species traits NMS: Here we calculated and plotted species in trait space using NMS

  • Subplot traits NMS: Here we calculated and plotted sublots in community weighted mean trait space using NMS

data files

  • 2020 biomass by subplot: biomass collected from 20x20cm subplots in 2020

    • plot_quad: represents the individual subplots from which biomass was collected.

      • First 4 letters refer to vegetation code (ARAR = Artemisia arbuscula = low sage steppe; ARRI = artimisia rigida= scabflat; SEEP = wet meadow)
      • Number after vegetation type references the block number (1*5)
      • "CC." indicates that these were all community plots with ventenata seed added
      • Final number indicates the number of the subplot at each plot (1*7)

    • resident_g: resident biomass (excluding ventenata) collected from within the subplots in 2020 measured in g/subplot

    • vedu_g: ventenata biomass collected from within the subplots in 2020 measured in g/subplot


  • nearest species distances: nearest species distance calculated from the "Calculating nearest species metrics" script

    • plot_quad: see above

    • finerootV.totalrootV (fine root Volume to total root volume ratio)

    • height (cm)

    • N.perc (% leaf N)

    • root.shoot (root to shoot mass ratio )

    • RootAvgDiam.mm (root average diameter (mm))

    • RootLength.cm (cm)

    • SLA (specific leaf area; leaf area/ leaf mass)

    • community weighted means: Community weighted means output from the "Calculating community weighted means" script

    • plot_quad: see above

    • finerootV.totalrootV (fine root Volume to total root volume ratio)

    • height (cm)

    • N.perc (% leaf N)

    • root.shoot (root to shoot mass ratio )

    • RootAvgDiam.mm (root average diameter (mm))

    • RootLength.cm (cm)

    • SLA (specific leaf area; leaf area/ leaf mass)

    • vegtype: (ARAR = Artemisia arbuscula = low sage steppe; ARRI = artimisia rigida= scabflat; SEEP = wet meadow)

    • plot: vegetation type with block number

    • plotno: just the block number

  • 2020_USDAspeciesCover_4traits_seededCommunitySubplots: Species foliar cover (values) for each subplot (row names) for each species represented by their USDA plants code (column names)

  • allTraits_averagedBySpecies_USDAcodes: Averaged raw trait values for each species (row names) for each trait (column names)