Realised I am essentially working on two different analyses. Split script into two: (1) LHvsGBR.R is original idea from Erika's thesis, (2) GBRSourcePool.R generates species pools for Andrew's GBR data.
Added new script - see previous comment