compare aed variables to field?

[hdugan]

Is there a way to calculate RMSE between field and modeled variables, as opposed to using a lakeAnalyzer function?
Example) Compare field data of oxygen to OXY_oxy?

[jordansread]

I think this is worth adding, and my efforts here https://github.com/USGS-R/glmtools/pull/179/files should make it simpler.

But short answer, no this doesn't exist, but I think I can get something workable for you quickly

[jordansread]

@hdugan are you converting your field Oxy data to use the same units as AED output?

[hdugan]

I've been using a homemade resample_to_field function that allows me to input a unit conversion

[jordansread]

ok. It would be nice if resample_to_field wasn't so locked into water temps

[jordansread]

as per convos in our last glmtools/GLMr workshop, it might be worth trying a function that adds the conversion as a new netcdf variable, so all of the _var functions would work.

[jordansread]

with #205 we are most of the way there, but

compare_to_field(nc_file, field_file, nml_file, metric, as_value = FALSE, na.rm = TRUE, ...)

doesn't do anything generic with metric, instead it requires it to be a function that is similar to rLakeAnalyzer or a function in that package. It you give the name of the var as metric, it would be nice to have this just do a straight up pass through of the variable for as_value or RMSE in the units of the variable. The way it is now, there is water.temperature that spoofs a rLA function and just returns the data.frame. Would be nice to do this in the same way for anything that is a variable in the netcdf that also matches the field file.

[lawinslow]

We could default metric to value or something like that that's just a straight pass-through? water.temperature is clearly weird.

[jordansread]

I like using "value". It fits too w/ as_value too