nemModelSelection            package:nem            R Documentation

_M_o_d_e_l _s_e_l_e_c_t_i_o_n _f_o_r _n_e_s_t_e_d _e_f_f_e_c_t _m_o_d_e_l_s

_D_e_s_c_r_i_p_t_i_o_n:

     Infers models with different regularization constants, compares
     them via the BIC or AIC criterion and returns the highest scoring
     one

_U_s_a_g_e:

     nemModelSelection(lambdas,D,inference="nem.greedy",models=NULL,control=set.default.parameters(unique(colnames(D))),verbose=TRUE,...)

_A_r_g_u_m_e_n_t_s:

 lambdas: vector of regularization constants

       D: data matrix with experiments in the columns (binary or
          continious)

inference: 'search' to use exhaustive enumeration, 'triples' for
          triple-based inference, 'pairwise' for the pairwise
          heuristic, 'ModuleNetwork' for the module based inference,
          'nem.greedy' for greedy hillclimbing, 'nem.greedyMAP' for
          alternating MAP optimization using log odds or log p-value
          densities

  models: a list of adjacency matrices for model search. If NULL, an
          exhaustive enumeration of all possible models is performed.

 control: list of parameters: see 'set.default.parameters'

 verbose: do you want to see progression statements? Default: TRUE

     ...: other arguments to pass to function 'nem' or 'network.AIC'

_D_e_t_a_i_l_s:

     'nemModelSelection' internally calls 'nem' to infer a model with a
     given regularization constant. The comparison between models is
     based on the BIC or AIC criterion, depending on the parameters
     passed to 'network.AIC'.

_V_a_l_u_e:

     nem object

_A_u_t_h_o_r(_s):

     Holger Froehlich

_S_e_e _A_l_s_o:

     'set.default.parameters', 'nem', 'network.AIC'

_E_x_a_m_p_l_e_s:

        data("BoutrosRNAi2002")
        D <- BoutrosRNAiDiscrete[,9:16]   
        hyper = set.default.parameters(unique(colnames(D)), para=c(0.13, 0.05), Pm=diag(4))
        res <- nemModelSelection(c(0.1,1,10), D, control=hyper)      
        
        plot.nem(res,main="highest scoring model")      

