Я запустил следующий скрипт для анализа непараметрической квантильной регрессии с использованием пакета quantreg
:
library (quantreg)
model.q = rq(Y_var ~ X1_var + X2_var + X3_var,
data = wideformat_data,
tau = 0.5)
summary(model.q)
model.nullq = rq(Y_var ~ 1,
data = wideformat_data,
tau = 0.5)
anova(model.q, model.nullq)
Все переменные Y и X были непрерывными и стандартизированными. В результате от anova ( )
появилось предупреждающее сообщение:
Warning:
In summary.rq(x, se = se, covariance = TRUE) : 2 non-positive fis
Кто-нибудь знает, почему это произошло и как с этим справиться?