Можем ли мы использовать функцию "предсказать" для модели Gee?

Question

Можем ли мы использовать функцию "предсказать" для модели Gee? Я знаю из модели gee, что мы можем использовать "fit.values", чтобы получить прогнозируемые значения для данных, используемых для построения модели. Тем не менее, я хочу использовать результаты регрессии Джи, чтобы генерировать прогнозные значения для тестовых данных.

data(warpbreaks)

smp_size <- floor(0.75 * nrow(warpbreaks))
set.seed(123)
train_ind <- sample(seq_len(nrow(warpbreaks)), size = smp_size)
train <- warpbreaks[train_ind, ]
test <- warpbreaks[-train_ind, ]

gee_model = gee(breaks ~ tension, id=wool, data=train, 
corstr="exchangeable")

predict(gee_model, train)
predict(gee_model, test)

Последние две строки приводят к следующей ошибке.

"Error in seq_len(p) : argument must be coercible to non-negative integer
In addition: Warning messages: