* for estimation of models using first wave of 2 wave panel data from BSA
* March 2 2008
* sex 1 = male 2 = female
* cohort is age group (19/25=1)(26/35=2)(36/45=3)(46/55=4)(56/65=5)(65/max=6)
* Y is opinion of premarital sex
* always  wrong, mostly  wrong, sometimes wrong, rarely wrong, not wrong

* this is mixture model with diff location and heterog thresholds 

lat 1
man 3
dim 2 2 5 6 
lab X S Y C 
mod X {X} 
    Y|XCS cum (a) {eff(Y,4,X,b) eff(C,5) eff(S,1) spe(X,1a) }  nin
iter 1000000
nec nfr
dum 1 1 5 1
des [
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
0 0 0 1 0
0 1 0 0 0 0
0 0 1 0 0 0
0 0 0 1 0 0
0 0 0 0 1 0
0 0 0 0 0 1
0 1
]
sta eff(Y) [1 2 3 4 ]
dat [
1	0	3	4	7	6
0	0	3	7	4	6
2	9	8	2	10	4
3	4	5	5	2	2
26	34	26	19	13	10
1	2	4	11	8	25
0	2	3	10	16	14
8	7	11	13	7	10
6	5	3	3	2	1
20	34	30	16	8	6
]