Skip to main content
Mathematica

You are not logged in. Your edit will be placed in a queue until it is peer reviewed.

We welcome edits that make the post easier to understand and more valuable for readers. Because community members review edits, please try to make the post substantially better than how you found it, for example, by fixing grammar or adding additional resources and hyperlinks.

Required fields*

4
  • 4
    $\begingroup$ Maybe x simply cannot be factored this way? Moreover, it is more realistic to condsider a relative error measure. E.g., Norm[w.h - x, "Frobenius"]/Norm[x, "Frobenius"] returns 0.00326206 which is not that bad... With MaxSteps -> 10000, one can get down to 0.00075928 or so. $\endgroup$ Commented Mar 6, 2019 at 20:41
  • $\begingroup$ If you create x = xL.xR then it for sure can be expressed as w.h, and there is still significant error in the Norm. But maybe you are right, the error is small compared to the size of x. $\endgroup$ Commented Mar 6, 2019 at 20:48
  • 1
    $\begingroup$ @HenrikSchumacher beat me to it! (BTW, the automatic precision goal is 4.) $\endgroup$ Commented Mar 6, 2019 at 20:58
  • $\begingroup$ "This is not a built-in function in Mathematica, but there is a package that implements it [...]" -- see the implementation and documentation "NonNegativeMatrixFactorization" published 12 days ago at Wolfram Function Repository. $\endgroup$ Commented Jan 1, 2020 at 16:20