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.