So, in short, when I use "ordered = c(NAMES OF OUTCOME VARIABLES)" the output does change when I am not using survey weights with the lavaan.survey package. But the output is exactly the same whether I include the ordered= argument or not in the fit variable when using lavaan.survey. I need to add survey weights and also want to use ordered regression. Does anyone know how to do this?
#######################################################################################
Model<- '
# latent variables
Dev1 =~ dev11 + dev12
Dev2 =~ dev21 + dev22
Fam =~ fam1 + fam2
# regressions
mental ~ Fam + SEX + income + education + Dev1 + Dev2 + interference
physical ~ Fam + SEX + income + education + Dev1 + Dev2 + interference
'
fit<- sem(Model, data = gss, ordered = c("mental", "physical "))
svy.df<-svydesign(id=RECID, weights=WGHT_PER, data=gss)
survey.fit <- lavaan.survey(fit, svy.df)
summary(survey.fit, fit.measures=TRUE, standardized=TRUE, rsquare=TRUE)
##########################################################################################
I can't seem to get the ordered= argument to do anything when I'm using survey weights with lavaan.survey. Any ideas and advice are appreciated.
Thank you