悬赏300论坛币解释面板分位数回归命令

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

这是R网站面板分位数回归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 modeldata(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 回帖

关键词：面板分位数回归 面板分位数 分位数回归 0论坛币 分位数 网站

首先，检验是固定效应还是随机效应。cre是随机效应的例子

初衷21 发表于 2016-10-22 16:47
首先，检验是固定效应还是随机效应。cre是随机效应的例子

固定效应啊，方便加QQ指导么？

showree 发表于 2016-10-23 15:30
固定效应啊，方便加QQ指导么？

我也不懂其实，它有固定效应的例子，可以出结果

fit <- rqpd(y ~ x | s, panel(lambda = 5))
sfit <- summary(fit)
固定效应模型  S是个体变量

初衷21 发表于 2016-10-23 16:43
fit

具体操作会出现error

showree 发表于 2016-10-23 17:24
具体操作会出现error

看帮助文件，按照帮助文件可以修改一下