- seq_len(m)[ apply( abs(params$phi[j,,]) > thresh, 1, any ) ]
+ #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 ) ] ) ))