gam0_allprobs: Compute allprobs
In simonecollier/lizardHMM: Fitting Lizard Movement Data with HMMs

This function computes the probability of each observation arising from each state.

gam0_allprobs(
  num_states,
  num_variables,
  num_subjects,
  num_time,
  x,
  pn,
  state_dep_dist_pooled = FALSE
)

`num_states`	The number of states in the desired HMM.
`num_variables`	The number of variables in the data.
`num_subjects`	The number of subjects/trials that generated the data.
`num_time`	A value indicating the length of the data.
`x`	The data to be fit with an HMM in the form of a 3D array. The first index (row) corresponds to time, the second (column) to the variable number, and the third (matrix number) to the subject number.
`pn`	A list containing the natural parameters for the gamma HMM.
`state_dep_dist_pooled`	A logical variable indicating whether the state dependent distribution parameters `alpha` and `theta` should be treated as equal for all subjects.