Data are included in a book Šmilauer & Lepš (2014) in Case study 3, and originally were published by Pyšek & Lepš (1991)

The experiment was focused on diversity and species richness of weed communities in barley fields after application of industrial fertilizers. The experiment was originally established as a part of hydrobiological research, and for practical reasons it was not established using properly randomized design - individual plots are practically pseudoreplications, but for the purpose of this exercise we ignore it.

All together 122 plots of the size 1 m² were fertilized by three different fertilizers (for the purpose of this study the types of fertilizers are not distinguished), and on each plot species composition of weed community was recorded (with species covers estimated using ordinal scale from 1 to 7). Fertilizing simultaneously influences also growth and cover of barley, and the cover of barley can in turn also influence species composition of weeds; for this reason, the barley cover in the plot was also recorded.

Figure 1: Examples of weeds occurring together with barley in the field (photographs from www.botanickafotogalerie.cz, credit go to authors).

Close to the village Samšín in South Bohemia, elevation 500 m a.s.l. For further details, see Pyšek & Lepš (1991).

Here we are using simplified dataset without distinguishing the types of fertilizers - we distinguish only the dose.

fertil.spe <- read.delim ('https://raw.githubusercontent.com/zdealveindy/anadat-r/master/data/fertil.spe.txt', row.names = 1)
fertil.env <- read.delim ('https://raw.githubusercontent.com/zdealveindy/anadat-r/master/data/fertil.env.txt', row.names = 1)
 

• Pyšek P. & Lepš J. (1991): Response of weed community to nitrogen fertilization: a multivariate analysis. Journal of Vegetation Science, 2: 237-244 pdf
• Šmilauer, P. & Lepš, J. (2014) Multivariate Analysis of Ecological Data using Canoco 5. Second Edition. Cambridge University Press, Cambridge, UK.online supplementary data - Chapter 14 (Case Study 3)