similar to: Psych package: fa.diagram, how to re-arrange layout so numbers do not over-write each other

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

Hello, I searched the archives but found no answers. How to modify the hisgram color of function pairs.panels of Psych package ? I tried col() but it was the line color modified. Thanks. Elaine [[alternative HTML version deleted]]

Dear list A strange bug in the psych package is due to the behavior of the diag function: It gives the expected values for 1, a vector (-1,1), but not for -1 Is this a known feature? > diag(1) [,1] [1,] 1 > diag(c(-1,1)) [,1] [,2] [1,] -1 0 [2,] 0 1 > diag(-1) Error in diag(-1) : invalid 'nrow' value (< 0) Bill William Revelle

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 =

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

Hello! I would like to draw a graph like the following: http://www.optics.rochester.edu/workgroups/cml/opt307/spr04/pavel/plot_small.jpg Aim is to plot confidence intervals of treatments for X(=response1) and Y(=response2) axis simultaneously to visualize aggreement of confidence interval for two responses. Can anyone please provide me some direction to start with? Thanks! -- Kim.

Hi All, I am using fleiss kappa for inter rater agreement. Are there any know issues with Fleiss kappa calculation in R? Even when I supply mock data with total agreement among the raters I do not get a kappa value of 1. instead I am getting negative values. I am using the irr package version 0.70 Any help is much appreciated. Thanks and Regards M [[alternative HTML version deleted]]

On 26.04.2015 13:23, William Revelle wrote: > This still seems to be the case. > > I tried uploading the most recent version of psych and got as far as the Step 3 page which says that it has sent out an email. > However, although normally this step takes agout 1 minute, nothing has happened for 16 hours. > > Should we just use the old system of uploading to the ftp site and

Andrie noticed that first, and I can confirm: from our end, it looks as if the backend to http://cran.r-project.org/submit.html is currently down. Dirk -- http://dirk.eddelbuettel.com | @eddelbuettel | edd at debian.org

Hi! I am trying to analyse with R a database that I have previously analysed with SPSS. Steps with SPSS: Factorial analysis Extraction options : I select = Principal component analysis Rotation: varimax Steps with R: I have tried it with varimax function with factanal or with princomp...and the results are different of what I have with SPSS. I think that varimax function is incorporated in

Dear fellow developers, Is there an easy way to determine how big my packages (psych and psychTools) will be on various versions of CRAN? I have been running into the dread 'you are bigger than 5 MB" message for some installations of R on CRAN but not others. The particular problem seems to be some of the mac versions (specifically r-oldrel-macos-arm64 and r-release-macos-X64 ) When

Dear friends, How do I stop partial matching of list names? e.g., x <- list(AAAA="aaaaa", BBBBB="bbbbb") is.null(x$A) #returns FALSE even though there is no element A. if(is.null(x$A)) {result <- x$BBBB} else {result <- x$A} result #is aaaa even though there is no x$A element x <- list(CCCC="aaaaa", BBBBB="bbbbb") if(is.null(x$A))

Hello, I need to do a principal component analysis with EQUAMAX-rotation. Unfortunately the function principal() I use normally for PCA does not offer this rotation specification. I could find out that this might be possible somehow with the package GPArotation but until now I could not figure out how to use this in the principal component analysis. Maybe someone can give an example on how to do

hi, I have following data and code; cov <- c (1.670028 ,-1.197685 ,-2.931445,-1.197685,1.765646,3.883839,-2.931445,3.883839,12.050816) cov.matrix <- matrix(cov, 3, 3, dimnames=list(c("y1","x1","x2"), c("y1","x1","x2"))) path.model <- specify.model() x1 -> y1, x1-y1 x2 <-> x1, x2-x1 x2 <->

Dear R developers, Is there a way to pretest a package on the Solaris-sparc and solaris-x86 systems equivalent to the win-builder check? My psych package (1.1.10 and 1.1.11) passes all checks for the Mac on my system, on the win-builder checking system supported by Uwe, and then passes Kurt's tests to install on CRAN. But it then fails when being built for the solaris systems. Rather

Dear R list, I have discovered a seemingly peculiar feature when using a matrix to index itself (yes, this is strange code, which I have now modified to be more reasonable). #this makes sense s <- matrix(1:3,nrow=1) s[s] #all three elements are shown #but when I try s <- matrix(1:2,nrow=1) s[1] #fine, the first element is shown s[2] #fine, the second element is

Hello R users, and perhaps William Revelle in particular, I'm curious as to how ICC deals with missing data, so for example you are sampling individuals over set periods in time and one individual is missing or was not recaptured at that given time point - leading to NA in the dataset. My thought was that it should then omit data by individual, but I'm not convinced that that is what it

hello, subsequently to a NMDS analysis (performed with metaMDS or isoMDS) is it possible to rotate the axis through a varimax-rotation? Thanks in advance. Bernd Panassiti

hi, How can I compute a reliability score of a scale consisting only of dichotomous items? thanks for any help!

Hello,   In SPSS, a multiple regression can be conducted by inputting the means, standard deviations, sample size, and correlation matrix without actually using the raw dataset. Is it possible to do the same in R?   Thanks in advance for your assistance.   Linda [[alternative HTML version deleted]]