Given a set of mutational signatures, inferring the activity of each mutational signature within individual tumors was performed by modifying the NMF multiplicative update process15,58. We randomly initialize a starting activity matrix H0 and update Hi to Hi+1 via the multiplicative update rule given for H15. W remains fixed and so is not updated. We keep track of the Frobenius norm of the error given by V-WHi at each iteration. We terminated the update process when error vs. iterations demonstrated horizontal asymptomatic behavior, which we defined as when 1/20 of the difference of the error between V-WHi-20 and V-WHi is below a given threshold (0.0001). The resulting matrix Hi is a good representation of the activity of the signatures in W.
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this article to respond.
Tips for asking effective questions
+ Description
Write a detailed description. Include all information that will help others answer your question including experimental processes, conditions, and relevant images.