Estou fazendo um programa de simulação cujo objetivo é comparar Teoria de Resposta ao Item com outros modelos de análise de variáveis latentes, como análise fatorial confirmatória com correlações de Pearson ou policórica. Em cada simulação, um array Exit é gerada e atribuído a um array tr (traço) que guarda todos os arrays "Exit" para posterior cálculo de média e desvio-padrão dos resultados (que serão, por sua vez, colocados no array Final). O problema é que, embora o array "Exit" seja gerado, os arrays tr e Final não o são. Não consigo entender o porquê....O código segue abaixo (retirei toda a parte de análise por TRI e CFA porque se não ficaria absurdametne grande)

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

}

Desde já, obrigado,

Felipe