cd "c:\mijn documenten\projecten\chapter_UH\haifa
use gss.dta, clear
desc
sum

sum sei
gen z_sei = (sei - r(mean)) / r(sd)

fre degree

gen ed01 = degree > 0 & degree < .
gen ed12 = degree > 1 if  ed01==1 & degree < .
tab ed01 ed12, missing
gen ed23 = degree > 2 if  ed12==1 & degree < .

gen baseline = 1

logit ed01 south sibs black z_sei baseline, or nocons
logit ed12 south sibs black z_sei baseline, or nocons
logit ed23 south sibs black z_sei baseline, or nocons

help seqlogit


seqlogit degree south sibs black baseline, nocons        ///
         tree(0 : 1 2 3, 1 : 2 3, 2 : 3) ///
		 ofinterest(z_sei) or 
est store null 
		 
seqlogit, or

seqsens [#1]_b[z_sei] using sens, replace sd(1, 2, 3, 4) rho(0)		 
seqsens [#2]_b[z_sei] using sens, append  	
seqsens [#3]_b[z_sei] using sens, append  

help test
test [#1=#2=#3]:z_sei

seqsens r(p) using sens, append test(test [#1=#2=#3]:z_sei)

seqsens using sens, describe

help estimates 

estimates use sens, number(4)
estimates esample:

seqlogit

uhdesc

/// sugestions:

est restore null
seqsens [#3]_b[z_sei] using sens, append rho(-.4 -.2 0 .2 .4)
seqsens [#3]_b[z_sei] using sens, append pr(.5 .5)
seqsens [#3]_b[z_sei] using sens, append pr(.2 .2 .2 .2 .2)
seqsens [#3]_b[z_sei] using sens, append delta(1.5)
seqsens [#3]_b[z_sei] using sens, append sd(1 2 2, 1 2 3)