 |
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]
st: permutations
Hi,
I just realized that stata is capable of doing permutations.
However, I find it quite difficult to figure out how:-)
Within the scope of a genetic analysis I'm using a command like:
xi:logistic affected_stata i.gene age ethno_x sex_x
where gene is coded as 1/2/3, referring to genotypes AA/AG/GG (for
instance).
The result is corrected for confounders age, ehtnicy and sex.
I played around a bit. I was trying this command line:
permute affected_stata "xi:logistic affected_stata i.MCP1_7656en age
ethno_x sex_x" e(chi2)
and I got a result which is probably wrong (or at least I don't
understand it):
Monte Carlo permutation statistics Number of obs
= 4356
Replications
= 100
------------------------------------------------------------------------------
T | T(obs) c n p=c/n SE(p) [95% Conf.
Interval]
-------------
+----------------------------------------------------------------
_pm_1 | 76.09804 0 100 0.0000 0.0000 0 .
0362167
------------------------------------------------------------------------------
Well. Any help would be appreciated:-)
Many thanks in advance,
Chris
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/