请问有人知道在r中那个包可以做面板数据的半参数时变系数，我需要系数岁时间变化的图

sun1018

28909

收藏 2012-04-09

请问有人知道在r中那个包可以做面板数据的半参数时变系数，我需要系数岁时间变化的图

扫码加我拉你入群

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

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

全部回复

ywh19860616

2012-4-9 16:11:43

rqpd_0.5.zip
大小:(173.59 KB)

马上下载

plm 包也许可以满足你的要求
在google输入 plm，R
会出来很多有用信息

扫码加我拉你入群

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

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

sun1018

2012-4-9 16:36:35

ywh19860616 发表于 2012-4-9 16:11
plm 包也许可以满足你的要求
在google输入 plm，R
会出来很多有用信息

这个包我有，的确是做面板的，不过没有我说的模型，这个模型估计要总非参数方法，在非参数方法里倒是有个包，但是那个无法估计时变系数，郁闷

扫码加我拉你入群

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

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

sun1018

2012-4-9 16:40:50

我说的那个包是Adaptive Semiparametic Regression，在里面有个例子好像是面板时变，可惜没有给出程序

扫码加我拉你入群

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

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

ywh19860616

2012-4-9 16:51:37

喔，没有看到你是实现半参数模型的
你可以看下np packages，这个是非参数的，但是不是面板的
现在面板的非参数应该较少

扫码加我拉你入群

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

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

ywh19860616

2012-4-9 16:57:52

sun1018 发表于 2012-4-9 16:40
我说的那个包是Adaptive Semiparametic Regression，在里面有个例子好像是面板时变，可惜没有给出程序

可以联系该作者试试

扫码加我拉你入群

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

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

点击查看更多内容…

sun1018

2012-4-9 17:27:49

ywh19860616 发表于 2012-4-9 16:57
可以联系该作者试试

的确是啊，好像实证的论文只找到两篇，但是都没有给出程序

扫码加我拉你入群

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

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

sun1018

2012-4-9 17:30:52

ywh19860616 发表于 2012-4-9 16:57
可以联系该作者试试

对了，有一个问题请教您，前些日子看到你们在讨论面板数据分位数回归，我看了一下那个Index of /~roger/research/panel，里面的那个程序，为什么在我的电脑里运行出现问题：
错误于cbind(deparse.level, ...) : args have differing numbers of rows
此外: 警告信息：
In as.matrix.csr(x, nrow = nrow, ncol = ncol, eps = eps) :
ncol*nrow indivisable by length(x)
你能直接欲行吗？

扫码加我拉你入群

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

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

ywh19860616

2012-4-10 08:18:13

sun1018 发表于 2012-4-9 17:30
对了，有一个问题请教您，前些日子看到你们在讨论面板数据分位数回归，我看了一下那个Index of /~roger/r ...

请安装rqpd packages，Koenker教授已经把那个函数发展成为了R packages，更为稳定些

扫码加我拉你入群

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

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

ywh19860616

2012-4-10 08:18:57

sun1018 发表于 2012-4-9 17:30
对了，有一个问题请教您，前些日子看到你们在讨论面板数据分位数回归，我看了一下那个Index of /~roger/r ...

这个应该是你数据问题

扫码加我拉你入群

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

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

sun1018

2012-4-10 13:01:38

ywh19860616 发表于 2012-4-10 08:18
这个应该是你数据问题

问题已经解决了，你说的那个package我一直按装不了，我的说版本没有，可是我是最新版本啊

扫码加我拉你入群

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

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

ywh19860616

2012-4-11 18:15:21

sun1018 发表于 2012-4-10 13:01
问题已经解决了，你说的那个package我一直按装不了，我的说版本没有，可是我是最新版本啊

rqpd.rar
大小:(172.08 KB)

马上下载

这个是我电脑上安装好的文件。如果还是不能正常运行，你试下载个R12.2
我在这个版本可以正常使用。

扫码加我拉你入群

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

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

epoh

2012-4-11 21:27:35

sun1018 发表于 2012-4-9 16:36
这个包我有，的确是做面板的，不过没有我说的模型，这个模型估计要总非参数方法，在非参数方法里倒是有个 ...

希望这篇文献对你有帮助.
NONPARAMETRIC AND SEMIPARAMETRIC METHODS IN R
page 28/36 3.12. Panel Data Models
###########
library(plm)
library(Ecdat)
data(Airline)
library(np)
attach(Airline)
lcost <- as.numeric(log(cost))
loutput <- as.numeric(log(output))
lpf <- as.numeric(log(pf))
lf <- as.numeric(lf)
bw <- npregbw(lcost~loutput + lpf + lf + ordered(year) + factor(airline), regtype="ll", bwmethod="cv.aic", ukertype="liracine", okertype="liracine")
summary(bw)

扫码加我拉你入群

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

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

sun1018

2012-5-5 20:20:11

ywh19860616 发表于 2012-4-11 18:15
这个是我电脑上安装好的文件。如果还是不能正常运行，你试下载个R12.2
我在这个版本可以正常使用。
...

