similar to: parametric bootstrap and computing factor scores

I am sure there is a simple way to do the following, but i haven't been able to find it. I am hoping a merciful soul on R-help could point me in the right direction. I am doing a factor analysis on survey data with missing values. to do this, I run: FA1<-factanal(na.omit(DATA), factors = X, rotation = 'oblimin', scores = 'regression') Now that I have my factors and

Hi, I wonder why factor scores produced by factanal are correlated, and I'd appreciate any hints from people that may help me to get a deeper understanding why that's the case. By the way: I'm a psychologist used to SPSS, so that question my sound a little silly to your ears. Here's my minimal example: *********************************************** v1 <-

I'm working with a data collected through complex survey design. My goal is to conduct a factor analysis to extract two a priori, known factors, and to get factor scores for these factors. Unfortunately, the "svyfactanal" procedure from the Survey package does not allow for the calculation of either Thompson regression scores or Bartlett scores. So, I found several sources that say

Dear expeRts, sorry for such a newbie question - in PCA/factor analysis e.g. in SPSS it is possible to save scores from the factors. Is it analogously possible to "save" the implied scores from the latent variables in a measurement model or structural model e.g. using the sem or lavaan packages, to use in further analyses? Best wishes Steve Powell www.promente.org | skype

Dear friends, I'm using R to produce the following Factor Analysis: > matriz.cor<-hetcor(matrix(as.factor(data), ncol=variables, byrow=T))$correlations > factanal(x=data, factors=2, covmat=matriz.cor, scores='regression') Then the screen output shows the following message: Error en factanal(x = data, factors = 2, covmat = matrix, : requested scores without

I was wondering if there is a way of predicting factor scores of new data for factor analysis in R (similar to "predict" in S-plus). So far I have not been able to find it, nor found reference to it. Thanks, Chris Fonnesbeck -- ~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~* Christopher J. Fonnesbeck Ph.D. Student Georgia Cooperative Wildlife Unit University of Georgia Athens, GA 30602

Greetings to you all, I am performing a semi parametric bootstrap in R on a Gamma Distributed data and a Binomial distributed data. The main challenge am facing is the fact that the residual variance depends on the mean (if I am correct). I strongly feel that the script below may be wrong due to mean-variance relationship #####R code####### fit1s

Hi I am trying to find a parametric bootstrap confidence interval and when I used the boot function I get zero bias and zero st.error? What could be my mistake? Thank you and take care. Laila [[alternative HTML version deleted]]

In R factanal output, I can't find a function to give me residuals e. I mannually got it by using x -lamda1*f1 -lamda2*f2 - ... -lamdan*fn, but the e I got are not uncorrelated with all the f's. What did I do wrong? Please help. Yijun ____________________________________________________________________________________ Be a better friend, newshound, and

Hi, I had a look at the MASS4 scripts in the MASS package, in Ch 11.3 Factor Analysis, there is a section of codes like: data(ability.cov) ability.FA <- factanal(covmat = ability.cov, factors = 1) ability.FA (ability.FA <- update(ability.FA, factors = 2)) #summary(ability.FA) round(loadings(ability.FA) %*% t(loadings(ability.FA)) + diag(ability.FA$uniq), 3)

--=-hiYzUeWcRJ/+kx41aPIZ Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: 8bit Hi, on March 10 I filed a wishlist bug report asking for the inclusion of some changes to factanal() and the associated print method. The changes were originally proposed by John Fox in 2005; they make print.factanal() display factor correlations if factanal() is called with rotation =

I am very new to factor analysis as well as R. I am trying to run a factor analysis on the residual returns on common stock (residual to some model) and trying to determine if there are any strong factors remaining. After running factanal, I can obtain the factor loadings but how do I get the values of the factor returns themselves? In other words if the relationship is r = lambda * f I

Dear all, I have a series of variables that looks roughly like the sample data below and I'm trying to conduct a factor analysis. I've omitted cases with missing values for the factor analysis, but now I'd like to use the scores on each component as new variables in the *original* data set for analysis. That is, I'd like to take the scores on each of the two factors and see how

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

I've read through many postings about principle component analysis in the R-help archives, but haven't been able to piece together the information I need. I'd like to recreate an SPSS-like experience of factor analysis using R. Here's what SPSS produces: 1. Scatterplots of all possible variable pairs, with regression lines. xyplot(my.dataframe) is perfect but for the lack of

Dear all, i wonder if there?s a command to obtain the actual values of a principal component or a factor (not as.factor, but factanal) . test=princomp(USArrests, cor = TRUE) summary(test) just outputs, standard deviation, Prop of Variance and cumulative proportion of variance. test$loadings offers yet another proportion of variance scheme. why is that? Apart from that: Is there a

I am having difficulty obtaining the scores from my principal component analysis. I have used this method before and have had no problems. The data set that I am using this time is similar to what I have used in the past. What do I need to do to my dataset in order for me to obtain these scores? R screen says the following message Error in factanal(covmat = pasa.cov, factors = 4) :

Hi, I was wondering whether this query was addressed on how to perform validation through boostrapping. I am currently trying to implement a boostrapping approach to validation but don't know where to start. Help please. Thank you and Regards, Vivienne Ozohili Risk Model Validation Manager Group Risk Independent Control Unit A5, Bank of Ireland Head Office Lower Baggot Street Dublin 2 Tel:

> Does factanal() force the user to use the formula interface if they > wish to specify an na.action? Short answer: yes. Long answer: The handling of na.action is a built in part of the formula processing functions, so it's automatic when dealing with a formula. There are also downstream effects on predict() and resid() that are worked out for the formula case, but aren't clear

hi, Does factanal() force the user to use the formula interface if they wish to specify an na.action? v1 <- c(1,1,1,1,1,1,1,1,NA,1,3,3,3,3,3,4,5,6) v2 <- c(1,2,1,1,1,1,2,1,2,1,3,NA,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,NA,3,1,1,2,1,1,1,1,2,NA,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