#Form a gauged pooling group and plot the diagnostics with gauged = TRUE
Pool.96001 <- Pool(GetCDs(96001))
DiagPlots(Pool.96001, gauged = TRUE)
#Form an ugauged pooling group and plot the diagnostics
Pool.96001 <- Pool(GetCDs(96001), exclude = 96001)
DiagPlots(Pool.96001)
Run the code above in your browser using DataLab