Могу ли я получить верхние / нижние границы при запуске квантильной регрессии на панели с использованием пакета rqpd?
Могу ли я получить верхние / нижние границы для моих коэффициентов при запуске квантильной регрессии на панели с использованием пакета rqpd?
Ниже приведена настройка (будьте добры, я новичок в R):
s = as.factor(rep(1:30,rep(63,30)))
s = as.matrix(s)
fit = rqpd(formula = y ~ x | s, panel = panel(lambda = 1, taus=c(0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95), tauw=c(0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95)), data = pdata, control = list(tmpmax = 500))
capture.output(summary(fit), file = "QR.txt")
Код работает хорошо, и вывод предоставляет мне обычные, коэффициенты, t-значения, p-значения и так далее. Но я хочу выяснить верхнюю и нижнюю границы для моей SE для 95% значимости.
Используя rq из пакета quantreg, я могу извлечь верхние / нижние границы, захватив выходные данные из:
fit = rq(y~z, tau = taeus)
plot.fit = summary(fit)
capture.output(plot.fit, file = "BOUNDS.txt")
Однако это не применимо для панельных данных, которые я сейчас просматриваю. У вас, ребята, есть указания?