assessStationarity bug

[spedygiorgio]

I have a question, however, about the assessStationarity package. It works well if there are no zeros in the frequency matrix but throws an NaN when zeros exist.

Below I provide two sequences vectors. They both come from Brian O’Conner’s article (https://people.ok.ubc.ca/brioconn/sequential/sequential.html). The first one is the original sequence and the second one includes a few more codes that eliminates the zero cells. Any chance that assessStationarity routine might include a switch that allows estimation of stationarity when zeros exist? Thanks very much!

Best,

Joe

oc <- c(2,5,2,4,2,5,2,4,2,4,2,2,5,2,4,2,5,3,4,3,4,1,4,4,1,2,5,2,1,5,4,2,5, 1,4,4,2,3,6,6,2,2,6,5,5,5,6,3,6,3,3,6,3,2,5,3,5,2,6,2,2,2,5,2,4,4, 5,2,5,2,5,5,4,2,5,5,5,2,5,2,5,2,5,2,2,5,2,5,2,2,5,5,2,2,2,5,2,5,2, 5,2,5,2,2,5,2,5,3,5,2,5,5,2,2,5,2,2,5,2,6,2,5,4,2,5,2,4,4,2,1,2,5, 2,5,2,5,5,5,3,5,2,5,5,2,5,5,2,5,2,5,2,4,2,4,2,5,4,2,2,5,2,5,5,5,2, 5,2,5,5,2,5,2,5,5,2,1,2,5,1,5,1,5,5,4,2,2,2,3,6,3,6,3,6,3,6,3,6,3, 3,6,6,6,5,1,2,5,2,5,5,2,4,5,5,2,5,5,2,5,2,5,2,5,1,5,4,2,2,5,2,5,2, 5,4,1,4,4,2,4,4,2,4,2,3,5,4,2,2,5,2,6,1,4,1,4,5,5,5,4,5,2,4,5,5,2, 5,2,2,2,5,2,2,5,2,2,4,2,2,2,6,2,4,3,3,3,3,2,3,6,3,5,3,3,5,3,2,2,2, 2,6,3,6,3,6,2,6,1,6,3,3,3,3,5,3,5,2,5,3,3,6,3,1,4,1,5,1,5,6,3,3,6, 3,6,3,6,3,6,3,6,3,3,6,3,6,3,2,5,2,5,2,5,2,5,2,2,2,3,5,1,6,1,6,6,6, 1,6,6,1,6,1,6,2,5,1,6,6,6,2,2,6,6,6,2,6,2,6,2,5,5,5,5,5,5,6)

gets rid of structural zeros in freq matrix

oc.alt <- c(2,5,2,4,2,5,2,4,2,4,2,2,5,2,4,2,5,3,4,3,4,1,4,4,1,2,5,2,1,5,4,2,5, 1,4,4,2,3,6,6,2,2,6,5,5,5,6,3,6,3,3,6,3,2,5,3,5,2,6,2,2,2,5,2,4,4, 5,2,5,2,5,5,4,2,5,5,5,2,5,2,5,2,5,2,2,5,2,5,2,2,5,5,2,2,2,5,2,5,2, 5,2,5,2,2,5,2,5,3,5,2,5,5,2,2,5,2,2,5,2,6,2,5,4,2,5,2,4,4,2,1,2,5, 2,5,2,5,5,5,3,5,2,5,5,2,5,5,2,5,2,5,2,4,2,4,2,5,4,2,2,5,2,5,5,5,2, 5,2,5,5,2,5,2,5,5,2,1,2,5,1,5,1,5,5,4,2,2,2,3,6,3,6,3,6,3,6,3,6,3, 3,6,6,6,5,1,2,5,2,5,5,2,4,5,5,2,5,5,2,5,2,5,2,5,1,5,4,2,2,5,2,5,2, 5,4,1,4,4,2,4,4,2,4,2,3,5,4,2,2,5,2,6,1,4,1,4,5,5,5,4,5,2,4,5,5,2, 5,2,2,2,5,2,2,5,2,2,4,2,2,2,6,2,4,3,3,3,3,2,3,6,3,5,3,3,5,3,2,2,2, 2,6,3,6,3,6,2,6,1,6,3,3,3,3,5,3,5,2,5,3,3,6,3,1,4,1,5,1,5,6,3,3,6, 3,6,3,6,3,6,3,6,3,3,6,3,6,3,2,5,2,5,2,5,2,5,2,2,2,3,5,1,6,1,6,6,6, 1,6,6,1,6,1,6,2,5,1,6,6,6,2,2,6,6,6,2,6,2,6,2,5,5,5,5,5,5,6,1,1,3, 4,6,4)