similar to: reliability, scale scores in the psych package

Hello I have 2 data frames DF1 and DF2 where DF2 is a subset of DF1: DF1= data.frame(V1=1:6, V2= letters[1:6]) DF2= data.frame(V1=1:3, V2= letters[1:3]) How do I create a new data frame of the difference between DF1 and DF2 newDF=data.frame(V1=4:6, V2= letters[4:6]) In my real data, the rows are not in order as in the example I provided. Thanks much Joseph [[alternative HTML version

Hello, I am using function princomp from the package psych. I have my principle component object mypc: mypc <- princomp(covmat=mycor) plot(mypc) # shows me a screeplot Question: how could I actually see the values displayed in the screeplot. I don't mean on the graph - I just want to know the actual value for each component (e.g., 10, 3.2, 1.8, etc.) I need to know how much variance,

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 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

--Multipart_Thu_Oct__9_10:41:03_1997-1 Content-Type: text/plain; charset=US-ASCII In case you did not realize how much this is related to R : --Multipart_Thu_Oct__9_10:41:03_1997-1 Content-Type: message/rfc822 Return-Path: s-sender@utstat.toronto.edu From: "Steven M. Boker" <sboker@calliope.psych.nd.edu> Date: Wed, 8 Oct 97 16:37:05 -0500 To: s-news@utstat.toronto.edu Subject:

This is about my searchable archive. The function RSiteSearch() now searches this archive. I'm considering the following changes. If you have comments, please write me. Try to avoid cc'ing the list. The Rhelp02a directory, which now contains all list mail from 2002 (about 100 MB), is getting larger and larger. This probably cannot go on forever, and performance might even improve if

Dear List, I sent a related message yesterday, but did not receive it through the mailing list; on the R-help archives it reads as "An embedded and charset-unspecified text was scrubbed...". So here it is again, with the (little) progress I have made since then. Any help would be greatly appreciated! I am working with a relatively large correlation matrix (~1600*1600), which I am

Hello: I s there a way to use xtable with objects from the psych package, particularly principal()? Is there a difference between princomp and principal? xtable seems to play better with princomp. Thank you. Yours, Simon Kiss ********************************* Simon J. Kiss, PhD Assistant Professor, Wilfrid Laurier University 73 George Street Brantford, Ontario, Canada N3T 2C9

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! I am running principle components analysis using princomp function in pacakge psych. mypc <- princomp(mydataforpc, cor=TRUE) Question: I'd like to use pairwise deletion of missing cases when correlations are calculated. I.e., I'd like to have a correlation between any 2 variables to be based on all cases that have valid values on both variables. What should my na.action be in

Hello everyone, Does anyone know what would cause the skewness() function (from e1071), as well as skew() from psych, to return a value of NaN? I have a vector of positively-skewed data (https://docs.google.com/file/d/0B6-m45Jvl3ZmYzlHRVRHRURzbVk/edit?usp=sharing) which these functions return a value for like normal: > skewness( data ) # returns 1.400405 but when I instead give those

My search site, http://finzi.psych.upenn.edu, has had several problems recently, all my fault, for which I apologize. But it now seems to be running reliably, on a new computer that is much faster than the old one. It uses htdig to permit search of the Rhelp mailing list, R documents, R functions, and various combinations of these. Search has several options, including Boolean search (with AND,

Hola a tod@s, ayer mismo gracias a CarlosJ. Gil Bellosta me he enterado de la existencia de la lista R en España. Excelente idea!! Aprovecho este mensaje para comentaros un problema en la utilización del paquete psych y Rgraphviz que nunca había utilizado antes. Despues de cambiar a la versión R.2.8.2 y haberla reinstalado varias veces con vista. Cuando ejecuto test.psych() aparece un mensaje

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 there, I wanted to use RedCloth with my RubyOnRails 3 application, but it seems that RoR uses 1.9.2''s Psych by default as the YAML parser if available, and Psych fails to load the LaTeX entities file: ==================================================================== #ruby -v: ruby 1.9.2p136 (2010-12-25 revision 30365) [x86_64-linux] irb(main):001:0> Dir.pwd =>

I have two hosts which have their respective keys loaded into each others .ssh/authorized-keys2 file for root. I want to move a directory tree from one host to the other via rsync to maintain a shadow structure of the application provided on the main host. This is the relevant part of the command: rsync -avz --rsh=ssh --delete-after /var/data/pas-redmine root at

Dear all, I like the way the Rcmdr package computes reliability. E.g reliability(cov(d[,c("q1", "q2", "q3", "q4", "q5", "q6")], use="complete.obs")) will not only give me the alpha score, but also for each variable, alpha.score if deleted. However - when writing scripts it's very tiresome to load the whole Rcmdr GUI just

Dear R help I have conducted a fa() analysis, and I want to use fa.diagram to assess the extent to which the 11 latent factors predict the 37 items in a psychological battery. However, the display on the screen has very large font size for the coefficients of the relationship between the 11 factors and the 37 items, so the numbers overlap and are therefore illegible. When I output this to a

Dear R users, I am trying to simulate some multitrait-multimethod models using the packages sem and psych but whatever I do to deal with models which do not converge I always get stuck and get error messages such as these: "Error in summary.sem(M1) : coefficient covariances cannot be computed" "Error in solve.default(res$hessian) : System ist f?r den Rechner singul?r: reziproke

Dear Professor Fox, Apologies for my oversight relating to the polychor command and thank you for your advice. I turned to the polychor command when trying to find an equivalent for the polychoric command found in the psych package (I am following a procedure outlined in Gadermann, Guhn & Zumbo, 2012 that uses this command rather than the polycor package) . The polychoric command is returning