# Analysis of community ecology data in R

David Zelený

### Site Tools

en:expl_var_exercise

# Differences

This shows you the differences between two versions of the page.

 en:expl_var_exercise [2019/02/10 16:07]David Zelený en:expl_var_exercise [2019/02/10 16:33]David Zelený [Explained variation and Monte Carlo permutation test (constrained ordination)] Both sides previous revision Previous revision 2019/02/26 23:35 David Zelený 2019/02/10 16:33 David Zelený [Explained variation and Monte Carlo permutation test (constrained ordination)] 2019/02/10 16:07 David Zelený 2019/01/26 20:04 David Zelený [Explained variance (constrained ordination)] 2019/01/26 19:51 David Zelený 2017/10/11 20:36 external edit2017/05/02 00:43 David Zelený created 2019/02/26 23:35 David Zelený 2019/02/10 16:33 David Zelený [Explained variation and Monte Carlo permutation test (constrained ordination)] 2019/02/10 16:07 David Zelený 2019/01/26 20:04 David Zelený [Explained variance (constrained ordination)] 2019/01/26 19:51 David Zelený 2017/10/11 20:36 external edit2017/05/02 00:43 David Zelený created Last revision Both sides next revision Line 6: Line 6: [[{|width: 7em; background-color:​ white; color: navy}expl_var_examples|Examples]] [[{|width: 7em; background-color:​ white; color: navy}expl_var_examples|Examples]] [[{|width: 7em; background-color:​ light; color: firebrick}expl_var_exercise|**Exercise** {{::​lock-icon.png?​nolink|}}]] [[{|width: 7em; background-color:​ light; color: firebrick}expl_var_exercise|**Exercise** {{::​lock-icon.png?​nolink|}}]] + + ===== Exercise 1 ===== + (This exercise is based on //Case study 4// from Šmilauer & Lepš 2014, and the coding was partly inspired by the [[http://​www.fromthebottomoftheheap.net/​2014/​11/​03/​randomized-complete-block-designs-and-vegan/​|blog post of Gavin Simpson]], analysing this data - I recommend it for checking on more details.)​ + + Use vegetation data from an [[en:​data:​removal_experiment|experimental study on wet meadow]], which tried to answer the question how is abundance and species composition of seedlings of herb species on wet meadow influenced by experimental removal of dominant clonal plant species (//Nardus stricta//), plant litter (remnants of plant biomass from the previous season in case that the meadow was not mown) and removal of litter+mosses (by raking). Before analysis, please get familiar with the data you will be using. Data are from completely randomized block design with four blocks, each with four plots (one per treatment, including control). The block design will need to be taken into consideration when testing the significance of explained variation using Monte Carlo permutation test. + + - Import datasets ''​seedlings.spe'' ​ and ''​seedlings.env''​ - the first containing species composition of seedlings per plot, the second contain assignment of the plot to the block and experimental treatment. + - First, calculate DCA on compositional data (no transformation) to decide whether to use linear (RDA) or unimodal (CCA) constrained ordination method. + - According to the result of previous analysis, use either RDA or CCA, where you will explain the species composition by treatment. Since the experiment design is a block design, you need to specify the variable ''​seedlings.env\$blocks''​ as covariable in the model. + - Draw the resulting ordination diagram. + - Test the significance of the variation explained by treatment (use model based permutation test). + - Test the significance of indiviudal ordination axes. Which of them is the most important, and to which level of the treatment it corresponds?​ + + + + - ''​read.delim''; ​ + - check the length of the first axis in ''​decorana''​ results. If over 4, use CCA, if less then 3, use RDA, 3-4 -both are fine; + - the formula will look like this: ''​seedlings.spe ~ treatment + Condition (block), data = seedlings.env'';​ + - ''​ordiplot''​ + - ''​anova''​ without any other arguments; since the ''​blocks''​ are defined as covariable, he will test the significance of variation after removal of their variation; + - ''​anova''​ with argument ''​by = "​axis"''​. + ​ + [[en:​solutions:​monte_carlo_solution_ex1|Solution]]