Question

我在r中使用fastICA包。在这个包中，我使用fastICA函数，它有一些参数。如果我将n.comp设置为2，那么工作正常，但如果我在此函数中将此参数设置为3或更多：

 ica<-fastICA(datalist,n.comp=3)

datalist在这里是一个包含20行和4列的矩阵：

     [,1]    [,2]     [,3]   [,4]
 [1,] 567.00 324.225 281.0889 538.25
 [2,] 557.75 317.500 269.5556 529.15
 [3,] 543.75 309.900 264.5778 515.95
 [4,] 557.00 316.225 265.0889 528.25
 [5,] 538.25 307.750 266.6667 510.95
 [6,] 531.25 301.025 250.0222 503.70
 [7,] 545.00 311.800 270.9333 517.40
 [8,] 550.00 316.925 284.3778 522.65
 [9,] 514.75 290.300 235.6000 487.75
[10,] 518.00 293.800 245.1556 491.20
[11,] 553.75 318.125 281.6667 526.00
[12,] 563.50 325.925 297.2667 535.75
[13,] 540.00 303.300 241.1556 511.40
[14,] 546.00 310.350 261.6444 517.90
[15,] 567.25 324.425 281.4889 538.50
[16,] 577.75 330.125 285.2222 548.40
[17,] 560.75 317.425 262.3778 531.60
[18,] 570.00 323.925 272.8222 540.65
[19,] 569.00 324.700 278.8444 540.00
[20,] 565.50 324.150 284.1333 537.00

我收到此错误：

Error in solve.default(w %*% t(w)) : 
  system is computationally singular: reciprocal condition number = 1.16873e-16

你可以告诉我为什么我会收到这个错误，我该如何解决？

Answer 1

在solve()中，使用较小的容差，例如solve(..., tol = 1e-17)。这应该没问题，因为你得到reciprocal condition number = 1.16873e-16。 help file和this related question中的更多信息。

系统是计算奇异误差

1 个答案: