qoiqpwqr 发表于 2012-12-27 08:52 
可以看下面这个例子
> set.seed(1)
而且算 rank 也是没有用的。
比如:> set.seed(1)
> a <- matrix(rnorm(100),10,
10)
> a[lower.tri(a)] <-0
> diag(a) <-rep(0.00001,
10)
> qr(a)$rank
[1]
10
> det(a)
[1] 1e-50
> solve(a)
错误于solve.default(a) : 系统计算上是奇异的: 倒条件数=9.93383e-50