similar to: How to get the descriptive statistic of the whole dataframe?

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

Dear list readers, I want to: 1. Get a table of basic descriptive statistics for my variables with the variable names one below the other like SPSS descriptive statistics: Varname N Min Max Mean SD xxxx x x x x x xxx x x x x x .... 2. Delete some variables from a data frame or exclude variables from beeing analyzed. 3. Create a text file / redirect the terminal output to a

The data is in the attachment. What I wanna get is: , , Sex = Male Eye Hair Brown Blue Hazel Green Black 32 11 10 3 Brown 53 50 25 15 Red 10 10 7 7 Blond 3 30 5 8 , , Sex = Female Eye Hair Brown Blue Hazel Green Black 36 9 5 2 Brown 66 34 29 14 Red 16 7 7 7 Blond 4

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 =

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

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

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

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

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

Hi all: There's a question about glht function. My data:data_ori,which inclue CD4, GROUP,time. f_GROUP<-factor(data_ori$GROUP) f_GROUP is a factor of 3 levels(0,1,2,3) result <- lme(sqrt(CD4) ~ f_GROUP*time ,random = ~time|ID,data=data_ori) glht(result, linfct = mcp(f_GROUP="Tukey") ) Error in `[.data.frame`(mf, nhypo[checknm]) : undefined columns selected I can't

I'm wanting to create a series of near-identical matrices via the addition of "white noise" to my starting matrix. Is there a function within R which will allow me to do this? Thank you Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596 fax: +44 113 343 6716 mail: laura at env.leeds.ac.uk

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!