Compare classifications between two analysis results
Examples
# \donttest{
results_2010 <- run_builtin_analysis(scale = "county", year = 2010)
#>
#> ======================================================================
#> BATCH MANURESHED ANALYSIS
#> ======================================================================
#> Year: 2010
#> Scale: county
#> Nutrients: nitrogen, phosphorus
#> ----------------------------------------------------------------------
#> Checking data availability...
#> Data availability confirmed
#> Available scales:county, huc8, huc2
#> Available years forcounty:1987-2016
#> Built-in WWTP data:Available (2007-2016)
#>
#> Loading built-in NuGIS data...
#> Using cached version of nugis_county_data
#> Loaded NuGIS county data for year 2010
#> Number of spatial units: 3064
#> Loading built-in spatial boundaries...
#> Using cached version of county_boundaries
#> Loaded county boundaries
#> Number of spatial units: 3112
#> Calculating cropland threshold...
#> Cropland threshold:1235.53acres
#>
#> Processing agricultural classifications...
#> Starting complete agricultural classification for county scale...
#> Processed NuGIS data for county scale:
#> Spatial units: 3064
#> Converted P2O5 to P using factor: 0.436
#> Nitrogen classification summary:
#> Excluded: 145 units
#> Sink_Deficit: 2362 units
#> Sink_Fertilizer: 287 units
#> Source: 83 units
#> Within_County: 187 units
#> Phosphorus classification summary:
#> Excluded: 145 units
#> Sink_Deficit: 1760 units
#> Sink_Fertilizer: 469 units
#> Source: 343 units
#> Within_County: 347 units
#> Agricultural classification complete!
#> Applied threshold: 1235.53 acres
#> Agricultural classification complete
#> Spatial units processed:3112
#> Nitrogen classes:Excluded ( 193 ), Sink_Deficit ( 2362 ), Sink_Fertilizer ( 287 ), Source ( 83 ), Within_County ( 187 )
#> Phosphorus classes:Excluded ( 193 ), Sink_Deficit ( 1760 ), Sink_Fertilizer ( 469 ), Source ( 343 ), Within_County ( 347 )
#>
#> Processing WWTP data...
#> Nutrients:nitrogen, phosphorus
#> WWTP year:2010
#> Load units:kg
#> Data source:Built-in (2010)
#> Loading built-in nitrogen WWTP data for2010...
#> Using cached version of wwtp_nitrogen_combined
#> Loaded WWTP nitrogen data for year 2010
#> Number of facilities: 24727
#> Loading built-in phosphorus WWTP data for2010...
#> Downloading wwtp_phosphorus_combined from OSF...
#> File ID: 689c97404cae9f37ea13c731
#> Downloaded successfully (3.75 MB)
#> Loaded WWTP phosphorus data for year 2010
#> Number of facilities: 7548
#> Processing nitrogen WWTP facilities...
#> Filtered for positive nitrogen loads:
#> Original: 24727 facilities
#> With positive loads: 24727 facilities
#> WWTP nitrogen source classification:
#> Minor Source: 20327 facilities
#> Small Source: 2445 facilities
#> Medium Source: 987 facilities
#> Large Source: 745 facilities
#> Very Large Source: 223 facilities
#> Created spatial WWTP data with 24727 facilities
#> Looking for boundary ID column: 'FIPS'
#> Available columns in boundaries: FIPS, County, State_Name, geometry
#> Aggregating WWTP nitrogen loads by spatial boundaries...
#> Aggregation complete:
#> WWTP facilities: 24727
#> Spatial units with facilities: 2370
#> Total nitrogen load: 2010631.56 tons/year
#> Processing phosphorus WWTP facilities...
#> Filtered for positive phosphorus loads:
#> Original: 7548 facilities
#> With positive loads: 7548 facilities
#> WWTP phosphorus source classification:
#> Minor Source: 4626 facilities
#> Small Source: 1572 facilities
#> Medium Source: 669 facilities
#> Large Source: 553 facilities
#> Very Large Source: 128 facilities
#> Created spatial WWTP data with 7548 facilities
#> Looking for boundary ID column: 'FIPS'
#> Available columns in boundaries: FIPS, County, State_Name, geometry
#> Aggregating WWTP phosphorus loads by spatial boundaries...
#> Aggregation complete:
#> WWTP facilities: 7548
#> Spatial units with facilities: 1539
#> Total phosphorus load: 99331.55 tons/year
#> WWTP data processing complete
#> nitrogen:24727facilities in2370spatial units
#> phosphorus:7548facilities in1539spatial units
#>
#> Integrating WWTP and agricultural data...
#> Integrating WWTP nitrogen data with agricultural classifications...
#> Using agricultural ID column: FIPS
#> Combined nitrogen classification summary:
#> Excluded: 193 units
#> Sink_Deficit: 2243 units
#> Sink_Fertilizer: 266 units
#> Source: 180 units
#> Within_County: 230 units
#> Integrating WWTP phosphorus data with agricultural classifications...
#> Using agricultural ID column: FIPS
#> Combined phosphorus classification summary:
#> Excluded: 193 units
#> Sink_Deficit: 1685 units
#> Sink_Fertilizer: 435 units
#> Source: 442 units
#> Within_County: 357 units
#> Integration complete
#>
#> Saving results...
#> Saved spatial data to: /tmp/RtmpJm2r1A/county_agricultural_2010.rds
#> File size: 1.04 MB
#> Rows: 3112, Columns: 17
#> Geometry type: POLYGON
#> CRS: EPSG:5070
#> Saved spatial data to: /tmp/RtmpJm2r1A/county_nitrogen_integrated_2010.rds
#> File size: 1.1 MB
#> Rows: 3112, Columns: 23
#> Geometry type: POLYGON
#> CRS: EPSG:5070
#> Saved centroid data to: /tmp/RtmpJm2r1A/county_nitrogen_centroids_2010.csv
#> File size: 752.62 KB
#> Rows: 3112, Columns: 24
#> Longitude range: [-124.158, -67.637]
#> Latitude range: [25.49, 48.826]
#> Saved spatial data to: /tmp/RtmpJm2r1A/county_phosphorus_integrated_2010.rds
#> File size: 1.09 MB
#> Rows: 3112, Columns: 23
#> Geometry type: POLYGON
#> CRS: EPSG:5070
#> Saved centroid data to: /tmp/RtmpJm2r1A/county_phosphorus_centroids_2010.csv
#> File size: 720.6 KB
#> Rows: 3112, Columns: 24
#> Longitude range: [-124.158, -67.637]
#> Latitude range: [25.49, 48.826]
#> Saved analysis summary to: /tmp/RtmpJm2r1A/analysis_summary_2010.rds
#> Format: RDS
#> File size: 0.77 KB
#> Results saved to:/tmp/RtmpJm2r1A
#> Files created:6
#>
#> ======================================================================
#> ANALYSIS COMPLETE
#> ======================================================================
#> Processing time:0.09minutes
#> Scale:county
#> Year:2010
#> Nutrients analyzed:nitrogen, phosphorus
#> Spatial units:3112
#> WWTP facilities:32275
#> Output directory:/tmp/RtmpJm2r1A
#> ======================================================================
results_2016 <- run_builtin_analysis(scale = "county", year = 2016)
#>
#> ======================================================================
#> BATCH MANURESHED ANALYSIS
#> ======================================================================
#> Year: 2016
#> Scale: county
#> Nutrients: nitrogen, phosphorus
#> ----------------------------------------------------------------------
#> Checking data availability...
#> Data availability confirmed
#> Available scales:county, huc8, huc2
#> Available years forcounty:1987-2016
#> Built-in WWTP data:Available (2007-2016)
#>
#> Loading built-in NuGIS data...
#> Using cached version of nugis_county_data
#> Loaded NuGIS county data for year 2016
#> Number of spatial units: 3058
#> Loading built-in spatial boundaries...
#> Using cached version of county_boundaries
#> Loaded county boundaries
#> Number of spatial units: 3112
#> Calculating cropland threshold...
#> Cropland threshold:1235.53acres
#>
#> Processing agricultural classifications...
#> Starting complete agricultural classification for county scale...
#> Processed NuGIS data for county scale:
#> Spatial units: 3058
#> Converted P2O5 to P using factor: 0.436
#> Nitrogen classification summary:
#> Excluded: 158 units
#> Sink_Deficit: 2450 units
#> Sink_Fertilizer: 214 units
#> Source: 83 units
#> Within_County: 153 units
#> Phosphorus classification summary:
#> Excluded: 158 units
#> Sink_Deficit: 1891 units
#> Sink_Fertilizer: 372 units
#> Source: 317 units
#> Within_County: 320 units
#> Agricultural classification complete!
#> Applied threshold: 1235.53 acres
#> Agricultural classification complete
#> Spatial units processed:3112
#> Nitrogen classes:Excluded ( 212 ), Sink_Deficit ( 2450 ), Sink_Fertilizer ( 214 ), Source ( 83 ), Within_County ( 153 )
#> Phosphorus classes:Excluded ( 212 ), Sink_Deficit ( 1891 ), Sink_Fertilizer ( 372 ), Source ( 317 ), Within_County ( 320 )
#>
#> Processing WWTP data...
#> Nutrients:nitrogen, phosphorus
#> WWTP year:2016
#> Load units:kg
#> Data source:Built-in (2016)
#> Loading built-in nitrogen WWTP data for2016...
#> Using cached version of wwtp_nitrogen_combined
#> Loaded WWTP nitrogen data for year 2016
#> Number of facilities: 20846
#> Loading built-in phosphorus WWTP data for2016...
#> Using cached version of wwtp_phosphorus_combined
#> Loaded WWTP phosphorus data for year 2016
#> Number of facilities: 10148
#> Processing nitrogen WWTP facilities...
#> Filtered for positive nitrogen loads:
#> Original: 20846 facilities
#> With positive loads: 20846 facilities
#> WWTP nitrogen source classification:
#> Minor Source: 17801 facilities
#> Small Source: 1799 facilities
#> Medium Source: 729 facilities
#> Large Source: 417 facilities
#> Very Large Source: 100 facilities
#> Created spatial WWTP data with 20846 facilities
#> Looking for boundary ID column: 'FIPS'
#> Available columns in boundaries: FIPS, County, State_Name, geometry
#> Aggregating WWTP nitrogen loads by spatial boundaries...
#> Aggregation complete:
#> WWTP facilities: 20846
#> Spatial units with facilities: 2581
#> Total nitrogen load: 582738.6 tons/year
#> Processing phosphorus WWTP facilities...
#> Filtered for positive phosphorus loads:
#> Original: 10148 facilities
#> With positive loads: 10148 facilities
#> WWTP phosphorus source classification:
#> Minor Source: 6846 facilities
#> Small Source: 1839 facilities
#> Medium Source: 795 facilities
#> Large Source: 561 facilities
#> Very Large Source: 107 facilities
#> Created spatial WWTP data with 10148 facilities
#> Looking for boundary ID column: 'FIPS'
#> Available columns in boundaries: FIPS, County, State_Name, geometry
#> Aggregating WWTP phosphorus loads by spatial boundaries...
#> Aggregation complete:
#> WWTP facilities: 10148
#> Spatial units with facilities: 1966
#> Total phosphorus load: 177064.66 tons/year
#> WWTP data processing complete
#> nitrogen:20846facilities in2581spatial units
#> phosphorus:10148facilities in1966spatial units
#>
#> Integrating WWTP and agricultural data...
#> Integrating WWTP nitrogen data with agricultural classifications...
#> Using agricultural ID column: FIPS
#> Combined nitrogen classification summary:
#> Excluded: 212 units
#> Sink_Deficit: 2358 units
#> Sink_Fertilizer: 185 units
#> Source: 157 units
#> Within_County: 200 units
#> Integrating WWTP phosphorus data with agricultural classifications...
#> Using agricultural ID column: FIPS
#> Combined phosphorus classification summary:
#> Excluded: 212 units
#> Sink_Deficit: 1819 units
#> Sink_Fertilizer: 344 units
#> Source: 389 units
#> Within_County: 348 units
#> Integration complete
#>
#> Saving results...
#> Saved spatial data to: /tmp/RtmpJm2r1A/county_agricultural_2016.rds
#> File size: 1.04 MB
#> Rows: 3112, Columns: 17
#> Geometry type: POLYGON
#> CRS: EPSG:5070
#> Saved spatial data to: /tmp/RtmpJm2r1A/county_nitrogen_integrated_2016.rds
#> File size: 1.1 MB
#> Rows: 3112, Columns: 23
#> Geometry type: POLYGON
#> CRS: EPSG:5070
#> Saved centroid data to: /tmp/RtmpJm2r1A/county_nitrogen_centroids_2016.csv
#> File size: 760.52 KB
#> Rows: 3112, Columns: 24
#> Longitude range: [-124.158, -67.637]
#> Latitude range: [25.49, 48.826]
#> Saved spatial data to: /tmp/RtmpJm2r1A/county_phosphorus_integrated_2016.rds
#> File size: 1.09 MB
#> Rows: 3112, Columns: 23
#> Geometry type: POLYGON
#> CRS: EPSG:5070
#> Saved centroid data to: /tmp/RtmpJm2r1A/county_phosphorus_centroids_2016.csv
#> File size: 736.84 KB
#> Rows: 3112, Columns: 24
#> Longitude range: [-124.158, -67.637]
#> Latitude range: [25.49, 48.826]
#> Saved analysis summary to: /tmp/RtmpJm2r1A/analysis_summary_2016.rds
#> Format: RDS
#> File size: 0.77 KB
#> Results saved to:/tmp/RtmpJm2r1A
#> Files created:6
#>
#> ======================================================================
#> ANALYSIS COMPLETE
#> ======================================================================
#> Processing time:0.07minutes
#> Scale:county
#> Year:2016
#> Nutrients analyzed:nitrogen, phosphorus
#> Spatial units:3112
#> WWTP facilities:30994
#> Output directory:/tmp/RtmpJm2r1A
#> ======================================================================
comparison <- compare_analyses(results_2010, results_2016, "nitrogen")
#>
#> Comparison ofNitrogenClassifications
#> ==================================================
#> Year 1:2010
#> Year 2:2016
#>
#> Classification 2010_Count 2016_Count Absolute_Change
#> Excluded Excluded 193 212 19
#> Sink_Deficit Sink_Deficit 2362 2450 88
#> Sink_Fertilizer Sink_Fertilizer 287 214 -73
#> Source Source 83 83 0
#> Within_County Within_County 187 153 -34
#> Percent_Change
#> Excluded 9.8
#> Sink_Deficit 3.7
#> Sink_Fertilizer -25.4
#> Source 0.0
#> Within_County -18.2
# }