Multiple instances of rw1 and rw2 in formula

[helske]

Currently, all time-varying covariates should be included in one rw1 and rw2 calls, with identical priors. This can be problematic in some cases. It would be relatively straightforward to parse multiple rw1 calls during the formula manipulation and then combine them before sending stuff to Stan side.