similar to: how to calculate factor Scoring Coefficients

Hi I have used fa() to perform a factor analysis of a psychological battery which is thought to have 11 factors. I can identify which factors the loadings relate to easily enough because I can see which items are loading onto each of the columns in the $loading output. However, how can I identify which items or loadings are being used to create each of the columns in the $scores output? I have

Hi everyone, I am trying to get the factor score for each individual case from a principal component analysis, as I understand, both princomp() and prcomp() can not produce this factor score, the principal() in psych package has this option: scores=T, but after running the code, I could not figure out how to show the factor score results. Here is my code, could anyone give me some advice please?

Dear Professor Revelle and R-helpers, This is a two-part question: the first part is general, and the second is specific to the psych package. First question: In the past I've constructed composite variables from questionnaire data by using rowMeans(), and then correlating items with the scale using cor() as an informal check for "bad" items. Over the weekend I decided to take a

Dear R users, When I tried to use princomp() from stats packages to do Principal Components Analysis, I am not very clear what is the "scale". And the scores are different from "PROC PRINCOMP" procedure from SAS. Using the example data from this package: restpc <- princomp(USArrests, cor = TRUE) > restpc$scale Murder Assault UrbanPop Rape 4.311735 82.500075

Hi, I am trying to calculate Principal Component Scores per id per year using the psych package. The following lines provide the scores per obeservation pca = data.frame(read.table(textConnection(" id year A B C D 1001 1972 64 56 14 23 1003 1972 60 55 62 111 1005 1972 57 51 10 47 1007 1972 59 49 7 10 1009 1972 65 50 9 32 1011 1972 52 58 3 5 1013

6.4.1 Estimation of fixed effects Heterogeneous team ability is a possible explanation for the result in Section 6.3. That result simply indicates that the more goals a team scores, the higher the probability that it will score more. However, teams that can score more goals also indicate teams with greater ability, or just greater scoring ability, than their opponents. As mentioned in Section

Dear all, 1)how can I easily get p value for the coefficients of factors in a multinomial model? 2)why the p values for "type III" test with Anova are not identical to that from SAS? for example: A,B and C are categorical variables,but the proportions of each level in each categorical variables are not balance. Y is a nominal variables(>=3 categories); To do

Hi, I am doing predictive modelling of Multivariate Time series Data of a Motor in R using various models such as Arima, H2O.Randomforest, glmnet, lm and few other models. I created a function to select a model of our choice and do prediction. Model1 <- function(){ .. return() } Model2 <- function(){ ... return() } Model3 <- function(){ ... return() } main <-

Hi All, I'm generating 10 different data sets with 1 and 0 in a matrix form and writing the output in separate files. Now I need to stack all these data sets in one vector and I know that stack only operates on list or data frame however I got these data sets by converting list to a matrix so can't go backwards now. Is there a way i can still use Stack? Please see the program:

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 R-List, I have carried out a factor analysis using fa (psych) with nfactors=2, rotation="oblimin" and fm="pa". Now I have to report both pattern AND structure matrix. As I have understood R-Documentation, this can be obtained by setting the "oblique.scores" argument TRUE (structure matrix) or FALSE (pattern matrix), respectively. However both produced the same

Hi All, I'm generating 10 different data sets with 1 and 0 in a matrix form and writing the output in separate files. Now I need to stack all these data sets in one vector and I know that stack only operates on list or data frame however I got these data sets by converting list to a matrix so can't go backwards now. Is there a way i can still use Stack? Please see the program:

The partial R-square (or coefficient of partial determination, or squared partial correlation coefficients) measures the marginal contribution of one explanatory variable when all others are already included in multiple linear regression model. The following link has very clear explanations on partial and semi-partial correlation: http://www.psy.jhu.edu/~ashelton/courses/stats315/week2.pdf In

Sorry for these multiple postings. I solved the problem using na.omit() to drop records with missing values for the time being. I will worry about imputation, etc. later. I calculated the sum of squared errors for 3 models, linear regression, neural networks, and support vector machines. This is the first run. Without doing any parameter tuning on the SVM or playing around with the number of

Hello forum, I have developed an IRT model for a set of items on a large reference sample. I now want to use the IRT parameters that were estimated in the reference sample to score new individuals, one at a time. What R package(s) would allow me to do this? Also, several covariates were associated with item thresholds in the reference sample, so I was wondering if any programs can take

Haj i try to perform a principal component analysis by using a tetrachoric correlation matrix as data input tetra <- tetrachoric (image_na, correct=TRUE) t_matrix <- tetra$rho pca.tetra <- principal(t_matrix, nfactors = 10, n.obs = nrow(image_na), rotate="varimax", scores=TRUE) the problem i have is to compute the individual factor scores from the pca. the code runs perfect

Dear all I'm unable to find an example of extracting the rotated scores of a principal components analysis. I can do this easily for the un-rotated version. data(mtcars) .PC <- princomp(~am+carb+cyl+disp+drat+gear+hp+mpg, cor=TRUE, data=mtcars) unclass(loadings(.PC)) # component loadings summary(.PC) # proportions of variance mtcars$PC1 <- .PC$scores[,1] # extract un-rotated scores of

I am trying to simulate a dataset using Parallel Latent CTT model and this is what i have done so far: (START) #Importing psych library for all the simulation related functions library(psych) # Settting the working directory path to C:/NCME path="C:/NCME" setwd(path) #Using the function to generate the data GenData <- congeneric.sim(N=500, loads =

Hi!   I have been reading lots of queries regarding logistic regression as well as Credit scoring model and the related replies, which I must admit are full of wisdom and help us to understand where do we stand. The contributions of Frank E Harrell Jr, etc are invaluable. Incidentally I am also trying to work on Credit scoring model and initially I tried to use the multiple regression etc. and now

Hi, I am using acts_as_ferret and have a problem with scoring. I would like to organize it in such way that, if any of the searched terms fits, I get 1.0 score as a result. I will explain it on the example. I have in index: a) "one two three four" b) "one two three" c) "one two" d) "one" When I search for "one" I would like to get 1.0 score for