有人会用R做面板分位数回归吗？会做例子不会套用，这个是例子，大家一起探讨

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

有没有人会用R做面板分位数回归吗？我已经安装上rqpd安装包了，也看见R中rqpd的实例了，做例子能运行出来，但是不太能看懂，m、n代表什么，套用不上自己的
Examples

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)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

好像就用后边这两句就行，把y和x换成你自己的变量名称，然后data也用自己的
fit <- rqpd(y ~ x | s, panel(lambda = 5))
sfit <- summary(fit)
我也是最近才整明白点，不知道对不对

请问rqpd怎么安装上去的，我一直安装不上，894664184，方便加qq讨论一下嘛？