[R] Protein profile by NMF

Marcal Plans m@rc@lpl@n@ @ending from gm@il@com
Thu Dec 6 07:07:01 CET 2018

Hi community,
I have an interesting challenge that i think that it can be resolved using
I have a total of 52 product with different protein profiles. I would like
to make two protein premixes (two unique profiles) to create all the SKU
using different concentrations of them (mixing). I have tried
the  following code and it partially works. The main problem is that the
levels that i have in the spec are the minimum protein concentrations
[minimum guarantee].
The solution that I got is the average that minimize the error, but with
this profile of proteins and concentrations I have profiles that do not
reach the minimum requirements. We can overdose, but we cannot under dose.


Data=read.table("clipboard", header=T, sep="\t", dec=".")
res=nmf(Data[,-1],2, nrun=30, seed=12345)

#write.csv(S%*%C, "Test1.csv")

Is there any way to add some restrictions to the NMF? for example W%*%H>
Is there a better algorithm to factorize this data set?

Thank you so much for your help!

More information about the R-help mailing list