Hi,

I want to do logistic regression based on a complex sample design. I used
package survey, but when I ran svyglm, an error message came out:
Error in onestrat(x[index, , drop = FALSE], clusters[index],
nPSU[index][1],  :
  Stratum (16) has only one PSU at stage 1


My code is below:

a.design<-svydesign(id = ~CASENUM ,strata = ~STRATUM ,data = a ,weights =
~SIZAGYWT )
summary(logistic1 <- svyglm(ANYCONTR ~ CHAIN+OWN+HPPAT, family =
binomial(link = "logit"), design=a.design))


How to solve this issue? Thank you.

	[[alternative HTML version deleted]]


