Matrix of count data; axis 0 is the number of time intervals per day and axis 1 is the number of days in the data.

List of different samples at all time periods

List with parameter values of prior distributions

Number of iterations over which to calcuate likelihood.

Parameter sharing controls <- c(S1, S2): S1 <- force sharing of delta (day effect) among days, S2 <- force sharing of eta (time of day) among days, Values: 1 (all days share), 2 (weekdays/weekends), 3 (none)