hmm                package:VanillaICE                R Documentation

_F_i_t_s _t_h_e _h_i_d_d_e_n _M_a_r_k_o_v _m_o_d_e_l

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

     Fits the hidden Markov model

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

     hmm(object, params)

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

  object: An object of class 'HmmOptions'

  params: An object of class 'HmmParameter'

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

     None yet.

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

     An object of class 'HmmPredict'

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

     R. Scharpf

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

     RB Scharpf et al. (2008) Hidden Markov Models for the assessment
     of chromosomal alterations using high-throughput SNP arrays,
     Annals of Applied Statistics

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

     'HmmParameter-class' 'HmmPredict-class' 'SnpLevelSet-class'

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

     data(chromosome1)
     chromosome1 <- chromosome1[1:500, ]
     options <- new("HmmOptions",
                    states=c("D", "N", "L", "A"),
                    snpset=chromosome1,
                    copyNumber.location=c(1, 2, 2, 3),
                    probHomCall=c(0.99, 0.7, 0.99, 0.7))
     validObject(options)
     params <- new("HmmParameter", 
                   states=states(options),
                   initialStateProbability=0.99)
     cn.emission <- copyNumber.emission(options)
     gt.emission <- calls.emission(options)
     emission(params) <- cn.emission + gt.emission ##log scale
     genomicDistance(params) <- exp(-2 *calculateDistance(options)/(100*1e6))
     ##no scaling
     transitionScale(params) <- matrix(1, length(states(options)), length(states(options)))
     if(validObject(params)) fit <- hmm(options, params)

