estimerrors              package:rHVDM              R Documentation

_c_o_m_p_u_t_e_s _t_h_e _s_t_a_n_d_a_r_d _d_e_v_i_a_t_i_o_n _o_f _t_h_e _m_e_a_s_u_r_e_m_e_n_t _e_r_r_o_r _u_s_i_n_g _p_r_e-_c_a_l_c_u_l_a_t_e_d _t_a_b_l_e_s _s_p_e_c_i_f_i_c _t_o _t_h_e _p_l_a_t_t_f_o_r_m _o_r _u_s_e_r-_d_e_f_i_n_e_d _t_a_b_l_e

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

     This method takes as main input and ExpressionSet object, and
     using pre-calculated tables contained in rHVDM, it returns the
     same object with an extra slot that contains the standard
     deviation of the measurement error which is indispensable to then
     run rHVDM proper. The table can be user-defined. To check which
     tables are stored in rHVDM, just run the command without arguments
     ('estimerrors()'), it will return a list of the supported rHVDM
     plattforms. All parameters of this function are optional, see
     below for details on the effect of omitting one or more of them.
     The method assumes the the amount of error is dependent both on
     the particular microarray and the signal intensity. For more
     details on how it is computed, please refer to the  textual
     description of rHVDM.

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

     estimerrors(eset,plattid,refchips,errtable)

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

    eset: an object of class ExpressionSet (Biobase), every time this
          parameter is omitted the function returns a list of supported
          plattforms (identifier+plattform description). 

 plattid: an optional argument (character format) allowing to specify
          the plattform identifiers. These identifiers are rHVDM
          specific. Any identifier will be overriden by a table given
          as input of the 'errtable' argument. In future versions of
          rHVDM, if this argument and the 'errtable' are missing, a
          search will be done among the supported plattforms to find a
          matching plattform, based on the individual genes
          identifiers. For now, omission of both arguments will be
          conductive to outputting a list of supported plattforms. 

refchips: a vector of names or column indexes compatible with the
          ExpressionSet useds (note that this compatibility is not
          verified by the function, it is up to the user to supply
          compatible array names). Some microarray might be known to be
          of not very good quality and although the measurement error
          will be computed for these chips it is better to leave them
          out of early stages of the computation. The list to be input
          should typically contain the least noisy microarrays. This
          parameter is optional and if omitted, all microarrays are
          used. 

errtable: a nx2 table in matrix format. The first column of the matrix
          should contain reference log of signals in ascending order
          and the second the variance corresponding to the signals. 

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

     Local interpolation and individual array normalisation are used to
     estimate the standard deviation of the measurement error for each
     individual transcript on each microarray. The precise method is
     laid out in the paper cited below, with some minor modifications.

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

     Normally an updated eset is returned. In case some crucial element
     is missing, the original eset is returned. If the latter is
     missing too a list of supported plattform is returned (thus, when
     using this function, it is better to be careful of what lies at
     the LHS of the assignement arrow).

_N_o_t_e:

     The HTML report is generated in the working directory.

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

     Martino Barenco

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

     M. Barenco, D. Tomescu, D. Brewer, R. Callard, J. Stark, M. Hubank
     (2006) Ranked predictions of p53 targets using Hidden Variable
     Dynamic Modelling. _Genome Biology_, *V7(3)*, R25.

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

     'HVDMcheck','training','screening','fitgene'

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

     data(HVDMexample)
     data(plattforms)
     fiveGybis<-estimerrors(eset=fiveGyMAS5,plattid="affy_HGU133A",refchips=leastnoisychips)
     #leastnoisy chips is a list of three chips identifiers stored in the HVDMexample data bundle

