Optimize omega in step 1
This function optimizes omega, in fact the convex mixing weights that govern prior over the reference cancer profile.
ISOpureS1.model_optimize.opt_omega( tumordata, model, NUM_ITERATIONS_RMINIMIZE, iter, NUM_GRID_SEARCH_ITERATIONS )
tumordata
: a GxD matrix representing gene expression profiles of tumour samplesmodel
: list containing all the parameters to be optimizedNUM_ITERATIONS_RMINIMIZE
: minimum number of iteration that the minimization algorithm runsiter
: the iteration numberNUM_GRID_SEARCH_ITERATIONS
: number of times to try restarting with different initial valuesThe model with the omega_weights and omega parameters updated
Gerald Quon, Catalina Anghel, Francis Nguyen
Useful links