0

我正在使用 R 使用自举标准误差进行分位数回归，以测试一个变量在分布的第 5、第 50 和第 95 个百分位数是否高于第二个变量。输出不包括 t stat 的自由度。我该如何计算？

summary(rq(data$var1~data$var2, tau=.05), se="boot")
summary(rq(data$var1~data$var2, tau=.5), se="boot")
4

1 回答 1

0

假设您使用了 library quantreg，如果您自己调用rq()，您将获得自由度。

看起来您对 SO 很陌生；欢迎来到社区！如果您想快速得到很好的答案，最好让您的问题可重现。这包括您使用的库或样本数据，例如来自 dput(head(dataObject))). 检查一下：制作 R 可重现的问题。

在这种情况下，捕捉自由度应该相对容易。

事实上，观察次数减去观察次数就是总自由度。剩余自由度是观察数减去公式中的变量数。

每个 t 统计量的自由度是代表该 t 统计量的变量数（通常为一个）。

如果您直接调用回归（而不是嵌套在汇总函数中），它也会为您提供有关自由度的信息。话虽如此，如果您不独立运行模型，则更难以测试数据必须满足分析的假设。最后，在这种形式下，您也无法测试模型是否过拟合。

library(quantreg)

data(mtcars)

(fit <- rq(mpg ~ wt, data = mtcars, tau = .05))
# Call:
# rq(formula = mpg ~ wt, tau = 0.05, data = mtcars)
# 
# Coefficients:
# (Intercept)          wt 
#   37.561538   -6.515837 
# 
# Degrees of freedom: 32 total; 30 residual 

(fit2 <- rq(mpg ~ wt, data = mtcars, tau = .5))
# Call:
# rq(formula = mpg ~ wt, tau = 0.5, data = mtcars)
# 
# Coefficients:
# (Intercept)          wt 
#   34.232237   -4.539474 
# 
# Degrees of freedom: 32 total; 30 residual 


summary(fit, se = "boot")
# 
# Call: rq(formula = mpg ~ wt, tau = 0.05, data = mtcars)
# 
# tau: [1] 0.05
# 
# Coefficients:
#             Value    Std. Error t value  Pr(>|t|)
# (Intercept) 37.56154  5.30762    7.07690  0.00000
# wt          -6.51584  1.58456   -4.11208  0.00028 

summary(fit2, se = "boot")
# 
# Call: rq(formula = mpg ~ wt, tau = 0.5, data = mtcars)
# 
# tau: [1] 0.5
# 
# Coefficients:
#             Value    Std. Error t value  Pr(>|t|)
# (Intercept) 34.23224  3.20718   10.67362  0.00000
# wt          -4.53947  1.04645   -4.33798  0.00015

我想指出这se = "boot"似乎没有做任何事情。此外，您可以tau在同一模型中运行这两种设置。该Quantreg软件包有几个工具可以在一起运行时比较模型。

于 2022-01-31T04:29:49.957 回答