Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: permutations


From   Christopher Intemann <intemann@gmail.com>
To   statalist@hsphsun2.harvard.edu
Subject   st: permutations
Date   Thu, 13 Mar 2008 10:46:43 +0100

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/




© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index