clear all
set obs 101
generate x = (_n-51)/10
generate pr = invlogit(x)
generate sd = sqrt(p*(1-p))
line pr x, name(g1)

// estimated as a blocked logit
generate pos = round(pr*100) // rounds exact values
generate pop = 100
glm pos x, family(binomial pop) link(logit)
estimates store orig

replace x = (_n-51)/20
replace pr = invlogit(x)
line pr x, name(g2)
replace pos = round(pr*100) // rounds exact values
*generate pop = 100
glm pos x, family(binomial pop) link(logit)
estimates store half

estimates table orig ., se

replace x = (_n-51)/50
replace pr = invlogit(x)
line pr x, name(g3)
replace pos = round(pr*100) // rounds exact values
*generate pop = 100
glm pos x, family(binomial pop) link(logit)
estimates store fifth
graph combine g1 g2 g3, xcommon ycommon

estimates table orig half fifth, se