Regresión por cuantiles en R

​

-----

 

# install.packages("quantreg")
library(quantreg)


mydata<- read.csv("C:/Users/minombre/mifolder/miarchivo.csv")
attach(mydata)

​

# Define variables
names(mydata)

​

Y <- cbind(Y)
X <- cbind(X1, X2, X3, X4)

​

# Estadísticas descriptivas
summary(Y)
summary(X)

​

# Regresión OLS
ols<- lm(Y ~ X, data=mydata)
summary(ols)

​

# Regresión para cada cuantil
qreg25 <- rq(Y ~ X, data=mydata, tau=0.25)
summary(qreg25)
qreg50 <- rq(Y ~ X, data=mydata, tau=0.5)
summary(qreg50)
qreg75 <- rq(Y ~ X, data=mydata, tau=0.75)
summary(qreg75)

​

# Regresión simultánea para todos los cuantiles
qregall <- rq(Y ~ X, data=mydata, tau=c(0.25, 0.50, 0.75))
summary(qregall)

​

# Prueba de diferencias en los coeficientes con ANOVA
anova(qreg25, qreg50, qreg75)

​

# Gráficos
qreg.all <- rq(Y ~ X, tau = seq(0.05, 0.95, by = 0.05), data=mydata)
qreg.plot <- summary(qreg.all)
plot(qreg.plot)

​

-----

 

Material de apoyo: