When simulation method is "constant" across qE, should do burn-in for each
rBatt opened this issue · comments
Right now the burn in is only at the beginning of the time series. If you want to supply several values for qE, and then have the model simulate dynamics at each constant value of qE, then you'd want it to do a burn in for each value ... but right now it only does it for the first.
This extra burn in is NOT desired for the linear case (? is it ?), because there would be a full burn in for each time step.
https://github.com/rBatt/edmShift/blob/master/R/functions/simFuns/FWsim.wrap.R#L37