Extracting factors from "factanal"

-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
project.org] On Behalf Of Virgile Capo-Chichi
Sent: Wednesday, January 02, 2013 9:33 AM
To: R help
Subject: [R] Extracting factors from "factanal"

Dear R users
Happy New year to all for a start. Below is some data that I ran a
factor analysis on. Using $score prints the scores for each of the
three factors.
However, I would like to access those factors as variable for new
computations. How do I do that? In SPSS we just call fact1_1, fact2_1
etc..
Thanks for your suggestions. V

============
v1 <- c(1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,4,5,6)
v2 <- c(1,2,1,1,1,1,2,1,2,1,3,4,3,3,3,4,6,5)
v3 <- c(3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,5,4,6)
v4 <- c(3,3,4,3,3,1,1,2,1,1,1,1,2,1,1,5,6,4)
v5 <- c(1,1,1,1,1,3,3,3,3,3,1,1,1,1,1,6,4,5)
v6 <- c(1,1,1,2,1,3,3,3,4,3,1,1,1,2,1,6,5,4)
m1 <- cbind(v1,v2,v3,v4,v5,v6)
factanal(m1, factors = 3,scores = "Bartlett")$scores

myscores[,1]

myscores[,"Factor1"]

	[[alternative HTML version deleted]]

Hi
I am not at all an expert in factor analysis but I would say that

myscores[1,1]

is a value of factor1 for the first row of your data frame. How would you like the link to be performed. Easiest way seems to me cbind your data frame with myscores.

m2 <- cbind(m1, myscores)

However grater expert than myself shall confirm that my assumptions are correct.

Regards

Petr

From: Virgile Capo-Chichi [mailto:vcapochichi at gmail.com]
Sent: Wednesday, January 02, 2013 11:20 AM
To: PIKAL Petr
Cc: R help
Subject: Re: [R] Extracting factors from "factanal"

Thanks Petr, very useful. Still, I need to link the scores to individual observations in my original data frame. Is that possible? V
2013/1/2 PIKAL Petr <petr.pikal at precheza.cz<mailto:petr.pikal at precheza.cz>>
Hi

-----Original Message-----
From: r-help-bounces at r-project.org<mailto:r-help-bounces at r-project.org> [mailto:r-help-bounces at r-<mailto:r-help-bounces at r->
project.org<http://project.org>] On Behalf Of Virgile Capo-Chichi
Sent: Wednesday, January 02, 2013 9:33 AM
To: R help
Subject: [R] Extracting factors from "factanal"

Dear R users
Happy New year to all for a start. Below is some data that I ran a
factor analysis on. Using $score prints the scores for each of the
three factors.
However, I would like to access those factors as variable for new
computations. How do I do that? In SPSS we just call fact1_1, fact2_1
etc..
Thanks for your suggestions. V

============
v1 <- c(1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,4,5,6)
v2 <- c(1,2,1,1,1,1,2,1,2,1,3,4,3,3,3,4,6,5)
v3 <- c(3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,5,4,6)
v4 <- c(3,3,4,3,3,1,1,2,1,1,1,1,2,1,1,5,6,4)
v5 <- c(1,1,1,1,1,3,3,3,3,3,1,1,1,1,1,6,4,5)
v6 <- c(1,1,1,2,1,3,3,3,4,3,1,1,1,2,1,6,5,4)
m1 <- cbind(v1,v2,v3,v4,v5,v6)
factanal(m1, factors = 3,scores = "Bartlett")$scores

myscores <- factanal(m1, factors = 3,scores = "Bartlett")$scores

myscores[,1]

  [1] -0.9039949 -0.8685952 -0.9082818 -1.0021975 -0.9039949 -0.7452711
  [7] -0.7098714 -0.7495580 -0.8080740 -0.7452711  0.9272282  0.9626279
[13]  0.9229413  0.8290256  0.9272282  0.4224366  1.4713902  1.8822320

or

myscores[,"Factor1"]

  [1] -0.9039949 -0.8685952 -0.9082818 -1.0021975 -0.9039949 -0.7452711
  [7] -0.7098714 -0.7495580 -0.8080740 -0.7452711  0.9272282  0.9626279
[13]  0.9229413  0.8290256  0.9272282  0.4224366  1.4713902  1.8822320

selects a factor score. myscore is a matrix so you cannot use $ notation. However you can convert it to data frame.

Regards
Petr

       [[alternative HTML version deleted]]

	[[alternative HTML version deleted]]