®️ rlinalg

以R语言实现的Python线性代数例程。

🗺️ 概述

The R Project for Statistical Computing provides an environment for using and developing statistical methods. Most of the array manipulation and linear algebra routines are implemented using LAPACK, which can be accessed in Python using SciPy and NumPy.

然而，当尝试将 R 语言中的代码移植到 Python 并重新生成时，人们会注意到几个例程的实现存在差异，尤其是在启用了转置的 QR 分解 方面。

> mat <- t(matrix(seq_len(9), nrow=3))
> qr.Q(mat)
           [,1]       [,2]       [,3]
[1,] -0.1230915  0.9045340  0.4082483
[2,] -0.4923660  0.3015113 -0.8164966
[3,] -0.8616404 -0.3015113  0.4082483

>>> mat = numpy.arange(1, 10).reshape(3, 3)
>>> scipy.linalg.qr(mat, pivoting = True)[0]
array([[-0.2672612  0.8728716  0.4082483]
       [-0.5345225  0.2182179 -0.8164966]
       [-0.8017837 -0.4364358  0.4082483]])

这里的罪魁祸首是 R 语言中的 qr 函数默认不使用 LAPACK 的 dgeqp3，而是使用一个修改过的、针对 R 语言的特有版本 LINPACK 例程 dqrdc（dqrdc2），该例程优化了转置策略。这意味着使用 R 中的 qr 函数的代码将不同于等效的 Python 使用 LAPACK 的代码，并且（直到现在）没有方法来重现 R 的行为。

rlinalg 库通过使用 Fortran 源代码提供 R 语言的线性代数例程，以允许可重复性。它提供了一个类似于 scipy 接口的 API（对于类似函数如 qr、cond、lstsq），可用于获得与 R 相同的结果。

>>> mat = numpy.arange(1, 10).reshape(3, 3)
>>> rlinalg.qr(mat).Q.round(7)
array([[-0.1230915  0.904534   0.4082483]
       [-0.492366   0.3015113 -0.8164966]
       [-0.8616404 -0.3015113  0.4082483]])

此库依赖于 NumPy 和系统上可用的 BLAS 库。它适用于所有现代 Python 版本（3.7+）。使用 Meson 进行构建，并且从源代码编译时需要 Fortran 编译器。

🔧 安装

直接从 PyPi 安装 rlinalg 包，该网站提供通用轮子，可以使用 pip 进行安装。

$ pip install rlinalg

💭 反馈

⚠️ 问题追踪器

发现了一个错误？有增强请求吗？如果您需要报告或询问某些内容，请访问 GitHub 问题追踪器。如果您正在提交错误报告，请尽可能提供关于问题的详细信息，并尝试在简单、易于复现的情况下重现相同的错误。

📋 更新日志

此项目遵循 语义版本控制，并在 保持更新日志 格式的 更新日志 中提供更新。

⚖️ 许可证

此库根据 GNU 通用公共许可证 v3.0 或更高版本提供。它包含一些来自 R 语言的代码，这些代码根据 GNU 通用公共许可证 v2.0 或更高版本进行许可。一些测试是根据 SciPy 修改的，SciPy 在 BSD-3-clause 许可证下开发。

本项目与 R 项目 无任何关联、赞助或任何形式的认可。它是由 Martin Larralde 在 莱顿大学医学中心 的 Zeller 实验室 进行博士项目期间开发的。