calc_markov_suff_stats: Calculate Markov chain sufficient statistics.
In a-parton/CTStepTurn: Fitting the continuous-time step-and-turn movement model

Description Usage Arguments Value See Also Examples

calc_markov_suff_stats calculates the sufficient statistics (time in each state and the number of transitions between state pairs) for a Markov chain.

1	calc_markov_suff_stats(num_states, switch_times, switch_states)

`num_states`	Numeric, number of behavioural states.
`switch_times`	Vector, times of behavioural switches.
`switch_states`	Vector, states at switching times.

List with components:

time_in_state: Vector (length num_states) with total time spent in each behavioural state.
num_trans: Matrix (square, size num_states with NA diagonal elements) with number of transitions from each state (row) to every other (column).

Other Behaviour parameters: calc_markov_post_param, check_valid_markov, update_behav_param

times <- c(0, 12, 15, 20, 26, 43, 50)
states <- c(1, 2, 1, 2, 1, 2, 2)

calc_markov_suff_stats(2, times, states)
# $time_in_state
# [1] 34 16
#
# $num_trans
# [,1] [,2]
# [1,]   NA    3
# [2,]    2   NA

times <- c(100, 140, 167, 189, 241, 264, 279, 304, 365, 397, 400)
states <- c(1, 2, 1, 3, 2, 3, 1, 2, 3, 2, 2)
calc_markov_suff_stats(3, times, states)
# time_in_state
# [1]  87 114  99
#
# $num_trans
# [,1] [,2] [,3]
# [1,]   NA    2    1
# [2,]    1   NA    2
# [3,]    1    2   NA