
Hi all,

Has any of you implemented code for non-negative matrix factorization to solve

Y=T P' +E; dim(Y)=n,p ; dim(T)=n,nc; dim (P)=(p,nc); dim(E)=n,p

where T and P must be non-negative and E either Gaussian or Poisson noise. 

I'm looking for two variants:

1. Easy (I think), T is known (that is we just want to solve the general inverse problem)

2. Harder (?), T is unknown (under some restrictions) [as an intermediate, we may want to fix nc]

Solutions/suggestions welcome

Have a nice day,


Christian Ritter
Senior Statistical Consultant
Analytical and Statistical Services
Shell Coordination Centre S.A.
Monnet Centre International Laboratory, Avenue Jean Monnet 1, B-1348 Louvain-La-Neuve, Belgium

Tel: +32 10 477  349 Fax: +32 10 477 219
Email: christian.ritter@shell.com
Internet: http://www.analyticalsolutions4u.eu


	[[alternative HTML version deleted]]

