Any comments are appreciated!
*/-----------START-----------------
input _study _status _count _genotype_d
1 1 23 0
1 1 118 1
1 0 13 0
1 0 126 1
2 1 4 0
2 1 49 1
2 0 8 0
2 0 53 1
3 1 10 0
3 1 101 1
3 0 9 0
3 0 91 1
4 1 6 0
4 1 70 1
4 0 4 0
4 0 49 1
end
*/--------END-----------------------
local studies = 4
xi: logistic _status i._study i._genotype_d [freq=_count]
qui esttab, se
matrix TEST = r(coefs)
scalar OR_genotype = exp(TEST[`studies',1] )
scalar SE_genotype = TEST[`studies',2 ]
scalar p_genotype = TEST[`studies',3 ]
scalar expected_z_statistic = ln(OR_genotype)/SE_genotype
scalar expected_p_two_sided_p_value = normprob(-abs(expected_z_statistic))*2
scalar list
dis as err "Are the standard errors from my -logistic- output wrong?"
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
