The documentation doesn't specify what values of k
are admissible, and k = 1
seems like a reasonable user choice, but results in an error, when using the example data used by the package.
nsY <- naiveRandRUV(Y, cIdx, nuCoeff = 0, k = 1) Error in svdYc$u[, 1:k] %*% diag(svdYc$d[1:k]) : non-conformable arguments
Surprisingly, a value such as k = 2.5 does not produce any error.
nsY <- naiveRandRUV(Y, cIdx, nuCoeff = 0, k = 2.5)
There's also no description in the documentation how naiveRandRUV
handles NA
values.
> Y[2, 2] <- NA > nsY <- naiveRandRUV(Y, cIdx, nuCoeff=0, k=2) >
It does not produce any error, but are the results correct ? Are there plans to improve the user documentation and data input checking in a future version ?