问题标签 [quantreg]
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.
r - 在分位数估计之后构建累积分布函数
我有这个包含两列(Y 和 X)的数据框。
使用quantreg
包我可以估计以 X 为条件的 Y 的分位数。
我无法以 X 为条件建立 Y 的累积密度函数。有人可以帮助我吗?
估计分位数:
在此之后,如何生成累积分布函数 (Fx)?
这是我的data.frame:
r - 从分位数回归中访问 R 中多个分位数的 p 值
我想使用多个分位数访问分位数回归的 p 值。
给出以下str(fit2_summ)
结果
我无法访问第 50 和第 95 分位数的 p 值,fit2_summ$coefficients[,4]
因为汇总模型存在 2 个分位数而不是单个分位数,在这种情况下,上述命令会起作用。
r - quantreg 包:predict.rq 只接受一个 tau
我正在使用quantreg
包来计算分位数回归。我正在拟合一些 taus 的回归,并想计算预测值和 95% 置信区间。问题是如果添加置信限,该函数predict.qr
似乎不允许超过一个分位数。到目前为止,我使用循环解决了这个问题,但这使我的代码非常长。最小的例子:
1)一个分位数(0.5)工作正常
2)不止一个分位数且只有预测值(无 ci)也可以正常工作
3)当超过一个分位数(0.5,0.75)和ci时,它不再起作用
因此我的问题是:是否可以仅使用 predict.rq 并避免循环来获得多个分位数回归的预测值和置信区间?
r - quantreg 编译导致`libgfortran.a: can not be used when making a shared object`
这是我尝试quantreg
在 R 中安装时出现的问题,当我看到安装错误时,我最初尝试安装 scde quantreg
,所以我尝试先安装quantreg
。
以前我有一个与 lgfortran 类似的问题,我按照上面的帖子说明通过这个链接 https://askubuntu.com/a/680100/586038解决了 我当前使用的 gcc 版本在
现在
给了我道路
然后我做了
目前我有文件 libgfortran.so 在
但不在
我该怎么做才能解决这个问题
r - 使用 rqpd 包在面板上运行分位数回归时,我可以获得上限/下限吗?
使用 rqpd-package 在面板上运行分位数回归时,我可以获得系数的上限/下限吗?
以下是设置(很好,我是 R 的初学者):
代码运行良好,输出为我提供了通常的系数、t 值、p 值等。但我希望找出我的 SE 的上限和下限,以获得 95% 的显着性。
使用 quantreg 包中的 rq 我可以通过捕获以下输出来提取上限/下限:
但是,这不适用于我目前正在查看的面板数据。各位有什么指点吗?
r - 无法安装 quantreg 包和更新 R 版本
我在 RStudio 界面和 Windows 10 操作系统上使用 R 版本 3.3.1。我试图为分位数回归安装“quantreg”包,但是我收到了错误消息
install.packages 中的警告:InternetOpenUrl 失败:“操作超时”
install.packages 中的警告:下载长度 180224 != 报告长度 1521162
install.packages 中的警告:包 'quantreg' 不可用(对于 R 版本 3.3.1)
之后我重复了这个过程几次,但一直得到上述错误的最后一行。认为这是由于使用了较旧的 R 版本,我尝试通过安装“installr”包进行更新,但这也给了我与上述相同的错误。我尝试从工具菜单更新 R,但它一直告诉我当前的 R 版本是最新的。这也是我第一次来这里,如果我搞砸了块引用,请道歉。
r - 使用多项式的 rq 函数计算 R 中分位数回归的 95% 置信区间
这是这个问题的扩展,关于分位数回归的 95% 置信区间,使用rquant
:
在这里,目标是确定多项式拟合的分位数回归的 95% 置信区间。
数据:
尝试使用上述问题中的方法:
但这当然独立地确定了每个系数的 95% CI,并且似乎高估了区间(见下图)。
我对一种方法有另一个想法,即简单地从数据的引导样本(即rq(y~x+I(x^2))
数千个 y 和 x 样本上)确定系数,但想看看是否有东西构建到包中。
r - 如何从分位数回归 rq() 中绘制交互效应
我想在 R 中运行分位数回归,其中包括一个三向交互项,其中位数为三个分类预测变量(tau = 0.5)。为了解释效果,我想创建一个情节。对于线性模型,我使用“效果”包中的函数 effect()。不幸的是,这不适用于“quantreg”包中的 rq()。您对如何显示交互效果有什么建议吗?
r - Quantile regression split plots
I have few years of daily rainfall data for a particular region. To get an insight of extreme rainfall events,I used quantile regression (quantreg package) in R. The plot for entire days is shown below. What I want is to split the regression line in the middle (or some other point) and fit for first and second half of the data separately to see the difference.
Here is how I used quantreg:
r - 在单个页面中按变量绘制分位数回归
我分别对几个自变量(相同的依赖)运行分位数回归。我只想在单个图中绘制每个变量的几个分位数的斜率估计值。
这是一个玩具数据:
我只想绘制分位数的斜率估计。因此,我放弃parm=2
了plot
。
现在,我想将这两个图合并到一个页面中。
到目前为止我已经尝试过什么;
- 我尝试设置
par(mfrow=c(2,2))
和绘制它们。但它正在产生一个空白页。 - 我试过使用
gridExtra
和 gridGraphics 没有成功。尝试将基础图转换为 Grob 对象,如此处所述 - 尝试使用本文档中的函数
layout
函数 - 我正在尝试查看
plot.rqs
. 但我无法理解它是如何绘制置信带的(我只能绘制分位数上的系数)或在mfrow
那里更改参数。
谁能指出我哪里出错了?我应该查看源代码plot.rqs
并更改那里的任何参数吗?