-
Notifications
You must be signed in to change notification settings - Fork 38
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
assessStationarity bug #63
Labels
Comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
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)
The text was updated successfully, but these errors were encountered: