R fixDependence 检测一个矩阵与另一个矩阵的线性相关性

说明

标识矩阵 X2 的列，这些列线性依赖于矩阵 X1 的列。主要用于为嵌套 GAM 设置可识别性约束。

用法

fixDependence(X1,X2,tol=.Machine$double.eps^.5,rank.def=0,strict=FALSE)

参数

X1	一个矩阵。
X2	一个矩阵，其列可能部分线性依赖于 X1 的列。
tol	评估线性相关性时使用的容差。
rank.def	如果在给定 X1 的情况下，X2 中的排名不足程度已知，则可以在此处提供，然后忽略 tol。除非为正且不大于 X2 中的列数，否则不使用。
strict	如果TRUE，则仅检测到单独依赖于X1的列，如果FALSE，则检测到足够的列以使减少的X2满秩并且独立于X1。

细节

该算法使用基于 QR 分解的简单方法：有关详细信息，请参阅 Wood(2017 年，第 5.6.3 节)。

值

X2 的列向量，线性依赖于 X1 的列(或者如果 strict==FALSE 则需要删除以实现独立性和满秩)。 NULL 如果两个矩阵独立。

例子

library(mgcv)
n<-20;c1<-4;c2<-7
X1<-matrix(runif(n*c1),n,c1)
X2<-matrix(runif(n*c2),n,c2)
X2[,3]<-X1[,2]+X2[,4]*.1
X2[,5]<-X1[,1]*.2+X1[,2]*.04
fixDependence(X1,X2)
fixDependence(X1,X2,strict=TRUE)

作者

Simon N. Wood simon.wood@r-project.org

参考

Wood S.N. (2017) Generalized Additive Models: An Introduction with R (2nd edition). Chapman and Hall/CRC Press.