findSegments           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 _t_o _a _s_e_q_u_e_n_c_e _o_f _n_u_m_b_e_r_s

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

     The function fits a piecewise constant curve to a sequence of
     numbers using a simple least squares cost function and the dynamic
     programming algorithm of Picard et al. (full citation see below).

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

     findSegments(x, maxcp, maxk, verbose=0)

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

       x: Numeric (real) vector.

   maxcp: Integer (length 1): maximum number of segments (= 1 + maximum
          number of change points).

    maxk: Integer (length 1): maximum length of a segment.

 verbose: Integer (length 1): if this parameter has a positive value,
          various diagnostic output is printed.

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

     A list with two elements 'J' and 'th'. See the vignette, and the
     paper cited below for details.

_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  
     \verb+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 = findSegments(x, maxcp=6, maxk=15)

