问题标签 [quantreg]

包quantreg安装失败，因为它没有链接到 lapack：

[剪辑]

[剪辑]

奇怪的是，从头开始创建新用户时，一切正常（并且日志中有一个“-llapack”）。

我怀疑某些环境变量对此负责，但无法确定是哪一个：取消设置 LIBRARY_PATH 和 LD_LIBRARY_PATH 时，仍然会发生这种情况。感谢您的任何想法！

我在quantregwin7 64bit 和 linux Mint Maya 64bit 上都使用库，我意识到速度差异很大。

我正在使用用 Fortran 编写的 Frisch-Newton 方法，并且我只测量 .Fortran() 调用的速度。

我正在分析长度为 1600 的向量。

在 Windows 上大约需要 2.5 秒，在 Linux 上大约需要 22 秒。而且我真的不明白为什么..

在此之后，我尝试通过 F2C 将 Fortran 代码翻译成 C 代码，并使用 clapack 函数将其编译成动态库。我通过 .C() 调用了这个函数，我再次测量了时间效率：

在 Windows 上需要 8 秒，在 Linux 上需要 22.5 秒。我真的很困惑。

一切都在 I7-3770、16GB RAM 上完成。

有人可以向我解释为什么它不起作用吗？

这是代码：

Fortran 代码：

C代码：

从 R 调用：

动态库的编译：

我正在为我使用的分位数回归绘制回归摘要quantreg。显然plot.summary.rqs这里使用了这个方法。问题是使用了很多解释变量，每个解释变量都显示在图中。大多数系数的表现与 OLS 没有显着差异，所以我只想挑选并显示其中的一些。

如何选择需要显示的图？我正在使用knitr我的报告，但不想显示几十个变量（你可以使用假人快速到达那里）。有办法摘樱桃吗？

我正在使用该quantreg包在 R 中运行以下分位数回归：

并通过以下方式提取系数和汇总统计量

我得到的结果是

但是，我想知道 R^2/调整后的 R^2 -summary()命令似乎为简单的 OLS 回归提供了它，但在分位数回归的情况下却没有。

有谁知道，如何提取它们？

我正在尝试使用包rq()中的函数运行一些分位数回归quantreg。我已经设置了我的数据并创建了模型等（按照另一个站点上的说明），一切看起来都很好。见下文; 是连续变量，RE定义如下。WTPCArmored

但是，当我调用模型时summary()，rq()我得到了这个

这是很好的信息，但我想要标准误差、p 值等，就像正常的线性回归一样。我检查了文档和其他示例，看起来这是执行此操作的方法。QR25建立summary(QR25)一个模型（

这是我第一次尝试分位数回归，所以我可能错过了前面的步骤或需要指定其他内容。有什么想法吗？我在 RStudio 版本 0.98.501 中运行它。

我想使用quantreg包从分位数回归中提取系数和上限和下限。这是帮助文件中的一个示例。

我知道我可以coefficients()用来获取系数。

但我不知道如何获得出现在summary(). 我看了看str(sf)，但我没有看到如何提取。

最终，我想将 taus、系数和上限/下限放入数据框中以进行进一步处理。

我目前正在使用 Roger Koenker 的quantreg包在 R 中运行分位数回归模型。我有一个有五个级别的有序分类变量和三个无序分类变量（星期几、季节以及是否是夜间），分别有七个、四个和两个类别。除了有序的分类变量（从 1 到 5 有五个级别）之外，所有最初都作为字符串存储在数据框中。这是我将前面提到的变量编码为因子的方式。请注意，之所以调用我的变量，是df_prime因为它是带有一些修剪行的原始数据框：

当我运行回归时，它完全从摘要和我的 acuity_id 的一个级别中删除了 Friday，这是我的有序变量。这可能与 R 如何创建虚拟变量有关吗？我注意到无论出于何种原因，seasonWinter 的 p 值为 1。

我的输出如下，尽管间距不同，因为我的工作场所在 VM 上有 R，我无法复制和粘贴：

我试图在quantreg包中查看函数“crq.fit.pen”的源代码。出于这个原因，我尝试做以下工作。一开始我只是输入了函数的名称，结果如下，

正如您在上面看到的，此函数的主要工作是由另一个函数完成的，该函数由以下函数引用：

我现在的问题是如何查看crqfnbFortran 基本函数？

之后，我做了以下任务并获得了结果，但我看不到 function 的完整代码crqfnb。

请您帮我看看如何查看函数的完整代码crqfnb？

我正在使用该quantreg包计算 R 中的分位数回归。我使用以下命令调用我的 QR。

现在，如果我查看返回值，变量 flightHours 和 flightHoursType 的系数为零/NULL。

添加时se = "nid"，我还会收到一条警告消息：

警告信息：在summary.rq(quantGsReg15, se = "nid") : 7 non-positive fis

我之前计算的 OLS 回归表明情况很可能并非如此。我尝试了不同的 QR 方法（method = "fn"）但没有成功。

对 0.5 分位数调用相同的命令会返回类似的输出。但是，我收到一条不同的警告消息：

警告消息：在 rq.fit.br(x, y, tau = tau, ci = TRUE, ...) 中：解决方案可能不唯一

我究竟做错了什么？我现在用谷歌搜索了 30 分钟，但没有成功。

我注意到，每当我尝试使用正常 OLS 系数及其 CI 绘制系数图及其置信区间 (CI) 时，每当我强制回归原点时都会出现错误。

因此，如果我使用此代码（engel 是 R 中分位数回归示例的数据）：

我没有问题，我的 coefficeint 图被绘制出来了。但如果我使用这个：

我有一个问题，因为拦截经过原点。如何在没有截距的情况下获得分位数回归的第一个代码中的图。