RJaCGH (version 1.1.1)

plot.Q.NH: Plot transition probabilities

Description

Plot the probabilities of staying in the same state for a non-homogenous hidden Markov model.

Usage

plot.Q.NH(x, beta, q=-beta, col = NULL, ...)

Arguments

x

Vector of distances between observations

beta

beta parameter of the transition matrix. Must be a square matrix with the same size as the number of hidden states

q

q parameter of the transition matrix. Must be a square matrix with the same size as the number of hidden states

col

vector of colors for each state. Must be of the same size as the number of hidden states

...

aditional arguments passed to plot

Value

• A plot is produced showing the probability of staying in the same hidden state versus distance between adjacent genes, for every state.

Details

Please note that RJaCGH model imposes that q is -beta, and distances are normalized to lay between 0 and 1.

References

Oscar Oscar M. Rueda and Ramon Diaz Uriarte. A flexible, accurate and extensible statistical method for detecting genomic copy-number changes. http://biostats.bepress.com/cobra/ps/art9/. {http://biostats.bepress.com/cobra/ps/art9/}.

Examples

## Model with two hidden states
## Note that RJaCGH normalizes distances to be between 0 and 1
x <- rexp(99)
x <- x/ max(x)
beta <- matrix(c(0, 1, 3, 0), 2, 2)
plot.Q.NH(x=x, beta=beta, q=-beta, col=c(1,2))