Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: RE: RE: Using command svy glm to obtain risk ratios


From   Steve Samuels <[email protected]>
To   [email protected]
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
[email protected]

On Feb 17, 2014, at 9:49 PM, jose maria pacheco de souza <[email protected]> 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
> [email protected]
> 
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/


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index