volcano               package:maanova               R Documentation

_V_o_l_c_a_n_o _p_l_o_t _f_o_r _F _t_e_s_t _r_e_s_u_l_t_s

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

     This function generates a volcano-like plot given the F test
     results.

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

     volcano(matestobj, threshold=c(0.001,0.05),
             method=c("unadj","unadj"), title="Volcano Plot",
             highlight.flag=TRUE, onScreen=TRUE)

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

matestobj: An object of class 'matest'.

threshold: A vector of three double values to indicate the thresholds
          for three F tests. The values should be between 0 and 1. Note
          that you need to put three values here even if you don't have
          all three F tests in matestobj.

  method: A flag indicates to use which P values to generate the plot
          and select genes. This is a vector with three elements, which
          corresponds to three F tests. Each element should be one of
          the following five selections:

        "_u_n_a_d_j" Unadjusted tabulated P values.

        "_n_o_m_i_n_a_l" Nominal permutation P values.

        "_f_w_e_r" FWER one-step adjusted P values.

        "_f_d_r" FDR adjusted tabulated P values.

        "_f_d_r_p_e_r_m" FDR adjusted nominal permutation P values.

          Default value is c("unadj", "unadj") which means to use
          tabulated P values for all tests.

          Note that you need to put three values here even if you don't
          have all three F tests in matestobj.

   title: Figure title. Default is "Volcano Plot".

highlight.flag: A logical value to indicate whether to highlight the
          genes with bad spots or not.

onScreen: A logical value to represent whether to display the plots on
          screen or not. If TRUE, the figure will be plotted on the
          screen. Otherwise, it will plot the figure on the current
          device. Default is TRUE.

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

     This function allows one to visualize the results from the F or T 
     tests. The figure looks like an erupting volcano. There will be
     one plot For F-test result and multiple plots for T-test result,
     each plot corresponds to one T-test. You must have F1 test result
     in the input object in order to do volcano plot. 

     On the plot, blue dots are the genes selected by the F1 test. The
     y-axis value is -log10(P-value) for the F1 test and x-axis value
     is proportional to the fold changes.  A horizontal line represents
     the significance threshold of the F1 test. The red dots are the
     genes selected by the Fs test (if there's Fs test result).  If
     there is flag information in the data and the user wants to
     highlight the flagged genes, the genes with any bad spots will be
     circled by a black circle.

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

     For F-test volcano plot, it returns an object which is a list of
     the following four fields:  

  idx.F1: The significant genes selected by F1 test.

  idx.Fs: The significant genes selected by Fs test.

 idx.all: The significant genes selected by all four F tests.


     For T-test volcano plot, it returns an array of the above object.
     Each element in the array corresponds to one T-test.

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

     Hao Wu

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

     ## Not run: 
     data(abf1)
     fit.full.mix <- fitmaanova(abf1, formula = ~Strain+Sample, 
         random = ~Sample)
     ftest.all = matest(abf1, fit.full.mix, test.method=c(1,1),
         shuffle.method="sample", term="Strain", n.perm= 100)
     volcano(ftest.all)
     ## End(Not run)

