O programa não recebe dados, mas faz uma simulação. Ele gera uma população de questões, com diferentes graus de dificuldade (que compõem um teste) e uma população de indivíduos, com uma abilidade unidimensional normalmente distribuída. Esses indivíduos "respondem" às questões e as respostas são analisadas por diferentes técnicas: TRI, análise fatorial confirmatória com correlação de Pearson e com correlação policórica.
O array Exit recebe os resultados dessas várias análises. A cada iteração, a matriz Exit gerada é armazenada no array tr.
Eu consigo ver a matriz Exit.
Arrumei o código. Espero que esta postagem fique mais clara.
NSim<-5 #Number of simulations
#Generates the trace matrix. trace[,,i] contains the i-th simulation of the Exit matrix
tr<-array(rep(NA,7*14*NSim),dim=c(7,14,NSim))
for(count in 1:NSim)
{
#%%%%%%%%%%%%%%%%%%%%%%%%%%
#THE ENTIRE CODE COMES HERE
#%%%%%%%%%%%%%%%%%%%%%%%%%%
#O código faz uma simulação de TRI.
#O código gera questões de um teste
#e respostas dadas por vários indivíduos com diferentes habilidades
#a essas questões. Essas respostas serão analisadas por TRI e análise fatorial confirmatória.
#Os resultados destas análises são colocados nas várias células do array Exit.
#Este código é rodado NSim vezes e os vários arrays Exit são guardados no array tr
Exit<-array(rep(NA,7*14),dim=c(7,14))
#O array Exit recebe os resultados de cálculos anteriores.
#Por exemplo, as linhas abaixo colocam na célula 1,1 do array o alpha de Cronbach de um subteste chamado E1.
Exit[1,1]<-alphaE1
#O array Exit existe. Aí, eu tenho passá-lo para o array tr.
tr[,,count]<-Exit
}
#Porém, ao pedir o array tr, recebo apenas NA's.
#Curiosamente, ao pedir o array Exit,um resultado não-NA aparece.