View source: R/conditional_probabilities.R
conditional_probabilities | R Documentation |
conditional_probabilities
returns a matrix with the conditional
probabilities of a categorical time series
conditional_probabilities(series, lag = 1)
series |
An object of type |
lag |
The considered lag (default is 1). |
Given a CTS of length T
with range \mathcal{V}=\{1, 2, \ldots, r\}
,
\overline{X}_t=\{\overline{X}_1,\ldots, \overline{X}_T\}
, the function computes the
matrix \widehat{\boldsymbol P}^c(l) = \big(\widehat{p}^c_{ij}(l)\big)_{1 \le i, j \le r}
,
with \widehat{p}^c_{ij}(l)=\frac{TN_{ij}(l)}{(T-l)N_i}
, where
N_i
is the number of elements equal to i
in the realization \overline{X}_t
and N_{ij}(l)
is the number
of pairs (\overline{X}_t, \overline{X}_{t-l})=(i,j)
in the realization \overline{X}_t
.
A matrix with the conditional probabilities.
Ángel López-Oriona, José A. Vilar
weiss2008measuringctsfeatures
sequence_1 <- GeneticSequences[which(GeneticSequences$Series==1),]
matrix_cp <- conditional_probabilities(series = sequence_1) # Computing the matrix of
# joint probabilities for the first series in dataset GeneticSequences
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.