############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Loading required package: RMySQL Loading required package: DBI Loading required package: fields Loading required package: spam Loading required package: dotCall64 Loading required package: grid Spam version 2.6-0 (2020-12-14) is loaded. Type 'help( Spam)' or 'demo( spam)' for a short introduction and overview of this package. Help for individual functions is also obtained by adding the suffix '.spam' to the function name, e.g. 'help( chol.spam)'. Attaching package: ‘spam’ The following objects are masked from ‘package:base’: backsolve, forwardsolve See https://github.com/NCAR/Fields for an extensive vignette, other supplements and source code [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.PM25_TOT_ob, d.PM25_TOT_mod ,d.POCode,s.state from ARDEQ_JUL_TRIAL as d, site_metadata as s WHERE d.PM25_TOT_ob is not NULL and d.network='IMPROVE' and s.stat_id=d.stat_id and d.ob_dates BETWEEN 20160701 and 20160731 and d.ob_datee BETWEEN 20160701 and 20160731 and ob_hour between 0 and 23 ORDER BY stat_id,ob_dates,ob_hour" There were 11 warnings (use warnings() to see them) null device 1 > proc.time() user system elapsed 3.205 0.228 2.527