Error: solve.default(CtC[vars, vars, drop = FALSE]) : 'a' is 0-diml

[Nicoziana]

Hello,

when I run nmf on a document-term matrix of dim 550 x 15000: nmf(tdm, method = "snmf/r", seed = "nndsvd" , rank = 60 , nrun = 1)

I get the following error:

Error in solve.default(CtC[vars, vars, drop = FALSE]) : 'a' is 0-diml

The error only appears for rank >= 60. Any idea of what the problem could be?