Boa tarde listeiros (@s),


Estou aplicando a análise fatorial e queria a ajuda para utilizar o pacote 'sem'.

De acordo com o Quick-r (http://www.statmethods.net/advstats/factor.html).


Segundo o script  apresentado no site os fatores F1 e F2 são descritos em função das variáveis X1, X2,...e as cargas lam1,lam2,...ect.


Bem minha dúvida é: essas posições das cargas (lam1,lam2,...) são fixas e  como posso colocar na equação do modelo as mais significativas, tipo (lam3, lam7).

E como posso retirar a equação usando meus dados. Segue minha rotina.


af=factanal(dad, factors=8, scores=c('Bartlett'), rotation='varimax')
> af

Call:
factanal(x = dad, factors = 8, scores = c("Bartlett"), rotation = "varimax")

Uniquenesses:
ï..s3    s7   s11   s15   s16   ad1   ad2   ad3   ad4   ad5   ad9  ad13  ad14  ad15  ad19
0.341 0.005 0.484 0.079 0.005 0.103 0.308 0.098 0.005 0.090 0.108 0.015 0.200 0.005 0.436

Loadings:
      Factor1 Factor2 Factor3 Factor4 Factor5 Factor6 Factor7 Factor8
ï..s3  0.295  -0.250   0.202   0.140           0.266           0.612
s7            -0.407   0.865  -0.135           0.163  -0.153         
s11           -0.426   0.541                   0.118                 
s15    0.930                          -0.164                  -0.139
s16    0.247   0.561  -0.103           0.771                         
ad1            0.910  -0.178           0.158                         
ad2            0.781  -0.229   0.101          -0.101                 
ad3            0.875  -0.202           0.134           0.248         
ad4   -0.134   0.538  -0.267                           0.777         
ad5    0.887  -0.156                                           0.295
ad9            0.128           0.926                                 
ad13   0.894   0.258          -0.227   0.235                         
ad14   0.813           0.183           0.103   0.256  -0.151         
ad15   0.389  -0.155   0.397                   0.744           0.305
ad19   0.622   0.132  -0.123   0.138   0.149                   0.312

               Factor1 Factor2 Factor3 Factor4 Factor5 Factor6 Factor7 Factor8
SS loadings      3.832   3.365   1.507   0.997   0.772   0.771   0.751   0.722
Proportion Var   0.255   0.224   0.100   0.066   0.051   0.051   0.050   0.048
Cumulative Var   0.255   0.480   0.580   0.647   0.698   0.750   0.800   0.848

Test of the hypothesis that 8 factors are sufficient.
The chi square statistic is 12.7 on 13 degrees of freedom.
The p-value is 0.471 

Muito obrigado a todos.


Bruce