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.