breaks              package:VanillaICE              R Documentation

_I_d_e_n_t_i_f_y _b_r_e_a_k_p_o_i_n_t_s _f_r_o_m _t_h_e _h_i_d_d_e_n _M_a_r_k_o_v _m_o_d_e_l _p_r_e_d_i_c_t_i_o_n_s

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

     Identify breakpoints: physical position of breaks, number of SNPs
     in region, and the called hidden state.

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

     breaks(x, states, position, chromosome, sampleNames, chromosomeAnnotation = NULL, verbose = FALSE)

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

       x: Locus X sample matrix of hidden states where the hidden
          states are represented as integers

  states: Labels for the hidden states

position: Physical position of loci

chromosome: integer indicating chromosome (23=X)

sampleNames: sample labels

chromosomeAnnotation: chromosome annotation. see details

 verbose: verbose output

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

     One may provide their own chromosome annotation with centromere
     start and stop sites.  The format must be the same as the
     chromosomeAnnotation dataset in the R package SNPchip.

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

     data.frame

  sample: sample label

     chr: chromosome (23 = X)

   start: starting physical position of segment

     end: last physical position of segment

  nbases: number of bases in segment

 nprobes: number of probes in segment

   state: label for the state of the segment

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

     R. Scharpf

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

     x <- matrix(rep(c(1, 2, 3, 1, 2), each=50), ncol=1)
     breaks(x, states=c("0", "1", "2"), position=1:nrow(x), chromosome=1,
            sampleNames="A")

