tol参数在QR分解中的作用
我是R中的新用户。您能告诉我或介绍一些关于计算R中QR分解的tol参数在QR分解中的作用,r,R,我是R中的新用户。您能告诉我或介绍一些关于计算R中QR分解的tol参数的参考资料吗 例如,这两行的区别是什么: qr(A, tol=1e-07) #Doesn't work qr(A, tol=1e-20) #Works 为什么我用这么小的tol值得到我想要的结果,而不是用更大的值?参数tol控制qr是否会根据列是否被判断为线性相关来返回列的值。我认为将tol值降低到1e-16以下会破坏该检查的目的。(这几乎是双精度数学中对零的实用定义。) 首先查看qr.default,然后找到FORTRA
tol
参数的参考资料吗
例如,这两行的区别是什么:
qr(A, tol=1e-07) #Doesn't work
qr(A, tol=1e-20) #Works
为什么我用这么小的
tol
值得到我想要的结果,而不是用更大的值?参数tol
控制qr
是否会根据列是否被判断为线性相关来返回列的值。我认为将tol
值降低到1e-16以下会破坏该检查的目的。(这几乎是双精度数学中对零的实用定义。)
首先查看qr.default
,然后找到FORTRAN代码:
这是FORTRAN例程中描述逻辑的注释:
c cycle the columns from l to p left-to-right until one
c with non-negligible norm is located. a column is considered
c to have become negligible if its norm has fallen below
c tol times its original norm. the check for l .le. k
c avoids infinite cycling.
您看过
qr()
上的文档了吗?是的,但我没有读过。这就是我问这个问题的原因。