segment             package:tilingArray             R Documentation

_F_i_t _a _p_i_e_c_e_w_i_s_e _c_o_n_s_t_a_n_t _c_u_r_v_e: _s_e_g_m_e_n_t_a_t_i_o_n _b_y _d_y_n_a_m_i_c
_p_r_o_g_r_a_m_m_i_n_g

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

     The function fits a piecewise constant curve to one or multiple
     sequences of measurements, using a least squares cost function and
     an O(n) dynamic programming algorithm (see references).

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

     segment(y, maxseg, maxk)

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

       y: Numeric matrix. Rows correspond to the 'x'-variable, columns
          to replicate measurements at the same value of 'x'.
          Breakpoints are fitted along the 'x'-axis. For example, the
          'x'-variable can be genomic coordinates or time. The
          segmentation will be along the rows of 'y'.

  maxseg: integer of length 1, maximum number of segments (= 1 +
          maximum number of change points).

    maxk: integer of length 1, maximum length of a single segment.

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

     The complexity of the algorithm is  'length(x)*maxk' in memory and
     'length(x)*maxk*maxcp' in time.

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

     An object of class 'segmentation'.

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

     W. Huber huber@ebi.ac.uk

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

     [1] Transcript mapping with high-density oligonucleotide tiling
     arrays. Huber W, Toedling J, Steinmetz, L. Bioinformatics 22,
     1963-1970 (2006).

     [2] A statistical approach for CGH microarray data analysis.
     Franck Picard, Stephane Robin, Marc Lavielle, Christian Vaisse,
     Gilles Celeux, Jean-Jacques Daudin. BMC Bioinformatics. 2005 Feb
     11; 6:27.

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

       x = rep( sin((0:4)/2*pi), each=3) + rnorm(3*5, sd=0.1)
       res = segment(x, maxseg=6, maxk=15)

