r语言能实现截面数据的分位数回归吗

刚开始学，想问一下截面数据能做分位数回归吗？用r能实现吗？还是说要用stata？能实现的话代码难不难啊？

用quantreg包中的rq()函数实现

# 加载 quantreg 包
library(quantreg)

# 创建示例数据
set.seed(123)
data <- data.frame(
  y = rnorm(100),
  x1 = rnorm(100),
  x2 = rnorm(100),
  x3 = rnorm(100)
)

# 分位数回归，取第50%分位数
fit <- rq(y ~ x1 + x2 + x3, data = data, tau = 0.5)

# 输出回归结果
summary(fit)