Сингулярные матричные ошибки в квантильных регрессиях
Я пытаюсь запустить некоторые квантильные регрессии в R, Переменная ответа является непрерывной. Предсказывающие переменные являются категориальными, принимая значения [0, 1, 2, 3, 4],
Каждый раз, когда я пытаюсь выполнить какой-либо анализ, я получаю сообщения об ошибках (примеры команд и ошибок приведены ниже).
summary(rq(HoursFlex ~ EnvtBenefits, tau = taus), se = "nid")
Error in base::backsolve(r, x, k = k, upper.tri = upper.tri, transpose = transpose, :
singular matrix in 'backsolve'. First zero in diagonal [5]
ggplot(data = s1, aes(EnvtBenefits, HoursFlex)) + geom_point() + geom_quantile(quantiles = taus, col = "gray") + geom_quantile(quantiles = 0.5, col = "blue") + geom_smooth(method="lm", col = 2)
Smoothing formula not specified. Using: y ~ x
Warning messages: 1: Computation failed in stat_quantile(): Singular design matrix
plot(summary(rq(HoursFlex ~ Milieu, tau = taus)))
Error in plot.window(...) : infinite axis extents [GEPretty(-inf,inf,5)]
И в результате, он даже не позволяет мне проверять наличие двухсторонних или трехсторонних взаимодействий между предикторами. Кто-нибудь знает, в чем может быть проблема? Я совсем новичок в R и не знаю, что происходит!!
Редактировать: я включаю небольшое подмножество данных; предикторы кодируются как категориальные с использованием функции as.factor
Непрерывная переменная ответа и выбор категориальных предикторов