Renamed psuAbund to psuDensity and changed output variable name
Changed avg.abun variable to dbar
Minor changes in rvcData internal var names
Changed vars to agg.by to match other functions
Added includes.protected parameter to stratData function
Allows user to include protected substrata, if available