#####
#####
#####read in appropriate libraries
library(AICcmodavg)
library(lme4)
library(lattice)
library(mgcv)
library(gamm4)
library(car)
library(cvTools)
library(influence.ME)
library(sp)
library(rgdal)
library(raster)
library(maptools)
library(ape)
library(geosphere)
library(geoR)
library(gstat)

#read in data
Ponded <- read.csv("BinomialTotalExtremeChangePresent3.csv", header = TRUE)
head(Ponded)
nrow(Ponded)

#####
#####
#####set nominal variables as factors
Ponded$Wetland = factor(Ponded$Wetland)
Ponded$Wettype = factor(Ponded$Wettype)
Ponded$Year = factor(Ponded$Year)
Ponded$agtype = factor(Ponded$agtype)

#####
#####
#####make predictions for ponded area manually
#with original data
Ponded$Pred.Log.Meters.ManualFixedOnly <- NA
for(i in 1:nrow(Ponded)){
	Ponded$Pred.Log.Meters.ManualFixedOnly[[i]] = 8.25906+
	+ 1.32752 * Ponded$SemiPerm[[i]]+
	- 0.76114 * Ponded$Temporary[[i]]+
	- 0.28251 * Ponded$Pivots[[i]]+
	+ 0.36406 * Ponded$Dryland[[i]]+
	- 0.56592 * Ponded$Gravity[[i]]+
	- ((0.40144 * Ponded$log.NewfootPA[[i]] - 0.40144 * -4.23480) / 0.53077)+
	+ ((0.07812 * Ponded$espro2PostChange[[i]] - 0.07812 * 0.28077) / 0.47154)+
	+ ((0.30366 * Ponded$log.lsprsPostChange[[i]] - 0.30366 * 5.45610) / 0.29539)+
	+ ((0.26013 * Ponded$log.wprsPostChange[[i]] - 0.26013 * 4.31325) / 0.76437)+
	+ ((0.77571 * Ponded$wtmaxsfPostChange[[i]] - 0.77571 * 26.01538) / 4.40183)
}

#####
#####
#####calculate the inverse log for predicted values
Ponded$Pred.Meters.ManualFixedOnly <- exp(Ponded$Pred.Log.Meters.ManualFixedOnly)

#obtain mean flooded area
mean(Ponded$Pred.Meters.ManualFixedOnly)

#obtain total flooded area
sum(Ponded$Pred.Meters.ManualFixedOnly)