R中的prcomp错误

Question

我正在使用R。我想在一个矩阵上运行prcomp。在一个Linux机器上安装R时，代码运行得很好，但在另一个相同(或者我认为是这样)的Linux机器上安装R时，代码就会中断。代码是

dataf = read.table("~/data/testdata.txt")
pca = prcomp(dataf)

坏实例上的错误消息是

> dataf = read.table("~/data/testdata.txt")
> pca = prcomp(dataf)
Error in La.svd(x, nu, nv) :
  BLAS/LAPACK routine 'DGESDD' gave error code -12

R的两个实例都有R version 2.9.2 (2009-08-24)，而且据我所知，所有的R库和环境变量也是以相同的方式配置的。

那么，有没有人对可能的错误有什么建议呢？这个错误代码是什么意思？(我在网上搜索了一下，没有找到任何有用的东西...)提前谢谢你！

Answer 1

您的第二个盒子可能有损坏的Blas / Lapack安装。不幸的是，我们不能说你的问题是不可重现的。

这是对Blas / Lapack的另一个简单调用--这对您有效吗？

R> crossprod(matrix(1:4, ncol=2))
     [,1] [,2]
[1,]    5   11
[2,]   11   25
R>