A survey of 41 plots in natural stands of pure or mixed Nothofagus forest at the foot of the Andes was conducted to investigate variation in stand structure. Each plot was classified by the dominant tree species into one of three stand types: Coigue (13 plots), Rauli (9 plots), or Roble (19 plots). For each plot, stand density (number of trees per hectare, SD) and mean quadratic diameter (in cm, QD) were recorded. To facilitate analysis, three dummy variables (d1, d2, d3) were created to represent the stand type factor. The main aim of the study was to model stand density as a function of quadratic diameter and to compare this relationship among the three types of Nothofagus stands.
forest2A data frame with 7 variables: DPlot, Type, SD, QD, d1, d2, d3.
Factor. Unique identifier for each forest plot.
Factor. Forest stand type classified by the dominant Nothofagus species: "Coigue", "Rauli", or "Roble".
Integer. Stand density, as the number of trees per hectare in each plot.
Numeric. Mean quadratic diameter (in centimeters) of trees in the plot.
Integer. Dummy variable for stand type Coigue (1 if Type is "Coigue", 0 otherwise).
Integer. Dummy variable for stand type Rauli (1 if Type is "Rauli", 0 otherwise).
Integer. Dummy variable for stand type Roble (1 if Type is "Roble", 0 otherwise).
anova(lm(log(SD) ~ log(QD) * (d2 + d3), data = forest2))
Run the code above in your browser using DataLab