问题描述
我正在尝试使用模型预测，但遇到了一个问题。
模型
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一对一的预测，然后再组合这些数据帧进行绘图，将花费大量代码和时间。
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）