Create functions to set/get formula attribute

Many twidlr methods involve setting and getting a formula as an attribute. It's a bit messy in base code. Eg...

```
attr(object, "formula")
```
Instead, create wrapper functions like `set_formula` and `get_formula`. Eg...

```
get_formula <- function(object) attr(object, "formula")
```