data(goats)
str(goats)
summary(goats)
goats$exp.HLI <- exp(goats$HLI)
goats$sin.SLOPE <- sin(pi * goats$SLOPE / 180)
goats$ELEVATION <- scale(goats$ELEVATION)
goats$ET <- scale(goats$ET)
goats$TASP <- scale(goats$TASP)
summary(m1 <- rspf(STATUS ~ TASP + sin.SLOPE + ELEVATION, goats, B = 99))
summary(m2 <- rspf(STATUS ~ TASP + ELEVATION, goats, B = 99))
AIC(m1, m2)
plot(m1)Run the code above in your browser using DataLab