applying if-else condition on 4 arrays containing 19 matrices of dim 50*50

Question

I have 4 50*50 matrices and I want to use if-else to get 4 sets of binary matrices. Now I if I want to apply a sequence of thresholds then am having problem.

Please see the following code:

X <- mvrnorm(n=50, mu = rep(0, 50), Sigma = matrix(1,50,50)+diag(50))
p1<- nrow(X)
threshold <- seq(from = 5, to = 95, by = 5)

array1<-array(0,c(p1,p1,length(threshold)))
array2<-array(0,c(p1,p1,length(threshold)))
array3<-array(0,c(p1,p1,length(threshold)))
array4<-array(0,c(p1,p1,length(threshold)))

criteria1 <- t(X)%*%X
criteria2 <- t(X)+X
criteria3 <- sqrt(t(X)%*%X)
criteria4 <- (t(X)+X)/2

for (i in 1:length(threshold)) {
  array1[,,i] <- ifelse( criteria1 <threshold[i],0,1)
  array2[,,i] <- ifelse( criteria2 <threshold[i],0,1)
  array3[,,i] <- ifelse( criteria3 <threshold[i],0,1)
  array4[,,i] <- ifelse( criteria4 <threshold[i],0,1)
}

Error in array4[, , i] <- ifelse(criteria2 < threshold[i], 0, 1) : 
  incorrect number of subscripts