Protein profile by NMF

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view

Protein profile by NMF

Marcal Plans
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!
[hidden email] mailing list -- To UNSUBSCRIBE and more, see
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.