Added murase(). Needs nnls package