如何编写具有不同行数的newdata参数？

我正在尝试使用模型预测，但遇到了一个问题。

模型

brm(postacute_therapy ~ (year | town))

我想对所有9年和15个城镇进行预测。

编码newdata参数

newdata = data.frame(year = c(2009,2010,2011,2012,2013,2014,2015,2016,2017),town= c("A","B","C","D","E","F","G","H","I","J","K","L","M","N","F"))

data.frame中的错误（year_centered = c（2009，2010，2011，2012，2013，2014，2015，2016，2017））：参数表示行数不同：9，15

难以理解的是9和15给出的行数不相等。但是，是否有一种简单的方法可以通过更聪明的方式编码 newdata参数来绕过该问题？进行9 * 15一对一的预测，然后再组合这些数据帧进行绘图，将花费大量代码和时间。