R语言 summary()函数 结果

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求助各位一个问题，本人在使用面板数据分位数回归（rqpd packages）的过程中，使用了summary（）函数来输出结果（这是rqpd: Regression Quantiles for panel data (longitudinal data) in rqpd: Regression Quantiles for Panel Data (rdrr.io)
网站上给出的example的代码）

1. library(rqpd)
   
2. 
   
3. set.seed(10)
   
4. m <- 3
   
5. n <- 10
   
6. s <- as.factor(rep(1:n,rep(m,n)))
   
7. x <- exp(rnorm(n*m))
   
8. u <- x*rnorm(m*n) + (1-x)*rf(m*n,3,3)
   
9. a <- rep(rnorm(n),rep(m,n))
   
10. y <- a + u
    
11. fit <- rqpd(y ~ x | s, panel=panel(lambda = 5))
    
12. **summary输出结果

复制代码

根据网页的说明和之前在论坛里找到的例子看，summary之后输出的结果应该是以下的样子：包含回归结果，标准误，t值和p-value,（图片来源于论坛）
但是，我运行这个程序或者运行本人自己的数据与程序后，进行summary（fit）之后的结果均为

与一般回归的结果在summary（）应用后的返回值不同，请问出现这个问题的原因是什么？我应该怎么样才能看到结果的标准误，t /p 数值。
谢谢大家

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Summary Summa Mary R语言 Mar

那个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内容

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

summary(fit)

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

yangqianmo 发表于 2022-3-13 20:01
不好意思，我换用了我同学的电脑和程序，但是仍然无法按照常见格式输出系数结果，后来在example中我们更改 ...

我有20万行数据的样子，请问您是怎么跑出来的呀，我感觉跑了好久都没出来，超级感谢！

巴卫w 发表于 2022-4-14 16:47
我有20万行数据的样子，请问您是怎么跑出来的呀，我感觉跑了好久都没出来，超级感谢！

数据结果是借了学校的工作站跑出来的， 时间也比较久，我自己的电脑也是跑不出的。可以问问看身边有没有配置比较高的电脑~