- #seq_len(m)[ apply( abs(params$phi[j,,]) > thresh, 1, any ) ]
- c(seq_len(m)[ apply( abs(params$phi[j,,]) > thresh, 1, any ) ],
- rep(0, m-length(seq_len(m)[ apply( abs(params$phi[j,,]) > thresh, 1, any ) ] ) ))
+ seq_len(m)[ apply( abs(params$phi[j,,]) > thresh, 1, any ) ]