similar to: New package: polycor

A bit late, but you might like to look at http://www.qimr.edu.au/davidD/polyr.R Regarding the original posters queries: You can analyse polychoric correlations as if they were Pearson correlations using standard software (eg sem), and this usually doesn't do too badly, or go to AWLS (Browne) in LISREL etc, or ML analysis of the full multidimensional contingency table using programs such as

I have been using the hetcor function in the polycor package. When I don't specify the use option everything runs smoothly. However, when I specify use either as "pairwise.complete.obs" or "complete.obs" I get this error Error in optim(rho, f, control = control, hessian = TRUE, method = "BFGS") : non-finite value supplied by optim Is this an error in

Dear R-help list members, I forgot to copy my reply to the r-help list. Here's most of it. John -------------------------------- John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -------------------------------- -----Original Message----- From: John Fox [mailto:jfox at mcmaster.ca] Sent: Sunday,

I am currently working with R's polycor package and I have encountered a problem. I tried to follow the steps as outlined in the sem.pdf file where a CFA model is run using polychoric correlations. Every time I run the command sem(model, data, N=.), I get the following warning message: Warning message: In if (orthogonal) { : the condition has length > 1 and only the first element will be

Dear Roman, I've attached the latex file for the copy-edited version of your manuscript -- correcting English grammar and style, and making the paper conform more closely to the R News style. Can you check this over to confirm that I've not altered the meaning of your text and to approve the copy-edits? Thanks, John -------------------------------- John Fox Department of Sociology

Hello R Users, While I have attempted to dig into the R help files and I have not identified the answer to these questions, I apologize in advance if my questions were answered in the past. I also recognize that one of my questions unfortunately verges on statistical rather than code territory. I have two rather unrelated questions about using the sem and polycor packages for a relatively

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

Hi, I created a matrix of correlation coefficients and their respective scatterplots using the commands hetcor() and panel.cor(). Can I edit this matrix? I would like to have the scale values display at the top, as opposed to the alternating top-bottom scale approach. Also, I would like to increase the font size of the graphic for print purposes. Any help on how to accomplish this would be

Hello, I am running polychoric correlations on a dataset composed of 12 ordinal and binary variables (N =384), using the polycor package. One of the association (between 2 dichotomous variables) is very high using the 2-step estimate (0.933 when polychoric run only between the two variables; but 0.801 when polychoric run on the 12 variables). The same correlation run with ML estimate returns a

Dear R-users Could anyone tell me which library contains a function to compute polycoric correlations? I wonder the same question was asked a while ago, but I could not locate the mail in the R-help archives. Sorry for bothering you. Sincerely ------------------------ Hiroto Miyoshi ???? h_m_ at po.harenet.ne.jp

Hello, I have a data set containing categorical and ordinal factors, as well as sampling weights (i.e., survey weights reflecting unequal probabilities of selection). I want to fit a structural equation model with sem(). I have run sem() on weighted covariance matrices using advice from John Fox (see <http://tolstoy.newcastle.edu.au/R/e5/help/08/12/8773.html> and

Dear r-devel list members, I've prepared a version of the Rcmdr package that uses the new tools for internationalization and localization in R 2.1.x. For the present, this is the development version of the Rcmdr package (1.1-0), available at <http://socserv.socsci.mcmaster.ca/jfox/Misc/Rcmdr/index.html>. There is currently a French translation provided with the package, courtesy of

Dear R-devel list members, For several years, I've maintained for my private use and for demonstrations a set of functions [e.g., for indexing, with(), summary(), etc.] that employ the RODBC package and a database to provide a "pseudo" data frame. Row "indexes" can be specified as SQL queries. I know that the idea isn't original, and I haven't done anything with

Dear list members, I've just uploaded version 1.0-0 of the Rcmdr package to CRAN. For people who haven't seen the package before, the "R Commander" provides a basic-statistics graphical user interface to R, based on the tcltk package. The new version incorporates a number of improvements to the R Commander interface, as documented in the CHANGES file distributed with the

Dear list members, I've just uploaded version 1.0-0 of the Rcmdr package to CRAN. For people who haven't seen the package before, the "R Commander" provides a basic-statistics graphical user interface to R, based on the tcltk package. The new version incorporates a number of improvements to the R Commander interface, as documented in the CHANGES file distributed with the

Hi list members I have a vector of strings x=c("w","ex","ee") And I want to get a logical vector showing the positions where my search string "e" matches the elements partially, i.e. is at least the left-hand part of the target strings, i.e. I want to get a vector FALSE TRUE TRUE. Any ideas? Thanks Steve Powell proMENTE social research research |

Hi there, Am I correct to believe that tetrachoric correlation is a special case of polychoric correlation when there are only two levels to the ordered factor? Thus it should be okay to use hetcor from the polycor package to build a matrix of correlations for binary variables? If this is true, how can one estimate 95% confidence intervals for the correlations? My guess would be mat =

Hello, It has been found that SEM analysis using polychoric correlations + maximum likelihood estimator produces incorrect test statistics and standard errors (e.g., Flora, D. B., & Curran, P. J. (2004). An Empirical Evaluation of Alternative Methods of Estimation for Con?rmatory Factor Analysis With Ordinal Data. Psychological Methods, 9(4), 466-491). Standard errors can be dealt with by

Dear r-help members, Can lmer() in the lme4 package fit models that have a specified within-group correlation structure, as provided, for example, by the correlation argument to lme() in the nlme package? Thanks, John -------------------------------- John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Dear list members, I think that it's a pity that there's no Windows binary for the rgl package on CRAN. I understand that rgl doesn't compile with the standard Windows tools and that there's a Windows binary on Daniel Adler's web site. But the lack of an rgl binary on CRAN provides presents a potential obstacle to at least some Windows users, such as many students in