Hi > -----Original Message----- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@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 <- 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]] > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting- > guide.html > and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
