Dear statalist
I want to estimate the risk for an exposed twin to develop a symptom
relative an unexposed twin. I have 169 pairs of monozygotic twins that are
discordant on my exposure variable. The basic design is a case (exposed) -
control (unexposed) study. The outcome variable (symptom) is binary coded.
The most obvious method should be -clogit- however, I'm not satisfied with
the statistics.
Like McNemar, it excludes the ties between cases and controls and the odds
ratio returned describes the ratio of cases "winning" over controls.
However, I would like to estimate the ratio of a case to HAVE the symptom
(not "winning") compared to the control.
I have been criticized for using logistic regression with the cluster option
and robust variance estimator. My OR is said to be invalid.
Actually, the criticism is harder than that: My analyses is said to be
"completely meaningless".
. logistic dv iv , cluster(twinpairid)
Do you have any advice?
I do not know your supervisor. However, if you want a risk ratio instead of
an odds ratio, and you want to cluster the standard errors by twin pair,
then you should use -glm- to fit a model with log link and binomial error.
I would type