全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2012-5-9 16:11:16
sun1018 发表于 2012-5-9 15:50
ywh19860616 老师,向您请教个问题,我今天看了几篇用面板数据分位数回归做的文章,比如文章有n个截面,作 ...
较早的做法都是借鉴一般面板数据模型中的reg+dummy,就是说pooled数据+虚拟变量等价于FE模型。
个人觉得比较正统的还是Quantile regression for longitudinal data一文的penalized quantile regression,因为他同时推到了渐进性,目前国外基本都是这样做。
二维码

扫码加我 拉你入群

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

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

2012-5-9 16:17:31
ywh19860616 发表于 2012-5-9 16:11
较早的做法都是借鉴一般面板数据模型中的reg+dummy,就是说pooled数据+虚拟变量等价于FE模型。
个人觉得 ...
对于面板数据估计方法的话,加虚拟变量是没有问题的,但是对于分位数回归的损失函数(加权绝对离差最小),用虚拟变量估计跟Quantile regression for longitudinal data提出的完全不一样啊,为什么还有人用呢?
二维码

扫码加我 拉你入群

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

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

2012-5-9 16:28:48

如果

sun1018 发表于 2012-5-9 16:17
对于面板数据估计方法的话,加虚拟变量是没有问题的,但是对于分位数回归的损失函数(加权绝对离差最小) ...
呵呵,你看别的面板分位的文章,应该很少看到别人写过程吧,所以作者让读者猜到底是用哪种方法?其实用了哪种我们也不太清楚,不知道他们是如何处理的。

ps:如果要用,建议你用penalized quantile,不必去深思别人的方法。其实方法都是一步一步发展的,没有绝对的正确。比如别人如果不知道有penalized方法(2004年才有),或者2004年以前的文章,他们利用reg+dummy,那么也就不足为奇。
二维码

扫码加我 拉你入群

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

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

2012-5-9 16:52:56
ywh19860616 发表于 2012-5-9 16:28
呵呵,你看别的面板分位的文章,应该很少看到别人写过程吧,所以作者让读者猜到底是用哪种方法?其实用了 ...
倒也是,呵呵。谢谢
二维码

扫码加我 拉你入群

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

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

2012-6-8 15:54:22
你好啊,前辈,我现在也在做跟你一个情况的论文,请问你这方面的资料能给我下吗?看了好多没头绪啊,谢谢
二维码

扫码加我 拉你入群

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

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

2012-6-9 15:56:26
K线人 发表于 2012-6-8 15:54
你好啊,前辈,我现在也在做跟你一个情况的论文,请问你这方面的资料能给我下吗?看了好多没头绪啊,谢谢
哪方面啊》你也没说清楚啊
二维码

扫码加我 拉你入群

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

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

2012-6-9 19:18:18
友情帮顶
二维码

扫码加我 拉你入群

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

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

2012-6-9 19:58:58
sun1018 发表于 2012-6-9 15:56
哪方面啊》你也没说清楚啊
嗯,不好意思。是这样的,我这里要处理这样一个问题,1个因变量Y,5个变量X,30个省市,20年的数据。想看这5个变量20年里系数的变化,求教处理方法,感觉跟你提问的内容相近,谢谢啊
二维码

扫码加我 拉你入群

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

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

2012-6-17 14:04:13
K线人 发表于 2012-6-9 19:58
嗯,不好意思。是这样的,我这里要处理这样一个问题,1个因变量Y,5个变量X,30个省市,20年的数据。想看 ...
你看看前面的就知道了,我们最终也没有实现啊,因为没程序啊
二维码

扫码加我 拉你入群

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

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

2012-6-17 14:05:34
ywh19860616 发表于 2012-5-9 16:11
较早的做法都是借鉴一般面板数据模型中的reg+dummy,就是说pooled数据+虚拟变量等价于FE模型。
个人觉得 ...
前辈请教您一下,面板分位数回归的伪R2是怎么实现的啊?
二维码

扫码加我 拉你入群

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

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

2012-9-13 17:25:07
真是感谢好心人
二维码

扫码加我 拉你入群

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

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

2012-9-29 17:06:48
收藏。
二维码

扫码加我 拉你入群

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

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

2012-10-21 19:20:34
学习了
二维码

扫码加我 拉你入群

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

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

2012-10-22 10:29:28
> 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'
>
为什么会有错误呢?
二维码

