Extracting factors from "factanal"
On 02/01/2013 10:34, PIKAL Petr wrote:
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.
If you give your observations names, the names would (AFAIK) be used for the scores.
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]]
______________________________________________ R-help at r-project.org<mailto:R-help at 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.
[[alternative HTML version deleted]]
______________________________________________ R-help at 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.
Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595