Извлечение формулы из модели в R

Я создаю функцию для многих типов моделей, которая должна извлечь формулу, используемую для создания модели. Есть ли гибкий способ сделать это? Например:

x <- rnorm(10)
y <- rnorm(10)
z <- rnorm(10)
equation <- z ~ x + y
model <- lm(equation)

Мне нужно извлечь формулу объекта «уравнение» после передачи модели.