scripts/sam_models/mod2 contains the latest antecedent model version. The code should run once gaps in the Q (verona), solar radiation, and air temp data are filled. The model now utilizes the variable past_topped to account for the number of days in the past 30 when there was inundation. In the input data, the vectors pA, pB, and pC indicate the number of days in each time step. For Qant, there are currently 5 timesteps of 1, 3, 3, 7, and 7 days, which can be experimented with. Also the current model does not include two-way interactions, but those can be added easily.
Because the chlorophyll dataset now includes 10 sites, I added a intercept term to account for site random effects. If there is a strict need for the REs to sum to zero, I can provide a code snipped to that effect.
scripts/sam_models/mod2contains the latest antecedent model version. The code should run once gaps in the Q (verona), solar radiation, and air temp data are filled. The model now utilizes the variablepast_toppedto account for the number of days in the past 30 when there was inundation. In the input data, the vectors pA, pB, and pC indicate the number of days in each time step. For Qant, there are currently 5 timesteps of 1, 3, 3, 7, and 7 days, which can be experimented with. Also the current model does not include two-way interactions, but those can be added easily.Because the chlorophyll dataset now includes 10 sites, I added a intercept term to account for site random effects. If there is a strict need for the REs to sum to zero, I can provide a code snipped to that effect.