Мне нужно проводить регрессионный анализ с непрерывными и категорированными данными (так как последние распространены в моей области).До сих пор я делал это следующим образом:
q <- lm(ep ~ cut2(exposure, g=5), data=data)
Когда я использую функции из пакета rms, это, похоже, не работает:
q <- ols(ep ~ cut2(exposure, g=5), data=data)
иприводит к сообщению об ошибке (индекс за пределами).Это работает нормально, если я предварительно определяю квантили, но это вызывает проблемы в сценариях, которые я использую.
Есть ли альтернатива?