library(rqpd)
set.seed(10)
m <- 3
n <- 10
s <- as.factor(rep(1:n,rep(m,n)))
x <- exp(rnorm(n*m))
u <- x*rnorm(m*n) + (1-x)*rf(m*n,3,3)
a <- rep(rnorm(n),rep(m,n))
y <- a + u
fit <- rqpd(y ~ x | s, panel(lambda = 5))
sfit <- summary(fit)

# A CRE model
data(bwd)

cre.form <- dbirwt ~ smoke + dmage + agesq +
   novisit + pretri2 + pretri3 | momid3 | smoke +
   dmage + agesq

# CRE-M type fit:
crem.fit <- rqpd(cre.form, panel(method="cre"), data=bwd)

# AD type fit:
ad.fit <- rqpd(cre.form, panel(method="cre", cre="ad"), data=bwd,
  subset=bwd$idx %in% 1:2)

你好，这个你能运行出结果吗？我运行的是结果是
Error in validObject(.Object) :
  invalid class “dsparseModelMatrix” object: superclass "dCsparseMatrix" not defined in the environment of the object's class

can you show me your code? the example I gave you, I have no problem with running it.

你好，我运行结果是这样