I am fitting a 3 level cross-classified model: level 1 is patient ID (ca.40000), the other 2 levels are physician ID (ca.7000) and region ID where patients live in (ca.700). Outcome variable is binary so this is a logistic model. Below is the code:

Code: Select all

```
runmlwin y sex age cons, level3(region:cons) level2(physician:cons) level1(id) discrete(distribution(binomial) link(logit) denominator(cons)) or nopause
runmlwin y sex age cons, level3(region:cons) level2(physician:cons) level1(id) discrete(distribution(binomial) link(logit) denominator(cons)) or mcmc(on) initsprevious nopause nogroup
```

Code: Select all

```
runmlwin y sex age cons, level3(physician:cons) level2(region:cons) level1(id) discrete(distribution(binomial) link(logit) denominator(cons)) or nopause
runmlwin y sex age cons, level3(physician:cons) level2(region:cons) level1(id) discrete(distribution(binomial) link(logit) denominator(cons)) or mcmc(on) initsprevious nopause nogroup
```

It would be great if anyone could help. Thank you indeed!

Best,

wj