René's URL Explorer Experiment

{"@context":"https://schema.org","@type":"DiscussionForumPosting","headline":"app3 ","articleBody":"In the rush to get rloadest released, some applications did not really get completed. App3 is a case in point--the model is not really well resolved. Below are my suggestions for improving the model.\r\n\r\nFigure 6 should be better commented and refined to better show the relation between flow and the timing of residuals. I suggest the code for figure 6 should be modified to:\r\n\r\n\\# setSweave is required for the vignette.\r\nsetSweave(\"app3_06\", 5, 5)\r\nAA.pl \u003c- timePlot(app3.calib$Date, residuals(app3.lra),\r\n         Plot=list(what=\"points\"), ytitle=\"Residuals\")\r\nrefLine(horizontal=0)\r\n\\# show water years\r\nrefLine(vertical=as.Date(\"2001-10-01\") + years(0:9))\r\n\\# show low flow values in red\r\npicks \u003c- app3.calib$Flow \u003c median(app3.calib$Flow)\r\naddXY(app3.calib$Date[picks], residuals(app3.lra)[picks], Plot=list(what=\"points\", color=\"red\"))\r\ngraphics.off()\r\n\r\nThe first sentence in the describing paragraph should be expanded from:\r\nThe residuals for water-year 2002 are mostly greater than 0; those for watery-year 2003 are\r\ntrending down and those for water-year 2004 are mostly less than 0.\r\nto:\r\nThe residuals for water-year 2002 are mostly greater than 0 and lower than normal flow; those for water-year 2003 are trending down and those for water-year 2004 are mostly less than 0 and represent greater than normal flow.\r\n\r\nThe diagnostic plots should for the anomaly model should include HFV, which indicates that a second order fit is suggested for HFV. The improved model is:\r\napp3.lra \u003c- loadReg(OrthoPhosphate.PO4 ~ a1yr + quadratic(HFV) + dectime(Date)\r\n                    + fourier(Date),\r\n                    data = app3.calib, \r\n                    flow = \"Flow\", dates = \"Date\",\r\n                    station=\"Potomac River at Chain Bridge, at Washington, DC\")\r\n\r\nI got Bp down to about 13% and a much better E value.\r\n","author":{"url":"https://github.com/DLLORENZ","@type":"Person","name":"DLLORENZ"},"datePublished":"2019-06-15T15:13:32.000Z","interactionStatistic":{"@type":"InteractionCounter","interactionType":"https://schema.org/CommentAction","userInteractionCount":0},"url":"https://github.com/29/rloadest/issues/29"}