// Shift cutpoint but not x
clear all
set obs 5000
generate ru = runiform()
generate x = runiform((-1),1)
generate x1 = x
generate y = .
postfile logitvar shift beta0 se0 beta1 se1 using "logitvar.dta", replace
forvalues i = -1(0.05)1 {
	di "shift= `i'"
	replace x1 = x+`i'
	replace y = ru < invlogit(x1)
	quietly logit y x
	post logitvar (`i') (_b[_cons]) (_se[_cons]) (_b[x]) (_se[x])
	}

postclose logitvar

use logitvar.dta, clear
line se0 beta0, xline(0) title("Constant") name(constant)
line se1 beta1, xline(1) xlabel(0.75(0.125)1.25) title("Discrimination") name(beta)
graph combine constant beta, title("Shifted Cutpoint") ///
	note("at constant x")