Alison Appling
Alison Appling
Hi Tzu-Yao, are you looking at this text?  If so, I think the needed clarification may be that we fixed process error (`err_proc_iid_sigma`) in the _simulations_ that generated DO...
See this [issue + answer](https://github.com/USGS-R/streamMetabolizer/issues/372#issuecomment-420134665) for a way to fix K to a single number. Basically, you should pick a distribution for K600_daily_meanlog that is very narrow (tiny meanlog_sdlog) and...
Hi Sam, it's a bit tricker to set K to multiple numbers, but it can be done. If you have distinct Q on each day or are willing to spoof...
By multi-node relationship I meant a piecewise linear relationship between K and a variable named Q or discharge, which is what's provided by a `b_Kb*_...` model (`pool_K600='binned*'` in `mm_name()`). See...
Thanks for sharing this code, Sam. I'm only seeing 9 values in sb4_lnK.csv, 8 of which get read by read.csv as non-header rows, but K600_lnQ_nodes_centers has length 10...is it possible...
Also, are you sure K is as high as you say, and in the right units? Those DO predictions are abysmal, and I'm surprised the model couldn't at least find...
Summary of the renaming I just did: * `frac_GPP` is now `light_mult_GPP` * `frac_ER` is now `const_mult_ER` * ~~`coef_GPP` is now `pp_mult_GPP`~~ it was briefly, but now it doesn't have...
All this code is in #376 , which resolves this issue unless we notice something amiss in the next few days.
may want to plot both DO_mod (the state, which is DO.obs for pi models) and DO_pure (aka DO_det for 'deterministic'?). And what about DO_mod_partial, which is deterministic based on the...
This function could also use updating for GPP-process-error models (`pp`, `err_proc_GPP=TRUE`).