Description Usage Arguments Value References See Also Examples

Perform inference of coupled hidden markov models.

1 2 3 |

`X` |
a matrix of observations. Columns correspond to series (individuals). |

`nb.states` |
a integer specifying the numbers of states. |

`S` |
a matrix of similarity between individuals. |

`omega.list` |
a vector of omega values. |

`var.equal` |
a logical variable indicating whether to treat the variances as being equal (var.equal = TRUE). |

`exact` |
a logical variable indicating whether to use VEM (exact = FALSE) or EM (exact = TRUE) algorithm for the inference of the model. |

`meth.init` |
a string specifying the initialization method ("mclust" or "kmeans") for the (V)-EM algorithm. The default method is "mclust". |

`viterbi` |
a logical variable indicating whether to use Maximum A Posteriori method (FALSE) or Viterbi algorithm (TRUE, by default) for recovering the most likely path. |

`itmax` |
an integer specifying the maximal number of iterations for the CHMM_(V)EM algorithm. |

`threshold` |
a value for the threshold used for the stopping criteria for the CHMM_(V)EM algorithm. |

A list of 4 objets.

`omega`

an integer corresponding to the selected value among the omega.list.

`model`

a list corresponding to the output of the

`CHMM-EM`

or`CHMM-VEM`

function for the selected model.`status`

a matrix with status associated to each series in column and each position in row.

`RSS.omega`

a dataframe with omega values and the associated Residuals Sum of Squares.

Wang, X., Lebarbier, E., Aubert, J. and Robin, S., Variational inference for coupled Hidden Markov Models applied to the joint detection of copy number variations.

1 2 3 4 5 6 7 | ```
data(toyexample)
# Variational inference of a coupled hidden Markov Chains
resCHMM <- coupledHMM(X = toydata, nb.states = 3, S = cor(toystatus),
omega.list = c(0.3, 0.5, 0.7, 0.9))
# Breakpoints positions and status of segments
info <- clusterseg(resCHMM$status)
# head(info)
``` |

```
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.