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,0.05,0.05),
             method=c("unadj","unadj","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 four double values to indicate the thresholds
          for four F tests. The values should be between 0 and 1. Note
          that you need to put four values here even if you don't have
          all four F tests in matestobj.

  method: A flag indicating to use which P values to generate the plot
          and select genes. This is a vector with four elements, which
          corresponds to four 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", "unadj", "unadj") which
          means to use tabulated P values for all four tests.

          Note that you need to put four values here even if you don't
          have all four 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 crresponds to one T-test. You must have F1 test result
     in the input object in order to do volcano plot. 

     On the plot, the y-axis value is -log10(P-value) for the F1 test.
     The x-axis value is propotional to the fold changes.  A horizontal
     line represents the significance threshold of the F1 test. The red
     dots are the genes selected by the F2 test (if there's F2 test
     result).  The green dots are the genes selected by the F3 test (if
     there's F3 test result). The orange 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.F2: The significant genes selected by F2 test.

  idx.F3: The significant genes selected by F3 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:

     data(paigen)
     paigen <- createData(paigen.raw, n.rep=2)

     # make model without interaction
     model.noint.fix <- makeModel(data=paigen, formula=~Array+Dye+Spot+Strain+Diet)

     # F-test strain effect
     ## Not run: 
     test.strain.fix <- matest(paigen, model.noint.fix, term="Strain", n.perm=500,
             shuffle.method="resid", test.method=rep(1,4))
     # volcano plot
     idx.strain.fix <- volcano(test.strain.fix, title="Strain test - fixed model")

     # T-test all pairwise comparison on strain
     C <- matrix(c(1,-1,0,1,0,-1, 0,1,-1), nrow=3, byrow=TRUE)
     ttest.strain.fix <- matest(paigen, model.noint.fix, term="Strain",
              Contrast=C, n.perm=500, test.method=rep(1,4))
     # volcano plot
     volcano(ttest.strain.fix)
     ## End(Not run)

