Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
| From | Steve Samuels <sjsamuels@gmail.com> |
| To | statalist@hsphsun2.harvard.edu |
| Subject | Re: st: RE: RE: Using command svy glm to obtain risk ratios |
| Date | Tue, 18 Feb 2014 01:41:23 -0500 |
I replied to Jose Maria privately with copies of the do files. The paper was: Waltermaurer, E, M Butsashvili, N Avaliani, S Samuels, and LA McNutt. 2013. An examination of domestic partner violence and its justification in the Republic of Georgia. BMC Womens Health 13, no. 1: 44. and is at: http://www.biomedcentral.com/1472-6874/13/44 It contained only the -svy: glm- analyses. To correct an error: In some instances of -svy: poisson-, the CI for a predicted *risk*, not the estimated RR, exceeded 1.0. log relative risk regression was the choice of my co-authors. My personal choice would have been -svy: logistic-, followed by -margins- for interpretation of risks. Steve sjsamuels@gmail.com On Feb 17, 2014, at 9:49 PM, jose maria pacheco de souza <jmpsouza@usp.br> wrote: Em 17/02/2014 22:40, Steve Samuels escreveu: > Amena: > > I recently co-authored a paper for which I first ran > -svy: poisson- > > and then > -svy: glm, link(log) > > which failed in about 2 of 20 models. > > I found four or five instances in which the relative risks differed so much that we had to abandon the idea of the Poisson approach. In others, the CIs for the RR exceeded 1.0 We were able to get -svy: glm, link(log) to work by slight modification of the models. > > You still have not shown us the actual commands you ran and the problematic -svy: glm- output. > Please do so. Right now we have no actual evidence of the problem.. The only difference that I can see in the output of -glm- and -svy: glm, with the "eform" option is that -glm- heads the coefficient table "Risk Ratio" whereas -svy: glm- has a heading "exp(b)", which is the same thing. > > > Steve Samuels > sjsamuels@gmail.com > Dear Steve: could you give the reference of your paper and the complete Stata code command? I am very interested in this subject. Thank you. josé maria -- Jose Maria Pacheco de Souza Professor Titular (aposentado), Professor Senior Departamento de Epidemiologia/Faculdade de Saude Publica, USP Av. Dr. Arnaldo, 715 01246-904 - S. Paulo/SP - Brasil fone (11)3768-8612 www.fsp.usp.br/~jmpsouza * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/ * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/