如何对 rqpd 对象的标准错误进行聚类?该功能summary.rqpd
似乎不允许“集群”选项。
最小工作示例:
# install.packages("rqpd", repos="http://R-Forge.R-project.org")
library(rqpd)
x <- 1:50
y <- c(x[1:48] + rnorm(48, 0, 5), rnorm(2, 150, 5))
cluster <- rep(letters[1:5], 10)
QR <- rqpd(y ~ x|cluster, panel(lambda=5))
#The function does not have a "cluster" option
summary.rqpd(QR, se='boot')