这个也不能使用啊，里面应该有个NAMESPACE文件的吧，你给我发的没有，安装后无法使用，您能把您R语言library目录下的程序包给我发一个吗？我下载了三个R的版本都不行，郁闷了，麻烦您在帮我看看，谢谢了，不胜感激额

扫码加我拉你入群

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

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

ywh19860616

2012-5-5 22:05:35

sun1018 发表于 2012-5-5 20:20
这个也不能使用啊，里面应该有个NAMESPACE文件的吧，你给我发的没有，安装后无法使用，您能把您R语言lib ...

噢，那抱歉了，这个我也不知道问题在哪
以前上传的rqpd就是我从安装好的library中解压出来的，在我这里也没有NAMESPACE文件，但是可以正常使用。
和其他包不一样，这个包并未正式发布的，所以是没有NAMESPACE文件。
建议你可以联系下开发包的作者，应该可以解决。

扫码加我拉你入群

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

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

epoh

2012-5-6 09:04:49

sun1018 发表于 2012-5-5 20:20
这个也不能使用啊，里面应该有个NAMESPACE文件的吧，你给我发的没有，安装后无法使用，您能把您R语言lib ...

from R 2.14.0 all packages have a namespace
package 'rqpd': Built: R 2.13.0; 2011-05-06 ; windows
所以你要安装 R version 2.13.0
ywh兄提供的版本可以使用
安装顺序
R version 2.13.0
package 'quantreg'
package  rqpd_0.5.zip
package 'MatrixModels'
package 'Formula'
######
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)
crem.fit

Call:
rqpd(formula = cre.form, panel = panel(method = "cre"), data = bwd)

Coefficients:
(Intercept)[0.25]    smoke[0.25]    dmage[0.25]    agesq[0.25]
   2440.2336143    -114.2131886       -5.5970239       0.4268880
novisit[0.25]    pretri2[0.25]    pretri3[0.25]    m.smoke[0.25]
   -380.3047734    -39.3595390       -5.3371753    -167.8668046
m.dmage[0.25]    m.agesq[0.25]  (Intercept)[0.5]       smoke[0.5]
   55.1567177       -1.1890590    2943.0862058    -118.9682812
   dmage[0.5]       agesq[0.5]    novisit[0.5]    pretri2[0.5]
   -37.6414221       0.9771566    -177.3260213    -23.9477189
   pretri3[0.5]    m.smoke[0.5]    m.dmage[0.5]    m.agesq[0.5]
   -18.4187786    -144.4995484       71.7491954       -1.4448411
(Intercept)[0.75]    smoke[0.75]    dmage[0.75]    agesq[0.75]
   3367.5109878    -120.6370785    -29.1507203       0.8348498
novisit[0.75]    pretri2[0.75]    pretri3[0.75]    m.smoke[0.75]
   -21.2166704       23.6103748    -39.4958273    -138.3875937
m.dmage[0.75]    m.agesq[0.75]
   54.8864519       -1.1462282

Degrees of freedom: 13154 total; 13124 residual

扫码加我拉你入群

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

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

ywh19860616

2012-5-6 09:11:24

epoh 发表于 2012-5-6 09:04
from R 2.14.0 all packages have a namespace
package 'rqpd': Built: R 2.13.0; 2011-05-06 ; windows ...

epoh老师，您好，我还有一个疑问。
这个rqpd包应该是未正式发布的，还在developing，按照作者给出的安装方式：
install.packages("rqpd", repos="http://R-Forge.R-project.org")

