Originally posted by Horea Feier
View Post
Code:
version 15.1
clear *
set seed `=strreverse("1436437")'
quietly set obs 500
generate double pid_u = rnormal()
generate int pid = _n
quietly expand 5
generate byte outcome = rbinomial(1, invlogit(pid_u))
melogit outcome || pid: , nolog
exit

Comment