问题标签 [quantile-regression]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
0 回答
62 浏览

r - 是否有一个函数可以模拟来自稳健线性回归或分位数回归模型的数据?

是否有任何函数可以模拟来自稳健线性模型或分位数回归模型的响应变量,例如线性模型(即stats::simulate.lm)?

如果没有,有没有办法调整代码来为任一模型执行此操作?

这是我正在处理的数据和模型的示例:

0 投票
0 回答
117 浏览

r - 如何运行分位数回归,在每个循环中添加额外四分之一的数据

我想计算一系列样本外预测。我需要运行一个分位数回归,它使用前 20 个季度的数据并生成预测,存储该系列的最后一个观察值,即预测,然后使用 21 个季度的数据重新运行分位数回归以产生新的预测下一季度,然后是 22 个季度的数据,依此类推。

我尝试使用循环来扩展 ts 数据的“窗口”,但我无法让循​​环工作。

然后我尝试使用 xts 数据,但我遇到了同样的问题

我希望 dyhat_qt_oos 的行不同,但它们都相同

0 投票
2 回答
354 浏览

r - 为时间序列的每个月运行回归 - apply.monthly 不起作用

我想对一年时间序列的每个月应用分位数回归(一个因变量和一个自变量),因此我将收到 12 个系数。我的数据集由给出,return_2000_xts并且rq()是分位数回归的函数。我的数据集是xts格式的,包括银行股票的每日回报。

我尝试使用apply.monthly()

不幸的是,我收到以下错误消息:

获取错误(as.character(FUN),模式=“功能”,环境=环境):找不到模式“功能”的对象“乐趣”

我的代码可能有什么问题?

0 投票
0 回答
1189 浏览

r - 如何估计 R 中的 Quantile-on-Quantile 回归?

我需要使用 Sim 和 Zhou (2015) [Sim, N., Zhou, H., 2015. 石油价格、美国股票收益和他们的分位数。银行与金融杂志 55, 1-8]。

我有两个变量 X 和 Y,它们代表增长率。我估计了分位数回归,但未能估计出 QQ 回归。尝试使用两个包:quantreg 和 quantreg.nonpar。

分位数回归通过 估计quantreg,但不适用于quantreg.nonpar(执行piv.bsp <- npqr(formula = form.par, data = data.growth,...导致错误:)Error in rq.fit.br(x, y, tau = tau, ...) : Singular design matrix。另外,我不确定如何指定分位数回归...

0 投票
1 回答
455 浏览

python - ValueError:在乘法运算中使用占位符变量时不支持无值

我在张量流的网络图中使用以下代码:

但是,由于乘法“不支持无值”,我得到一个错误。

谁能告诉我如何使用我的占位符正确执行乘法?

0 投票
0 回答
112 浏览

r - 使用引导、分类变量的分位数回归的置信区间

我正在尝试获取分位数回归的引导置信区间。我的回归响应因子主要是分类变量。当我运行代码 boot.rq 时,出现以下错误;

boot.rq.xy(x, y, U, tau) 中 coercionError 引入的 NA:外部函数调用中的 NA/NaN/Inf (arg 7)

下面是一个例子,我想计算类型的置信区间;

插图来自该论坛之前的讨论;

使用 bootstrap 进行分位数回归的置信区间

0 投票
1 回答
192 浏览

r - ggplot中带有虚拟变量的分位数回归误差

所以我试图计算分位数回归并绘制结果由于ggplot 某种原因在绘制结果时由于某种原因无法显示虚拟变量ggplot

示例mtcars数据集上的代码,如果如下所示:

当使用一切工作正常绘制结果时plot(summary(QR.2, se = "boot")),但由于某种原因使用 ggplot 显示错误。

0 投票
0 回答
370 浏览

r - 如何在 R 中运行分段分位数回归?

我正在尝试使用 R 中的分位数回归将“约束线”拟合到双变量散点图。我的响应变量 ( WaterEr) 代表年度土壤流失,我的自变量代表 ( RainfallEr) 降雨侵蚀性(见下图)。 在此处输入图像描述

我想要达到的结果是基于Hao 等人的一系列出版物。2017 年,其中 x 轴上的值均分为 100 个部分(或箱)以创建 100 列(见下文)。然后,他们使用分段分位数回归来定义与列数一样多的边界点来拟合约束线。有关信息,他们使用 Origin 9(美国 OriginLabs)来运行这些计算。 在此处输入图像描述

由于响应似乎不是线性的,我尝试使用包中的nlrq()函数应用非线性分位数回归quantreg::。这是我应用的代码示例:

quantreg::用包 制作的情节:在此处输入图像描述

nlrq()正如您可能看到的,我测试过的函数的回归线与发布的图表完全不同。

你有什么建议吗?

0 投票
1 回答
169 浏览

r - RMS::orm 中的分位数图

我不确定这个问题是否与 CrossValidated 或 StackOverflow 更相关。如有必要,我很乐意迁移它。

我试图从Roger Koenker分位数回归小插图中重现这些图,但使用Frank Harrell序数回归模型

这是 Koenker 在第 8 页上的图(我最感兴趣的是分位数与预测器的图): Koenker 的阴谋 第 8 页

代码

我可以制作类似的图,但不完全相同。例如:

我相信 SAS 在proc quantreg.

0 投票
1 回答
542 浏览

r - 如何使用 stargazer 创建分位数回归表?

quantreg我使用包计算了以下分位数回归

我正在尝试使用stargazer. 我试图运行的代码如下:

但是,我收到以下错误消息

我假设此错误消息与标准错误函数rq.se有关stargazer,例如summary(qr_10, se = "iid")可以正常工作。

有人有解决这个问题的方法吗?谢谢你。