楼主: aalavender
2934 16

[程序分享] 1000金币求助R语言进行面板分位数回归 [推广有奖]

  • 0关注
  • 0粉丝

高中生

25%

还不是VIP/贵宾

-

威望
0
论坛币
163 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
254 点
帖子
16
精华
0
在线时间
36 小时
注册时间
2005-4-21
最后登录
2018-4-17

楼主
aalavender 发表于 2017-8-29 19:24:44 来自手机 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
1000金币求助R语言进行面板分位数回归
二维码

扫码加我 拉你入群

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

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

关键词:面板分位数回归 面板分位数 求助R语言 分位数回归 分位数

沙发
deem 学生认证  发表于 2017-8-29 19:37:09
  1. install.packages("rqpd", repos="http://R-Forge.R-project.org")
复制代码
已有 1 人评分论坛币 收起 理由
admin_kefu + 10 热心帮助其他会员

总评分: 论坛币 + 10   查看全部评分

藤椅
aalavender 发表于 2017-8-29 19:43:42
deem 发表于 2017-8-29 19:37
你好,这个代码我知道,我找见r里的例子了,就是不知道每个字母的意思,例子的代码能运行出来,就是自己的数据带进去就不行

板凳
aalavender 发表于 2017-8-29 19:51:01
set.seed(5)
m<-5
n<-5
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)
cre.form<-cons~f+m+incom+assert+novisit+pretri2+pretri3|momid3|f+m+incom+assert
crem.fit<-rqpd(cre.form,panel(method="cre"),data=d1)
crem.fit
我想知道m,n,seed括号里的数字都是什么意思,cre.form那个模型该怎么建立

报纸
丁兆海4 在职认证  学生认证  发表于 2017-8-30 10:17:05
我一般不跟别人说的,就看在金币的份上:
链接一:http://www.xiaowanxue.com/up_files/2012121819040.html#_ftn1
链接二:http://www.cnblogs.com/shangfr/p/5736738.html
别忘了,金币!!!
已有 1 人评分经验 论坛币 收起 理由
admin_kefu + 20 + 10 热心帮助其他会员

总评分: 经验 + 20  论坛币 + 10   查看全部评分

地板
deem 学生认证  发表于 2017-8-30 10:50:48
aalavender 发表于 2017-8-29 19:51
set.seed(5)
m
就是simulate出来一个面板数据,产生n个个体,每个个体有m个观测。rqpd之前的都是产生数据的过程,set.seed就是产生一个随机数种子,让别人可以reproduce随机模拟的结果,数字可以任意正整数。cre.form看文档呀

7
aalavender 发表于 2017-9-5 09:32:13
deem 发表于 2017-8-30 10:50
就是simulate出来一个面板数据,产生n个个体,每个个体有m个观测。rqpd之前的都是产生数据的过程,set.se ...
大神,那我已经整理好面板数据了,怎么能直接把我的数据带到代码里运行

8
aalavender 发表于 2017-9-5 09:33:41
丁兆海4 发表于 2017-8-30 10:17
我一般不跟别人说的,就看在金币的份上:
链接一:http://www.xiaowanxue.com/up_files/2012121819040.htm ...
你好,第一个链接那个我有,第二个链接打不开啊

9
deem 学生认证  发表于 2017-9-5 09:52:31
aalavender 发表于 2017-9-5 09:32
大神,那我已经整理好面板数据了,怎么能直接把我的数据带到代码里运行
文档都有:
CRE:        y ~ x1 + x2 + ... | s | z1 + z2 + ...

s:  id variable
z: CRE component. These are possibly endogenous variables (in the sense that they are affected by the fixed effects) and must be time-varying

10
丁兆海4 在职认证  学生认证  发表于 2017-9-5 13:00:52
aalavender 发表于 2017-9-5 09:33
你好,第一个链接那个我有,第二个链接打不开啊
但是,我可以打开。

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-22 03:38