es_discover_perturbation_sample

ES_DISCOVER_PERTURBATION_SAMPLE

 [ris, es_result_list] = es_discover_perturbation_sample(N,W,ext_ind,constraints,...
  options,c_ratio,v_ratio,expansion,expansion_order,nrun,u_sigma_prior,s_sigma_prior)

 Compute the most likely perturbation of enzyme expression (u) and external metabolites (s)
 that would give rise to a given differential profile of balanced concentrations c and fluxes j
 do this for a number of model variants
 with sampled elasticities and present a statistics of the results

 N, W, ext_ind: structural information about the network
 x_ratio:       expression ratio column vector (or matrix, for several samples)
 nrun:          number of samples
 constraints, options: options for elasticity sampling
 expansion:     {'logarithmic','non-logarithmic'}
                -> depending on the expansion type, the output refers to logarithmic 
                   or non-logarithmic values; 
 expansion_order:       1 or 2 for first- or second-expansion_order expansion

 Outputs:
 ris:            contains statistics over predicted effects
 es_result_list: contains the results of all sampling runs