similar to: Help with SEM package: Error message

Hi all, I am trying to set up a simple path analysis in the SEM package, but I am having some trouble. I keep getting the following error message or something similar with my model, and I'm not sure what I'm doing wrong: Error in solve.default(C) : system is computationally singular: reciprocal condition number = 2.2449e-20 In addition: Warning message: In sem.default(ram = ram, S = S,

I have been trying to use lavaan (version 0.4-7) for a simple path model, but the program seems to be computing far less degrees of freedom for my model then it should have. I have 7 variables, which should give (7)(8)/2 = 28 covariances, and hence 28 DF. The model seems to only think I have 13 DF. The code to reproduce the problem is below. Have I done something wrong, or is this something I

Hello, I've been unsuccessfully trying to reproduce a sem from Grace et al. (2010) published in Ecological Monographs: http://www.esajournals.org/doi/pdf/10.1890/09-0464.1 The model in question is presented in Figure 8, page 81. The errors that I've been getting are: 1. Using a correlation matrix: res.grace <- sem(grace.model, S = grace, N = 190) Warning message: In sem.default(ram

Hello all I'm trying to run some path analysis in either sem or lavaan (preferably lavaan because I find its interface easier to use). Most of my variables are continuously distributed and fairly well-behaved but I have a single exogenous variable (sex) which is not continuously distributed. Preliminary model fitting suggests that there aren't any sex by (anything else) interactions. The

Hi all, I am attempting to learn my way through the sem package by constructing a simple structural model for some of my data on bird diversity, abundance, and primary productivity. I have constructed a covariance matrix between these variables as per the following: >S_matrix = matrix(c( >+ 0.003083259, 0, 0, >+ 0.143870284, 89.7648490, 0, >+ 0.276950919,

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

Hi, we are trying to do structural equation modelling on R. However, one of our predictor variables is categorical (smoker/nonsmoker). Now, if we want to run the sem() command (from the sem library), we need to specify a covariance matrix (cov). However, Pearson's correlation does not work on the dichotomous variable, so instead we produced a covariance matrix using the Spearman's (or

Hi, New to both R and SEM, so this may be a very simple question. I am trying to run a very simple path analysis using the sem package. There are 2 exogenous (FARSCH, LOCUS10) and 2 endogenous (T_ATTENT, RMTEST) observed variables in the model. The idea is that T_ATTENT mediates the effect of FARSCH and LOCUS10 on RMTEST. The RAM specification I used is FARSCH -> T_ATTENT, y1x1, NA

Hi there, as a student sociology, I'm starting to learn about SEM. The course I follow is based on LISREL, but I want to use the SEM-package on R parallel to it. Using LISREL, I found it to be very usable to be able to see the total direct and total indirect effects (standardized and unstandardized) in the output. Can I create these effects using R? I know how to calculate them

Dear All, I am new for R and SEM. I try to fit the model with Y (ordinal outcome), X (4 categorical data), M1-M3 (continuous), and 2 covariates (Age&sex) as a diagram. library(polycor) model.ly <-specify.model() 1: x -> m1, gam11, NA 2: x -> m2, gam12, NA 3: x -> m3, gam13, NA 4: age -> m1, gam14, NA 5: age -> m2, gam15, NA 6: age -> m3, gam16, NA 7: sex -> m1,

Dear Helplist, I want to produce forecasts from a structural equation model. With the SEM package the model setup and its estimation is possible. However, I have not figured out how to obtain ex ante forecasts, i.e. applying the Gauss-Seidel algorithm to the estimated structural equations for provided values of the exogenous variables (i.e.: y_t = -inv(A)*B*x_t). Does anyone know if the there is

Hello all, I have an issue where I am generating data and trying to confirm the estimates using a sem. I keep getting an error about the degree of freedom being negative "Error in sem.default(ram, S = S, N = N, raw = raw, data = data, pattern.number = pattern.number, : The model has negative degrees of freedom = -6" Can someone explain this error or tell me what is wrong with my

Dear all, How do I get the standardized covariance (the correlation) between two latent variables? 'standardized.coefficients' gives standardized path coefficients, but not covariances. The covariance estimates are easily obtained from fit$coeff or 'summary', but EQS reports both the covariance and the correlation, how can I get that? best wishes, Mike

Can somebody help me with the following issue (SEM in R), please:   When I run the model (includes second order models) in R, it gives me the following:   1)       In sem.default(ram = ram, S = S, N = N, param.names = pars, var.names = vars,  :   Could not compute QR decomposition of Hessian. Optimization probably did not converge.   2)       I have aliased parameters and NaNS   or sometimes when

Dear R users, Jarrett Byrnes and I would like to announce version 2.0-0 of the sem package for fitting observed- and latent-variable structural equation models. This is a general reworking of the original sem package (which is still available on R-Forge as package sem1). Some highlights of sem 2.0-0 include: o More convenient and compact model specification, including the default automatic

Dear R users, Jarrett Byrnes and I would like to announce version 2.0-0 of the sem package for fitting observed- and latent-variable structural equation models. This is a general reworking of the original sem package (which is still available on R-Forge as package sem1). Some highlights of sem 2.0-0 include: o More convenient and compact model specification, including the default automatic

Hi, I am doing a hierarchical CFA using the sem package. I have 20 items, and I have 2 factors (F3 and F4), and also F1 and F2 are nested within F3. Here is the code that I have, but it is giving me an error message "Warning message: In eval(expr, envir, enclos) : Negative parameter variances. Model may be underidentified." and a further error "Error in summary.objectiveML(cfa,

Dear Spencer and Prof. Fox, Thank you for your replies. I'll very appreciate, if you have any ideas concerning the problem described below. First, I'd like to describe the model in brief. In general I consider a model with three equations. First one is for annual GRP growth - in general it looks like: 1) GRP growth per capita = G(investment, migration, initial GRP per

Greetings, I am working in the sem package on a model with 3 exogenous variables (2 are nominal-categorical), and 4 endogenous, continuous variables. To use sem with the nominal variables, I created dummy variables. Now, in my sem output I have estimates for path coefficients for the relationship between each level of the nominal variables and the endogenous variables they are associated

Hi, I am using tsls function from sem package to estimate a model which includes large number of data. Among its predictors, it has a nominal data which has about 10 possible values. So I expand this parameter into 9-binary-value predictors with the coefficient of base value equals 0. I also have another continuous predictor. The problem is that, whenever I run the tsls, I will get 'System