r - 制作一个全排序的矩阵？

Question

如何将常规矩阵转换为 R 中的全排序矩阵？有没有可用的方法？

我有一个可能具有线性相关列的矩阵，我需要将它传递给一个函数，该函数需要它的参数是一个具有满秩的矩阵。由于无论如何对线性相关的列都不感兴趣，因此我正在寻找一个删除此类列的函数，直到矩阵满秩。当然，可能有几种解决方案，但其中任何一种都应该没问题。

现在我只是逐列构建矩阵，并且只在其结果矩阵仍然是满秩的情况下添加一列，但感觉应该有更好的方法来做到这一点。

Answer 1

另一种方法是最小化 |y - Ax| ² + c |x| ²，通过将单位矩阵附加到 A 和零到 y。参数 c (aka λ) 权衡拟合 y - Ax，并保持 |x| 小的。然后用 x 的 r 个最大分量进行第二次拟合，r = rank(A)（或任何你喜欢的数字）。