Quick sanity-check for thresholded and survey-weighted models

Hi all,

I'm trying to run a set of models with a binary outcome variable that also incorporates some survey weights. I'm pretty sure that I'm misspecifying something somewhere (my standard errors are madness), but I'm not 100% sure where. 

In building things out, I'm getting a few worrying error messages: 

-when I try to run the umxThresholdMatrix, I'm getting an "Error in rep(TRUE, nThreshThisVar) : invalid 'times' argument"

-when I try to run the AutoStart, I'm getting an "Error: In model 'PMTMod' column 'unsecureOrd' is not an ordered factor. Use mxFactor() on this column."

-when I run the refModels, I'm getting an "Error: In model 'Saturated PMTMod' column 'unsecureOrd' is not an ordered factor. Use mxFactor() on this column."

Any chance one of you all could look over the code below to help me identify the issue?

Thanks!

******

pmt.std<-na.omit(pmt.std) %>%
 zap_labels()

pmt.std$unsecureOrd<-mxFactor(pmt.std$unsecure, levels=c(0,1))

unsecureThreshold<-umx::umxThresholdMatrix(pmt.std, fullVarNames="unsecureOrd", method="Mehta")

PMTMod<-mxModel(model="PMTMod", type = "RAM",
               manifestVars = c(ta.manifests, ca.manifests, pm.manifests, "unsecureOrd"),
               latentVars = c("ThreatAppraisal", "CopingAppraisal", "ProtectionMotivation"),
               mxPath(from = "ThreatAppraisal", to = ta.manifests, arrows = 1, free = T, values = .1),
               mxPath(from = "ThreatAppraisal", to = ta.manifests[1], arrows = 1, free = F, values = 1),
               mxPath(from = "CopingAppraisal", to = ca.manifests, arrows = 1, free = T, values = .1),
               mxPath(from = "CopingAppraisal", to = ca.manifests[1], arrows = 1, free = F, values = 1),
               mxPath(from = "ProtectionMotivation", to = pm.manifests, arrows = 1, free = T, values = .1),
               mxPath(from = "ProtectionMotivation", to = pm.manifests[1], arrows = 1, free = F, values = 1),
               
               mxPath(from = "ThreatAppraisal", to = "ProtectionMotivation", arrows = 1, free = T, values = .1, labels = "ThA"),
               mxPath(from = "CopingAppraisal", to = "ProtectionMotivation", arrows = 1, free = T, values = .1, labels = "CoA"),
               mxPath(from = "ThreatAppraisal", to = "CopingAppraisal", arrows = 2, free = T, values = .5, labels= "TCr"),
               
               mxPath(from = "ProtectionMotivation", to = "unsecureOrd", arrows = 1, free = T, values = .1, labels = "B"),
               mxPath(from = "ThreatAppraisal", to = "unsecureOrd", arrows = 1, free = T, values = .1, labels = "ThC"),
               mxPath(from = "CopingAppraisal", to = "unsecureOrd", arrows = 1, free = T, values = .1, labels = "CoC"),
               
               mxThreshold(vars="unsecureOrd", nThresh=1, free=T, values=0),
               
               mxPath(from='one', to="unsecureOrd", free=FALSE, arrows=1, values=0),
               mxPath(from="unsecureOrd", free=FALSE, arrows=2, values=1),
               
               mxAlgebra(ThA * B, name = "ThAB"),
               mxAlgebra(CoA * B, name = "CoAB"),
               mxAlgebra(ThA * B + CoA * B, name = "ToTAB"),
               
               mxPath(from=c(ta.manifests, ca.manifests, pm.manifests, "ThreatAppraisal", "CopingAppraisal", "ProtectionMotivation"), arrows=2, free=TRUE, values=.8, lbound=.0001),
               mxPath(from='one', to=c(ta.manifests, ca.manifests, pm.manifests), arrows=1, free=TRUE, values=0),
               mxData(observed = pmt.std, type = "raw", weight = "weight"),
               mxCI(c("ThAB", "CoAB", "TCr", "B", "ThC", "CoC", "ToTAB"))
               )

PMTAuto<-mxAutoStart(PMTMod)
PMTOut<-mxTryHard(PMTAuto, intervals = T)
PMTRef<-mxRefModels(PMTOut, run = T)
PMTFull<-summary(PMTOut, refModels=PMTRef)