系统是计算奇异误差

Kaj*_*aja 12 r

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

 ica<-fastICA(datalist,n.comp=3)
Run Code Online (Sandbox Code Playgroud)

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
Run Code Online (Sandbox Code Playgroud)

我收到此错误:

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

你能告诉我为什么我会收到这个错误,我该如何解决?

Kon*_*nos 10

solve(),使用较小的公差,如solve(..., tol = 1e-17).这应该没问题,因为你得到了reciprocal condition number = 1.16873e-16.帮助文件此相关问题中的更多信息.

  • .@pidosaurus - 如果“倒数条件数”根据大输入数据集不断变化怎么办?我在“R”中使用“mahalanobis()”时遇到这个问题,因为它还调用“solve()”来计算距离。我已在您提到的问题中提出评论。只是想,你对此可能有不同的看法。 (2认同)