Chapter 21 Arne data analysis

21.1 Load packages and data

library(tidyverse)
library(mgcv)
library(aqm)
library(sf)

In order to write up the assignment you can use maps generated in the previous chapters to illustrate how GIS techniques were applied. You can then use the data that have been processed using the GIS components of R to conduct further analysis. It might be speculated that the density of regenerating pines (numbers of pines per square meter) in the heathland area may be influenced by the moisture of the soil, the amount of insolation the microsite receives, the height above sea level and/or the distance to the nearest seed source. Are any of these elements clearly related to pine density? You can analyse these data using the simpler non-aptial techniques that you have seen in the fist section of the unit.

data("arne_assignment_2022")
d %>% st_drop_geometry() ->d
d$Year<-as.factor(d$Year)
str(d)
## 'data.frame':    184 obs. of  10 variables:
##  $ Year        : Factor w/ 3 levels "2017","2018",..: 3 3 3 3 3 3 3 3 3 3 ...
##  $ pine_density: num  0.875 0.557 0.955 0.239 0.318 0.637 0.398 0.637 0.557 0.477 ...
##  $ dtm         : num  12.6 12.8 13 13.3 13.2 ...
##  $ insol       : num  1.83 1.88 1.81 1.78 1.84 ...
##  $ twi         : num  4.24 3.96 3.64 4.14 4.26 ...
##  $ slope       : num  0.0224 0.0338 0.0339 0.015 0.0213 ...
##  $ aspect      : num  1.57 2.15 1.98 1.38 2.13 ...
##  $ path_dist   : num  44.4 36.56 21.48 7.22 4.4 ...
##  $ min_dist    : num  16.93 14.38 8.09 4.27 24.98 ...
##  $ n50         : int  5 5 7 10 5 8 8 5 10 8 ...

Year = Year data was collected Pine density = numbers of pines per meter squared
twi = Topographic wetness index (no units)
sol = Direct beam insolation
slope = slope in degrees
aspect = aspect in degrees
pat_dist = Distance to nearest path in meters min_dist = minimum distance to nearest seed source.
n50 = Number of large trees within 50 meters

Are there any clear relationships between pine density and these variables? Analyse the data using techniques that you have learned on the course.

If you do not find relationships then this finding would still be useful. It would imply that the pine regeneration essentially takes place quite randomly over the area within which the data were collected.