r中的rqpd命令用于运行面板分位数

Question

我想使用rqpd包在r中运行面板分位数来获得固定效果。我试过运行示例

set.seed(10)
m <- 3
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 <- rep(1:n,rep(m,n)) + u 
fit <- rqpd(y ~ x | s, panel(lambda = 5))

但我得到错误：

  

validObject（.Object）中的错误：无效的类“dsparseModelMatrix”   对象：超类＆＃34; mMatrix＆＃34;没有在环境中定义   对象的课程

Answer 1

您缺少为该示例定义n。原作完全没问题。要运行此程序包，我使用R版本2.15.3

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)

Answer 2

您安装了MatrixModels软件包吗？在运行rqpd之前，您需要这样做。