扫码加我 拉你入群

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

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

2013-4-24 10:29:45
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老师您好,我用的是R3.0.0,按照上面的顺序安装但是一调用就出现“程序包是在R版本的3.0.0建立的,您的重新安装”但我没有找到R.2.14等低版本的,老师能不能发一个给我,我的邮箱liu_xinwen@126.com,由于初学,还不懂,非常感谢!
二维码

扫码加我 拉你入群

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

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

2013-4-24 10:32:08
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”是在R版本的3.0.0建立的,您的重新安装”但我没有找到R.2.14等低版本的,谢谢老师!
二维码

扫码加我 拉你入群

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

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

2014-7-27 13:48:18
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 ...
请教一个问题:我安装了rqpd包,但是运行同样的语句,出现了下面提示:
Error in UseMethod("determinant") :
  no applicable method for 'determinant' applied to an object of class "matrix.csr"
请问这是怎么回事呢?语句肯定没错的,我直接复制的,而且运行别的也是一样。能否帮忙解答一下,非常感激。
二维码

扫码加我 拉你入群

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

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

2014-7-29 11:00:40
多谢,学习了。
二维码

扫码加我 拉你入群

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

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

2014-8-3 11:35:04
glxfgsh 发表于 2014-7-27 13:48
请教一个问题:我安装了rqpd包,但是运行同样的语句,出现了下面提示:
Error in UseMethod("determinan ...
我也遇到了同样的问题,试了各种数据都是这个错误,不知道你的解决了么?
二维码

扫码加我 拉你入群

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

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

2014-8-3 16:25:56
yangge744 发表于 2014-8-3 11:35
我也遇到了同样的问题,试了各种数据都是这个错误,不知道你的解决了么?
没有,如果你解决了,麻烦也跟我说一声,谢谢
二维码

扫码加我 拉你入群

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

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

2014-8-6 09:49:35
glxfgsh 发表于 2014-8-3 16:25
没有,如果你解决了,麻烦也跟我说一声,谢谢
出现这个问题是因为安装的是source包,就是那个tar.gz的,不管是windows还是Ubuntu下都是这个错误。换成前面ywh19860616给出的rqpd.zip就好了,泪奔啊~~o(>_<)o ~~
二维码

扫码加我 拉你入群

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

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

2014-8-6 17:30:09
glxfgsh 发表于 2014-8-3 16:25
没有,如果你解决了,麻烦也跟我说一声,谢谢
发现ywh19860616提供的包只能做cre,做不了固定效应的,会报错:
错误于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)
后来看了代码发现是有一处cbind应该改为as.matrix(对比了Roger Koenker的2004年源程序),这样就可以做固定效应了,另外det函数那段,会出之前det(),matrix.csr的错误,不知道怎么改合适,我把它注释掉了,那是个判断语句,只要保证数据没问题应该就好了~
二维码

扫码加我 拉你入群

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

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

2014-8-8 00:12:12
yangge744 发表于 2014-8-6 17:30
发现ywh19860616提供的包只能做cre,做不了固定效应的,会报错:
错误于cbind(deparse.level, ...) : ar ...
非常感谢回复,我试一下看看
二维码

扫码加我 拉你入群

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

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

2014-10-10 20:05:25
给顶起来。。。。。
二维码

扫码加我 拉你入群

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

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

2014-10-12 22:33:39
学习了~~~有用
二维码

扫码加我 拉你入群

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

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

2014-10-13 12:18:58
非常感谢
二维码

扫码加我 拉你入群

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

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

2015-3-2 15:12:30
大神的光辉照得我一愣一愣{:3_41:}
二维码

扫码加我 拉你入群

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

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

2015-12-20 18:09:15
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 ...
thanks赶紧下按顺序下载试试
二维码

扫码加我 拉你入群

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

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

2016-1-1 11:46:57
yangge744 发表于 2014-8-6 17:30
发现ywh19860616提供的包只能做cre,做不了固定效应的,会报错:
错误于cbind(deparse.level, ...) : ar ...
谢谢分享,这方法能做固定效应真赞
二维码

扫码加我 拉你入群

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

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

2016-2-27 19:40:42
请问您做出来了吗,我最近也在纠结这个问题,不会求函数型的系数。求帮忙,必有重谢
二维码

扫码加我 拉你入群

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

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

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群