> 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))
错误于validObject(.Object) :
invalid class “dsparseModelMatrix” object: superclass "mMatrix" not defined in the environment of the object's class
> sfit <- summary(fit)
错误于summary(fit) :
在为'summary'函数选择方法时评估'object'参数出了错: 错误: 找不到对象'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)
错误于validObject(.Object) :
invalid class “dsparseModelMatrix” object: superclass "mMatrix" not defined in the environment of the object's class
> crem.fit
错误: 找不到对象'crem.fit'
>
为什么会有错误呢?
|