inverseVST               package:lumi               R Documentation

_I_n_v_e_r_s_e _V_S_T _t_r_a_n_s_f_o_r_m

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

     Inverse transform of VST (variance stabilizing transform), see
     'vst'.

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

     inverseVST(x, fun = c('asinh', 'log'), parameter)

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

       x: a VST transformed LumiBatch object or a numeric matrix or
          vector

     fun: function used in VST transform 

parameter: parameter of VST function 

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

     Recover the raw data from VST transformed data returned by 'vst'.
     This function can be directly applied to the VST transformed or
     VST + RSN normalized LumiBatch object to reverse transform the
     data to the original scale.

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

     Return the raw data before VST transform

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

     Pan Du

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

     Lin, S.M., Du, P., Kibbe, W.A.,  "Model-based Variance-stabilizing
     Transformation for Illumina Mi-croarray Data", submitted

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

     'vst'

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

     ## load example data
     data(example.lumi)

     ## get the gene expression mean for one chip
     u <- exprs(example.lumi)[,1]
     ## get the gene standard deviation for one chip
     std <- se.exprs(example.lumi)[,1]

     ## do variance stabilizing transform
     transformedU <- vst(u, std)

     ## do inverse transform and recover the raw data
     parameter <- attr(transformedU, 'parameter')
     transformFun <- attr(transformedU, 'transformFun')
     recoveredU <- inverseVST(transformedU, fun=transformFun, parameter=parameter)

     ## compare with the raw data
     print(u[1:5])
     print(recoveredU[1:5])

     ## do inverse transform of the VST + RSN processed data
     lumi.N <- lumiExpresso(example.lumi[,1:2])
     ## Inverse transform. 
     ## Note: as the normalization is involved, the processed data will be different from the raw data.
     lumi.N.raw <- inverseVST(lumi.N)

