Я генерирую Q-матрицу для CDM. Каждый раз, когда я запускаю код, у меня появляется эта ошибка: «Ошибка в Q [1: J,]: индекс за пределами границ»
N = 1000
K = 4
J = 100
Q = matrix(rep(diag(K),2),2*K,K,byrow=TRUE)
for(mm in 2:K){
temp = combn(1:K,m=mm)
tempmat = matrix(0,ncol(temp),K)
for(j in 1:ncol(temp)) tempmat[j,temp[,j]] = 1
Q = rbind(Q,tempmat)
}
Q = Q[1:J,]
Also the output should be 2^k, in this example it should be 16 output, but I got 19!! I dont know what is the error?
Спасибо в авансовый