dataWO <- data %>%
    filter(Species == "QUGA") %>%
    filter(!is.na(BAI_avg))

# unstandardized
dataWO_n <- dataWO

# standardize variables
dataWO$DBH <- arm::rescale(dataWO$DBH)
dataWO$H_14m_all <- arm::rescale(dataWO$H_14m_all)
dataWO$H_14m_DF <- arm::rescale(dataWO$H_14m_DF)
dataWO$H_14m_nDF <- arm::rescale(dataWO$H_14m_nDF)
dataWO$H_12m_all <- arm::rescale(dataWO$H_12m_all)
dataWO$H_12m_DF <- arm::rescale(dataWO$H_12m_DF)
dataWO$H_12m_nDF <- arm::rescale(dataWO$H_12m_nDF)
dataWO$DF_shade <- arm::rescale(dataWO$DF_shade)
dataWO$nDF_shade <- arm::rescale(dataWO$nDF_shade)

M3_B_WO.17 <- lme(fixed = log(BAI_avg) ~ DBH + H_14m_nDF + DF_shade, data=dataWO, random = ~ 1|Plot_No)
AICc(M3_B_WO.17)

## [1] 174.8177

M3_B_WO.17_n <- lme(fixed = log(BAI_avg) ~ DBH + H_14m_nDF + DF_shade, data=dataWO_n, random = ~ 1|Plot_No)
AICc(M3_B_WO.17_n)

## [1] 184.1404

M3_B_WO.31 <- lme(fixed = log(BAI_avg) ~ DBH + H_12m_DF + H_14m_nDF, data=dataWO, random = ~ 1|Plot_No)
AICc(M3_B_WO.31)

## [1] 178.2103

M3_B_WO.31_n <- lme(fixed = log(BAI_avg) ~ DBH + H_12m_DF + H_14m_nDF, data=dataWO_n, random = ~ 1|Plot_No)
AICc(M3_B_WO.31_n)

## [1] 191.7922

M3_B_WO.71 <- lme(fixed = log(BAI_avg) ~ DBH + H_14m_all, data=dataWO, random = ~ 1|Plot_No)
AICc(M3_B_WO.71)

## [1] 179.2528

M3_B_WO.71_n <- lme(fixed = log(BAI_avg) ~ DBH + H_14m_all, data=dataWO_n, random = ~ 1|Plot_No)
AICc(M3_B_WO.71_n)

## [1] 190.5091