Using function overall_JMV() I get "Error in Y[, i] : subscript out of bounds". It seems the error comes from the group_data_gen() function inside the overall_JMV() function, but I don't know why. Here is the reproducible example:
- Code: Select all
library(R2ucare)
ch <- matrix(scan(text = "
0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0
0 0 0 0 0 0 0 0 3 0 4 0 0 0 0 0 0 4 0
0 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 0
0 3 0 0 0 0 0 0 0 0 4 0 0 0 4 0 4 0 0
0 0 3 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 0 0 0 0 4 0 0 0 4 0 0 0 0 0 4 0 4 4
1 1 1 0 0 0 0 0 0 3 3 0 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 0 0 0
0 3 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 4 4 0 4 0 4 0 0 0 4 4 4 4 0 4 0 4 0
0 4 4 0 0 0 0 0 0 0 0 0 4 0 4 4 4 0 0
0 3 3 3 0 0 0 0 0 3 3 0 0 4 4 4 0 4 4
0 3 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 4 0
0 0 0 0 0 0 0 0 0 1 0 1 0 0 2 0 0 0 0
0 4 0 0 4 0 0 0 0 0 4 0 0 0 0 0 0 0 0
0 1 0 1 1 0 0 0 0 2 2 0 0 0 0 3 0 3 0
0 4 0 0 0 0 0 0 0 0 0 0 4 0 4 0 0 4 4
0 3 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 3 0 3 0
0 0 3 0 0 0 0 0 4 0 0 4 0 0 0 0 0 4 0
0 0 4 4 0 0 0 0 0 4 0 0 0 0 0 0 0 4 0
0 0 3 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0
0 0 3 0 0 3 0 0 0 0 0 4 0 0 3 0 0 0 0
0 0 0 4 0 0 0 0 0 0 4 0 0 4 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0
0 0 0 0 0 3 0 0 0 0 0 0 4 0 0 0 0 0 0
0 0 0 0 0 0 3 0 0 3 3 3 3 3 3 3 4 0 4
0 0 0 0 0 0 0 0 4 0 0 0 0 0 4 4 0 4 0
0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 1 0 0 0 0 2 0 0 0 3 0
0 0 0 0 0 0 0 0 0 3 0 0 3 3 0 4 4 0 0
0 0 0 0 0 0 0 0 0 3 3 0 0 0 3 0 0 4 4
0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 3 0 0 0 4 0
0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 4 0
0 0 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 3 0 0 0 3 0 0 0 0
0 0 0 0 0 0 0 0 0 0 3 3 3 4 4 4 0 4 0
0 0 0 0 0 0 0 0 0 0 0 3 3 3 3 0 0 4 0
0 0 0 0 0 0 0 0 0 0 0 4 0 0 4 4 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 3 0 3 3 0 4 4
0 0 0 0 0 0 0 0 0 0 0 0 0 3 3 3 0 3 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 4 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 3 0 0 3
0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 2 0 3 0
0 0 0 0 0 0 0 0 0 0 0 3 0 0 3 4 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 4 0 0
4 4 0 0 4 0 0 0 0 0 4 0 0 0 0 0 0 4 4
4 4 4 0 4 0 0 4 0 0 4 0 0 4 0 0 0 0 4"),nrow = 50,byrow = TRUE)
freq=rep(1,nrow(ch))
overall_JMV(ch,freq)
Many thanks in advance for your hints,
Simone