similar to: estimate the reliability of a scale with dichotomous items

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

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

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

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

Hello, I was wondering if someone could tell me how to output, to file, the residuals from a REML model-fit. The type of residuals I am interested in are the simple "original raw values - model fit" type. Thanks in advance, Josh B. [[alternative HTML version deleted]]

Dear List, I use Rgraphviz for display of graphs in some packages. Since Rgraphviz is no longer on CRAN it needs to be installed from Bioconductor and that is fine, but I have trouble figureing out the following: I create a plot method which - if Rgraphviz is installed - uses Rgraphviz for displaying and otherwise does nothing. This is implemented as: if (!("package:Rgraphviz" %in%

Hi R, Let, x=1:80 I want to sum up first 8 elements of x, then again next 8 elements of x, then again another 8 elements..... So, my new vector should look like: c(36,100,164,228,292,356,420,484,548,612) I used: aggregate(x,list(rep(1:10,each=8)),sum)[-1] or rowsum(x,group=rep(1:10,each=8)) But without grouping, can I achieve the required? Any other ways of doing

Dear Users, I ran factor analysis using R and SAS. However, I had different outputs from R and SAS. Why they provide different outputs? Especially, the factor loadings are different. I did real dataset(n=264), however, I had an extremely different from R and SAS. Why this things happened? Which software is correct on? Thanks in advance, - TY #R code with example data # A little

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 Kevin, In addition to the advice you've received, take a look at the survey package. It's not quite what you're asking for, but in fact it's probably more useful, in that it provides correct statistical inference for data collected in complex surveys. The package is described in an article, T. Lumley (2004), Analysis of complex survey samples, Journal of Statistical Software

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.

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

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

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

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