training                package:rHVDM                R Documentation

_P_e_r_f_o_r_m_s _t_h_e _H_V_D_M _t_r_a_i_n_i_n_g _s_t_e_p _a_n_d _r_e_t_u_r_n_s _a _l_i_s_t _c_o_n_t_a_i_n_i_n_g _t_h_e _r_e_s_u_l_t_s

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

     This method performs the training step of the HVDM algorithm. It
     returns a list that will then be used in the subsequent screening
     step.

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

     training(eset,genes,transforms,degrate,actname,pdata,forcetransforms)

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

    eset: an ExpressionSet object (Biobase) 

   genes: a vector containg the gene identifiers of the training genes 

transforms: a vector containing the kintetic parameter identifiers that
          have to be transformed during optimisation (optional) 

 degrate: value of the anchoring gene degradation rate, expressed in
          inverse unit time (optional) 

   pdata: a data frame, phenoData to be used for the training
          (optional)

 actname: name of the transcription factor of interest (optional)

forcetransforms: Boolean, whether the transformation in argument
          transforms have to be applied 

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

     The first entry in the genes vector is the anchoring gene. This
     means that the sensitivity (Sj) for this genes is set at 1.0 by
     default and that if a degradation rate is supplied it applies to
     that gene.

     An exponential transform is set by default for both the basal (Bj)
     and degradation (Dj) rates (through the transforms argument). This
     forces the values for both these parameters to be positive. It
     also helps to reach a better fit. To turn this off, set the
     'forcetransforms' switch to FALSE. Even in this case the
     degradation rate will not be allowed to take non-positive values
     as it causes problems with the differential operator used
     internally. The value in the vector indicates the parameter to be
     transformed: "Bj": basal rate of transcription, "Sj": sensitivity,
     "Dj": degrdation rate. The entry label indicates the transform to
     be applied; presently, only log-tranforms are implemented (ie
     "exp").

     The 'degrate' argument is optional, but it is recommended to
     provide the algorithm with an externally measured degradation
     rate, as this greatly improves the accuracy and robustness of the
     outcome.

     The 'pdata' argument is also optional. By default the method will
     use the phenoData contained in the expression set. This argument
     can be used for excluding a time point, or an entire replicate. To
     extract the phenoData from the expression set,  use
     'dataframe<-pData(eset)'. The 'dataframe' object obtained can then
     be manipulated as desired.

     The default name of the transcription factor is "trfact1".

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

     a list containing the results (see documentation for more
     details).

_N_o_t_e:

     It is recommended to run the HVDMcheck method before running this
     command.

_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','screening','fitgene','HVDMreport'

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

     data(HVDMexample)
     tHVDMp53<-training(eset=fiveGyMAS5,genes=p53traingenes,degrate=0.8,actname="p53")

