Error with sem function df = -6

Dear Immanuel,

In the absence of more information, your question is unanswerable, since there's no way of knowing which variables in your model are observed and which are latent. It was my guess that Y1, Z1, and Z2 (along with UNIT) are observed and CF is latent, but this is clearly wrong given the df calculation.

Here, as the posting guide suggests, is a reproducible example, using your model and call to sem(), but clearly defined data:

---- snip ---

set.seed(12345)
data <- matrix(rnorm(500*3), 500, 3)
S2 <- rawMoments(cbind(1, data))
rownames(S2) <- colnames(S2) <- c("UNIT", "Z1", "Z2", "Y1")

sem.m1<-sem(model=model.ram1, S=S2, N=500, fixed.x="UNIT", raw=T)

Error in sem.default(ram, S = S, N = N, raw = raw, data = data, pattern.number = pattern.number,  : 
  The model has negative degrees of freedom = -2

---- snip ---

This model is underidentified, but at least it's possible to figure out why: There are 4*5/2 = 10 unique observed moments; your model has 11 parameters and also has to account for the second moment of the fixed exogenous UNIT. Thus 10 - (11 + 1) = -2, and the model is underidentified.

Some additional comments:

(1) Aside from UNIT there are no exogenous variables in this model.

(2) It's usually easier to use specifyEquations() than specifyModel() to define the model.

(3) If CF is really meant to be a latent variable, then there's a missing normalizing constraint, underidentifying the model apart from the gross comparison of moments and parameters.

(4) If, as here, you have the data, you don't have to compute the raw-moment matrix in a separate step, but could call sem() with the data argument.

Best,
 John

On Mon, 22 Jul 2013 18:55:31 -0400

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 model?
Thank you.
Here is the code:




model.ram1 <- specifyModel()
UNIT -> Y1, ty,0.3
UNIT -> Z1, tz1,-0.1
UNIT -> Z2, tz2,0.1
CF -> Y1, lamy,0.5
CF -> Z1, lamz1,0.85
CF -> Z2, lamz2,0.2
UNIT -> CF, k
Y1 <-> Y1, psi3, NA
Z1 <-> Z1, psi1, NA
Z2 <-> Z2, psi2, NA
CF <-> CF,vCF1,NA

sem.m1<-sem(model=model.ram1,S=S2,N=500,fixed.x="UNIT",raw=T)

-- 

I.J. Williams

Ph.D. Student in Education Measurement and Statistics
Statistics MS

Mathematics BS

Rutgers University

	[[alternative HTML version deleted]]

______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

------------------------------------------------
John Fox
Sen. William McMaster Prof. of Social Statistics
Department of Sociology
McMaster University
Hamilton, Ontario, Canada
http://socserv.mcmaster.ca/jfox/