fom                 package:maanova                 R Documentation

_F_i_g_u_r_e _o_f _M_e_r_i_t

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

     K-means clustering needs a given number of groups, which is
     difficult to guess in most of the cases. This function calculates
     the Figure of Merit values for different number of groups and
     generates the FOM plot (FOM value versus number of groups). Lower
     FOM value means better grouping. User can decide the number of
     groups in kmeans cluster based on that result.

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

     fom(anovaobj, idx.gene, term, ngroups)

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

anovaobj: An object of class 'maanova'.

idx.gene: The index of genes to be clustered.

    term: The factor (in formula) used in clustering. The expression
          level for this term will be used in clustering. This term has
          to  correspond to the gene list, e.g, idx.gene in this
          function. The gene list should be the significant hits in
          testing this term.

 ngroups: The number of groups for K-means cluster. This could be a
          vector or an integer.

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

     A vector of FOM values for the given number of groups

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

     Hao Wu

_R_e_f_e_r_e_n_c_e_s:

     Yeung, K.Y., D.R. Haynor, and W.L.Ruzzo (2001). Validating
     clustering for gene expression data. _Bioinformatics_,
     *17:309-318*.

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

     'macluster', 'consensus', 'kmeans'

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

     # load in data
     data(abf1)
     # fit the anova model
     ## Not run: 
     fit.fix = fitmaanova(abf1,formula = ~Strain)
     # test Strain effect 
     test.fix = matest(abf1, fit.fix, term="Strain",n.perm= 1000)
     # pick significant genes - pick the genes selected by Fs test
     idx <- volcano(test.fix)$idx.Fs
     # generate FOM
     m <- fom(fit.fix, idx, "Strain", 10)
     ## End(Not run)