我利用这个命令安装成功过，但是现在在R2.14.2安装会出现错误：
警告: 无法在貯藏處http://R-Forge.R-project.org/bin/windows/contrib/2.14中读写索引
警告信息：
In getDependencies(pkgs, dependencies, available, lib) :
package ‘rqpd’ is not available (for R version 2.14.2

似乎这种错误也很常见，您知道原因吗？

扫码加我拉你入群

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

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

epoh

2012-5-6 09:29:46

ywh19860616 发表于 2012-5-6 09:11
epoh老师，您好，我还有一个疑问。
这个rqpd包应该是未正式发布的，还在developing，按照作者给出的安装 ...

R version 2.13.0
也发生一样错误
Warning: unable to access index for repository http://R-Forge.R-project.org/bin/windows/contrib/2.13
Warning message:
In getDependencies(pkgs, dependencies, available, lib) :
package ‘rqpd’ is not available (for R version 2.13.0)

我觉得跟project administrators有关
以rqpd为例
Request to join project rqpd
You can request to join a project by clicking the submit button.
The administrators will be emailed to approve or deny your request.

扫码加我拉你入群

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

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

ywh19860616

2012-5-6 09:31:53

epoh 发表于 2012-5-6 09:29
R version 2.13.0
也发生一样错误
Warning: unable to access index for repository http://R-Forge.R- ...

谢谢epoh老师，那还是最好下载到本地安装。

扫码加我拉你入群

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

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

epoh

2012-5-6 09:42:59

ywh19860616 发表于 2012-5-6 09:11
epoh老师，您好，我还有一个疑问。
这个rqpd包应该是未正式发布的，还在developing，按照作者给出的安装 ...

rqpd : Regression Quantiles for Panel Data
      Build status: Failed to build
Build status codes

0 - Current: the package is available for download. The corresponding package
   passed checks on the Linux and Windows platform without ERRORs.
1 - Scheduled for build: the package has been recognized by the build system
   and provided in the staging area.
2 - Building: the package has been sent to the build machines. It will be built
   and checked using the latest patched version of R. Note that it is included
   in a batch of several packages. Thus, this process will take some time to finish.
3 - Failed to build: the package failed to build or did not pass the checks on the
   Linux and/or Windows platform. It is not made available since it does not meet the policies.
4 - Conflicts: two or more packages of the same name exist. None of them will be built.
   Maintainers are asked to negotiate further actions.
5 - Offline: the package is not available. The build system may be offline or
the package maintainer did not trigger a rebuild (done e.g., via committing to the package repository).

扫码加我拉你入群

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

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

ywh19860616

2012-5-6 09:56:42

epoh 发表于 2012-5-6 09:42
rqpd : Regression Quantiles for Panel Data
Build status: Failed to build
Build statu ...

呵呵，epoh老师，是不是说明这个包未通过审核意思？
5 - Offline: the package is not available. The build system may be offline or
the package maintainer did not trigger a rebuild (done e.g., via committing to the package repository)

这个状态5应该是最新的吧，the package maintainer did not trigger a rebuild
可能Koenker教授没有再修正和更新了。

谢谢epoh老师了。

扫码加我拉你入群

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

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

sun1018

2012-5-6 14:22:47

epoh 发表于 2012-5-6 09:04
from R 2.14.0 all packages have a namespace
package 'rqpd': Built: R 2.13.0; 2011-05-06 ; windows ...

呵呵，epoh老师真的是无所不能啊，谢谢啦，已经可以了，非常感谢

扫码加我拉你入群

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

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

sun1018

2012-5-6 14:23:25

ywh19860616 发表于 2012-5-6 09:31
谢谢epoh老师，那还是最好下载到本地安装。

本地安装可以了，谢谢

扫码加我拉你入群

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

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

sun1018

2012-5-6 15:46:12

epoh 发表于 2012-5-6 09:04
from R 2.14.0 all packages have a namespace
package 'rqpd': Built: R 2.13.0; 2011-05-06 ; windows ...

对啦，epoh老师，您知道R里面有面板数据的BMA变量选择程序吗？，好像BMA只能对时间序列模型变量进行选择

扫码加我拉你入群

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

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

epoh

2012-5-6 17:07:37

sun1018 发表于 2012-5-6 15:46
对啦，epoh老师，您知道R里面有面板数据的BMA变量选择程序吗？，好像BMA只能对时间序列模型变量进行选择

try package 'BMS'

BMS and the Fixed Effects Estimator - A Tutorial

This tutorial illustrates how to use Bayesian Model Averaging (BMA)
with panel data using the R package BMS

扫码加我拉你入群

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

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

sun1018

2012-5-6 20:54:59

epoh 发表于 2012-5-6 17:07
try package 'BMS'

BMS and the Fixed Effects Estimator - A Tutorial

我刚看了一下，是否将各个截面设置为不同的虚拟变量，然后再估计的时候将虚拟变量设为fixed.reg，不知道这样是否正确，请epho老师指教

扫码加我拉你入群

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

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

epoh

2012-5-6 21:09:58

sun1018 发表于 2012-5-6 20:54
我刚看了一下，是否将各个截面设置为不同的虚拟变量，然后再估计的时候将虚拟变量设为fixed.reg，不知道 ...

这里有个范例,看过没?
http://bms.zeugner.eu/blog/2011-05-26/panel-estimation

扫码加我拉你入群

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

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

sun1018

2012-5-6 21:29:19

epoh 发表于 2012-5-6 21:09
这里有个范例,看过没?
http://bms.zeugner.eu/blog/2011-05-26/panel-estimation

我看了一下，好像就是这样的在实际操作过程中，如果有N个截面，设置N-1个虚拟变量，就是例子中的countryDummies，然后BMS运算是将N-1个虚拟变量设为fixed.reg，不知道我的想法对不对啊，请指教

扫码加我拉你入群

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

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

sun1018

2012-5-6 22:08:48

epoh 发表于 2012-5-6 21:09
这里有个范例,看过没?
http://bms.zeugner.eu/blog/2011-05-26/panel-estimation

没看仔细，呵呵，我再看看，让你见笑了，好像不完全是

扫码加我拉你入群

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

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

sun1018

2012-5-9 15:50:22

ywh19860616 发表于 2012-4-10 08:18
这个应该是你数据问题

ywh19860616 老师，向您请教个问题，我今天看了几篇用面板数据分位数回归做的文章，比如文章有n个截面，作者就是添加n-1个虚拟变量，然后直接对数据用EVIEW或者STATA进行估计，再看了Quantile regression for longitudinal data后，觉得这样做是错误的，但是说不出来哪里有问题，希望您能指教

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群