Exit é inicialmente definida como NAs, mas depois uma série de valores são colocados nelas. O código todo desde a definição de Exit até a de trace é o seguinte. Exit<-array(rep(NA,7*14),dim=c(7,14)) #alpha1 Exit[1,1]<-alphaE1 Exit[2,1]<-alphaE1 Exit[3,1]<-alphaE1 Exit[4,1]<-alphaI1 Exit[5,1]<-alphaD1 Exit[6,1]<-alphaS1 Exit[7,1]<-alphaA1 #alpha2 Exit[1,2]<-alphaE2 Exit[2,2]<-alphaE2 Exit[3,2]<-alphaE2 Exit[4,2]<-alphaI2 Exit[5,2]<-alphaD2 Exit[6,2]<-alphaS2 Exit[7,2]<-alphaA2 #RMSEA1 Exit[1,3]<-CFA4[1,4,1] Exit[2,3]<-CFA4[1,4,1] Exit[3,3]<-CFA4[1,4,1] Exit[4,3]<-CFA4[2,4,1] Exit[5,3]<-CFA4[3,4,1] Exit[6,3]<-CFA4[4,4,1] Exit[7,3]<-CFA10[1,4,1] #RMSEA2 Exit[1,4]<-CFA4[1,4,2] Exit[2,4]<-CFA4[1,4,2] Exit[3,4]<-CFA4[1,4,2] Exit[4,4]<-CFA4[2,4,2] Exit[5,4]<-CFA4[3,4,2] Exit[6,4]<-CFA4[4,4,2] Exit[7,4]<-CFA10[1,4,2] #2nd greatest eigenvalue Exit[7,5]<-eigen2[1,1] Exit[7,5]<-eigen2[1,2] #R2 regression sumated scale Exit[1,6]<-summary(rEE)$r.squared Exit[2,6]<-summary(rED)$r.squared Exit[3,6]<-summary(rEI)$r.squared Exit[4,6]<-summary(rII)$r.squared Exit[5,6]<-summary(rDD)$r.squared Exit[6,6]<-summary(rSS)$r.squared #p regression sumated scale Exit[1,7]<-summary(rEE)$coef[7] Exit[2,7]<-summary(rED)$coef[7] Exit[3,7]<-summary(rEI)$coef[7] Exit[4,7]<-summary(rII)$coef[7] Exit[5,7]<-summary(rDD)$coef[7] Exit[6,7]<-summary(rSS)$coef[7] Exit[7,7]<-summary(rAA)$coef[7] #R2 regression IRT Exit[1,8]<-summary(rScEE)$r.squared Exit[2,8]<-summary(rScED)$r.squared Exit[3,8]<-summary(rScEI)$r.squared Exit[4,8]<-summary(rScII)$r.squared Exit[5,8]<-summary(rScDD)$r.squared Exit[6,8]<-summary(rScSS)$r.squared Exit[6,8]<-summary(rScAA)$r.squared #p regression IRT Exit[1,9]<-summary(rScEE)$coef[7] Exit[2,9]<-summary(rScED)$coef[7] Exit[3,9]<-summary(rScEI)$coef[7] Exit[4,9]<-summary(rScII)$coef[7] Exit[5,9]<-summary(rScDD)$coef[7] Exit[6,9]<-summary(rScSS)$coef[7] Exit[7,9]<-summary(rScAA)$coef[7] #--------- #SEM part #--------- #Easy-Easy Exit[1,10]<-SEM[1,1] #Chi2 Exit[1,11]<-SEM[1,2] #Chi2/df Exit[1,12]<-SEM[1,3] #CFI Exit[1,13]<-SEM[1,4] #RMSEA Exit[1,14]<-slink[1] #p #E-D Exit[2,10]<-SEM[5,1] #Chi2 Exit[2,11]<-SEM[5,2] #Chi2/df Exit[2,12]<-SEM[5,3] #CFI Exit[2,13]<-SEM[5,4] #RMSEA Exit[2,14]<-slink[5] #p #I-I Exit[4,10]<-SEM[2,1] #Chi2 Exit[4,11]<-SEM[2,2] #Chi2/df Exit[4,12]<-SEM[2,3] #CFI Exit[4,13]<-SEM[2,4] #RMSEA Exit[4,14]<-slink[2] #p #D-D Exit[5,10]<-SEM[3,1] #Chi2 Exit[5,11]<-SEM[3,2] #Chi2/df Exit[5,12]<-SEM[3,3] #CFI Exit[5,13]<-SEM[3,4] #RMSEA Exit[5,14]<-slink[3] #p #S-S Exit[6,10]<-SEM[4,1] #Chi2 Exit[6,11]<-SEM[4,2] #Chi2/df Exit[6,12]<-SEM[4,3] #CFI Exit[6,13]<-SEM[4,4] #RMSEA Exit[6,14]<-slink[4] #p tr[,,count]<-Exit mais uma vez, obrigado, Felipe 2011/12/8 Benilton Carvalho <beniltoncarvalho@gmail.com>

mas sua sua matriz 'tr' eh de NA e vc esta' substituindo por outra de NA, o q vc esperaria obter? (q pedaco de informacao eu perdi?) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

Desculpe minha pouca clareza na minha postagem inicial. 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. Mas, quando eu peço a matriz tr, recebo só NA's Arrumei o código. Espero que esta postagem fique mais clara. Mais uma vez obrigado, Felipe 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. 2011/12/8 Benilton Carvalho <beniltoncarvalho@gmail.com>

Felipe,

leia o guia de postagem (

http://www.leg.ufpr.br/doku.php/software:rbr/#guia_de_postagem_da_r-br_em_co... ) e tambem https://gist.github.com/1088208 .

Feito isso e produzido o CMR, volte a postar.

b _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.