全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
6821 4
2022-03-09
求助各位一个问题,本人在使用面板数据分位数回归(rqpd packages)的过程中,使用了summary()函数来输出结果(这是rqpd: Regression Quantiles for panel data (longitudinal data) in rqpd: Regression Quantiles for Panel Data (rdrr.io)
网站上给出的example的代码)
复制代码
根据网页的说明和之前在论坛里找到的例子看,summary之后输出的结果应该是以下的样子:包含回归结果,标准误,t值和p-value,(图片来源于论坛) 1.png
但是,我运行这个程序或者运行本人自己的数据与程序后,进行summary(fit)之后的结果均为
2.png
与一般回归的结果在summary()应用后的返回值不同,请问出现这个问题的原因是什么?我应该怎么样才能看到结果的标准误,t /p 数值。
谢谢大家

二维码

扫码加我 拉你入群

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

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

全部回复
2022-3-9 22:37:58
那个example的代码运行后,summary结果是正常的。

summary(fit)

Call: rqpd(formula = y ~ x | s, panel = panel(lambda = 5))

taus: [1] 0.25 0.50 0.75

tau weights: [1] 0.25 0.50 0.25

Coefficients:
                  Value    Std. Error t value  Pr(>|t|)
(Intercept)[0.25]  2.16738  1.20475    1.79902  0.08460
x[0.25]           -3.52922  2.09202   -1.68699  0.10456
(Intercept)[0.5]   3.22301  1.33314    2.41760  0.02358
x[0.5]            -3.10681  1.20002   -2.58897  0.01610
(Intercept)[0.75]  4.51892  1.54618    2.92263  0.00745
x[0.75]           -1.89637  1.60018   -1.18509  0.24758

说明summary()没问题。

为什么你换了数据或模型,summary出不来,原因不清楚。
但是你可以用fit\$coef试试fit的系数结果,或者直接点fit查看list内容
二维码

扫码加我 拉你入群

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

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

2022-3-13 20:01:19
llb_321 发表于 2022-3-9 22:37
那个example的代码运行后,summary结果是正常的。

summary(fit)
不好意思,我换用了我同学的电脑和程序,但是仍然无法按照常见格式输出系数结果,后来在example中我们更改命令为summary.rqpd(fit)解决了这个问题。
但是由于我自己的数据有430万,最终输出的fit$coefficient矩阵长度为16万,所以仍然没有办法通过更改命令获得我想要的结果,不知道您是否对数据量较大的情况下的输出有了解。   
二维码

扫码加我 拉你入群

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

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

2022-4-14 16:47:47
yangqianmo 发表于 2022-3-13 20:01
不好意思,我换用了我同学的电脑和程序,但是仍然无法按照常见格式输出系数结果,后来在example中我们更改 ...
我有20万行数据的样子,请问您是怎么跑出来的呀,我感觉跑了好久都没出来,超级感谢!
二维码

扫码加我 拉你入群

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

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

2022-4-18 20:24:25
巴卫w 发表于 2022-4-14 16:47
我有20万行数据的样子,请问您是怎么跑出来的呀,我感觉跑了好久都没出来,超级感谢!
数据结果是借了学校的工作站跑出来的, 时间也比较久,我自己的电脑也是跑不出的。可以问问看身边有没有配置比较高的电脑~
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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