我正在尝试在R中使用grplasso来分析一些遗传数据。然而，这是我第一次使用它，我真的很挣扎。我的X设计矩阵有问题。
该数据由5677个预测因子（SNP）组成，它们形成 136 组（基因）和 747 观察结果。
我在grplasso(x,y,index,lambda)函数中用作X的矩阵是一个矩阵，其中 747 行和 5677 +1 列。另外一列是矩阵的第一列，其中包含截距。
当我尝试运行grplasso时，我得到一个ER的“子矩阵”的错误， 对应于第一组（第2列至第26列）的，不是完整的 秩。
我试图对整个X（没有拦截列）进行正交化，以及对应于不同组的不同“子矩阵”;使用“远程包中的正交归一化函数。但是在整个X矩阵的情况下，我得到一个警告，列的数量太大。另一方面，在第二种情况下，我得到一个列之间存在共线性，并且无法完成正交归一化的错误
您对如何处理该问题有任何想法吗？我使用错误的设计矩阵吗？