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

st: permutations

From   Christopher Intemann <>
Subject   st: permutations
Date   Thu, 13 Mar 2008 10:46:43 +0100

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,
* For searches and help try:

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