[R-br] Por que o array tr não recebe a matriz Exit?

[Felipe Buchbinder]

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 em 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 em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111208/43ce1e43/attachment.html>