similar to: SEM modelling

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

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

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

Hi all, For computational reasons, I need to estimate an 18x18 polychoric correlation matrix two variables at a time (rather than trying to estimate them all simultaneously using ML). The resulting polychoric correlation matrix I am getting is non-positive definite, which is problematic because I'm using this matrix later on as if it were a legitimately estimated correlation matrix (in order

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

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

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 =

Hi, I would like to use the sem package to perform a path analysis (no latent variables) with a mixture of 2 nominal exogenous, 1 continuous exogenous, and 4 continuous endogenous variables. I seek advice as to how to calculate the appropriate covariance matrix for use with the sem package. I have read through the polycor package, and am confused as to the use of "numeric" for

Dear all, Has someone implemented in R (or any other language) Knol DL, ten Berge JMF. Least-squares approximation of an improper correlation matrix by a proper one. Psychometrika, 1989, 54, 53-61. or any other similar algorithm? Best regards Jens Oehlschl?gel Background: I want to factanal() matrices of polychoric correlations which have negative eigenvalue. I coded Highham 2002

Hello, I'd like to know if there is a way to fit a Structural equation model with dichotomous indicators (ex: problem with a phone solved/ or not) having effects on a ordinal variable. How I do that using R? Do you have an example with the code in R that you can send to me? Thanks a lot! Renata Estrella UFRJ, Brasil, Rio de Janeiro Renata Leite Estrella Assistente de

Hi I am looking for a structural equation modeling package in R which can be used for categorical data. Is anyone aware of the existence of such a package? Would appreciate any help on this. Thank you Upasna -- --------------------------------------------------------------------- Upasna Sharma Research Scholar Shailesh J. Mehta School of Management, Indian Institute of Technology, Bombay

Someone can help me? I tried several things and always don't converge # Model library(sem) dados40.cov <- cov(dados40,method="spearman") model.dados40 <- specify.model() F1 -> Item11, lam11, NA F1 -> Item31, lam31, NA F1 -> Item36, lam36, NA F1 -> Item54, lam54, NA F1 -> Item63, lam63, NA F1 -> Item65, lam55, NA F1 -> Item67, lam67, NA F1 ->

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

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

Hi, I'm running this for a simulation study, so many combinations of parameter produce many predictions that I need to correlate with data. The problem ---------------- I'm using rating data with 3 to 5 categories (e.g., too low, correct, too high). The underlying continuous scales should be normal, so I chose the polychoric correlation. I'm using library(polychor) in its

I am trying to run a confirmatory factor analysis using the SEM package. My data are ordinal. I have read http://socserv.mcmaster.ca/jfox/Misc/sem/SEM-paper.pdf. When I apply the hetcor function, I receive the following error: Error in checkmvArgs(lower = lower, upper = upper, mean = mean, corr = corr, : at least one element of 'lower' is larger than 'upper' Example:

Hi, I'm trying to import categorical data from SPSS to R using the script: xxx <-spss.get("xxx.por", use.value.labels=TRUE) but unfortunately am getting an error message 'error reading portable-file dictionary'. I have successfully imported data in the past. What could be the problem with this data? Thanks Simo

Dear R users, I am a psychology postgraduate student who is relatively new to using R. I am currently developing a psychometric scale and have run into a few problems when using R to calculate a polychoric correlation matrix for my dataset. I am trying to produce a polychoric correlation matrix for calculating ordinal reliability estimates (eg. Alpha, omega).The set consists of 439 observations

Sorry if this post should be long but I tried to give you a piece of my data to reproduce my error message using hetcor: Fehler in result$rho : $ operator is invalid for atomic vectors Zus?tzlich: Warning messages: 1: In polychor(x, y, ML = ML, std.err = std.err) : 1 row with zero marginal removed 2: In polychor(x, y, ML = ML, std.err = std.err) : the table has fewer than 2 rows Error in

I hope someone here knows the answer to this since it will save me from delving deep into documentation. Based on 22 pairs of vectors, I have noticed that tetrachoric correlation coefficients in stata are almost uniformly higher than those in R, sometimes dramatically so (TCC=.61 in stata, .51 in R; .51 in stata, .39 in R). Stata's estimate is higher than R's in 20 out of 22