clear all
set obs 200
forval i=1/4 {
	gen x`i'=rnormal()
}

egen ylin=rowtotal(x*)
gen y=runiform()<invlogit(ylin)

logit y x*
margins, post
estimates store full
logit y x1
margins, post
estimates table full ., se

logit y x*
margins, at(x1=(-2/2)) post
estimates store full
logit y x1
margins, at(x1=(-2/2)) post
estimates table full ., se

logit y x*
margins, dydx(x1) post
estimates store full
logit y x1
margins, dydx(x1) post
estimates table full ., se