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:

     A statistical approach for CGH microarray data analysis. Franck
     Picard, Stephane Robin, Marc Lavielle, Christian Vaisse, Gilles
     Celeux, Jean-Jacques Daudin, Rapport de recherche No. 5139, Mars
     2004, Institut National de Recherche en Informatique et en
     Automatique (INRIA), ISSN 0249-6399. The code of this function is
     based on the Matlab implementation presented at   <URL:
     http://www.inapg.fr/ens_rech/mathinfo/recherche/mathematique/outil
     .html>, but it has evolved.

_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)

